1.1 --- /dev/null Thu Jan 01 00:00:00 1970 +0000
1.2 +++ b/lai/04.lai_ensemble_1.9.ncl Mon Jan 26 22:08:20 2009 -0500
1.3 @@ -0,0 +1,46 @@
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/lai/ob/"
1.14 +
1.15 + diro = "/fis/cgd/cseg/people/jeff/clamp_data/lai/ob/"
1.16 + filo1 = "LAI_"+year_start+"-"+year_end+"_MONS_1.9.nc"
1.17 + filo2 = "LAI_"+year_start+"-"+year_end+"_ANN_1.9.nc"
1.18 + c = addfile(diro+filo1,"c")
1.19 + d = addfile(diro+filo2,"c")
1.20 +
1.21 + nyear = year_end - year_start + 1
1.22 +
1.23 + files = new((/nyear/),"string")
1.24 +
1.25 + do n = year_start,year_end
1.26 + i = n - year_start
1.27 + files(i)= diri + "LAI_" + n + "_monthly_1.9.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 + x@_FillValue = 1.e+36
1.36 +
1.37 + y = dim_avg_Wrap(x(time|:,lat|:,lon|:,case|:))
1.38 + printVarSummary (y)
1.39 + y@_FillValue = 1.e+36
1.40 +
1.41 + c->LAI = y
1.42 +
1.43 + z = dim_avg_Wrap(y(lat|:,lon|:,time|:))
1.44 + printVarSummary (z)
1.45 + z@_FillValue = 1.e+36
1.46 +
1.47 + d->LAI = z
1.48 +end
1.49 +