[ncl-talk] ncl-talk Digest, Vol 161, Issue 41(1) Pressure level (Karin Meier-Fleischer)

[WIND]

Thanks Karin,
   It's not the whole script, I choose the important part to show you. I just check VarSummary, 
   level variable is listed below:
   standard_name :  air_pressure
   long_name :   vertical level
   units :       hPa
   positive :    down
   axis :        Z
(0)      200
(1)      150
(2)      100
(3)       70
(4)       50
(5)       40
(6)       30
(7)       20
(8)       10
(9)        7
(10)       5
(11)       4
(12)       3
(13)       2
(14)       1


    and variable  merra's VarSummary is listed below:
    Dimensions and sizes:   [lev | 15] x [lat | 73]
    Coordinates:
            lev: [ 200..   1]
            lat: [  90.. -90]


But I still don't know how to show level axis above 10hPa 
via “gsn_csm_pres_hgt” to show Pressure log coordinate system.


Yan

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170427/2e439787/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: my example.jpg
Type: application/octet-stream
Size: 613107 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170427/2e439787/attachment-0001.obj