1.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
     1.2 +++ b/beta/01.co2_f.ncl	Mon Jan 26 22:08:20 2009 -0500
     1.3 @@ -0,0 +1,25 @@
     1.4 +; ****************************************************
     1.5 +; combine scatter, histogram, global and zonal plots
     1.6 +; compute all correlation coef and M score
     1.7 +; ****************************************************
     1.8 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl.test"
     1.9 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl.test"
    1.10 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
    1.11 +;************************************************
    1.12 +
    1.13 +; Main code.
    1.14 +begin
    1.15 +
    1.16 + dirm = "/fis/cgd/cseg/people/jeff/clamp_data/model/"
    1.17 + film = "surf_co2.clamp.T42.20070208.nc"
    1.18 + fm   = addfile (dirm+film,"r")
    1.19 +  
    1.20 + surf_co2  = fm->surf_co2
    1.21 + 
    1.22 + print (surf_co2(:,0,0)) 
    1.23 +
    1.24 + avg_co2_1990_2004 = avg(surf_co2(192:206,0,0))
    1.25 +
    1.26 + print (avg_co2_1990_2004)
    1.27 +
    1.28 +end