Time series of WRF rainfall - values

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 Tue Jan 20 2009 - 19:02:28 MST