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=65.  rzero=170. ; 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 
 bavq=0.6   bava=0.6  bavxz=0.6  nbavxz=6000 ; back-averaging
      ; back-averaging usually needed for equiltype='qsolver'
 equiltype='qsolver'
 ; 'qsolver' is far and away the least robust option but
 ; it is the one that allows you to specify q profile 
 ; it's also the only option which needs the runtoq.sc script
 modelp=2 modelf=2 modelbnd=2 modelq=2 
  ; modelx=2 means get arrays from dskeqdata file
  ; modelp=2  get pressure profile from dskeqdata
  ; for other modelp options see routine psetup.f
  ; modelf=2  get f and ff'  profile from dskeqdata--not needed for 
  ;                                                equiltype='qsolver'
  ; for other modelf options see routine sigset.f
  ; modelbnd=2 get boundary shape from dskeqdata file
  ; for other boundary options (modelbnd=1) see routine initbndry.f
  ; if modelbnd=1 then shape params determine boundary
  ; for dee shape with elongation=eshape and triangularity=sin(xshape)
  ;     ishape=4 eshape=2.0 xshape=0.5
  ; set ishape=0 to get a circle
  ; modelq=2 get q profile from dskeqdata file 
  ; for other q options see routine qsetup.f
  ; e.g., if modelq=3 use spline data as below
  ;nqspline=9
  ;qspline = 3.85 2.338 2.189 2.263 2.248 3.009 4.420 5.973
  ;xqspline = 0.0 .09375 .1406 .1875 .5938 .7969 .9063 .9531
  ;qpaxspline=-2. qplimspline=15.0
   jlimzero=0
 dchi=2.493e7 ; adjust to get desired baxis0 or totcur if equiltype='qsolver'
; baxis0=3.508e4 ; not used with  equiltype='qsolver'
; totcur=1.571e6   ; not used with  equiltype='qsolver' 
 fixcur=.f. 
 prteqdata=.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"
 ; the following specify pressure model if modelp=1
 ;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
 ; q profile if modelq=1
 ; qax=3.9  qaxp=-1.0  qlim=5.2  qlimp=20. ; alphaq=6.0
 ;  gammaq=0.9  iqval=0  jlimzero=1
 ; 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