Re: how to modify the code for using wrf output

From: Mary Haley <haley_at_nyahnyahspammersnyahnyah>
Date: Wed Aug 03 2011 - 15:17:08 MDT

Further info: you can use wrf_user_ll_to_ij to get the ij indexes of the lat/lon values you want.

http://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_ll_to_ij.shtml

--Mary

On Aug 3, 2011, at 12:05 PM, Dennis Shea wrote:

> You already have the vertical velocity (w, not omega).
> At each grid point, you have or can calculate the pressure levels
> (bottom_top_stag) ... call it p(Time,bottom_top_stag,south_north,west_east)
>
> You can not use the {...} syntax because the grid is curvilinear and not
> rectilinear, You must use standard indexing.
>
> t = 5
> n = (/32, 73/) ; lat index values of the points you want
> m = (/52, 105/) ; lon
>
> x = w(t,:,n,m)
> p = p(t,:,n,m)
>
> plt = gsn_csm_xy (wks,x,p,resxy)
>
>
> On 8/3/11 11:06 AM, wei huang wrote:
>> Hi, could anyone give me some help to modifying the code (omega_1.ncl
>> <http://www.ncl.ucar.edu/Applications/Scripts/omega_1.ncl>) for using
>> the WRF output. I’m interested to plot the vertical profile like this
>> omega plot. The wrf variable print is as follows:
>> thanks and regards
>>
>> Variable: w
>>
>> Type: float
>>
>> Total Size: 130572000 bytes
>>
>> 32643000 values
>>
>> Number of Dimensions: 4
>>
>> Dimensions and sizes:[Time | 48] x [bottom_top_stag | 27] x [south_north
>> | 100] x [west_east | 80]
>>
>> Coordinates:
>>
>> Number Of Attributes: 6
>>
>> coordinates : XLONG XLAT
>>
>> stagger :Z
>>
>> units :m s-1
>>
>> description : z-wind component
>>
>> MemoryOrder : XYZ
>>
>> FieldType :104
>>
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Aug 3 15:17:15 2011

This archive was generated by hypermail 2.1.8 : Fri Aug 05 2011 - 14:53:55 MDT