diff -r 000000000000 -r 0c6405ab2ff4 npp/02.read_binary_0.05deg.ncl
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/npp/02.read_binary_0.05deg.ncl	Mon Jan 26 22:08:20 2009 -0500
@@ -0,0 +1,79 @@
+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/"
+  fili  = "Npp_0.05deg_mean.int16"
+
+  diro  = "/fis/cgd/cseg/people/jeff/clamp/"
+  filo  = "Npp_0.05deg_mean_3.nc"
+  c = addfile(diro+filo,"c")
+
+; I think  (a) the file is "little endian"; read into short
+;          (b) the _FillValue is 32700
+
+  setfileoption("bin","ReadByteOrder","LittleEndian")
+  xShort= fbindirread(diri+fili,0, (/nlat,mlon/), "short")
+  xShort@_FillValue= inttoshort(32700)
+
+;;xShort@scale_factor = ?????
+;;xShort@add_offset   = ?????
+
+  x     = short2flt( xShort )
+
+;;print(xShort(:,1800)+"   "+x(:,1800))      ; look at the values
+  delete(xShort)
+
+  x@long_name  = "net primary production"
+  x@units      = "gC/m^2/year"
+; x@_FillValue = -17281                      ; orig fill    
+  x@_FillValue =  1e20                       ; 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
+
+  c->NPP  = x
+
+  exit
+  
+;************************************************
+; create default plot
+;************************************************
+  
+  setvalues NhlGetWorkspaceObjectId()
+    "wsMaximumSize" : 199999999
+  end setvalues
+
+  wks = gsn_open_wks("ps","Npp_rdBinPlt")                  ; open a ps file
+ ;gsn_define_colormap(wks,"wgne15")          ; choose colormap
+
+  res                     = True             ; Use plot options
+  res@cnFillOn            = True             ; Turn on color fill
+  res@cnFillMode          = "RasterFill"     ; Turn on raster color
+  res@lbLabelAutoStride   = True
+  res@cnLinesOn           = False            ; Turn off contourn lines
+ ;res@gsnSpreadColors     = True             ; use full colormap
+  res@mpFillOn            = False            ; Turn off map fill
+
+ ;res@cnLevelSelectionMode = "ManualLevels"    ; Manual contour invtervals
+ ;res@cnMinLevelValF       =                   ; Min level
+ ;res@cnMaxLevelValF       =                   ; Max level
+ ;res@cnLevelSpacingF      =                   ; interval
+  res@tiMainString         = fili
+
+  plot = gsn_csm_contour_map_ce(wks,x,res)    
+
+end
+