1 ; ***********************************************
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 ;************************************************
9 ;************************************************
10 ; read in observed data
11 ;************************************************
12 f = addfile ("Npp_T42_mean.nc","r")
14 ;************************************************
16 ;************************************************
17 g = addfile ("i01.03cn_1545-1569_ANN_climo.nc","r")
18 ;g = addfile ("i01.04casa_1605-1629_ANN_climo.nc","r")
23 y = where(ismissing(y).and.(ismissing(x0).or.(x0.lt.delta)),0.,y)
27 sec_to_year = 86400.*365.
37 bias = sum(((v-u)/u)^2)
39 M = 1.- sqrt(bias/dimsizes(u))