1 ; ****************************************************
2 ; combine scatter, histogram, global and zonal plots
3 ; compute all correlation coef and M score
4 ; ****************************************************
5 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl.test"
6 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl.test"
7 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
8 ;************************************************
13 dirm = "/fis/cgd/cseg/people/jeff/clamp_data/model/"
14 film = "surf_co2.clamp.T42.20070208.nc"
15 fm = addfile (dirm+film,"r")
17 surf_co2 = fm->surf_co2
19 print (surf_co2(:,0,0))
21 avg_co2_1990_2004 = avg(surf_co2(192:206,0,0))
23 print (avg_co2_1990_2004)