;************************************************

;    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