1.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
     1.2 +++ b/biomass/01.read.byte.ncl	Mon Jan 26 22:08:20 2009 -0500
     1.3 @@ -0,0 +1,71 @@
     1.4 +;************************************************
     1.5 +;    Read byte, Write nc
     1.6 +; input data is  :   13.8533N -> -21.1277S
     1.7 +;                   -82.7209 -> -33.5739                       
     1.8 +; output data is :  -21.1277S -> 13.8533N
     1.9 +;                    277.2791 -> 326.4261   
    1.10 +;************************************************
    1.11 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"  
    1.12 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"  
    1.13 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"   
    1.14 +;************************************************
    1.15 +begin
    1.16 +; final output
    1.17 +  diro = "/fis/cgd/cseg/people/jeff/clamp_data/biomass/ob/"
    1.18 +  filo =  "amazon_biomass.nc"
    1.19 +  c = addfile(diro+filo,"c")
    1.20 +
    1.21 +; input byte file
    1.22 +  diri = "/fis/cgd/cseg/people/jeff/clamp_data/biomass/ob/"
    1.23 +  fili =  "amazon_biomass_byte.nc"
    1.24 +  b = addfile(diri+fili,"r")
    1.25 +
    1.26 +  data = byte2flt(b->biomass)
    1.27 +
    1.28 +  dsizes_data = dimsizes(data)
    1.29 +  nlat = dsizes_data(0)
    1.30 +  nlon = dsizes_data(1)
    1.31 +; print (nlat)
    1.32 +; print (nlon)
    1.33 + 
    1.34 +  lat  = new((/nlat/),float)
    1.35 +  lon  = new((/nlon/),float)
    1.36 +
    1.37 +  latS = -21.1277
    1.38 +  latN =  13.8533
    1.39 +  lonL = -82.7209
    1.40 +  lonR = -33.5739
    1.41 +
    1.42 +  dy = abs((latN-latS)/(nlat-1))
    1.43 +  dx = abs((lonR-lonL)/(nlon-1))
    1.44 +
    1.45 +  do n = 0,nlat-1
    1.46 +     lat(n) = latS + n*dy
    1.47 +  end do
    1.48 +
    1.49 +  do n = 0,nlon-1
    1.50 +     lon(n) = lonL + n*dx + 360.
    1.51 +  end do 
    1.52 +                 
    1.53 +  lat!0  = "lat"
    1.54 +  lon!0  = "lon"
    1.55 +  lat@units      = "degrees_north"
    1.56 +  lat@long_name  = "Latitude"
    1.57 +  lon@units      = "degrees_east"
    1.58 +  lon@long_name  = "Longitude"
    1.59 +  print (lat)
    1.60 +  print (lon)
    1.61 +          
    1.62 +  data!0  = "lat"
    1.63 +  data!1  = "lon"
    1.64 +
    1.65 +  data@units      = "Mg C/ha"
    1.66 +  data@long_name  = "Amoazon Biomass"
    1.67 +  data@_FillValue = 1.e+36
    1.68 +
    1.69 +; print (data)
    1.70 +  
    1.71 +  c->lat       = lat
    1.72 +  c->lon       = lon
    1.73 +  c->BIOMASS   = data(::-1,:)  
    1.74 +end