LAT A=1.4--kappa=3.0000; delta=0.40000
 &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=67 	nthet=65	; mesh size-- psi & theta points
 ;  does substitutions on these values-- acceptable values are
 ;  npsi=19  nthet=17 
 ;  npsi=35  nthet=33 
 ;  npsi=67  nthet=65 
 ;  npsi=131  nthet=129 
  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.99 ; determines the relationship between psic--psi coordinate
           ;  and psiv--psi value=poloidal flux(*2pi)
           ; alpsi=0 makes psic = sqrt(psiv) normalized to 0-1
    ; for alpsi<0
    ; psi_norm(i)=(sin(psic(i)*alpsi*pi*0.5))**2/sin(-alpsi*pi*0.5)**2
    ; idea for alpsi~-1 is to put lots of surfaces near axis and 
    ; near edge
 modelbnd=1  ; this says to use one of the shape formulas below to get
              ; the outermost boundary      
  ; modelbnd=2 get boundary shape from dskeqdata file
  ; for dee shape with elongation=eshape and triangularity=sin(xshape)
  ;     ishape=4 eshape=2.0 xshape=0.5236 to get triangularity=0.5
  ; set ishape=0 to get a circle
  ishape=4   ashape=0.    eshape=3.0000	xshape=.5236
  equiltype='jdotb'   ; this says use j.b and p' to determine equilibrium
                      ; other equiltype values are 
                      ; 'qsolver' and 'ffprime'
  baxis0=2.00e4       ; vacuum b field in gauss
  fixcur=.false.      ; what happens if you set to true??
  ; following set of parameters determine j.b due to current drive "cd"
  ncdcoeff=7  ; cd will be a polynomial in x=psi_norm with 7 terms  
  cdcoeff=0.075 0.0 0.0 0.0 0.925 0.0 -1.00 ; here are the 7 coefficients
              ; 0.075+0.925*psi_norm^4-1.00*psi_norm^6  
  cdx=0.50  cddelx=0.04  ; cd polynomial multiplied by
                         ;( 1-tanh((psi_norm-cdx)/cddelx))/2
			 ; i.e., here it is shut off beyond ~ psi_norm~0.54
  cdnorm=1.2e7  ; cdnorm sets j.b on axis, i.e., raising cdnorm will
                ; lower q axis
  cdfreeze=.true.  ncdfreeze=12  ; don't ask and don't change
  ; following parameter is  what make this 'jdotb' run into a 
  ;   bootstrap current case
  jbs=.true.	
  jbtype='combined' ; says if bootstrap current exceeds jdotb on a given
                    ; flux surface--just use jdotb  
  zeff=2. gamboot=0.50000 ; bootstrap model params
                          ; gamboot=L_p/L_T the ratio
                          ; of pressure to temperature scale lenghts
  ifillboot=5             ; bootstrap used only beyond 5th mesh point
                          ; from axis 
  betafix='bt'	          ; will adjust magnitude of pressure to achieve
                          ; beta on axis value of betiso.
                          ; if betafix='none' pressure not modified
  betiso=0.75             ; desired beta on axis
  ieqdsk=0  ;if =2 use dskeqi to get initial guess
            ; =0 don't read dskeqi
            ; dskeqi is dskeqo generated from a previous run
  prtbal=.true.	  ; print out file dskbalnew for baloo
  prtboot=.true.  ; print out file dskboot for bootstrap current analysis
  prtgato=.false. ; print out file dskgato for gato
  prteqdata=.true. ; print out file dskeqdatao for iput to another toq run
                   ; perhaps after modifying p,p',f,f',q or boundary
                   ; to use in a later run, change dskeqdatao to dskeqdata
                   ; and set desired model parameter to 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
                   ; for other modelf options see routine sigset.f
                   ; modelbnd=2 get boundary shape from dskeqdata file
                   ; modelq=2 get q profile from dskeqdata file
  prtfixb=.true.   ; print out dskfixb file which can be used with
                   ; fixbnd.x to generate an eqdsk
  ; the following are control parameters on multi-grid algorithm "mg"
  bavmg=0.5	minmg=3	maxmg=6	tolmg=1.e-9	loopinmg=2     
  ; specification of pressure profile in this run uses a spline
  ; xppspline are psi_norm locations of knots
  modelp=  5  ; means use spline option for specifying pprime
	      ; see routine psetup.f for other options
  nppspline= 11     
  xppspline=  0.0000  0.1000  0.2000  0.3000  0.4000  
  0.5000  0.6000  0.7000  0.8000  0.9000    
  1.0000        
  ppspline=  2.5000E-02  2.5875E-02  3.1200E-02  4.3225E-02  6.1800E-02    
  8.4375E-02  1.0600E-01  1.1933E-01  1.1460E-01  7.9675E-02  
  0.0000E+00      
  nhior=15	nbndry=199 ; don't ask you don't want to know
  nbug(1)=1      ; controls some multigrid printout--set to zero to get
                 ; less output to stdout
  toleq=1.e-3	 ; convergence criterion--I almost never change this
  iteqmx=120	 ; will iterate (regrid)  120 times before it gives up      
  &end