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
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
|
from __future__ import print_function
import cmor
import numpy
error_flag = cmor.setup(inpath='Tables', netcdf_file_action=cmor.CMOR_REPLACE)
error_flag = cmor.dataset_json("Test/CMOR_input_example.json")
# creates 1 degree grid
nlat = 180
nlon = 360
alats = numpy.arange(180) - 89.5
bnds_lat = numpy.arange(181) - 90
alons = numpy.arange(360) + .5
bnds_lon = numpy.arange(361)
cmor.load_table("Tables/CMIP6_Amon.json")
# cmor.load_table("Test/IPCC_table_A1")
ilat = cmor.axis(
table_entry='latitude',
units='degrees_north',
length=nlat,
coord_vals=alats,
cell_bounds=bnds_lat)
ilon = cmor.axis(
table_entry='longitude',
length=nlon,
units='degrees_east',
coord_vals=alons,
cell_bounds=bnds_lon)
mlev_val = """
0.000000 0.000000 0.000000 0.000000 0.000000
0.000000 0.000000 0.000000 0.000000 0.000000
0.000000 0.000000 0.000000 0.000000 0.000000
0.000000 0.000000 0.000000 0.000000 0.000000
0.000200 0.001650 0.006050 0.014750 0.028650
0.048250 0.073700 0.104950 0.141700 0.183550
0.229950 0.280200 0.333650 0.389650 0.447450
0.506300 0.565500 0.624350 0.682100 0.738000
0.791300 0.841100 0.886350 0.925950 0.958600
0.982650 0.996150""".split()
levs = []
for l in mlev_val:
levs.append(float(l))
BS_bnds = """
0.000000 0.000000 0.000000 0.000000 0.000000
0.000000 0.000000 0.000000 0.000000 0.000000
0.000000 0.000000 0.000000 0.000000 0.000000
0.000000 0.000000 0.000000 0.000000 0.000000
0.000000 0.000400 0.002900 0.009200 0.020300
0.037000 0.059500 0.087900 0.122000 0.161400
0.205700 0.254200 0.306200 0.361100 0.418200
0.476700 0.535900 0.595100 0.653600 0.710600
0.765400 0.817200 0.865000 0.907700 0.944200
0.973000 0.992300 1.000000
""".split()
levs_bnds = []
for l in BS_bnds:
levs_bnds.append(float(l))
levs = numpy.array(levs)
levs_bnds = numpy.array(levs_bnds)
nlevs = len(levs)
ntimes = 12
itim = cmor.axis(
table_entry='time',
units='months since 2030-1-1',
length=ntimes,
interval='1 month')
zlevs = numpy.array((0.1, 0.3, 0.55, 0.7, 0.9))
zlev_bnds = numpy.array((0., .2, .42, .62, .8, 1.))
table_entry = 'alternate_hybrid_sigma'
# for i in range(nlevs):
# print i,levs_bnds[i],levs[i],levs_bnds[i+1]
# if not (levs_bnds[i]<=levs[i]<=levs_bnds[i+1]) :
# print 'Yikes'
ilev = cmor.axis(
table_entry=table_entry,
units='',
length=nlevs,
coord_vals=levs,
cell_bounds=levs_bnds)
|
