1.1 --- a/all/11.ameriflux.ncl Mon Jan 26 22:08:20 2009 -0500
1.2 +++ b/all/11.ameriflux.ncl Thu Mar 26 14:02:21 2009 -0400
1.3 @@ -557,7 +557,8 @@
1.4 ; compute annual cycle correlation coef and M score
1.5 ;************************************************************
1.6
1.7 - score_max = 1.
1.8 + ;score_max = 1.
1.9 + score_max = (/6., 0., 9., 9., 6., 0./)
1.10
1.11 ccr = new ((/nstation, nfield/),float)
1.12 M_score = new ((/nstation, nfield/),float)
1.13 @@ -566,7 +567,7 @@
1.14 do m=0,nfield-1
1.15 ccr(n,m) = esccr(data_ob_ann(m,:,n),data_mod_ann(m,:,n),0)
1.16 bias = sum(abs(data_mod_ann(m,:,n)-data_ob_ann(m,:,n))/(abs(data_mod_ann(m,:,n))+abs(data_ob_ann(m,:,n))))
1.17 - M_score(n,m) = (1. -(bias/nmonth)) * score_max
1.18 + M_score(n,m) = (1. -(bias/nmonth)) * score_max(m)
1.19 end do
1.20 end do
1.21