Re: average wind field

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Fri Jun 08 2012 - 07:15:58 MDT

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-talk
Received 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