diff -r 000000000000 -r 0c6405ab2ff4 beta/01.co2_f.ncl --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/beta/01.co2_f.ncl Mon Jan 26 22:08:20 2009 -0500 @@ -0,0 +1,25 @@ +; **************************************************** +; combine scatter, histogram, global and zonal plots +; compute all correlation coef and M score +; **************************************************** +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl.test" +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl.test" +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" +;************************************************ + +; Main code. +begin + + dirm = "/fis/cgd/cseg/people/jeff/clamp_data/model/" + film = "surf_co2.clamp.T42.20070208.nc" + fm = addfile (dirm+film,"r") + + surf_co2 = fm->surf_co2 + + print (surf_co2(:,0,0)) + + avg_co2_1990_2004 = avg(surf_co2(192:206,0,0)) + + print (avg_co2_1990_2004) + +end