1.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
     1.2 +++ b/npp/42.bias_normalized.ncl	Mon Jan 26 22:08:20 2009 -0500
     1.3 @@ -0,0 +1,60 @@
     1.4 +; ***********************************************
     1.5 +; xy_4.ncl
     1.6 +; ***********************************************
     1.7 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
     1.8 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
     1.9 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
    1.10 +;************************************************
    1.11 +begin
    1.12 +;************************************************
    1.13 +; read in data
    1.14 +;************************************************
    1.15 + g     = addfile ("data.81.nc","r")
    1.16 + a     = g->SITE_ID  
    1.17 + xo    = g->LONG_DD  
    1.18 + yo    = g->LAT_DD
    1.19 + c     = g->TNPP_C
    1.20 +;c     = g->ANPP_C
    1.21 +;c     = g->BNPP_C
    1.22 + cp    = g->PREC_ANN
    1.23 + do i= 0,dimsizes(cp)-1
    1.24 + if (ismissing(cp(i))) then
    1.25 +    cp (i) = 100
    1.26 + end if
    1.27 + end do  
    1.28 + d     = int2flt(c)/int2flt(cp)
    1.29 + print (d)
    1.30 +
    1.31 + nx = dimsizes(xo)
    1.32 + do i= 0,nx-1
    1.33 +    if (xo(i) .lt. 0.) then
    1.34 +        xo(i) = xo(i)+ 360.
    1.35 +    end if
    1.36 + end do
    1.37 +;print (xo)
    1.38 + 
    1.39 + f     = addfile ("i01.03cn_1545-1569_ANN_climo.nc","r")
    1.40 +;f     = addfile ("i01.04casa_1605-1629_ANN_climo.nc","r")
    1.41 + b     = f->NPP
    1.42 +;b     = f->AGNPP
    1.43 +;b     = f->BGNPP    
    1.44 + xi    = f->lon     
    1.45 + yi    = f->lat
    1.46 + bp    = f->RAIN      
    1.47 +
    1.48 + sec_to_year = 86400.*365.
    1.49 +
    1.50 + bo = linint2_points(xi,yi,b,True,xo,yo,0) * sec_to_year
    1.51 +;print (bo)
    1.52 + bq = linint2_points(xi,yi,bp,True,xo,yo,0) * sec_to_year 
    1.53 +
    1.54 + bo = bo/bq
    1.55 + 
    1.56 + ccr = esccr(bo,d,0)
    1.57 + print (ccr)
    1.58 + bias = sum(((bo(0,:)-d(:))/d(:))^2)
    1.59 + print (bias)
    1.60 + M    = 1. - sqrt(bias/nx)
    1.61 + print (M)
    1.62 +
    1.63 +end