diff -r 000000000000 -r 0c6405ab2ff4 time_series/10.write_ameriflux_casa.ncl
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/time_series/10.write_ameriflux_casa.ncl	Mon Jan 26 22:08:20 2009 -0500
@@ -0,0 +1,167 @@
+;************************************************
+;    Read ascii, Write nc
+;    output: lat: N->S     lon: -180W->180E                         
+;************************************************
+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
+
+; output
+
+  model_name = "i01.10casa"
+  
+  diro  = "/ptmp/jeff/"
+  filo  = model_name + "_ameriflux_1990-2004_monthly.nc"
+  c = addfile(diro+filo,"c")
+  filedimdef(c,"year",-1,True) 
+
+; input data
+
+  nlat = 64
+  nlon = 128
+
+  year_start = 1990
+  year_end   = 2004
+
+  nyear  = year_end - year_start + 1
+  nmonth = 12
+  
+  data1  = new((/nyear,nmonth,nlat,nlon/),float)
+  data2  = new((/nyear,nmonth,nlat,nlon/),float)
+  data3  = new((/nyear,nmonth,nlat,nlon/),float)
+  data4  = new((/nyear,nmonth,nlat,nlon/),float)
+  data5  = new((/nyear,nmonth,nlat,nlon/),float)
+  data6  = new((/nyear,nmonth,nlat,nlon/),float)
+  data7  = new((/nyear,nmonth,nlat,nlon/),float)
+  date   = new((/nyear,nmonth/),integer)
+
+;  name dimensions and assign coordinate variables
+
+  data1!0    ="year"
+  data1!1    ="month"
+  data1!2    = "lat"
+  data1!3    = "lon"
+
+  data2!0    ="year"
+  data2!1    ="month"
+  data2!2    = "lat"
+  data2!3    = "lon"
+
+  data3!0    ="year"
+  data3!1    ="month"
+  data3!2    = "lat"
+  data3!3    = "lon"
+
+  data4!0    ="year"
+  data4!1    ="month"
+  data4!2    = "lat"
+  data4!3    = "lon"
+
+  data5!0    ="year"
+  data5!1    ="month"
+  data5!2    = "lat"
+  data5!3    = "lon"
+
+  data6!0    ="year"
+  data6!1    ="month"
+  data6!2    = "lat"
+  data6!3    = "lon"
+
+  data7!0    ="year"
+  data7!1    ="month"
+  data7!2    = "lat"
+  data7!3    = "lon"
+
+  date!0    ="year"
+  date!1    ="month"
+  date@long_name  = "date: yyyymm"
+
+; input dir
+  diri  = "/ptmp/jeff/"+model_name+"/"
+
+  do m = 0,nyear-1
+     year = year_start + m
+  do n = 0,nmonth-1
+     nn = n+1
+     date(m,n) = year*100 +nn 
+     
+;    input file
+     if (nn .ge. 10) then
+        fili = model_name+".clm2.h0."+year+"-"+nn+".nc"
+     else
+        fili = model_name+".clm2.h0."+year+"-0"+nn+".nc"
+     end if
+
+     print (fili)
+
+     g = addfile (diri+fili,"r")
+
+     b = g->NEE
+     data1(m,n,:,:) = (/b(0,:,:)/)
+     data1@long_name  = b@long_name
+     data1@units      = b@units
+     data1@_FillValue = b@_FillValue
+     delete (b)
+
+     b = g->GPP
+     data2(m,n,:,:) = (/b(0,:,:)/)
+     data2@long_name  = b@long_name
+     data2@units      = b@units
+     data2@_FillValue = b@_FillValue
+     delete (b)
+
+     b = g->AR
+     data3(m,n,:,:) = (/b(0,:,:)/)
+     b = g->HR
+     data3(m,n,:,:) = data3(m,n,:,:)+(/b(0,:,:)/)
+     data3@long_name  = "total ecosystem respiration, autotrophic + heterotrophic"
+     data3@units      = b@units
+     data3@_FillValue = b@_FillValue
+     delete (b)
+
+     b = g->NETRAD
+     data4(m,n,:,:) = (/b(0,:,:)/)
+     data4@long_name  = b@long_name
+     data4@units      = b@units
+     data4@_FillValue = b@_FillValue
+     delete (b)
+
+     b = g->LATENT
+     data5(m,n,:,:) = (/b(0,:,:)/)
+     data5@long_name  = b@long_name
+     data5@units      = b@units
+     data5@_FillValue = b@_FillValue
+     delete (b)
+
+     b = g->FSH
+     data6(m,n,:,:) = (/b(0,:,:)/)
+     data6@long_name  = b@long_name
+     data6@units      = b@units
+     data6@_FillValue = b@_FillValue
+     delete (b)
+
+     b = g->FSDS
+     data7(m,n,:,:) = (/b(0,:,:)/)
+     data7@long_name  = b@long_name
+     data7@units      = b@units
+     data7@_FillValue = b@_FillValue
+     delete (b)
+  
+  end do
+  end do
+         
+  c->NEE = data1
+  c->GPP = data2
+  c->ER  = data3
+  c->NETRAD = data4
+  c->LATENT = data5
+  c->FSH    = data6
+  c->FSDS    = data7
+  c->date   = date
+  c->area   = g->area
+
+  print (date)
+
+end