diff -r 000000000000 -r 0c6405ab2ff4 class/02.class.ncl
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/class/02.class.ncl	Mon Jan 26 22:08:20 2009 -0500
@@ -0,0 +1,72 @@
+;*************************************************
+; ce_1.ncl
+;************************************************
+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
+ 
+;************************************************
+; read in model data
+;************************************************
+  diri  = "/fis/cgd/cseg/people/jeff/clamp/class/"
+  fili  = "T42.clamp.surface-data.060412.nc"
+  f     = addfile(diri+fili,"r")
+
+  y     = f->PCT_PFT
+  printVarSummary(y)
+
+  diri  = "/fis/cgd/cseg/people/jeff/surface_data/"
+  fili  = "lnd_T42.nc"   
+  g = addfile(diri+fili,"r")
+
+  landmask = g->landmask
+  lat      = g->lat
+  lon      = g->lon
+
+  nlat = dimsizes(lat)
+  nlon = dimsizes(lon)
+
+  x = y(0,:,:)
+
+  x!0   = "lat"
+  x&lat = lat            
+  x!1   = "lon"
+  x&lon = lon
+  x@_FillValue = 1.e36
+  x@long_name  = "Model PFT Classes"  
+  
+  do j= 0,nlat-1
+  do i= 0,nlon-1
+     x(j,i) = maxind(y(:,j,i))      
+  end do
+  end do 
+
+; print (x)  
+
+  x = where(landmask .lt. 1.,x@_FillValue,x)
+
+;************************************************
+; create default plot
+;************************************************
+
+  wks = gsn_open_wks("ps","class")              ; open a ps file
+  gsn_define_colormap(wks,"gui_default")     ; choose colormap
+
+  res                     = True             ; Use plot options
+  res@cnFillOn            = True             ; Turn on color fill
+  res@gsnSpreadColors     = True             ; use full colormap
+  res@cnLinesOn           = False            ; Turn off contourn lines
+  res@mpFillOn            = False            ; Turn off map fill
+  res@tiMainString        = "Class"
+
+  res@gsnSpreadColors      = True            ; use full colormap
+  res@cnLevelSelectionMode = "ManualLevels"  ; Manual contour invtervals
+  res@cnMinLevelValF       = 0.              ; Min level
+  res@cnMaxLevelValF       = 16.             ; Max level
+  res@cnLevelSpacingF      = 1.              ; interval
+
+  plot = gsn_csm_contour_map_ce(wks,x,res)   
+
+end
\ No newline at end of file