[ncl-talk] Pressure Vertical Axis
Soma Roy
somaroy892 at gmail.com
Wed Jul 25 13:25:16 MDT 2018
You can use pressure height vector in this case.
Hope this may help.
Thanks.
On Thu, Jul 26, 2018, 00:28 MELISSA KAZEMI RAD <
mk1369 at scarletmail.rutgers.edu> wrote:
> Hi,
>
> I am trying to create a vertical cross_section plot from a wrf output
> file. I have spent a lot of time looking through all example ncl codes,
> trying to figure out how to have the y-axis, representing the vertical
> height show pressure values in pa, rather than the number of eta-levels.
>
> Here is my NCL code:
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
> ;----------------------------------------------------------------------
> ; Main code
> ;----------------------------------------------------------------------
> begin
> a = addfile("wrfout_d03_2017-07-17_06_00_00.nc","r")
>
> ; type = "x11"
> ; type = "pdf"
> ; type = "ps"
> type = "png"
>
> wks = gsn_open_wks(type,"EQTemp_Press_Cross_166")
>
> time = 166
>
> p = wrf_user_getvar(a, "pressure", time) ; Total Pressure
> P = wrf_user_getvar(a, "P", time)
> PB = wrf_user_getvar(a, "PB",time)
> T = wrf_user_getvar(a,"T",time)
> U = wrf_user_getvar(a, "ua", time)
> V = wrf_user_getvar(a, "va", time)
> ; z = wrf_user_getvar(a,"z",time)
> qv = wrf_user_getvar(a, "QVAPOR", time)
>
>
> T_t = T + 290. ; potential temperature in K.
> P_t = P + PB
> tk = wrf_tk(P_t, T_t)
> eth = wrf_eth(qv, tk, P_t)
> plane = (/342, 1387/)
> t_plane = wrf_user_intrp3d(eth,p,"v",plane,50.,False)
> P_plane = wrf_user_intrp3d(P,p,"v",plane,50.,False)
>
> printVarSummary(t_plane)
>
> res = True
> res at MainTitle = "Vertical Cross-Section of Potential
> Temperature"
> res at InitTime = False
> res at TimePos = "Left"
> res at Footer = False
> res at lbTitleOn = True
> res at lbTitleString = "Temperature"
> res at tmYLMode = "Explicit"
> res at tmYLValues = (/1000, 850, 700, 500, 400, 300, 250,\
> 200, 150, 100, 50, 10/)
> res at tmYLLabels = "" + res at tmYLValues
> res at tmXBLabelFontHeightF = 0.015
> res at tmYLLabelFontHeightF = 0.015
>
>
> tres = res
> tres at cnFillOn = True
> tres at cnFillPalette = "BlWhRe"
> tres at cnSpreadcolor = True
>
>
> pres = res
> pres at cnFillOn = False
> pres at ContourParameters = (/ 2. /)
> pres at cnInfoLabelFontHeightF = 0.01
>
>
>
>
> T_contour = wrf_contour(a,wks,t_plane,tres)
> P_contour = wrf_contour(a,wks,P_plane,pres)
>
> pltres = True
> pltres at gsnMaximize = True
>
> plot = wrf_overlays(a,wks,(/T_contour, P_contour/), pltres)
>
> end
>
> I'd really appreciate your help.
>
> Bests,
> Melissa
>
>
>
>
> _______________________________________________
> 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/20180726/308bfd72/attachment.html>
More information about the ncl-talk
mailing list