clamp: beta/02.site.ncl@0c6405ab2ff4

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

     2 ; combine  scatter, histogram, global and zonal plots

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

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

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

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

     8 ; **********************************************************

    10 ; Main code.

    11 begin

    13 ;===================================================

    14 ; observed station name and location (lat,lon)

    15 ;===================================================

    17  beta_4_ob = 0.60

    19  station = (/"DukeFACE" \

    20             ,"AspenFACE" \

    21             ,"ORNL-FACE" \

    22             ,"POP-EUROFACE" \

    23             /)

    25  lat_ob = (/ 35.58,  45.40,  35.54, 42.22/)

    26  lon_ob = (/-79.05, -89.37, -84.20, 11.48/)

    27  lon_ob = where(lon_ob.lt.0.,lon_ob+360.,lon_ob)

    28 ;print (lon_ob)

    30  n_sta  = dimsizes(station)

    32 ;====================================================

    33 ; model data

    34 ;====================================================

    36  co2_i = 283.1878

    37  co2_f = 364.1252

    39 ;model_name_i = "i01.07cn"

    40 ;model_name_f = "i01.10cn"

    42  model_name_i = "i01.07casa"

    43  model_name_f = "i01.10casa"

    45  dirm = "/fis/cgd/cseg/people/jeff/clamp_data/model/"

    46  film_i = model_name_i + "_1990-2004_ANN_climo.nc"

    47  film_f = model_name_f + "_1990-2004_ANN_climo.nc"

    49  fm_i   = addfile (dirm+film_i,"r")

    50  fm_f   = addfile (dirm+film_f,"r")

    52  npp_i  = fm_i->NPP

    53  npp_f  = fm_f->NPP

    55  xm       = fm_i->lon

    56  ym       = fm_i->lat

    58  npp_i_4  =linint2_points(xm,ym,npp_i,True,lon_ob,lat_ob,0)

    60  npp_f_4  =linint2_points(xm,ym,npp_f,True,lon_ob,lat_ob,0)

    62 ; Units for these variables are:

    63 ; npp_i: g C/m^2/s

    65  print (npp_i_4)

    66  print (npp_f_4)

    67 ;============================

    68 ;compute beta

    69 ;============================

    71  beta_4 = new((/n_sta/),float)

    73  beta_4 = ((npp_f_4/npp_i_4) - 1.)/log(co2_f/co2_i)

    75  beta_4_avg = avg(beta_4)

    77  print (beta_4)

    78  print (beta_4_avg)

    80  M_beta = abs((beta_4_avg/beta_4_ob) - 1.)* 5.

    82  print (M_beta)

    83 end

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