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
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
|
"""
Test SPECS DRS structure.
NOTE: this structure is not finalised!
"""
from drslib.specs import SpecsFileSystem, SpecsDRS
from .drs_tree_shared import TestEg, TestListing
specs_fs = SpecsFileSystem('/specs')
def test_1():
path = '/SPECS/output/MPI-M/MPI-ESM-LR/decadal/S19610101/day/atmos/day/pr/r1i1p1/v20010112/pr_day_MPI-ESM-LR_decadal_S19610101_r1i1p1_19610101-19701231.nc'
drs = specs_fs.filepath_to_drs(path)
drs.activity = 'specs'
assert drs.activity == 'specs'
assert drs.product == 'output'
assert drs.institute == 'MPI-M'
assert drs.model == 'MPI-ESM-LR'
assert drs.experiment == 'decadal'
assert drs.start_date[:3] == (1961, 1, 1)
assert drs.frequency == 'day'
assert drs.table == 'day'
assert drs.realm == 'atmos'
assert drs.ensemble == (1,1,1)
assert drs.variable == 'pr'
assert drs.version == 20010112
assert drs.subset[0][:3] == (1961, 1, 1)
assert drs.subset[1][:3] == (1970, 12, 31)
def test_2():
drs_id = 'specs.output.MPI-M.MPI-ESM-LR.decadal.S19610101.day.atmos.day.pr.r1i1p1.v20010112'
drs = specs_fs.drs_cls.from_dataset_id(drs_id)
assert drs.activity == 'specs'
assert drs.product == 'output'
assert drs.institute == 'MPI-M'
assert drs.model == 'MPI-ESM-LR'
assert drs.experiment == 'decadal'
assert drs.start_date[:3] == (1961, 1, 1)
assert drs.frequency == 'day'
assert drs.table == 'day'
assert drs.realm == 'atmos'
assert drs.ensemble == (1,1,1)
assert drs.variable == 'pr'
assert drs.version == 20010112
assert drs.subset == None
def test_3():
filename = 'pr_day_MPI-ESM-LR_decadal_S19610101_r1i1p1_19610101-19701231.nc'
drs = specs_fs.filename_to_drs(filename)
assert drs.activity == 'specs'
assert drs.experiment == 'decadal'
assert drs.ensemble == (1,1,1)
assert drs.model == 'MPI-ESM-LR'
assert drs.start_date[:3] == (1961, 1, 1)
assert drs.table == 'day'
assert drs.variable == 'pr'
assert drs.subset[0][:3] == (1961, 1, 1)
assert drs.subset[1][:3] == (1970, 12, 31)
def test_4():
filename = 'msftmyzba_Omon_IPSL-CM5A-LR_decadal_S19620101_r3i1p1_196201-197112.nc'
drs = specs_fs.filename_to_drs(filename)
assert drs.activity == 'specs'
assert drs.variable == 'msftmyzba'
assert drs.table == 'Omon'
assert drs.model == 'IPSL-CM5A-LR'
assert drs.experiment == 'decadal'
assert drs.start_date[:3] == (1962, 1, 1)
assert drs.ensemble == (3,1,1)
assert drs.subset[0][:2] == (1962, 0o1)
assert drs.subset[1][:2] == (1971, 12)
def test_5():
filename = 'goro_fx_IPSL-CM5A-LR_decadal_S19820101_r0i0p0.nc'
drs = specs_fs.filename_to_drs(filename)
assert drs.activity == 'specs'
assert drs.variable == 'goro'
assert drs.table == 'fx'
assert drs.model == 'IPSL-CM5A-LR'
assert drs.ensemble == (0,0,0)
assert drs.start_date[:3] == (1982, 1, 1)
class TestSpecsListing1(TestListing):
__test__ = True
listing_file = 'specs_1.ls'
def setUp(self):
super(TestSpecsListing1, self).setUp()
# incoming is not tmpdir/output.
self.incoming = self.tmpdir
def _init_drs_fs(self):
self.drs_fs = SpecsFileSystem(self.tmpdir)
def test_1(self):
self.dt.discover(self.incoming, activity='specs',
product='output',
realm='ocean',
frequency='mon',
institute='IPSL',
)
print (len(self.dt.pub_trees))
assert len(self.dt.pub_trees) == 3
def test_2(self):
self.dt.discover(self.incoming, activity='specs',
product='output',
realm='ocean',
frequency='mon',
institute='IPSL',
)
print (len(self.dt.pub_trees))
assert len(self.dt.pub_trees) == 3
pt = list(self.dt.pub_trees.values())[0]
self._do_version(pt)
pt.deduce_state()
assert len(list(pt.versions.values())[0]) == 1
class TestSpecsListing2(TestSpecsListing1):
# Test avoiding .ftpaccess files.
listing_file = 'specs_2.ls'
class TestSpecsMapfile(TestListing):
__test__ = True
listing_file = 'specs_1.ls'
# Test mapfile generation
def setUp(self):
super(TestSpecsMapfile, self).setUp()
# incoming is not tmpdir/output.
self.incoming = self.tmpdir
self.dt.discover(self.incoming, activity='specs',
product='output',
realm='ocean',
frequency='mon',
institute='IPSL',
)
self.pt = list(self.dt.pub_trees.values())[0]
self._do_version(self.pt)
def _init_drs_fs(self):
self.drs_fs = SpecsFileSystem(self.tmpdir)
# def version_to_mapfile(self, version, fh=None, checksum_func=None):
def test_mapfile(self):
from io import StringIO
fh = StringIO()
self.pt.version_to_mapfile(list(self.pt.versions.keys())[0], fh)
assert 'specs.output.IPSL.IPSL-CM5A-LR.decadal.S19630101.mon.ocean.day.clt.r3i1p1' in fh.getvalue()
#-----------------------------------------------------------------------------
# drs_tool tests
from drslib.drs_command import main as drstool_main
def test_componentparse1():
# Test that non-string components are parsed correctly, e.g. start_date
drstool_main('drs-tool list --scheme=specs --component=start_date=S20100101 -R drs_tool_example'.split())
|
