diff -r 000000000000 -r 0c6405ab2ff4 biomass/01.read.byte.ncl
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/biomass/01.read.byte.ncl	Mon Jan 26 22:08:20 2009 -0500
@@ -0,0 +1,71 @@
+;************************************************
+;    Read byte, Write nc
+; input data is  :   13.8533N -> -21.1277S
+;                   -82.7209 -> -33.5739                       
+; output data is :  -21.1277S -> 13.8533N
+;                    277.2791 -> 326.4261   
+;************************************************
+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"   
+;************************************************
+begin
+; final output
+  diro = "/fis/cgd/cseg/people/jeff/clamp_data/biomass/ob/"
+  filo =  "amazon_biomass.nc"
+  c = addfile(diro+filo,"c")
+
+; input byte file
+  diri = "/fis/cgd/cseg/people/jeff/clamp_data/biomass/ob/"
+  fili =  "amazon_biomass_byte.nc"
+  b = addfile(diri+fili,"r")
+
+  data = byte2flt(b->biomass)
+
+  dsizes_data = dimsizes(data)
+  nlat = dsizes_data(0)
+  nlon = dsizes_data(1)
+; print (nlat)
+; print (nlon)
+ 
+  lat  = new((/nlat/),float)
+  lon  = new((/nlon/),float)
+
+  latS = -21.1277
+  latN =  13.8533
+  lonL = -82.7209
+  lonR = -33.5739
+
+  dy = abs((latN-latS)/(nlat-1))
+  dx = abs((lonR-lonL)/(nlon-1))
+
+  do n = 0,nlat-1
+     lat(n) = latS + n*dy
+  end do
+
+  do n = 0,nlon-1
+     lon(n) = lonL + n*dx + 360.
+  end do 
+                 
+  lat!0  = "lat"
+  lon!0  = "lon"
+  lat@units      = "degrees_north"
+  lat@long_name  = "Latitude"
+  lon@units      = "degrees_east"
+  lon@long_name  = "Longitude"
+  print (lat)
+  print (lon)
+          
+  data!0  = "lat"
+  data!1  = "lon"
+
+  data@units      = "Mg C/ha"
+  data@long_name  = "Amoazon Biomass"
+  data@_FillValue = 1.e+36
+
+; print (data)
+  
+  c->lat       = lat
+  c->lon       = lon
+  c->BIOMASS   = data(::-1,:)  
+end