Re: FillValue question :TRMM

From: Adam Phillips <asphilli_at_nyahnyahspammersnyahnyah>
Date: Thu Oct 15 2009 - 17:41:12 MDT

Hi Bold,
Your printVarSummary output shows that your prate_join array is a short.
So when you go to assign an integer _FillValue, you are getting an error
message. Change this:
prate_join = fall[:]->PRC ; retrieve the PRC from many TRMM.nc
prate_join@_FillValue = -9999
to this:
prate_join = short2flt(fall[:]->PRC) ; retrieve the PRC from many TRMM.nc
prate_join@_FillValue = -9999.
Good luck,
Adam

B N wrote:
> Hi. I used the script from the NCL example site to extract
> precipitation rate from TRMM files.
>
> Unfortunately, the resulting fillValues (_FillValue : 32767) are
> messing up all of my plots.
> How do I change this? I tried "prate_join@_FillValue = -9999" and I
> got an error.
> Thank you for the help.
> Sincerely,
> Bold
>
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
> fils = systemfunc("ls
> /Volumes/Data-Store_Terabyte2/DATA/Precipitation_Observations/TRMM3B42/3B*z.nc")
> fall = addfiles (fils, "r")
> prate_join = fall[:]->PRC ; retrieve the PRC from many TRMM.nc
> files (here, 240 files that = 1 month)
> printVarSummary(prate_join)
> prate_join@_FillValue = -9999
>
> $ ncl PRATE_TRMM_plot_daily_Avg_05deg.ncl
> Copyright (C) 1995-2009 - All Rights Reserved
> University Corporation for Atmospheric Research
> NCAR Command Language Version 5.1.1
> The use of this software is governed by a License Agreement.
> See http://www.ncl.ucar.edu/ for more details.
>
> Variable: prate_join
> Type: short
> Total Size: 276480000 bytes
> 138240000 values
> Number of Dimensions: 3
> Dimensions and sizes: [time | 240] x [lat | 400] x [lon | 1440]
> Coordinates:
> time: [84720..85437]
> lat: [-49.875..49.875]
> lon: [-179.875..179.875]
> Number Of Attributes: 9
> _FillValue : 32767
> long_name : TRMM_3B42 3-hourly average rainfall rate
> units : mm/hr
> hdf_name : precipitation
> add_offset : 0
> scale_factor : 0.001
> vMin_original_data : 0
> vMax_original_data : 49.88
> vRange : 49.88
>
>
> fatal:Type Mismatch: The type of missing value could not be converted
> to type of variable (prate_join)
> fatal:Execute: Error occurred at or near line 9 in file
> PRATE_TRMM_plot_daily_Avg_05deg.ncl
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
__________________________________________________
Adam Phillips 
asphilli@ucar.edu
National Center for Atmospheric Research   tel: (303) 497-1726
Climate and Global Dynamics Division         fax: (303) 497-1333
P.O. Box 3000				
Boulder, CO 80307-3000    http://www.cgd.ucar.edu/cas/asphilli
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu, 15 Oct 2009 17:41:12 -0600

This archive was generated by hypermail 2.1.8 : Thu Oct 22 2009 - 12:22:52 MDT