forrest@0: ; ****************************************************
forrest@0: ; combine scatter, histogram, global and zonal plots
forrest@0: ; compute all correlation coef and M score
forrest@0: ; ****************************************************
forrest@0: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl.test"
forrest@0: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl.test"
forrest@0: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
forrest@0: ;************************************************
forrest@0: 
forrest@0: ; Main code.
forrest@0: begin
forrest@0: 
forrest@0:  dirm = "/fis/cgd/cseg/people/jeff/clamp_data/model/"
forrest@0:  film = "surf_co2.clamp.T42.20070208.nc"
forrest@0:  fm   = addfile (dirm+film,"r")
forrest@0:   
forrest@0:  surf_co2  = fm->surf_co2
forrest@0:  
forrest@0:  print (surf_co2(:,0,0)) 
forrest@0: 
forrest@0:  avg_co2_1990_2004 = avg(surf_co2(192:206,0,0))
forrest@0: 
forrest@0:  print (avg_co2_1990_2004)
forrest@0: 
forrest@0: end