wrf_uvmet odd output values

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
Received on Tue Jan 27 2009 - 12:29:56 MST