&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. /