[ncl-talk] Pressure level
Karin Meier-Fleischer
meier-fleischer at dkrz.de
Wed Apr 26 07:57:28 MDT 2017
Hi,
if this is your full script comment the lines res at gsnDraw and
res at gsnFrame because you wouldn't get any plot.
Please, do a printVarSummary of your level variable (e.g. lev or level).
Is the unit attribute "hPa"? Do also a print of the level variable to
see which values are stored.
You don't have to add anything your variable should have all
information. Do also a printVarSummary(merra).
Bye,
Karin
Am 26.04.17 um 04:41 schrieb WIND:
> Hi all,
> I'm trying to draw the pressure v.s. latitude contour via
> "gsn_csm_pres_hgt", but the pressure axis just don't show pressure
> data above 10hPa (from 10hPa to 1hPa). (Besides, I want to emphasize
> the upper stratosphere data.)
> How should I add the whole pressure data?
> Here is my script
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> fq=addfile("/merra2/200-1/MERRA2.nc4","r")
> merra=fq->T
> wks = gsn_open_wks("eps" ,"t")
> res = True
> res at cnFillOn = True
> res at tmYRMode = "Automatic"
> res at gsnDraw = False
> res at gsnFrame = False
> res at tmXBMinorOn = False
> res at tiYAxisString="Pressure"
> plot = gsn_csm_pres_hgt(wks,merra,res)
>
> Any help would be grateful!
>
> Yan
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170426/da147c61/attachment.html
More information about the ncl-talk
mailing list