Mercurial Mercurial > repos > public_software > clamp / file revision
summary | shortlog | changelog | graph | tags | bookmarks | branches | files | changeset | file | latest | revisions | annotate | diff | comparison | raw | help
lai/04.lai_ensemble_T42.ncl.x
author Forrest Hoffman <forrest@climatemodeling.org>
Mon, 26 Jan 2009 22:08:20 -0500
changeset 0 0c6405ab2ff4
permissions -rw-r--r--
Initial commit of C-LAMP Diagnostics from Jeff Lee 
     1 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"  

     2 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"  

     3 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"  

     4 

     5 begin 

     6 ;----------------------------------------------------------

     7   year_start = 2000

     8   year_end   = 2005

     9 ;----------------------------------------------------------

    10   diri  = "/fis/cgd/cseg/people/jeff/clamp_data/"

    11   diro  = "/fis/cgd/cseg/people/jeff/clamp_data/"

    12   filo1 = "MOD15A2_LAI_"+year_start+"-"+year_end+"_ensemble.nc"

    13   filo2 = "MOD15A2_LAI_"+year_start+"-"+year_end+"_mean.nc"

    14   c = addfile(diro+filo1,"c")

    15   d = addfile(diro+filo2,"c")

    16 

    17   nyear = year_end - year_start + 1

    18 

    19   files  = new((/nyear/),"string")

    20 

    21   do n = year_start,year_end

    22      i = n - year_start

    23 ;    files(i)= diri + "MOD15A2_LAI_" + n + "_monthly.nc"

    24      files(i)= diri + "MOD15A2_LAI_" + n + "_monthly_T42.nc"

    25   end do

    26   print (files)  

    27 

    28   file_list = addfiles(files,"r")

    29   ListSetType(file_list,"join")

    30   x = addfiles_GetVar(file_list, files, "LAI")

    31   printVarSummary (x)

    32 

    33   y = addfiles_GetVar(file_list, files, "time")

    34   time = y(0,:)

    35   delete(time@case)

    36   delete(y)

    37   printVarSummary (time)

    38 

    39   y = addfiles_GetVar(file_list, files, "lat")

    40   lat = y(0,:)

    41   delete(lat@case)

    42   delete(y)

    43   printVarSummary (lat)

    44 

    45   y = addfiles_GetVar(file_list, files, "lon")

    46   lon = y(0,:)

    47   delete(lon@case)

    48   delete(y)

    49   printVarSummary (lon)

    50 

    51   x!0  = "case"

    52   x!1  = "time"

    53   x!2  = "lat"

    54   x!3  = "lon"

    55 

    56   y = dim_avg(x(time|:,lat|:,lon|:,case|:))

    57 

    58   y!0  = "time"

    59   y!1  = "lat"

    60   y!2  = "lon"

    61   y&time = time

    62   y&lat  = lat

    63   y&lon  = lon

    64   y@units     = "none"

    65   y@long_name = "Leaf Area Index"

    66   printVarSummary (y)

    67 

    68   c->LAI  = y

    69 

    70   z = dim_avg(y(lat|:,lon|:,time|:))

    71 

    72   z!0  = "lat"

    73   z!1  = "lon"

    74   z&lat  = lat

    75   z&lon  = lon

    76   z@units     = "none"

    77   z@long_name = "Leaf Area Index"

    78   printVarSummary (z)

    79   

    80   d->LAI  = z     

    81 end

    82
clamp