equal
deleted
inserted
replaced
|
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) |
|
16 |
|
17 nsec_per_year = 60*60*24*365 |
|
18 |
|
19 v = v * nsec_per_year |
|
20 |
|
21 v@long_name = "NPP (gC/m2/year)" |
|
22 ;************************************************ |
|
23 ; plotting parameters |
|
24 ;************************************************ |
|
25 wks = gsn_open_wks ("png","xy") |
|
26 |
|
27 res = True |
|
28 ;res@tiMainString = "Model i01.03cn" |
|
29 res@tiMainString = "Model i01.04casa" |
|
30 |
|
31 plot = gsn_csm_xy (wks,v&lat,v(0,:),res) |
|
32 end |