|
forrest@0
|
1 |
;************************************************
|
|
forrest@0
|
2 |
; Read ascii, Write nc
|
|
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 |
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
|
|
forrest@0
|
7 |
;************************************************
|
|
forrest@0
|
8 |
begin
|
|
forrest@0
|
9 |
diri = "/fis/cgd/cseg/people/jeff/clamp_data/co2/globalView-C4Mip/seas/"
|
|
forrest@0
|
10 |
b = (/"zot005_45D2_seas.co2", "zot015_45D2_seas.co2"/)
|
|
forrest@0
|
11 |
do n = 0,1
|
|
forrest@0
|
12 |
print (b(n))
|
|
forrest@0
|
13 |
ncol1 = 4
|
|
forrest@0
|
14 |
data1 = readAsciiTable(diri+b(n), ncol1, "float", (/10,17/))
|
|
forrest@0
|
15 |
print (data1)
|
|
forrest@0
|
16 |
|
|
forrest@0
|
17 |
ncol2 = 5
|
|
forrest@0
|
18 |
data2 = readAsciiTable(diri+b(n), ncol2, "float", 16)
|
|
forrest@0
|
19 |
print (data2)
|
|
forrest@0
|
20 |
end do
|
|
forrest@0
|
21 |
|
|
forrest@0
|
22 |
end
|