1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
|
from __future__ import print_function
import cmor
import numpy
breq = "100000. 80000. 80000. 68000. 68000. 56000. 56000. 44000. 44000. 31000. 31000. 18000. 18000. 0.".split()
bnds_req = []
for b in breq:
bnds_req.append(float(b))
bnds_req = numpy.array(bnds_req)
bnds_req.shape = (7, 2)
print(bnds_req[-2], bnds_req.shape)
levs = []
for b in bnds_req:
levs.append((b[0] + b[1]) / 2.)
levs = numpy.array(levs)
print(levs)
cmor.setup(inpath="Tables",
set_verbosity=cmor.CMOR_NORMAL,
netcdf_file_action=cmor.CMOR_REPLACE,
logfile=None)
cmor.dataset_json("Test/common_user_inputNOBOUNDS.json")
cmor.load_table("TestTables/python_test_table_A")
nlat = 90
dlat = 180 / nlat
nlon = 180
dlon = 360. / nlon
lats = numpy.arange(-90 + dlat / 2., 90, dlat)
lons = numpy.arange(0, 360., dlon)
ntime = 12
data = numpy.random.random((ntime, 7, nlat, nlon)) + 280.
itim = cmor.axis(
table_entry='time',
coord_vals=numpy.arange(
0,
ntime,
1),
units='month since 2008')
ilat = cmor.axis(
table_entry='latitude',
coord_vals=lats,
units='degrees_north')
ilon = cmor.axis(
table_entry='longitude',
coord_vals=lons,
units='degrees_east')
print('so far', itim, ilat, ilon)
print(bnds_req)
ilev = cmor.axis(
table_entry="pressure2",
coord_vals=levs,
cell_bounds=bnds_req,
units="Pa")
iv = cmor.variable(
table_entry='ta', axis_ids=numpy.array(
(itim, ilev, ilat, ilon)), units='K')
cmor.write(iv, data)
|
