this first line is a comment line and is required
 &all
 ; this namelist is read in routine init.f
 ; Note: all TOQ units are cgs--Sorry about that
 rmax=100.  rzero=140. ; minor and major radii 
 npsi=xx  nthet=xx ; mesh size-- set to xx here because script runtoq.sc
 ;  does substitutions on these values-- acceptable values are
 ;  npsi=19  nthet=17 
 ;  npsi=35  nthet=33 
 ;  npsi=67  nthet=65 
 ;  npsi=131  nthet=129 
 ishape=4 eshape=1.5 xshape=0.5
 bavq=0.6   bava=0.6  bavxz=0.6  nbavxz=6000 ; back-averaging
      ; back-averaging usually needed for equiltype='qsolver'
 equiltype='qsolver'
 ;modelq=3 
 ; nqspline=5
 ; qspline = 0.5 0.375 0.25 0.125 5.e-4
 ; xqspline = 0.0 0.25 0.5 0.75 1.0
 ; qpaxspline=-0.5 qplimspline=-0.5
   jlimzero=0
 dchi=2.493e7 ; adjust to get desired baxis0 or totcur if equiltype='qsolver'
 prteqdata=.t.
 prtmercy=.t.
 betafix="none"  ; to rescale pressure profile to get betiso on axis
                 ;    set betafix="bt"
 betiso=11.0e-2  ; beta on axis--not used if betafix="none"
  modelp=2
 ncent=1. nedge=0.0 nexpin=1. nexpout=5.00 ; density
 tcent=1. tedge=1.0 texpin=1. texpout=1.0 ; temperature
 rcent=0. redge=0. ; rotation
 ieqdsk=0 ;if =2 use dskeqi to get initial guess
          ; =0 don't read dskeqi
          ; dskeqi is dskeqo generated from a previous run
   modelq=1
  qax=0.9  qaxp=-0.9  qlim=5.e-3  qlimp=-0.9  alphaq=3.0
 iqval=0
 ; debugging parameters--ignore
 ;bug=0. 0. 0. 0. 0. 0. 0. 0. 0. 2.5 0.25 
 ;bug1=0.  bug2=0.    bug3=0. 
  imislo=0 ; sometimes data near axis is garbage--in that case set
           ; imislo=3 to use extrapolation to axis which should 
           ; clean it up
 alpsi=0.0 ; determines the relationship between psic--psi coordinate
           ;  and psiv--psi value=poloidal flux(*2pi)
           ; alpsi=0 makes psic = sqrt(psiv) normalized to 0-1
           ; set routine initpsi.f for other options
 ; I don't think you'll need to change following
 iteqmx=1000 ; will iterate (regrid)  1000 times before it gives up
             ; this is a lot
 toleq=1.e-3 ; convergence criterion--I almost never change this
 
 &end