clamp: npp/41.correlation_zonal.ncl@0c6405ab2ff4

     1 ; ***********************************************

     2 ; xy_4.ncl

     3 ; ***********************************************

     4 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"

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

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

     7 ;************************************************

     8 begin

     9 ;************************************************

    10 ; read in data

    11 ;************************************************

    12  g     = addfile ("Npp_T31_mean.nc","r")

    13  c     = g->NPP

    14  u     = zonalAve(c)

    16 ;f     = addfile ("i01.03cn_1545-1569_ANN_climo.nc","r")

    17  f     = addfile ("i01.04casa_1605-1629_ANN_climo.nc","r")

    18  b     = f->NPP

    20  sec_to_year = 86400.*365.

    22  b = b * sec_to_year

    23  v = zonalAve(b)

    25  b@long_name = "NPP (gC/m2/year)"

    27  ccr = esccr(u,v(0,:),0)

    29  print (ccr)

    31 end

author	Forrest Hoffman <forrest@climatemodeling.org>
	Mon, 26 Jan 2009 22:08:20 -0500
changeset 0	0c6405ab2ff4
permissions	-rw-r--r--