equal
deleted
inserted
replaced
|
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_T42_mean.nc","r") |
|
13 c = g->NPP |
|
14 u = ndtooned(c) |
|
15 |
|
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 |
|
19 v = ndtooned(b) |
|
20 |
|
21 sec_to_year = 86400.*365. |
|
22 v = v * sec_to_year |
|
23 |
|
24 good = ind(.not.ismissing(u) .and. .not.ismissing(v)) |
|
25 uu = u(good) |
|
26 nx = dimsizes(uu) |
|
27 ;print (nx) |
|
28 |
|
29 do i = 0,nx-1 |
|
30 if (uu(i) .lt. 1.) then |
|
31 uu(i) =10. |
|
32 end if |
|
33 end do |
|
34 |
|
35 ccr = esccr(uu,v(good),0) |
|
36 print (ccr) |
|
37 ;print (uu) |
|
38 bias = sum(((v(good)-uu)/uu)^2) |
|
39 print (bias) |
|
40 M = 1.- sqrt(bias/nx) |
|
41 print (M) |
|
42 |
|
43 end |