|
1 ; **************************************************** |
|
2 ; combine scatter, histogram, global and zonal plots |
|
3 ; ********************************************************* |
|
4 |
|
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 ; ********************************************************** |
|
9 |
|
10 ; Main code. |
|
11 begin |
|
12 |
|
13 ;=================================================== |
|
14 ; observed station name and location (lat,lon) |
|
15 ;=================================================== |
|
16 |
|
17 beta_4_ob = 0.60 |
|
18 |
|
19 station = (/"DukeFACE" \ |
|
20 ,"AspenFACE" \ |
|
21 ,"ORNL-FACE" \ |
|
22 ,"POP-EUROFACE" \ |
|
23 /) |
|
24 |
|
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) |
|
29 |
|
30 n_sta = dimsizes(station) |
|
31 |
|
32 ;==================================================== |
|
33 ; model data |
|
34 ;==================================================== |
|
35 |
|
36 co2_i = 283.1878 |
|
37 co2_f = 364.1252 |
|
38 |
|
39 ;model_name_i = "i01.07cn" |
|
40 ;model_name_f = "i01.10cn" |
|
41 |
|
42 model_name_i = "i01.07casa" |
|
43 model_name_f = "i01.10casa" |
|
44 |
|
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" |
|
48 |
|
49 fm_i = addfile (dirm+film_i,"r") |
|
50 fm_f = addfile (dirm+film_f,"r") |
|
51 |
|
52 npp_i = fm_i->NPP |
|
53 npp_f = fm_f->NPP |
|
54 |
|
55 xm = fm_i->lon |
|
56 ym = fm_i->lat |
|
57 |
|
58 npp_i_4 =linint2_points(xm,ym,npp_i,True,lon_ob,lat_ob,0) |
|
59 |
|
60 npp_f_4 =linint2_points(xm,ym,npp_f,True,lon_ob,lat_ob,0) |
|
61 |
|
62 ; Units for these variables are: |
|
63 ; npp_i: g C/m^2/s |
|
64 |
|
65 print (npp_i_4) |
|
66 print (npp_f_4) |
|
67 ;============================ |
|
68 ;compute beta |
|
69 ;============================ |
|
70 |
|
71 beta_4 = new((/n_sta/),float) |
|
72 |
|
73 beta_4 = ((npp_f_4/npp_i_4) - 1.)/log(co2_f/co2_i) |
|
74 |
|
75 beta_4_avg = avg(beta_4) |
|
76 |
|
77 print (beta_4) |
|
78 print (beta_4_avg) |
|
79 |
|
80 M_beta = abs((beta_4_avg/beta_4_ob) - 1.)* 5. |
|
81 |
|
82 print (M_beta) |
|
83 end |