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