Good afternoon,
I am trying to ingest a test.nc file (lat=59.125 to 59.375; lon=150.125
to 299.875):
netcdf test {
dimensions:
lat = 2 ;
lon = 600 ;
bin = 55 ;
variables:
short bin_cnt(lat, lon, bin) ;
bin_cnt:missing_value = -999s ;
float lat(lat) ;
lat:units = "degrees_north" ;
lat:long_name = "latitude" ;
lat:_FillValue = -999.f ;
float lon(lon) ;
lon:units = "degrees_east" ;
lon:long_name = "longitude" ;
lon:_FillValue = -999.f ;
I would like to plot a histogram containing the sum of the counts per
bin for a portion of the globe using the NCL code below:
begin
data = addfile(test.nc,"r")
bin_cnt = data->bin_cnt
bin_cnt@_FillValue = floattoshort(-999)
glb_hist = dim_sum_Wrap(bin_cnt(bin|:,lon|:,lat|:))
glb_hist1 = dim_sum_Wrap(glb_hist(bin|:,lon|:)) ;Contains a bin count
for the entire globe
end
For some strange reason, the glb_hist1 array contains negative values
when the input variable (bin_cnt) does not. The bin_cnt variable
contains positive (short format) numbers. I can provide the test.nc
file on request.
-- Scott Capps University of California, Irvine Department of Earth System Science Croul Hall Irvine, CA 92617
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Sep 25 2007 - 13:38:35 MDT
This archive was generated by hypermail 2.2.0 : Mon Oct 01 2007 - 07:42:01 MDT