James's memory constraint: wrf_user_getvar.

Hi Mary:

I solved his problem by this sample script:

Here you are:

  u4t = f4->U
  v4t = f4->V
  lat4 = f4->XLAT(:,:,:)
  lon4 = f4->XLONG(:,:,:)
  cen_long4 = f4_at_STAND_LON
cone = 0.0
printVarSummary(u4t)
uvmet4 = wrf_uvmet(u4t,v4t,lat4,lon4,cen_long4,cone)
printVarSummary(uvmet4)

Function wrf_uvmet gives the output below with more than one time slice and
no loop.

My uvmet4 looks like this:

Variable: uvmet4
Type: float
Total Size: 94832640 bytes
            23708160 values
Number of Dimensions: 5
Dimensions and sizes: [768] x [2] x [35] x [21] x [21]
Coordinates:
Number Of Attributes: 3
  _FillValue : -999
  units : m/s
  description : u,v met velocity

Please note lat4, lon4 need have time dimension as well.

Thanks,

Yi Wang, PhD
Pacific Northwest National Lab

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri May 23 2008 - 10:10:05 MDT