WRF post-processing questions should be sent to wrfhelp@ucar.edu
cc'd here
--- http://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml See Example 1 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl" a = addfile("wrfout_d01_2000-01-24_12:00:00.nc","r") ua = wrf_user_getvar(a,"ua",-1) ; u on mass points va = wrf_user_getvar(a,"va",-1) ; v on mass points printVarSummary(ua) wspd = ua ; 'trick' for meta data wspd = sqrt(ua^2 + va^2) ; (Time, bottom_top, south_north, west_east) wspd@description = "wind speed at mass points" printVarSummary(wspd) wspd_avg = dim_avg_n_Wrap(wspd, 0) ; (bottom_top,south_north,west_east) printVarSummary(wspd) On 6/8/12 2:12 AM, nlcheng wrote: > In wrf output files the U and V were defined below > float U(Time, bottom_top, south_north, west_east_stag) ; > > float V(Time, bottom_top, south_north_stag, west_east) ; > you can see that south_north_stag is not equal to south_north. I want to > caculate the average wind field in april,2011 > is there any example script on the NCL website? > thank you > > > > > _______________________________________________ > 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-talkReceived on Fri Jun 8 07:16:07 2012
This archive was generated by hypermail 2.1.8 : Tue Jun 12 2012 - 13:58:38 MDT