clamp: npp/51.zonavg.ncl@4be95183fbcd

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

     2 ; xy_1.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 ;f     = addfile ("i01.03cn_1545-1569_ANN_climo.nc","r")

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

    14  u     = f->NPP

    15  v     = zonalAve(u)

    17  nsec_per_year = 60*60*24*365

    19  v = v * nsec_per_year

    21  v@long_name = "NPP (gC/m2/year)"

    22 ;************************************************

    23 ; plotting parameters

    24 ;************************************************

    25  wks   = gsn_open_wks ("png","xy")

    27  res                  = True

    28 ;res@tiMainString     = "Model i01.03cn"

    29  res@tiMainString     = "Model i01.04casa"

    31  plot  = gsn_csm_xy (wks,v&lat,v(0,:),res)

    32 end

author	Forrest Hoffman <forrest@climatemodeling.org>
	Thu, 26 Mar 2009 14:02:21 -0400
changeset 1	4be95183fbcd
permissions	-rw-r--r--