Re: xy-plot using wrf output

Wei,

The min/max that you printed out is from a different variable (rain_tot) than what you are actually plotting (rain).

Please examine all your rain*** variables and make sure the missing values (if any) are what you expect.

If this doesn't reveal anything, then can you provide me with the wrfout_d01.nc file?

--Mary

On Nov 27, 2010, at 8:41 AM, wei huang wrote:

> 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
> <xyplot.pdf>_______________________________________________
> 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 Mon Nov 29 15:36:39 2010