Re: Time series of WRF rainfall - values

Erik,

Can you put the WRF file somewhere that I can download it? If not, email me
directly and we'll figure something out.

I'm wondering if there might be a missing value issue.

--Mary

On Tue, 20 Jan 2009, Erik Noble wrote:

> Dear ncl,
> To plot a time series of rainfall at pojnt, I am trying to follow the ncl
> WRF xy plot example at
> http://www.ncl.ucar.edu/Applications/Scripts/WRF_xy_1.ncl
>
> But I¹m using it for rainfall . Yet, when I print out of the rainfall
> variable at a point (5.5 N, 25 W), I get exponential answers even though the
> units are in millimeters; the plot looks unrealistic when I try to plot a
> time series.
> Is there another way to plot WRF rainfall as a timeseries?
> Sincerely,
> Erik Noble
>
>
> 107 squall:noble% ncl
> Copyright (C) 1995-2007 - All Rights Reserved
> University Corporation for Atmospheric Research
> NCAR Command Language Version 5.0.0
> The use of this software is governed by a License Agreement.
> See http://www.ncl.ucar.edu/ for more details.
> ncl 0> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> ncl 1> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> ncl 2> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> ncl 3> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
> ncl 4> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> ncl 5> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> ncl 6> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
> ncl 7> filename =
> "WRF_Output_gfs_3_20060901/wrfout_d01_2006-09-01_00:00:00.nc"
> ncl 8> f=addfile(filename,"r")
> ncl 9> ;************************************************
> ncl 10> ; Read all the times. Convert to numeric units for plotting.
> ncl 11> ;************************************************
> ncl 12> Times = f->Times ; Times is of type
> character
> ncl 13> Time = WRF_Times2Udunits_c(Times, 0) ; convert to "hours since"
> ncl 14> nc = f->RAINNC
> ncl 15> cc = f->RAINC
> ncl 16> rainfall = cc+nc
> ncl 17> copy_VarMeta(cc,rainfall); Copies all attributes, named dimensions
> and coordinate
> ncl 18> rainfall_at_description ="Total Accumulated Precipitation"
> ncl 19> loc = wrf_user_latlon_to_ij(f,5.5,-25)
> ncl 20> print("X/Y location is: " +loc)
> (0) X/Y location is: 56
> (1) X/Y location is: 20
> ncl 21> y = rainfall(:,56,20)
> ncl 22> printVarSummary(y)
> Variable: y
> Type: float
> Total Size: 68 bytes
> 17 values
> Number of Dimensions: 1
> Dimensions and sizes: [Time | 17]
> Coordinates:
> Number Of Attributes: 6
> coordinates : XLONG XLAT
> stagger :
> units : mm
> description : Total Accumulated Precipitation
> MemoryOrder : XY
> FieldType : 104
> ncl 23> print(max(y))
> (0) 2.187776e-10
> ncl 24> print(min(y))
> (0) 0
> ncl 25> print(y)
>
>
> Variable: y
> Type: float
> Total Size: 68 bytes
> 17 values
> Number of Dimensions: 1
> Dimensions and sizes: [Time | 17]
> Coordinates:
> Number Of Attributes: 6
> coordinates : XLONG XLAT
> stagger :
> units : mm
> description : Total Accumulated Precipitation
> MemoryOrder : XY
> FieldType : 104
> (0) 0
> (1) 0
> (2) 0
> (3) 0
> (4) 1.831067e-20
> (5) 5.492178e-14
> (6) 5.492178e-14
> (7) 5.492178e-14
> (8) 5.492178e-14
> (9) 5.492178e-14
> (10) 5.492178e-14
> (11) 2.187776e-10
> (12) 2.187776e-10
> (13) 2.187776e-10
> (14) 2.187776e-10
> (15) 2.187776e-10
> (16) 2.187776e-10
> (END)
>
> ncl 26> exit
>
> Why the high values? How can I just plot a simple time series?
> Thank you for your help.
> Sincerely,
> Erik Noble
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Dec 22 2008 - 16:42:45 MST