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