forrest@0: ; ***********************************************
forrest@0: ; xy_4.ncl
forrest@0: ; ***********************************************
forrest@0: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
forrest@0: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
forrest@0: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
forrest@0: ;************************************************
forrest@0: begin
forrest@0: ;************************************************
forrest@0: ; read in data
forrest@0: ;************************************************
forrest@0:  g     = addfile ("Npp_T42_mean.nc","r")
forrest@0:  c     = g->NPP  
forrest@0:  u     = ndtooned(c)
forrest@0: 
forrest@0: ;f     = addfile ("i01.03cn_1545-1569_ANN_climo.nc","r")
forrest@0:  f     = addfile ("i01.04casa_1605-1629_ANN_climo.nc","r")
forrest@0:  b     = f->NPP      
forrest@0:  v     = ndtooned(b)
forrest@0: 
forrest@0:  sec_to_year = 86400.*365.
forrest@0:  v = v * sec_to_year
forrest@0: 
forrest@0:  good = ind(.not.ismissing(u) .and. .not.ismissing(v))
forrest@0:  uu = u(good)
forrest@0:  nx = dimsizes(uu)
forrest@0: ;print (nx)
forrest@0: 
forrest@0:  do i = 0,nx-1
forrest@0:     if (uu(i) .lt. 1.) then
forrest@0:        uu(i) =10.
forrest@0:     end if
forrest@0:  end do
forrest@0: 
forrest@0:  ccr = esccr(uu,v(good),0)
forrest@0:  print (ccr)
forrest@0: ;print (uu)
forrest@0:  bias = sum(((v(good)-uu)/uu)^2)
forrest@0:  print (bias)
forrest@0:  M    = 1.- sqrt(bias/nx)
forrest@0:  print (M)
forrest@0: 
forrest@0: end