diff -r 000000000000 -r 0c6405ab2ff4 lai/02.read_byte_LAI.ncl.x
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/lai/02.read_byte_LAI.ncl.x	Mon Jan 26 22:08:20 2009 -0500
@@ -0,0 +1,54 @@
+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
+  nlat  = 3600
+  mlon  = 7200
+
+  diri  = "/fis/cgd/cseg/people/jeff/clamp_data/"
+  fili  = "MOD15A2_GEO_2000001.nc"
+  a  = addfile(diri+fili,"r")
+
+  y = a->Lai_0_05deg    ; byte
+  y@_FillValue = inttobyte(253)
+  y@_FillValue = inttobyte(255)
+
+  x = byte2flt(y)
+  
+  x@_FillValue =  1.e+36                     ; new fill
+
+  lat  = latGlobeFo(nlat, "lat", "latitude", "degrees_north")
+  lat  = (/ lat(::-1) /)                     ; make N->S
+  lon  = lonGlobeFo(mlon, "lon", "longitude", "degrees_east")
+  lon  = (/ lon - 180. /) ; subtract 180 from all values 
+  lon&lon = lon           ; update coordinates
+
+  x!0  = "lat"
+  x!1  = "lon"
+  x&lat=  lat
+  x&lon=  lon
+
+  y!0  = "lat"
+  y!1  = "lon"
+  y&lat=  lat
+  y&lon=  lon
+
+  iy   = new (dimsizes(y), "integer", -999)
+  iy   = y
+
+  iy!0  = "lat"
+  iy!1  = "lon"
+  iy&lat=  lat
+  iy&lon=  lon
+
+  print(y(:,{130})+"   "+ iy(:,{130})+"   "+x(:,{130}))
+
+
+
+; diro  = "/fis/cgd/cseg/people/jeff/clamp_data/"
+; filo  = "LAI_2000001_2.nc"
+; c = addfile(diro+filo,"c")
+; c->LAI  = x 
+end
+