&testgridgen
 source="eik.out",
 nperiodout=1,  ! output nperiod
 nthetaout=32,  ! output ntheta
 nlambdaout=30, ! half of ntheta, plus 6 (and it will find an optimal
! nthetaout=16,  ! output ntheta
! nlambdaout=14, ! half of ntheta, plus 6 (and it will find an optimal
                ! number of pitch angles mu/E to get bounce points right).
 !npadd=4,       ! # of points added locally when hunting for local wells
 npadd=0,       ! # of points added locally when hunting for local wells
                ! npadd=1 or 2 may work for me.
 alknob=0.0,   ! how lambda's distributed, 
 epsknob=1e-5,  ! figure out later
 smoothknob=0.,  ! something to do with B(theta) smoothing

! additional parameters in namelist in rungridgen.f90, with their defaults:
! extrknob = 10.0
! nfinegrid=200
! gingrid = "gingrid"

! Increase theta resolution by a factor of 1+deltaw, in a region with a 
! Lorentzian width of +-widthw, centered around +thetamax, and another
! centered around -thetamax (see tfunc in gridgen4.f90):

 deltaw=0.,     ! 0 to give uniform weights, and (almost) uniform theta grid
 widthw=1.0,    ! width of Lorentzians
 thetamax=0.00, ! usually 0 to focus resolution there.

 tension=1.0,  ! for splines, leave at 1

 screenout=.false.,  ! .true. to print more info to screen

 auto_width=.false.,  ! use the above parameters to optimize grid
!auto_width=.true.,   ! use following 3? parameters to optimize grid
                      ! (ignores the above parameters)

 cv_fraction=0.6,   ! fraction of theta grid points you want to be in the bad
                    ! curvature region

 delth_max=0.5,   ! estimate of maximum delta(theta) allowed.

 ! trying to balance between putting lots of points in the bad curvature
 ! region, and not let theta grid points get too spread out in good
 ! curvature region

 max_autoiter=4,  ! max # of iterations to get auto width.

/
