clamp: npp/42.bias_global.nc@4be95183fbcd

     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)

    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)

    21  sec_to_year = 86400.*365.

    22  v = v * sec_to_year

    24  good = ind(.not.ismissing(u) .and. .not.ismissing(v))

    25  uu = u(good)

    26  nx = dimsizes(uu)

    27 ;print (nx)

    29  do i = 0,nx-1

    30     if (uu(i) .lt. 1.) then

    31        uu(i) =10.

    32     end if

    33  end do

    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)

    43 end

author	Forrest Hoffman <forrest@climatemodeling.org>
	Thu, 26 Mar 2009 14:02:21 -0400
changeset 1	4be95183fbcd
permissions	-rw-r--r--