author | Forrest Hoffman <forrest@climatemodeling.org> |
Mon, 26 Jan 2009 22:08:20 -0500 | |
changeset 0 | 0c6405ab2ff4 |
permissions | -rw-r--r-- |
forrest@0 | 1 |
;************************************************* |
forrest@0 | 2 |
; ce_1.ncl |
forrest@0 | 3 |
;************************************************ |
forrest@0 | 4 |
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" |
forrest@0 | 5 |
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" |
forrest@0 | 6 |
;************************************************ |
forrest@0 | 7 |
begin |
forrest@0 | 8 |
|
forrest@0 | 9 |
;************************************************ |
forrest@0 | 10 |
; read in model data |
forrest@0 | 11 |
;************************************************ |
forrest@0 | 12 |
diri = "/fis/cgd/cseg/people/jeff/clamp/class/" |
forrest@0 | 13 |
fili = "T42.clamp.surface-data.060412.nc" |
forrest@0 | 14 |
f = addfile(diri+fili,"r") |
forrest@0 | 15 |
|
forrest@0 | 16 |
y = f->PCT_PFT |
forrest@0 | 17 |
printVarSummary(y) |
forrest@0 | 18 |
|
forrest@0 | 19 |
x = y(0,:,:) |
forrest@0 | 20 |
|
forrest@0 | 21 |
dimx = dimsizes(y) |
forrest@0 | 22 |
nlat = dimx(1) |
forrest@0 | 23 |
nlon = dimx(2) |
forrest@0 | 24 |
|
forrest@0 | 25 |
do j= 0,nlat-1 |
forrest@0 | 26 |
do i= 0,nlon-1 |
forrest@0 | 27 |
x(j,i) = sum(y(:,j,i)) |
forrest@0 | 28 |
end do |
forrest@0 | 29 |
end do |
forrest@0 | 30 |
|
forrest@0 | 31 |
print (x) |
forrest@0 | 32 |
|
forrest@0 | 33 |
end |