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
This archive was generated by hypermail 2.1.8 : Tue Nov 30 2010 - 08:26:08 MST