diff -r 0c6405ab2ff4 -r 4be95183fbcd all/11.ameriflux.ncl --- a/all/11.ameriflux.ncl Mon Jan 26 22:08:20 2009 -0500 +++ b/all/11.ameriflux.ncl Thu Mar 26 14:02:21 2009 -0400 @@ -557,7 +557,8 @@ ; compute annual cycle correlation coef and M score ;************************************************************ - score_max = 1. + ;score_max = 1. + score_max = (/6., 0., 9., 9., 6., 0./) ccr = new ((/nstation, nfield/),float) M_score = new ((/nstation, nfield/),float) @@ -566,7 +567,7 @@ do m=0,nfield-1 ccr(n,m) = esccr(data_ob_ann(m,:,n),data_mod_ann(m,:,n),0) 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)))) - M_score(n,m) = (1. -(bias/nmonth)) * score_max + M_score(n,m) = (1. -(bias/nmonth)) * score_max(m) end do end do