file='rho=1.2.sav'
file='rho=1.4.sav'
file='rho=1.8.sav'
file='rho=2.2.sav'
;
;  produce z array
;
s=size(map)
nz=s[1]
Lz=.2
dz=Lz/nz
zvid=dz*findgen(nz)-.1+dz/2.
;
save,file=file,tt,ppm,rhom,TTm,tvid,map,zvid
END