1.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
     1.2 +++ b/lai/04.lai_ensemble_T42.ncl.x	Mon Jan 26 22:08:20 2009 -0500
     1.3 @@ -0,0 +1,82 @@
     1.4 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"  
     1.5 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"  
     1.6 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"  
     1.7 +
     1.8 +begin 
     1.9 +;----------------------------------------------------------
    1.10 +  year_start = 2000
    1.11 +  year_end   = 2005
    1.12 +;----------------------------------------------------------
    1.13 +  diri  = "/fis/cgd/cseg/people/jeff/clamp_data/"
    1.14 +  diro  = "/fis/cgd/cseg/people/jeff/clamp_data/"
    1.15 +  filo1 = "MOD15A2_LAI_"+year_start+"-"+year_end+"_ensemble.nc"
    1.16 +  filo2 = "MOD15A2_LAI_"+year_start+"-"+year_end+"_mean.nc"
    1.17 +  c = addfile(diro+filo1,"c")
    1.18 +  d = addfile(diro+filo2,"c")
    1.19 +
    1.20 +  nyear = year_end - year_start + 1
    1.21 +
    1.22 +  files  = new((/nyear/),"string")
    1.23 +
    1.24 +  do n = year_start,year_end
    1.25 +     i = n - year_start
    1.26 +;    files(i)= diri + "MOD15A2_LAI_" + n + "_monthly.nc"
    1.27 +     files(i)= diri + "MOD15A2_LAI_" + n + "_monthly_T42.nc"
    1.28 +  end do
    1.29 +  print (files)  
    1.30 +
    1.31 +  file_list = addfiles(files,"r")
    1.32 +  ListSetType(file_list,"join")
    1.33 +  x = addfiles_GetVar(file_list, files, "LAI")
    1.34 +  printVarSummary (x)
    1.35 +
    1.36 +  y = addfiles_GetVar(file_list, files, "time")
    1.37 +  time = y(0,:)
    1.38 +  delete(time@case)
    1.39 +  delete(y)
    1.40 +  printVarSummary (time)
    1.41 +
    1.42 +  y = addfiles_GetVar(file_list, files, "lat")
    1.43 +  lat = y(0,:)
    1.44 +  delete(lat@case)
    1.45 +  delete(y)
    1.46 +  printVarSummary (lat)
    1.47 +
    1.48 +  y = addfiles_GetVar(file_list, files, "lon")
    1.49 +  lon = y(0,:)
    1.50 +  delete(lon@case)
    1.51 +  delete(y)
    1.52 +  printVarSummary (lon)
    1.53 +
    1.54 +  x!0  = "case"
    1.55 +  x!1  = "time"
    1.56 +  x!2  = "lat"
    1.57 +  x!3  = "lon"
    1.58 +
    1.59 +  y = dim_avg(x(time|:,lat|:,lon|:,case|:))
    1.60 +
    1.61 +  y!0  = "time"
    1.62 +  y!1  = "lat"
    1.63 +  y!2  = "lon"
    1.64 +  y&time = time
    1.65 +  y&lat  = lat
    1.66 +  y&lon  = lon
    1.67 +  y@units     = "none"
    1.68 +  y@long_name = "Leaf Area Index"
    1.69 +  printVarSummary (y)
    1.70 +
    1.71 +  c->LAI  = y
    1.72 +
    1.73 +  z = dim_avg(y(lat|:,lon|:,time|:))
    1.74 +
    1.75 +  z!0  = "lat"
    1.76 +  z!1  = "lon"
    1.77 +  z&lat  = lat
    1.78 +  z&lon  = lon
    1.79 +  z@units     = "none"
    1.80 +  z@long_name = "Leaf Area Index"
    1.81 +  printVarSummary (z)
    1.82 +  
    1.83 +  d->LAI  = z     
    1.84 +end
    1.85 +