Modifications to scoring and graphics production for the final version of code for the C-LAMP paper in GCB.
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 ;************************************************
11 ;************************************************
12 g = addfile ("data.81.nc","r")
13 ;g = addfile ("data.933.nc","r")
24 if (xo(i) .lt. 0.) then
30 f = addfile ("i01.03cn_1545-1569_ANN_climo.nc","r")
31 ;f = addfile ("i01.04casa_1605-1629_ANN_climo.nc","r")
38 sec_to_year = 86400.*365.
40 bo = linint2_points(xi,yi,b,True,xo,yo,0) * sec_to_year
47 ;bias = sum(((bo(0,:)-c(:))/c(:))^2)
48 ;M = (1. - sqrt(bias/nx))*5.
51 bias = sum(abs(bo(0,:)-c(:))/(bo(0,:)+c(:)))
52 M = (1. - (bias/nx))*5.