[ncl-talk] Plotting vertical temperature from WRF

Dennis Shea shea at ucar.edu
Fri Mar 5 11:13:45 MST 2021


Sorry , I do not know. You will have to ask someone knowledgeable about WRF data. 

Sent from my iPhone

> On Mar 5, 2021, at 8:40 AM, ali mughal <mughalali655 at gmail.com> wrote:
> 
> Thank you for sharing this.
> 
> Is terrain height included in z and do I need to subtract the terrain height ?
> 
>> On Friday, March 5, 2021, Dennis Shea <shea at ucar.edu> wrote:
>> The simplest would be
>> 
>>     a = addfile("wrf...","r")
>>   
>>     tc  = wrf_user_getvar(a,"tc",-1)         ; T in C
>>     p   = wrf_user_getvar(a, "pressure",-1)  ; grid point pressure
>>     z   = wrf_user_getvar(a, "z",-1)         ; grid point height
>> 
>>     printVarSummary(tc)
>> 
>>     nt = 1
>>     nl = 50
>>     ml = 30
>>     TC = tc(nt,:,nl,ml)
>>     P  =  p(nt,:,nl,ml)
>> 
>>     wks   = gsn_open_wks ("png","WRF_VERT_PROF")       ; send graphics to PNG file
>> 
>>     res                   = True                       ; plot mods desired
>>     res at tiMainString      = " Sample Profile Plot"     ; add title
>>     res at trYReverse        = True                       ; reverse Y-axis
>>     plot  = gsn_csm_xy (wks,TC,P,res) ; create plot
>>    
>> 
>>> On Fri, Mar 5, 2021 at 3:47 AM ali mughal via ncl-talk <ncl-talk at mailman.ucar.edu> wrote:
>>> Dear all
>>> 
>>> Is there any example code available to plot vertical temperature profile from WRF ?
>>> 
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at mailman.ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210305/2f6df573/attachment.html>


More information about the ncl-talk mailing list