xy-plot using wrf output

From: wei huang <whua27_at_nyahnyahspammersnyahnyah>
Date: Sat Nov 27 2010 - 08:41:39 MST

Dear All,

When I plotted the xy-plot using the WRF output data, the figure is look
like different to me. Because when I got it print the total precipitation:
min=-7.09444e-08 max=230.888 and the figure output is different. However,
for your convenience the script and the figure are also attached herewith.
Could someone help me to fix the problem?
Thanks
wei

**************************************

load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"

   begin

   a = addfile("wrfout_d01.nc","r")

   ; Calculate total precipitation

    rain_exp = wrf_user_getvar(a,"RAINNC",-1)

    rain_con = wrf_user_getvar(a,"RAINC",-1)

    rain_tot = rain_exp + rain_con

    rain_tot@description = "Total Precipitation"

    rain_tot@units =rain_con@units

    printVarSummary(rain_tot)

    printMinMax(rain_tot, True)

    rain = wgt_areaave_Wrap(rain_tot,1.0,1.0,0)

    printVarSummary(rain)

 ;*****************************

  wks = gsn_open_wks("x11","xy") ; open a workstation

  res = True

;***********;

;************

  plot = gsn_csm_xy(wks,taus,rain, res)

  end

*******************************

Print:

Variable: rain_tot

Type: float

Total Size: 879660 bytes

            219915 values

Number of Dimensions: 3

Dimensions and sizes: [9] x [135] x [181]

Coordinates:

Number Of Attributes: 2

  units : mm

  description : Total Precipitation

(0)

(0) Total Precipitation: min=-7.09444e-08 max=230.888

 Variable: rain

Type: float

Total Size: 36 bytes

            9 values

Number of Dimensions: 1

Dimensions and sizes: [9]

Coordinates:

Number Of Attributes: 4

  _FillValue : -999

  units : mm

  description : Total Precipitation

  wgt_areaave_op_ncl : Area Average

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

Received on Sat Nov 27 08:41:46 2010

This archive was generated by hypermail 2.1.8 : Tue Nov 30 2010 - 08:26:08 MST