1 ;************************************************
3 ; input data is : 13.8533N -> -21.1277S
5 ; output data is : -21.1277S -> 13.8533N
7 ;************************************************
8 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
9 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
10 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
11 ;************************************************
14 diro = "/fis/cgd/cseg/people/jeff/clamp_data/biomass/ob/"
15 filo = "amazon_biomass.nc"
16 c = addfile(diro+filo,"c")
19 diri = "/fis/cgd/cseg/people/jeff/clamp_data/biomass/ob/"
20 fili = "amazon_biomass_byte.nc"
21 b = addfile(diri+fili,"r")
23 data = byte2flt(b->biomass)
25 dsizes_data = dimsizes(data)
31 lat = new((/nlat/),float)
32 lon = new((/nlon/),float)
39 dy = abs((latN-latS)/(nlat-1))
40 dx = abs((lonR-lonL)/(nlon-1))
47 lon(n) = lonL + n*dx + 360.
52 lat@units = "degrees_north"
53 lat@long_name = "Latitude"
54 lon@units = "degrees_east"
55 lon@long_name = "Longitude"
62 data@units = "Mg C/ha"
63 data@long_name = "Amoazon Biomass"
64 data@_FillValue = 1.e+36
70 c->BIOMASS = data(::-1,:)