FillValue question :TRMM

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
Received on Thu, 15 Oct 2009 18:58:37 -0400