pro zp
;
; zp = -ln(p/p0) -> -zp = ln(p/p0) -> p/p0 = exp(-zp) -> p(dynes/cm2) = p0*exp(-zp)
; pmb = p*.001
;
nlev=49
zp = [-17, -16.5, -16, -15.5, -15, -14.5, -14, -13.5, -13, -12.5, -12,$
    -11.5, -11, -10.5, -10, -9.5, -9, -8.5, -8, -7.5, -7, -6.5, -6, -5.5, -5,$
    -4.5, -4, -3.5, -3, -2.5, -2, -1.5, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3,$
    3.5, 4, 4.5, 5, 5.5, 6, 6.5, 7]
p0=5.e-4
for k=0,nlev-1 do begin
  pmb = p0*exp(-zp[k])*1.e-3  ; find pressure in mb
  lev = alog(p0/(pmb/1.e-3))  ; reverse to find zp level
  print,format="('k=',i3,' zp[k]=',f8.2,' pmb=',e12.4,' lev=',f8.2)",k,zp[k],pmb,lev
endfor

end