Re: wrf_uvmet odd output values

Hi Melissa,

The wrf_uvmet function has undergone a bit of an overhaul
for the V5.1.0 release coming soon, but don't think there
were any changes that might address the problem you mention.

You are right about the variable being returned is not the correct
advertised size. I think this was due to staggered versus unstaggered
grids.

I've CC'ed wrfhelp_at_ucar.edu to see if they have any ideas on this.

You have a rather old version of NCL, by the way. The current version
is 5.0.0 which was released over a year ago. Are you able to try
your script using V5.0.0?

--Mary

On Tue, 27 Jan 2009, Melissa Bukovsky wrote:

> Hi,
>
> As wrf_user_getvar can't handle arrays of files, I've been trying to rework
> parts of it to create the same variables but for many files. However, I'm
> receiving some unrealistic values from the part using wrf_uvmet. At least
> they seem that way (~ -150 mi/hr converted to -1500 mi/hr winds seems pretty
> extreme). The minimum U and V values that come out of wrf_uvmet are
> significantly lower that those that went in. I was wondering if anyone had
> the code for this function so that I could try to figure out why, or if
> someone knows specifically why this is. It doesn't seem like rotating them
> should create that big of a difference (even at the farthest extents of my
> domain). There is no problem when I use a smaller subset of my domain over
> the US; however, I'd like use the entire domain. Also note that the
> dimensions of the output from wrf_uvmet are not what the function description
> says they should be, so I am wondering if I am feeding it variables
> improperly and I just don't see it.
> The relevant part of my code is pasted below. The print output from it
> (including lat/lon to provide an idea of my domain extent) is also included.
> I am using ncl version 4.3.1 (with the updated WRFUserARW.ncl script).
> Thanks in advance for any help!
>
> ~Melissa
>
>
>
> ;YEAR LOOP
>
> nyr = dimsizes(year)
>
> do l = 0,nyr-1
>
> dir = rtdir+runname+"/"+year(l)+tag+"/"
>
> ;MONTH LOOP (since ncl can't handle over 4 months of files at once)
>
> nmon = dimsizes(mons)
>
> do n = 0,nmon-1
>
> smon = sprinti("%0.2i",mons(n))
>
> fils = systemfunc("ls "+dir+"wrfout_"+dom+"_"+year(l)+"-"+smon+"*")
> f = addfiles(fils+".nc","r")
> ListSetType (f,"cat") ;concatenate
>
> ;mass point lat and lon
>
> latitude = f[:]->XLAT longitude = f[:]->XLONG
> u = f[:]->U
> v = f[:]->V
> ua = uv_to_massGrid(u,"U") ;must send wrf_uvmet unstaggered grid,
> not working with staggered!
> va = uv_to_massGrid(v,"V")
>
> cone = wrf_calc_cone(f[0])
>
> if(isatt(f[0],"STAND_LON"))then
> cenlon = f[0]@STAND_LON
> else
> cenlon = f[0]@CEN_LON
> end if
>
> uvmet = wrf_uvmet(ua,va,latitude,longitude,cenlon,cone)
>
> printVarInfo(latitude,"lat")
> printVarInfo(longitude,"lon")
> printVarInfo(ua,"ua")
> printVarInfo(va,"va")
> printVarInfo(uvmet,"uvmet")
>
> ua = uvmet(:,0,:,:,:) ;output from wrf_uvmet is (Time, 2,
> bottom_top, south_north, west_east) (WHY!)
> va = uvmet(:,1,:,:,:) ;not (2, Time, bottom_top, south_north,
> west_east) as function description states!!
>
> printVarInfo(ua,"ua")
> printVarInfo(va,"va")
>
> PRINT OUTPUT...
>
>
>
> (0) ==========> printVarInfo: lat
>
>
> Variable: x (parameter)
> Type: float
> Total Size: 4097952 bytes
> 1024488 values
> Number of Dimensions: 3
> Dimensions and sizes: [248] x [51] x [81]
> Coordinates: (0) Minimum: 11.4338 Maximum: 55.8994
> (0) (0) (0) ==========> printVarInfo: lon
>
>
> Variable: x (parameter)
> Type: float
> Total Size: 4097952 bytes
> 1024488 values
> Number of Dimensions: 3
> Dimensions and sizes: [248] x [51] x [81]
> Coordinates: (0) Minimum: -138.778 Maximum: -40.2224
> (0) (0) (0) ==========> printVarInfo: ua
>
>
> Variable: x (parameter)
> Type: float
> Total Size: 110644704 bytes
> 27661176 values
> Number of Dimensions: 4
> Dimensions and sizes: [Time | 248] x [bottom_top | 27] x [south_north | 51]
> x [west_east | 81]
> Coordinates: Number Of Attributes: 2
> description : U Velocity
> units : m/s
> (0) Minimum: *-39.8189* Maximum: 83.9085
> (0) (0) (0) ==========> printVarInfo: va
>
>
> Variable: x (parameter)
> Type: float
> Total Size: 110644704 bytes
> 27661176 values
> Number of Dimensions: 4
> Dimensions and sizes: [Time | 248] x [bottom_top | 27] x [south_north | 51]
> x [west_east | 81]
> Coordinates: Number Of Attributes: 2
> description : V Velocity
> units : m/s
> (0) Minimum: *-64.4031* Maximum: 80.8059
> (0) (0) (0) ==========> printVarInfo: uvmet
>
>
> Variable: x (parameter)
> Type: float
> Total Size: 221289408 bytes
> 55322352 values
> Number of Dimensions: 5
> Dimensions and sizes: *[248] x [2] x [27] x [51] x [81]*
> Coordinates: Number Of Attributes: 3
> _FillValue : -999
> units : m/s
> description : u,v met velocity
> (0) Minimum: -651.28 Maximum: 84.4652
> (0) (0) ==========> printVarInfo: ua
>
>
> Variable: x (parameter)
> Type: float
> Total Size: 110644704 bytes
> 27661176 values
> Number of Dimensions: 4
> Dimensions and sizes: [248] x [27] x [51] x [81]
> Coordinates: Number Of Attributes: 3
> description : u,v met velocity
> units : m/s
> _FillValue : -999
> (0) Minimum: *-350.442* Maximum: 84.4652
> (0) (0) (0) ==========> printVarInfo: va
>
>
> Variable: x (parameter)
> Type: float
> Total Size: 110644704 bytes
> 27661176 values
> Number of Dimensions: 4
> Dimensions and sizes: [248] x [27] x [51] x [81]
> Coordinates: Number Of Attributes: 3
> description : u,v met velocity
> units : m/s
> _FillValue : -999
> (0) Minimum: *-651.28* Maximum: 64.0054
> (0)
>
>
> _______________________________________________
> 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 Tue Jan 27 2009 - 17:33:04 MST