<div dir="ltr"><div class="gmail_default" style="font-size:small">Hi Melissa,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">To expand on what Soma suggested, you could visit the WRF online tutorial that was developed by the WRF developers:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www2.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/NCL_examples.htm">http://www2.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/NCL_examples.htm</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">In particular, look at the "vertical cross section" section</div><div class="gmail_default"><br></div><div class="gmail_default">Further, and I don't know how helpful this example would be, but you can also go to our cross section page and look at "narr_5.ncl". This shows how to do different cross-sections of curvilinear data, which is what you have:</div><div class="gmail_default"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/transect.shtml">http://www.ncl.ucar.edu/Applications/transect.shtml</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">You can still use wrf_user_getvar to retrieve the required variables.</div><div class="gmail_default"><br></div><div class="gmail_default">HTH,</div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div><div class="gmail_default"><br></div><div class="gmail_default"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jul 25, 2018 at 1:25 PM, Soma Roy <span dir="ltr"><<a href="mailto:somaroy892@gmail.com" target="_blank">somaroy892@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="auto">You can use pressure height vector in this case.<div dir="auto"><br></div><div dir="auto">Hope this may help.</div><div dir="auto"><br></div><div dir="auto">Thanks.</div></div><br><div class="gmail_quote"><div><div class="h5"><div dir="ltr">On Thu, Jul 26, 2018, 00:28 MELISSA KAZEMI RAD <<a href="mailto:mk1369@scarletmail.rutgers.edu" target="_blank">mk1369@scarletmail.rutgers.<wbr>edu</a>> wrote:<br></div></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr"><div class="gmail_default" style="font-family:georgia,serif">Hi,<br><br></div><div class="gmail_default" style="font-family:georgia,serif">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.<br><br></div><div class="gmail_default" style="font-family:georgia,serif">Here is my NCL code:<br><br>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_code.ncl"<br> load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_csm.ncl"<br> load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/wrf/WRFUserARW.ncl"<br>;-----------------------------<wbr>------------------------------<wbr>-----------<br> ; Main code<br> ;----------------------------<wbr>------------------------------<wbr>------------<br> begin<br> a = addfile("<a href="http://wrfout_d03_2017-07-17_06_00_00.nc" rel="noreferrer" target="_blank">wrfout_d03_2017-07-<wbr>17_06_00_00.nc</a>","r")<br><br>; type = "x11"<br>; type = "pdf"<br>; type = "ps"<br> type = "png"<br><br> wks = gsn_open_wks(type,"EQTemp_<wbr>Press_Cross_166")<br><br> time = 166<br><br> p = wrf_user_getvar(a, "pressure", time) ; Total Pressure<br> P = wrf_user_getvar(a, "P", time)<br> PB = wrf_user_getvar(a, "PB",time)<br> T = wrf_user_getvar(a,"T",time)<br> U = wrf_user_getvar(a, "ua", time)<br> V = wrf_user_getvar(a, "va", time)<br>; z = wrf_user_getvar(a,"z",time)<br> qv = wrf_user_getvar(a, "QVAPOR", time)<br><br><br> T_t = T + 290. ; potential temperature in K.<br> P_t = P + PB<br> tk = wrf_tk(P_t, T_t)<br> eth = wrf_eth(qv, tk, P_t)<br> plane = (/342, 1387/)<br> t_plane = wrf_user_intrp3d(eth,p,"v",<wbr>plane,50.,False)<br> P_plane = wrf_user_intrp3d(P,p,"v",<wbr>plane,50.,False)<br><br> printVarSummary(t_plane)<br><br> res = True<br> res@MainTitle = "Vertical Cross-Section of Potential Temperature"<br> res@InitTime = False<br> res@TimePos = "Left"<br> res@Footer = False<br> res@lbTitleOn = True<br> res@lbTitleString = "Temperature"<br> res@tmYLMode = "Explicit"<br> res@tmYLValues = (/1000, 850, 700, 500, 400, 300, 250,\<br> <wbr> 200, 150, 100, 50, 10/)<br> res@tmYLLabels = "" + res@tmYLValues<br> res@tmXBLabelFontHeightF = 0.015<br> res@tmYLLabelFontHeightF = 0.015<br><br><br> tres = res<br> tres@cnFillOn = True<br> tres@cnFillPalette = "BlWhRe"<br> tres@cnSpreadcolor = True<br><br><br> pres = res<br> pres@cnFillOn = False<br> pres@ContourParameters = (/ 2. /)<br> pres@cnInfoLabelFontHeightF = 0.01<br><br><br><br><br> T_contour = wrf_contour(a,wks,t_plane,<wbr>tres)<br> P_contour = wrf_contour(a,wks,P_plane,<wbr>pres)<br><br> pltres = True<br> pltres@gsnMaximize = True<br><br> plot = wrf_overlays(a,wks,(/T_<wbr>contour, P_contour/), pltres)<br><br></div><div class="gmail_default" style="font-family:georgia,serif">end<br><br></div><div class="gmail_default" style="font-family:georgia,serif">I'd really appreciate your help.<br><br></div><div class="gmail_default" style="font-family:georgia,serif">Bests,<br></div><div class="gmail_default" style="font-family:georgia,serif">Melissa <wbr> <wbr> <wbr> <br></div><div class="gmail_default" style="font-family:georgia,serif"><br><br></div></div></div></div>
______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" rel="noreferrer" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
</blockquote></div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>