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
|
from __future__ import print_function
import cmor
import numpy
nlat = 10
dlat = 180. / nlat
nlon = 20
dlon = 360. / nlon
nlev = 17
ntimes = 1
lats = numpy.arange(90 - dlat / 2., -90, -dlat)
blats = numpy.arange(90, -90 - dlat, -dlat)
lats2 = numpy.arange(-90 + dlat / 2., 90, dlat)
blats2 = numpy.arange(-90, 90 + dlat, dlat)
lons = numpy.arange(0 + dlon / 2., 360., dlon)
blons = numpy.arange(0, 360. + dlon, dlon)
cmor.setup(inpath='.', netcdf_file_action=cmor.CMOR_REPLACE)
cmor.dataset_json("Test/CMOR_input_example.json")
table = 'Tables/CMIP6_Amon.json'
cmor.load_table(table)
data = lats[:, numpy.newaxis] * lons[numpy.newaxis, :]
data = (data + 29000) / 750. + 233.2
ilat = cmor.axis(
table_entry='latitude',
coord_vals=lats,
cell_bounds=blats,
units='degrees_north')
ilat2 = cmor.axis(
table_entry='latitude',
coord_vals=lats2,
cell_bounds=blats2,
units='degrees_north')
ilon = cmor.axis(
table_entry='longitude',
coord_vals=lons,
cell_bounds=blons,
units='degrees_east')
itim = cmor.axis(table_entry='time', units='months since 2010')
ivar = cmor.variable(
table_entry='tasmin',
units='K',
axis_ids=[
itim,
ilat,
ilon])
cmor.write(ivar, data, ntimes_passed=1, time_vals=[0.], time_bnds=[0., 1.])
fnm = cmor.close(ivar, file_name=True)
print('*******************************')
ivar2 = cmor.variable(
table_entry='tasmin',
units='K',
axis_ids=[
itim,
ilat2,
ilon])
cmor.write(ivar2, data, ntimes_passed=1, time_vals=[1.], time_bnds=[1., 2.])
fnm2 = cmor.close(ivar2, file_name=True)
## import cdms2,vcs
# x=vcs.init()
# x.portrait()
## import EzTemplate
# M=EzTemplate.Multi(columns=1,rows=3)
# print 'Getting templates'
# t1=M.get()
# t2=M.get()
# print 'Plotting'
# x.plot(data,t1)
# f=cdms2.open(fnm)
# s=f("tasmin")
# print 'Ok data read, shape:',s.shape
# x.plot(s,t2)
# f=cdms2.open(fnm2)
# s=f("tasmin")
# print 'Ok data read, shape:',s.shape
# x.plot(s,M.get())
# raw_input()
|
