mapv takes the refined equilibrium data eqb1 from Jsolver,
and output two files mpout1 and mapdsk for nova.
It is available through lmhd svn repository
svn+ssh://svnsrv.pppl.gov/svn/mhd/trunk/mapv.
The syntax is
mapvThe preallocated fixed array sizes are given as nq for theta grids and ny for radial grids in file clichpar.h :
prameter(ns=nmbrgrd,nt0=mnths0,nt=nt0+5, nq=526,ny=526,nsgrd=526)
The meaning of nq and ny are the same as defined in
mapck.
The actual array sizes are changed at runtime in namelist comap
&dfun / p &wdat nosuri=201, mthi=128, igrid=1 /or
&dfun / p &wdat nosuri=129, mthi=128, igrid=1 /and these changes should be matched in file gridparam:
parameter(nmbrgrd=201,mnths0=128,ndat=50)
or
parameter(nmbrgrd=129,mnths0=128,ndat=50)
Detail:
dmapb2.f
45 common ntitle(20),dat,nxx(5),axx(13),nxy(10),axy(10)
46 common/var/ p(ny),pp(ny),q(ny),qp(ny),
47 1 g(ny),gp(ny),fb(ny),fbp(ny),f(ny),fp(ny),
48 1 psival(ny),x(nq,ny),z(nq,ny),aj3(nq,ny),aj(nq,ny)
1. array passing from Jsolver
133 write(*,*)'nxx',nxx
nxx 128 129 21 40 0
2. array used to work on equilibrium data
157 nthe=nxx(1)
158 npsi=nxx(2)
3. array (comap) used to interpolate equilibrium data and passed to nova code
253 write ( *, wdat )
254 write(*,*) "WDAT NAMELIST"
255 read ( *, wdat )
256 write ( 16, wdat )
271 nosurf=nosuri
272 mth=mthi
273 nxx(3)=nosurf
274 nxx(4)=mth