gridgen4_1 Subroutine

public subroutine gridgen4_1(n, nbmag, thetain, bmagin, npadd, alknob, epsknob, bpknob, extrknob, thetamax, deltaw, widthw, tension, ntheta, nlambda, thetagrid, bmaggrid, alambdagrid)

FIXME : Add documentation

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: n
integer, intent(in) :: nbmag
real, intent(in), dimension (nbmag) :: thetain
real, intent(in), dimension (nbmag) :: bmagin
integer, intent(in) :: npadd
real, intent(in) :: alknob
real, intent(in) :: epsknob
real, intent(in) :: bpknob
real, intent(in) :: extrknob
real, intent(in) :: thetamax
real, intent(in) :: deltaw
real, intent(in) :: widthw
real, intent(in) :: tension
integer, intent(inout) :: ntheta
integer, intent(inout) :: nlambda
real, intent(out), dimension (ntheta+1) :: thetagrid
real, intent(out), dimension (ntheta+1) :: bmaggrid
real, intent(out), dimension (nlambda) :: alambdagrid

Contents

Source Code


Source Code

  subroutine gridgen4_1 (n,nbmag,thetain,bmagin, npadd, &
       alknob,epsknob,bpknob,extrknob,thetamax,deltaw,widthw,tension, &
       ntheta,nlambda,thetagrid,bmaggrid,alambdagrid)
    implicit none
    integer, intent (in) :: n
    integer, intent (in) :: nbmag
    real, dimension (nbmag), intent (in) :: thetain, bmagin
    integer, intent (in) :: npadd
    real, intent (in) :: alknob, epsknob, bpknob, extrknob
    real, intent (in) :: thetamax, deltaw, widthw, tension
    integer, intent (in out) :: ntheta, nlambda
    real, dimension (ntheta+1), intent (out) :: thetagrid, bmaggrid
    real, dimension (nlambda), intent (out) :: alambdagrid
    
    call gridgen4_2 (n,nbmag,thetain,bmagin, npadd, &
         alknob,epsknob,bpknob,extrknob,thetamax,deltaw,widthw,tension, &
         ntheta,nlambda,thetagrid,bmaggrid,alambdagrid)
    alambdagrid(1:nlambda) = 1.0/alambdagrid(1:nlambda)
  end subroutine gridgen4_1