run this with dskeqdata generated from eqdtoq.x
 &all
 rmax=65.  rzero=170. ; major and minor radii
 npsi=67  nthet=65 ; mesh size
 bavq=0.0   bava=0.0  bavxz=0.0  nbavxz=6000 ; back-averaging
 equiltype='ffprime'
 modelp=2 ; read p' from dskeqdata
 modelf=2 ; read ff' from dskeqdata
 modelbnd=2 ; read boundary points from dskeqdata
; baxis0=3.508e4 ; not used with  equiltype='qsolver'
; totcur=1.571e6   ; not used with  equiltype='qsolver' 
 fixcur=.f. ; don't fix total current, i.e., don't scale ff'
 prtboot=.f.
 prteqdata=.t. ; generate dskeqdatao
 prtfixb=.f. ; set to true if you want to generate dskfixb
             ; dskfixb can be read by fixbnd.x to generate eqdsk
 prtgato=.f.
 prtmercy=.f.
 prtbal=.f.
 betafix="none" 
 betiso=11.0e-2  ; not used since betafix="none"
 ieqdsk=0 ;if =2 use dskeqi to get initial guess
          ; =0 don't read dskeqi
          ; dskeqi is dskeqo generated from a previous run
 iteqmx=120	 ; will iterate (regrid)  120 times before it gives up      
 imislo=0 ; sometimes data near axis is garbage--in that case set
           ; imislo=3 to use extrapolation to axis which should 
           ; clean it up
 toleq=1.e-3	 ; convergence criterion--I almost never change this
 alpsi=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
  nhior=15	nbndry=199 ; don't ask you don't want to know

 &end