Context Navigation

source: flex_extract.git/Documentation/html/_modules/tools.html @ 0690a6c

ctbtodev

Last change on this file since 0690a6c was ba99230, checked in by Anne Philipp <anne.philipp@…>, 5 years ago
capitalize all directory names and adapt pathes in files
Property mode set to `100644`
File size: 95.7 KB

Line
1
2
3	<!DOCTYPE html>
4	<!--[if IE 8]><html class="no-js lt-ie9" lang="en" > <![endif]-->
5	<!--[if gt IE 8]><!--> <html class="no-js" lang="en" > <!--<![endif]-->
6	<head>
7	<meta charset="utf-8">
8
9	<meta name="viewport" content="width=device-width, initial-scale=1.0">
10
11	<title>tools — flex_extract 7.1 alpha documentation</title>
12
13
14
15
16
17
18
19
20	<script type="text/javascript" src="../_static/js/modernizr.min.js"></script>
21
22
23	<script type="text/javascript" id="documentation_options" data-url_root="../" src="../_static/documentation_options.js"></script>
24	<script type="text/javascript" src="../_static/jquery.js"></script>
25	<script type="text/javascript" src="../_static/underscore.js"></script>
26	<script type="text/javascript" src="../_static/doctools.js"></script>
27	<script type="text/javascript" src="../_static/language_data.js"></script>
28	<script async="async" type="text/javascript" src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/latest.js?config=TeX-AMS-MML_HTMLorMML"></script>
29
30	<script type="text/javascript" src="../_static/js/theme.js"></script>
31
32
33
34
35	<link rel="stylesheet" href="../_static/css/theme.css" type="text/css" />
36	<link rel="stylesheet" href="../_static/pygments.css" type="text/css" />
37	<link rel="stylesheet" href="../_static/css/custom.css" type="text/css" />
38	<link rel="stylesheet" href="../_static/css/theme_overrides.css" type="text/css" />
39	<link rel="index" title="Index" href="../genindex.html" />
40	<link rel="search" title="Search" href="../search.html" />
41	</head>
42
43	<body class="wy-body-for-nav">
44
45
46	<div class="wy-grid-for-nav">
47
48	<nav data-toggle="wy-nav-shift" class="wy-nav-side">
49	<div class="wy-side-scroll">
50	<div class="wy-side-nav-search" >
51
52
53
54	<a href="../index.html" class="icon icon-home"> flex_extract
55
56
57
58	</a>
59
60
61
62
63	<div class="version">
64	7.1
65	</div>
66
67
68
69
70	<div role="search">
71	<form id="rtd-search-form" class="wy-form" action="../search.html" method="get">
72	<input type="text" name="q" placeholder="Search docs" />
73	<input type="hidden" name="check_keywords" value="yes" />
74	<input type="hidden" name="area" value="default" />
75	</form>
76	</div>
77
78
79	</div>
80
81	<div class="wy-menu wy-menu-vertical" data-spy="affix" role="navigation" aria-label="main navigation">
82
83
84
85
86
87
88	<p class="caption"><span class="caption-text">Table of Contents:</span></p>
89	<ul>
90	<li class="toctree-l1"><a class="reference internal" href="../installation.html">Installation</a></li>
91	<li class="toctree-l1"><a class="reference internal" href="../quick_start.html">Quick Start</a></li>
92	<li class="toctree-l1"><a class="reference internal" href="../ecmwf_data.html">ECMWF Data</a></li>
93	<li class="toctree-l1"><a class="reference internal" href="../documentation.html">Documentation</a></li>
94	<li class="toctree-l1"><a class="reference internal" href="../evaluation.html">Evaluation</a></li>
95	<li class="toctree-l1"><a class="reference internal" href="../dev_guide.html">Developer Guide</a></li>
96	<li class="toctree-l1"><a class="reference internal" href="../changelog.html">Changelog</a></li>
97	<li class="toctree-l1"><a class="reference internal" href="../support.html">Support</a></li>
98	<li class="toctree-l1"><a class="reference internal" href="../Support/faq.html">FAQ - Frequently asked questions</a></li>
99	<li class="toctree-l1"><a class="reference internal" href="../authors.html">Developer Team</a></li>
100	</ul>
101
102
103
104	</div>
105	</div>
106	</nav>
107
108	<section data-toggle="wy-nav-shift" class="wy-nav-content-wrap">
109
110
111	<nav class="wy-nav-top" aria-label="top navigation">
112
113	<i data-toggle="wy-nav-top" class="fa fa-bars"></i>
114	<a href="../index.html">flex_extract</a>
115
116	</nav>
117
118
119	<div class="wy-nav-content">
120
121	<div class="rst-content">
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139	<div role="navigation" aria-label="breadcrumbs navigation">
140
141	<ul class="wy-breadcrumbs">
142
143	<li><a href="../index.html">Docs</a> »</li>
144
145	<li><a href="index.html">Module code</a> »</li>
146
147	<li>tools</li>
148
149
150	<li class="wy-breadcrumbs-aside">
151
152	</li>
153
154	</ul>
155
156
157	<hr/>
158	</div>
159	<div role="main" class="document" itemscope="itemscope" itemtype="http://schema.org/Article">
160	<div itemprop="articleBody">
161
162	<h1>Source code for tools</h1><div class="highlight"><pre>
163	<span></span><span class="ch">#!/usr/bin/env python</span>
164	<span class="c1"># -- coding: utf-8 --</span>
165	<span class="c1">#*******************************************************************************</span>
166	<span class="c1"># @Author: Anne Philipp (University of Vienna)</span>
167	<span class="c1">#</span>
168	<span class="c1"># @Date: May 2018</span>
169	<span class="c1">#</span>
170	<span class="c1"># @Change History:</span>
171	<span class="c1"># October 2014 - Anne Fouilloux (University of Oslo)</span>
172	<span class="c1"># - created functions silent_remove and product (taken from ECMWF)</span>
173	<span class="c1">#</span>
174	<span class="c1"># November 2015 - Leopold Haimberger (University of Vienna)</span>
175	<span class="c1"># - created functions: interpret_args_and_control, clean_up</span>
176	<span class="c1"># my_error, normal_exit, init128, to_param_id</span>
177	<span class="c1">#</span>
178	<span class="c1"># April - December 2018 - Anne Philipp (University of Vienna):</span>
179	<span class="c1"># - applied PEP8 style guide</span>
180	<span class="c1"># - added documentation</span>
181	<span class="c1"># - moved all non class methods from former file Flexparttools in here</span>
182	<span class="c1"># - seperated args and control interpretation</span>
183	<span class="c1"># - added functions get_list_as_string, read_ecenv, send_mail, make_dir,</span>
184	<span class="c1"># put_file_to_ecserver, submit_job_to_ecserver, get_informations,</span>
185	<span class="c1"># get_dimensions, execute_subprocess, none_or_int, none_or_str</span>
186	<span class="c1">#</span>
187	<span class="c1"># @License:</span>
188	<span class="c1"># (C) Copyright 2014-2019.</span>
189	<span class="c1"># Anne Philipp, Leopold Haimberger</span>
190	<span class="c1">#</span>
191	<span class="c1"># This work is licensed under the Creative Commons Attribution 4.0</span>
192	<span class="c1"># International License. To view a copy of this license, visit</span>
193	<span class="c1"># http://creativecommons.org/licenses/by/4.0/ or send a letter to</span>
194	<span class="c1"># Creative Commons, PO Box 1866, Mountain View, CA 94042, USA.</span>
195	<span class="c1">#</span>
196	<span class="c1"># @Methods:</span>
197	<span class="c1"># none_or_str</span>
198	<span class="c1"># none_or_int</span>
199	<span class="c1"># get_cmdline_args</span>
200	<span class="c1"># read_ecenv</span>
201	<span class="c1"># clean_up</span>
202	<span class="c1"># my_error</span>
203	<span class="c1"># send_mail</span>
204	<span class="c1"># normal_exit</span>
205	<span class="c1"># product</span>
206	<span class="c1"># silent_remove</span>
207	<span class="c1"># init128</span>
208	<span class="c1"># to_param_id</span>
209	<span class="c1"># get_list_as_string</span>
210	<span class="c1"># make_dir</span>
211	<span class="c1"># put_file_to_ecserver</span>
212	<span class="c1"># submit_job_to_ecserver</span>
213	<span class="c1"># get_informations</span>
214	<span class="c1"># get_dimensions</span>
215	<span class="c1"># execute_subprocess</span>
216	<span class="c1">#*******************************************************************************</span>
217	<span class="sd">'''This module contains a collection of diverse tasks within flex_extract.</span>
218	<span class="sd">'''</span>
219
220	<span class="c1"># ------------------------------------------------------------------------------</span>
221	<span class="c1"># MODULES</span>
222	<span class="c1"># ------------------------------------------------------------------------------</span>
223	<span class="kn">from</span> <span class="nn">__future__</span> <span class="k">import</span> <span class="n">print_function</span>
224
225	<span class="kn">import</span> <span class="nn">os</span>
226	<span class="kn">import</span> <span class="nn">errno</span>
227	<span class="kn">import</span> <span class="nn">sys</span>
228	<span class="kn">import</span> <span class="nn">glob</span>
229	<span class="kn">import</span> <span class="nn">subprocess</span>
230	<span class="kn">import</span> <span class="nn">traceback</span>
231	<span class="k">try</span><span class="p">:</span>
232	<span class="kn">import</span> <span class="nn">exceptions</span>
233	<span class="k">except</span> <span class="ne">ImportError</span><span class="p">:</span>
234	<span class="kn">import</span> <span class="nn">builtins</span> <span class="k">as</span> <span class="nn">exceptions</span>
235	<span class="kn">from</span> <span class="nn">datetime</span> <span class="k">import</span> <span class="n">datetime</span><span class="p">,</span> <span class="n">timedelta</span>
236	<span class="kn">from</span> <span class="nn">argparse</span> <span class="k">import</span> <span class="n">ArgumentParser</span><span class="p">,</span> <span class="n">ArgumentDefaultsHelpFormatter</span>
237
238	<span class="c1"># ------------------------------------------------------------------------------</span>
239	<span class="c1"># METHODS</span>
240	<span class="c1"># ------------------------------------------------------------------------------</span>
241
242	<div class="viewcode-block" id="setup_controldata"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.setup_controldata">[docs]</a><span class="k">def</span> <span class="nf">setup_controldata</span><span class="p">():</span>
243	<span class="sd">'''Collects, stores and checks controlling arguments from command line,</span>
244	<span class="sd"> CONTROL file and ECMWF_ENV file.</span>
245
246	<span class="sd"> Parameters</span>
247	<span class="sd"> ----------</span>
248
249	<span class="sd"> Return</span>
250	<span class="sd"> ------</span>
251	<span class="sd"> c : ControlFile</span>
252	<span class="sd"> Contains all the parameters of CONTROL file and</span>
253	<span class="sd"> command line.</span>
254
255	<span class="sd"> ppid : str</span>
256	<span class="sd"> Parent process id.</span>
257
258	<span class="sd"> queue : str</span>
259	<span class="sd"> Name of queue for submission to ECMWF (e.g. ecgate or cca )</span>
260
261	<span class="sd"> job_template : str</span>
262	<span class="sd"> Name of the job template file for submission to ECMWF server.</span>
263	<span class="sd"> '''</span>
264	<span class="kn">import</span> <span class="nn">_config</span>
265	<span class="kn">from</span> <span class="nn">Classes.ControlFile</span> <span class="k">import</span> <span class="n">ControlFile</span>
266
267	<span class="n">args</span> <span class="o">=</span> <span class="n">get_cmdline_args</span><span class="p">()</span>
268	<span class="n">c</span> <span class="o">=</span> <span class="n">ControlFile</span><span class="p">(</span><span class="n">args</span><span class="o">.</span><span class="n">controlfile</span><span class="p">)</span>
269	<span class="n">c</span><span class="o">.</span><span class="n">assign_args_to_control</span><span class="p">(</span><span class="n">args</span><span class="p">)</span>
270	<span class="k">if</span> <span class="n">os</span><span class="o">.</span><span class="n">path</span><span class="o">.</span><span class="n">isfile</span><span class="p">(</span><span class="n">_config</span><span class="o">.</span><span class="n">PATH_ECMWF_ENV</span><span class="p">):</span>
271	<span class="n">env_parameter</span> <span class="o">=</span> <span class="n">read_ecenv</span><span class="p">(</span><span class="n">_config</span><span class="o">.</span><span class="n">PATH_ECMWF_ENV</span><span class="p">)</span>
272	<span class="n">c</span><span class="o">.</span><span class="n">assign_envs_to_control</span><span class="p">(</span><span class="n">env_parameter</span><span class="p">)</span>
273	<span class="n">c</span><span class="o">.</span><span class="n">check_conditions</span><span class="p">(</span><span class="n">args</span><span class="o">.</span><span class="n">queue</span><span class="p">)</span>
274
275	<span class="k">return</span> <span class="n">c</span><span class="p">,</span> <span class="n">args</span><span class="o">.</span><span class="n">ppid</span><span class="p">,</span> <span class="n">args</span><span class="o">.</span><span class="n">queue</span><span class="p">,</span> <span class="n">args</span><span class="o">.</span><span class="n">job_template</span></div>
276
277	<div class="viewcode-block" id="none_or_str"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.none_or_str">[docs]</a><span class="k">def</span> <span class="nf">none_or_str</span><span class="p">(</span><span class="n">value</span><span class="p">):</span>
278	<span class="sd">'''Converts the input string into pythons None-type if the string</span>
279	<span class="sd"> contains string "None".</span>
280
281	<span class="sd"> Parameters</span>
282	<span class="sd"> ----------</span>
283	<span class="sd"> value : str</span>
284	<span class="sd"> String to be checked for the "None" word.</span>
285
286	<span class="sd"> Return</span>
287	<span class="sd"> ------</span>
288	<span class="sd"> None or value:</span>
289	<span class="sd"> Return depends on the content of the input value. If it was "None",</span>
290	<span class="sd"> then the python type None is returned. Otherwise the string itself.</span>
291	<span class="sd"> '''</span>
292	<span class="k">if</span> <span class="n">value</span> <span class="o">==</span> <span class="s1">'None'</span><span class="p">:</span>
293	<span class="k">return</span> <span class="kc">None</span>
294	<span class="k">return</span> <span class="n">value</span></div>
295
296	<div class="viewcode-block" id="none_or_int"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.none_or_int">[docs]</a><span class="k">def</span> <span class="nf">none_or_int</span><span class="p">(</span><span class="n">value</span><span class="p">):</span>
297	<span class="sd">'''Converts the input string into pythons None-type if the string</span>
298	<span class="sd"> contains string "None". Otherwise it is converted to an integer value.</span>
299
300	<span class="sd"> Parameters</span>
301	<span class="sd"> ----------</span>
302	<span class="sd"> value : str</span>
303	<span class="sd"> String to be checked for the "None" word.</span>
304
305	<span class="sd"> Return</span>
306	<span class="sd"> ------</span>
307	<span class="sd"> None or int(value):</span>
308	<span class="sd"> Return depends on the content of the input value. If it was "None",</span>
309	<span class="sd"> then the python type None is returned. Otherwise the string is</span>
310	<span class="sd"> converted into an integer value.</span>
311	<span class="sd"> '''</span>
312	<span class="k">if</span> <span class="n">value</span> <span class="o">==</span> <span class="s1">'None'</span><span class="p">:</span>
313	<span class="k">return</span> <span class="kc">None</span>
314	<span class="k">return</span> <span class="nb">int</span><span class="p">(</span><span class="n">value</span><span class="p">)</span></div>
315
316	<div class="viewcode-block" id="get_cmdline_args"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.get_cmdline_args">[docs]</a><span class="k">def</span> <span class="nf">get_cmdline_args</span><span class="p">():</span>
317	<span class="sd">'''Decomposes the command line arguments and assigns them to variables.</span>
318	<span class="sd"> Apply default values for non mentioned arguments.</span>
319
320	<span class="sd"> Parameters</span>
321	<span class="sd"> ----------</span>
322
323	<span class="sd"> Return</span>
324	<span class="sd"> ------</span>
325	<span class="sd"> args : Namespace</span>
326	<span class="sd"> Contains the commandline arguments from script/program call.</span>
327	<span class="sd"> '''</span>
328
329	<span class="n">parser</span> <span class="o">=</span> <span class="n">ArgumentParser</span><span class="p">(</span><span class="n">description</span><span class="o">=</span><span class="s1">'Retrieve FLEXPART input from </span><span class="se">\</span>
330	<span class="s1"> ECMWF MARS archive'</span><span class="p">,</span>
331	<span class="n">formatter_class</span><span class="o">=</span><span class="n">ArgumentDefaultsHelpFormatter</span><span class="p">)</span>
332
333	<span class="c1"># control parameters that override control file values</span>
334	<span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">"--start_date"</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s2">"start_date"</span><span class="p">,</span>
335	<span class="nb">type</span><span class="o">=</span><span class="n">none_or_str</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span>
336	<span class="n">help</span><span class="o">=</span><span class="s2">"start date YYYYMMDD"</span><span class="p">)</span>
337	<span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">"--end_date"</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s2">"end_date"</span><span class="p">,</span>
338	<span class="nb">type</span><span class="o">=</span><span class="n">none_or_str</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span>
339	<span class="n">help</span><span class="o">=</span><span class="s2">"end_date YYYYMMDD"</span><span class="p">)</span>
340	<span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">"--date_chunk"</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s2">"date_chunk"</span><span class="p">,</span>
341	<span class="nb">type</span><span class="o">=</span><span class="n">none_or_int</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span>
342	<span class="n">help</span><span class="o">=</span><span class="s2">"# of days to be retrieved at once"</span><span class="p">)</span>
343	<span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">"--job_chunk"</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s2">"job_chunk"</span><span class="p">,</span>
344	<span class="nb">type</span><span class="o">=</span><span class="n">none_or_int</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span>
345	<span class="n">help</span><span class="o">=</span><span class="s2">"# of days to be retrieved within a single job"</span><span class="p">)</span>
346	<span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">"--controlfile"</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s2">"controlfile"</span><span class="p">,</span>
347	<span class="nb">type</span><span class="o">=</span><span class="n">none_or_str</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="s1">'CONTROL_EA5'</span><span class="p">,</span>
348	<span class="n">help</span><span class="o">=</span><span class="s2">"The file with all CONTROL parameters."</span><span class="p">)</span>
349	<span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">"--basetime"</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s2">"basetime"</span><span class="p">,</span>
350	<span class="nb">type</span><span class="o">=</span><span class="n">none_or_int</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span>
351	<span class="n">help</span><span class="o">=</span><span class="s2">"base such as 0 or 12 (for half day retrievals)"</span><span class="p">)</span>
352	<span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">"--step"</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s2">"step"</span><span class="p">,</span>
353	<span class="nb">type</span><span class="o">=</span><span class="n">none_or_str</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span>
354	<span class="n">help</span><span class="o">=</span><span class="s2">"Forecast steps such as 00/to/48"</span><span class="p">)</span>
355	<span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">"--levelist"</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s2">"levelist"</span><span class="p">,</span>
356	<span class="nb">type</span><span class="o">=</span><span class="n">none_or_str</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span>
357	<span class="n">help</span><span class="o">=</span><span class="s2">"Vertical levels to be retrieved, e.g. 30/to/60"</span><span class="p">)</span>
358	<span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">"--area"</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s2">"area"</span><span class="p">,</span>
359	<span class="nb">type</span><span class="o">=</span><span class="n">none_or_str</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span>
360	<span class="n">help</span><span class="o">=</span><span class="s2">"area defined as north/west/south/east"</span><span class="p">)</span>
361
362	<span class="c1"># some switches</span>
363	<span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">"--debug"</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s2">"debug"</span><span class="p">,</span>
364	<span class="nb">type</span><span class="o">=</span><span class="n">none_or_int</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span>
365	<span class="n">help</span><span class="o">=</span><span class="s2">"debug mode - leave temporary files intact"</span><span class="p">)</span>
366	<span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">"--oper"</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s2">"oper"</span><span class="p">,</span>
367	<span class="nb">type</span><span class="o">=</span><span class="n">none_or_int</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span>
368	<span class="n">help</span><span class="o">=</span><span class="s1">'operational mode - prepares dates with '</span>
369	<span class="s1">'environment variables'</span><span class="p">)</span>
370	<span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">"--request"</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s2">"request"</span><span class="p">,</span>
371	<span class="nb">type</span><span class="o">=</span><span class="n">none_or_int</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span>
372	<span class="n">help</span><span class="o">=</span><span class="s2">"list all mars requests in file mars_requests.dat"</span><span class="p">)</span>
373	<span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">"--public"</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s2">"public"</span><span class="p">,</span>
374	<span class="nb">type</span><span class="o">=</span><span class="n">none_or_int</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span>
375	<span class="n">help</span><span class="o">=</span><span class="s2">"public mode - retrieves the public datasets"</span><span class="p">)</span>
376	<span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">"--rrint"</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s2">"rrint"</span><span class="p">,</span>
377	<span class="nb">type</span><span class="o">=</span><span class="n">none_or_int</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span>
378	<span class="n">help</span><span class="o">=</span><span class="s1">'Selection of old or new precipitation '</span>
379	<span class="s1">'interpolation:</span><span class="se">\n</span><span class="s1">'</span>
380	<span class="s1">' 0 - old method</span><span class="se">\n</span><span class="s1">'</span>
381	<span class="s1">' 1 - new method (additional subgrid points)'</span><span class="p">)</span>
382
383	<span class="c1"># set directories</span>
384	<span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">"--inputdir"</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s2">"inputdir"</span><span class="p">,</span>
385	<span class="nb">type</span><span class="o">=</span><span class="n">none_or_str</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span>
386	<span class="n">help</span><span class="o">=</span><span class="s1">'Path to the temporary directory for the '</span>
387	<span class="s1">'retrieval grib files and other processing files.'</span><span class="p">)</span>
388	<span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">"--outputdir"</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s2">"outputdir"</span><span class="p">,</span>
389	<span class="nb">type</span><span class="o">=</span><span class="n">none_or_str</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span>
390	<span class="n">help</span><span class="o">=</span><span class="s1">'Path to the final directory where the final '</span>
391	<span class="s1">'FLEXPART ready input files are stored.'</span><span class="p">)</span>
392
393	<span class="c1"># this is only used by prepare_flexpart.py to rerun a postprocessing step</span>
394	<span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">"--ppid"</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s2">"ppid"</span><span class="p">,</span>
395	<span class="nb">type</span><span class="o">=</span><span class="n">none_or_str</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span>
396	<span class="n">help</span><span class="o">=</span><span class="s1">'This is the specify parent process id of a '</span>
397	<span class="s1">'single flex_extract run to identify the files. '</span>
398	<span class="s1">'It is the second number in the GRIB files.'</span><span class="p">)</span>
399
400	<span class="c1"># arguments for job submission to ECMWF, only needed by submit.py</span>
401	<span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">"--job_template"</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s1">'job_template'</span><span class="p">,</span>
402	<span class="nb">type</span><span class="o">=</span><span class="n">none_or_str</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="s2">"job.temp"</span><span class="p">,</span>
403	<span class="n">help</span><span class="o">=</span><span class="s1">'The job template file which are adapted to be '</span>
404	<span class="s1">'submitted to the batch system on ECMWF server.'</span><span class="p">)</span>
405	<span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">"--queue"</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s2">"queue"</span><span class="p">,</span>
406	<span class="nb">type</span><span class="o">=</span><span class="n">none_or_str</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span>
407	<span class="n">help</span><span class="o">=</span><span class="s1">'The ECMWF server name for submission of the '</span>
408	<span class="s1">'job script to the batch system '</span>
409	<span class="s1">'(e.g. ecgate \| cca \| ccb)'</span><span class="p">)</span>
410
411	<span class="n">args</span> <span class="o">=</span> <span class="n">parser</span><span class="o">.</span><span class="n">parse_args</span><span class="p">()</span>
412
413	<span class="k">return</span> <span class="n">args</span></div>
414
415	<div class="viewcode-block" id="read_ecenv"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.read_ecenv">[docs]</a><span class="k">def</span> <span class="nf">read_ecenv</span><span class="p">(</span><span class="n">filepath</span><span class="p">):</span>
416	<span class="sd">'''Reads the file into a dictionary where the key values are the parameter</span>
417	<span class="sd"> names.</span>
418
419	<span class="sd"> Parameters</span>
420	<span class="sd"> ----------</span>
421	<span class="sd"> filepath : str</span>
422	<span class="sd"> Path to file where the ECMWF environment parameters are stored.</span>
423
424	<span class="sd"> Return</span>
425	<span class="sd"> ------</span>
426	<span class="sd"> envs : dict</span>
427	<span class="sd"> Contains the environment parameter ecuid, ecgid, gateway</span>
428	<span class="sd"> and destination for ECMWF server environments.</span>
429	<span class="sd"> '''</span>
430	<span class="n">envs</span><span class="o">=</span> <span class="p">{}</span>
431	<span class="k">try</span><span class="p">:</span>
432	<span class="k">with</span> <span class="nb">open</span><span class="p">(</span><span class="n">filepath</span><span class="p">,</span> <span class="s1">'r'</span><span class="p">)</span> <span class="k">as</span> <span class="n">f</span><span class="p">:</span>
433	<span class="k">for</span> <span class="n">line</span> <span class="ow">in</span> <span class="n">f</span><span class="p">:</span>
434	<span class="n">data</span> <span class="o">=</span> <span class="n">line</span><span class="o">.</span><span class="n">strip</span><span class="p">()</span><span class="o">.</span><span class="n">split</span><span class="p">()</span>
435	<span class="n">envs</span><span class="p">[</span><span class="nb">str</span><span class="p">(</span><span class="n">data</span><span class="p">[</span><span class="mi">0</span><span class="p">])]</span> <span class="o">=</span> <span class="nb">str</span><span class="p">(</span><span class="n">data</span><span class="p">[</span><span class="mi">1</span><span class="p">])</span>
436	<span class="k">except</span> <span class="ne">OSError</span> <span class="k">as</span> <span class="n">e</span><span class="p">:</span>
437	<span class="nb">print</span><span class="p">(</span><span class="s1">'... ERROR CODE: '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">e</span><span class="o">.</span><span class="n">errno</span><span class="p">))</span>
438	<span class="nb">print</span><span class="p">(</span><span class="s1">'... ERROR MESSAGE:</span><span class="se">\n</span><span class="s1"> </span><span class="se">\t</span><span class="s1"> '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">e</span><span class="o">.</span><span class="n">strerror</span><span class="p">))</span>
439
440	<span class="n">sys</span><span class="o">.</span><span class="n">exit</span><span class="p">(</span><span class="s1">'</span><span class="se">\n</span><span class="s1">... Error occured while trying to read ECMWF_ENV '</span>
441	<span class="s1">'file: '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">filepath</span><span class="p">))</span>
442
443	<span class="k">return</span> <span class="n">envs</span></div>
444
445	<div class="viewcode-block" id="clean_up"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.clean_up">[docs]</a><span class="k">def</span> <span class="nf">clean_up</span><span class="p">(</span><span class="n">c</span><span class="p">):</span>
446	<span class="sd">'''Remove files from the intermediate directory (inputdir).</span>
447
448	<span class="sd"> It keeps the final FLEXPART input files if program runs without</span>
449	<span class="sd"> ECMWF Api and keywords "ectrans" or "ecstorage" are set to "1".</span>
450
451	<span class="sd"> Parameters</span>
452	<span class="sd"> ----------</span>
453	<span class="sd"> c : ControlFile</span>
454	<span class="sd"> Contains all the parameters of CONTROL file and</span>
455	<span class="sd"> command line.</span>
456
457	<span class="sd"> Return</span>
458	<span class="sd"> ------</span>
459
460	<span class="sd"> '''</span>
461
462	<span class="nb">print</span><span class="p">(</span><span class="s2">"... clean inputdir!"</span><span class="p">)</span>
463
464	<span class="n">cleanlist</span> <span class="o">=</span> <span class="p">[</span><span class="n">file</span> <span class="k">for</span> <span class="n">file</span> <span class="ow">in</span> <span class="n">glob</span><span class="o">.</span><span class="n">glob</span><span class="p">(</span><span class="n">os</span><span class="o">.</span><span class="n">path</span><span class="o">.</span><span class="n">join</span><span class="p">(</span><span class="n">c</span><span class="o">.</span><span class="n">inputdir</span><span class="p">,</span> <span class="s2">"*"</span><span class="p">))</span>
465	<span class="k">if</span> <span class="ow">not</span> <span class="n">os</span><span class="o">.</span><span class="n">path</span><span class="o">.</span><span class="n">basename</span><span class="p">(</span><span class="n">file</span><span class="p">)</span><span class="o">.</span><span class="n">startswith</span><span class="p">(</span><span class="n">c</span><span class="o">.</span><span class="n">prefix</span><span class="p">)]</span>
466
467	<span class="k">if</span> <span class="n">cleanlist</span><span class="p">:</span>
468	<span class="k">for</span> <span class="n">element</span> <span class="ow">in</span> <span class="n">cleanlist</span><span class="p">:</span>
469	<span class="n">silent_remove</span><span class="p">(</span><span class="n">element</span><span class="p">)</span>
470	<span class="nb">print</span><span class="p">(</span><span class="s2">"... done!"</span><span class="p">)</span>
471	<span class="k">else</span><span class="p">:</span>
472	<span class="nb">print</span><span class="p">(</span><span class="s2">"... nothing to clean!"</span><span class="p">)</span>
473
474	<span class="k">return</span></div>
475
476
477	<div class="viewcode-block" id="my_error"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.my_error">[docs]</a><span class="k">def</span> <span class="nf">my_error</span><span class="p">(</span><span class="n">message</span><span class="o">=</span><span class="s1">'ERROR'</span><span class="p">):</span>
478	<span class="sd">'''Prints a specified error message which can be passed to the function</span>
479	<span class="sd"> before exiting the program.</span>
480
481	<span class="sd"> Parameters</span>
482	<span class="sd"> ----------</span>
483	<span class="sd"> message : str, optional</span>
484	<span class="sd"> Error message. Default value is "ERROR".</span>
485
486	<span class="sd"> Return</span>
487	<span class="sd"> ------</span>
488
489	<span class="sd"> '''</span>
490
491	<span class="n">trace</span> <span class="o">=</span> <span class="s1">'</span><span class="se">\n</span><span class="s1">'</span><span class="o">.</span><span class="n">join</span><span class="p">(</span><span class="n">traceback</span><span class="o">.</span><span class="n">format_stack</span><span class="p">())</span>
492	<span class="n">full_message</span> <span class="o">=</span> <span class="n">message</span> <span class="o">+</span> <span class="s1">'</span><span class="se">\n\n</span><span class="s1">'</span> <span class="o">+</span> <span class="n">trace</span>
493
494	<span class="nb">print</span><span class="p">(</span><span class="n">full_message</span><span class="p">)</span>
495
496	<span class="n">sys</span><span class="o">.</span><span class="n">exit</span><span class="p">(</span><span class="mi">1</span><span class="p">)</span>
497
498	<span class="k">return</span></div>
499
500
501	<div class="viewcode-block" id="send_mail"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.send_mail">[docs]</a><span class="k">def</span> <span class="nf">send_mail</span><span class="p">(</span><span class="n">users</span><span class="p">,</span> <span class="n">success_mode</span><span class="p">,</span> <span class="n">message</span><span class="p">):</span>
502	<span class="sd">'''Prints a specific exit message which can be passed to the function.</span>
503
504	<span class="sd"> Parameters</span>
505	<span class="sd"> ----------</span>
506	<span class="sd"> users : list of str</span>
507	<span class="sd"> Contains all email addresses which should be notified.</span>
508	<span class="sd"> It might also contain just the ecmwf user name which wil trigger</span>
509	<span class="sd"> mailing to the associated email address for this user.</span>
510
511	<span class="sd"> success_mode : str</span>
512	<span class="sd"> States the exit mode of the program to put into</span>
513	<span class="sd"> the mail subject line.</span>
514
515	<span class="sd"> message : str, optional</span>
516	<span class="sd"> Message for exiting program. Default value is "Done!".</span>
517
518	<span class="sd"> Return</span>
519	<span class="sd"> ------</span>
520
521	<span class="sd"> '''</span>
522
523	<span class="k">for</span> <span class="n">user</span> <span class="ow">in</span> <span class="n">users</span><span class="p">:</span>
524	<span class="k">if</span> <span class="s1">'$</span><span class="si">{USER}</span><span class="s1">'</span> <span class="ow">in</span> <span class="n">user</span><span class="p">:</span>
525	<span class="n">user</span> <span class="o">=</span> <span class="n">os</span><span class="o">.</span><span class="n">getenv</span><span class="p">(</span><span class="s1">'USER'</span><span class="p">)</span>
526	<span class="k">try</span><span class="p">:</span>
527	<span class="n">p</span> <span class="o">=</span> <span class="n">subprocess</span><span class="o">.</span><span class="n">Popen</span><span class="p">([</span><span class="s1">'mail'</span><span class="p">,</span> <span class="s1">'-s flex_extract_v7.1 '</span> <span class="o">+</span>
528	<span class="n">success_mode</span><span class="p">,</span> <span class="n">os</span><span class="o">.</span><span class="n">path</span><span class="o">.</span><span class="n">expandvars</span><span class="p">(</span><span class="n">user</span><span class="p">)],</span>
529	<span class="n">stdin</span><span class="o">=</span><span class="n">subprocess</span><span class="o">.</span><span class="n">PIPE</span><span class="p">,</span>
530	<span class="n">stdout</span><span class="o">=</span><span class="n">subprocess</span><span class="o">.</span><span class="n">PIPE</span><span class="p">,</span>
531	<span class="n">stderr</span><span class="o">=</span><span class="n">subprocess</span><span class="o">.</span><span class="n">PIPE</span><span class="p">,</span>
532	<span class="n">bufsize</span><span class="o">=</span><span class="mi">1</span><span class="p">)</span>
533	<span class="n">pout</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">communicate</span><span class="p">(</span><span class="nb">input</span><span class="o">=</span><span class="n">message</span><span class="o">.</span><span class="n">encode</span><span class="p">()</span> <span class="o">+</span> <span class="s1">'</span><span class="se">\n\n</span><span class="s1">'</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
534	<span class="k">except</span> <span class="ne">ValueError</span> <span class="k">as</span> <span class="n">e</span><span class="p">:</span>
535	<span class="nb">print</span><span class="p">(</span><span class="s1">'... ERROR: '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">e</span><span class="p">))</span>
536	<span class="n">sys</span><span class="o">.</span><span class="n">exit</span><span class="p">(</span><span class="s1">'... Email could not be sent!'</span><span class="p">)</span>
537	<span class="k">except</span> <span class="ne">OSError</span> <span class="k">as</span> <span class="n">e</span><span class="p">:</span>
538	<span class="nb">print</span><span class="p">(</span><span class="s1">'... ERROR CODE: '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">e</span><span class="o">.</span><span class="n">errno</span><span class="p">))</span>
539	<span class="nb">print</span><span class="p">(</span><span class="s1">'... ERROR MESSAGE:</span><span class="se">\n</span><span class="s1"> </span><span class="se">\t</span><span class="s1"> '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">e</span><span class="o">.</span><span class="n">strerror</span><span class="p">))</span>
540	<span class="n">sys</span><span class="o">.</span><span class="n">exit</span><span class="p">(</span><span class="s1">'... Email could not be sent!'</span><span class="p">)</span>
541	<span class="k">else</span><span class="p">:</span>
542	<span class="nb">print</span><span class="p">(</span><span class="s1">'Email sent to '</span> <span class="o">+</span> <span class="n">os</span><span class="o">.</span><span class="n">path</span><span class="o">.</span><span class="n">expandvars</span><span class="p">(</span><span class="n">user</span><span class="p">))</span>
543
544	<span class="k">return</span></div>
545
546
547	<div class="viewcode-block" id="normal_exit"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.normal_exit">[docs]</a><span class="k">def</span> <span class="nf">normal_exit</span><span class="p">(</span><span class="n">message</span><span class="o">=</span><span class="s1">'Done!'</span><span class="p">):</span>
548	<span class="sd">'''Prints a specific exit message which can be passed to the function.</span>
549
550	<span class="sd"> Parameters</span>
551	<span class="sd"> ----------</span>
552	<span class="sd"> message : str, optional</span>
553	<span class="sd"> Message for exiting program. Default value is "Done!".</span>
554
555	<span class="sd"> Return</span>
556	<span class="sd"> ------</span>
557
558	<span class="sd"> '''</span>
559
560	<span class="nb">print</span><span class="p">(</span><span class="nb">str</span><span class="p">(</span><span class="n">message</span><span class="p">))</span>
561
562	<span class="k">return</span></div>
563
564
565	<div class="viewcode-block" id="product"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.product">[docs]</a><span class="k">def</span> <span class="nf">product</span><span class="p">(</span><span class="o"></span><span class="n">args</span><span class="p">,</span> <span class="o">*</span><span class="n">kwds</span><span class="p">):</span>
566	<span class="sd">'''Creates combinations of all passed arguments.</span>
567
568	<span class="sd"> This method combines the single characters of the passed arguments</span>
569	<span class="sd"> with each other. So that each character of each argument value</span>
570	<span class="sd"> will be combined with each character of the other arguments as a tuple.</span>
571
572	<span class="sd"> Note</span>
573	<span class="sd"> ----</span>
574	<span class="sd"> This method is taken from an example at the ECMWF wiki website.</span>
575	<span class="sd"> https://software.ecmwf.int/wiki/display/GRIB/index.py; 2018-03-16</span>
576
577	<span class="sd"> Example</span>
578	<span class="sd"> -------</span>
579	<span class="sd"> product('ABCD', 'xy') --> Ax Ay Bx By Cx Cy Dx Dy</span>
580
581	<span class="sd"> product(range(2), repeat = 3) --> 000 001 010 011 100 101 110 111</span>
582
583	<span class="sd"> Parameters</span>
584	<span class="sd"> ----------</span>
585	<span class="sd"> \*args : list or str</span>
586	<span class="sd"> Positional arguments (arbitrary number).</span>
587
588	<span class="sd"> \\kwds : dict</span>
589	<span class="sd"> Contains all the keyword arguments from \*args.</span>
590
591	<span class="sd"> Return</span>
592	<span class="sd"> ------</span>
593	<span class="sd"> prod : :obj:`tuple`</span>
594	<span class="sd"> Return will be done with "yield". A tuple of combined arguments.</span>
595	<span class="sd"> See example in description above.</span>
596	<span class="sd"> '''</span>
597	<span class="k">try</span><span class="p">:</span>
598	<span class="n">pools</span> <span class="o">=</span> <span class="p">[</span><span class="nb">tuple</span><span class="p">(</span><span class="n">arg</span><span class="p">)</span> <span class="k">for</span> <span class="n">arg</span> <span class="ow">in</span> <span class="n">args</span><span class="p">]</span> <span class="o">*</span> <span class="n">kwds</span><span class="o">.</span><span class="n">get</span><span class="p">(</span><span class="s1">'repeat'</span><span class="p">,</span> <span class="mi">1</span><span class="p">)</span>
599	<span class="n">result</span> <span class="o">=</span> <span class="p">[[]]</span>
600	<span class="k">for</span> <span class="n">pool</span> <span class="ow">in</span> <span class="n">pools</span><span class="p">:</span>
601	<span class="n">result</span> <span class="o">=</span> <span class="p">[</span><span class="n">x</span> <span class="o">+</span> <span class="p">[</span><span class="n">y</span><span class="p">]</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">result</span> <span class="k">for</span> <span class="n">y</span> <span class="ow">in</span> <span class="n">pool</span><span class="p">]</span>
602	<span class="k">for</span> <span class="n">prod</span> <span class="ow">in</span> <span class="n">result</span><span class="p">:</span>
603	<span class="k">yield</span> <span class="nb">tuple</span><span class="p">(</span><span class="n">prod</span><span class="p">)</span>
604	<span class="k">except</span> <span class="ne">TypeError</span> <span class="k">as</span> <span class="n">e</span><span class="p">:</span>
605	<span class="n">sys</span><span class="o">.</span><span class="n">exit</span><span class="p">(</span><span class="s1">'... PRODUCT GENERATION FAILED!'</span><span class="p">)</span>
606
607	<span class="k">return</span></div>
608
609
610	<div class="viewcode-block" id="silent_remove"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.silent_remove">[docs]</a><span class="k">def</span> <span class="nf">silent_remove</span><span class="p">(</span><span class="n">filename</span><span class="p">):</span>
611	<span class="sd">'''Remove file if it exists.</span>
612	<span class="sd"> The function does not fail if the file does not exist.</span>
613
614	<span class="sd"> Parameters</span>
615	<span class="sd"> ----------</span>
616	<span class="sd"> filename : str</span>
617	<span class="sd"> The name of the file to be removed without notification.</span>
618
619	<span class="sd"> Return</span>
620	<span class="sd"> ------</span>
621
622	<span class="sd"> '''</span>
623	<span class="k">try</span><span class="p">:</span>
624	<span class="n">os</span><span class="o">.</span><span class="n">remove</span><span class="p">(</span><span class="n">filename</span><span class="p">)</span>
625	<span class="k">except</span> <span class="ne">OSError</span> <span class="k">as</span> <span class="n">e</span><span class="p">:</span>
626	<span class="c1"># errno.ENOENT = no such file or directory</span>
627	<span class="k">if</span> <span class="n">e</span><span class="o">.</span><span class="n">errno</span> <span class="o">==</span> <span class="n">errno</span><span class="o">.</span><span class="n">ENOENT</span><span class="p">:</span>
628	<span class="k">pass</span>
629	<span class="k">else</span><span class="p">:</span>
630	<span class="k">raise</span> <span class="c1"># re-raise exception if a different error occured</span>
631
632	<span class="k">return</span></div>
633
634
635	<div class="viewcode-block" id="init128"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.init128">[docs]</a><span class="k">def</span> <span class="nf">init128</span><span class="p">(</span><span class="n">filepath</span><span class="p">):</span>
636	<span class="sd">'''Opens and reads the grib file with table 128 information.</span>
637
638	<span class="sd"> Parameters</span>
639	<span class="sd"> ----------</span>
640	<span class="sd"> filepath : str</span>
641	<span class="sd"> Path to file of ECMWF grib table number 128.</span>
642
643	<span class="sd"> Return</span>
644	<span class="sd"> ------</span>
645	<span class="sd"> table128 : dict</span>
646	<span class="sd"> Contains the ECMWF grib table 128 information.</span>
647	<span class="sd"> The key is the parameter number and the value is the</span>
648	<span class="sd"> short name of the parameter.</span>
649	<span class="sd"> '''</span>
650	<span class="n">table128</span> <span class="o">=</span> <span class="nb">dict</span><span class="p">()</span>
651	<span class="k">try</span><span class="p">:</span>
652	<span class="k">with</span> <span class="nb">open</span><span class="p">(</span><span class="n">filepath</span><span class="p">)</span> <span class="k">as</span> <span class="n">f</span><span class="p">:</span>
653	<span class="n">fdata</span> <span class="o">=</span> <span class="n">f</span><span class="o">.</span><span class="n">read</span><span class="p">()</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="s1">'</span><span class="se">\n</span><span class="s1">'</span><span class="p">)</span>
654	<span class="k">except</span> <span class="ne">OSError</span> <span class="k">as</span> <span class="n">e</span><span class="p">:</span>
655	<span class="nb">print</span><span class="p">(</span><span class="s1">'... ERROR CODE: '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">e</span><span class="o">.</span><span class="n">errno</span><span class="p">))</span>
656	<span class="nb">print</span><span class="p">(</span><span class="s1">'... ERROR MESSAGE:</span><span class="se">\n</span><span class="s1"> </span><span class="se">\t</span><span class="s1"> '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">e</span><span class="o">.</span><span class="n">strerror</span><span class="p">))</span>
657
658	<span class="n">sys</span><span class="o">.</span><span class="n">exit</span><span class="p">(</span><span class="s1">'</span><span class="se">\n</span><span class="s1">... Error occured while trying to read parameter '</span>
659	<span class="s1">'table file: '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">filepath</span><span class="p">))</span>
660	<span class="k">else</span><span class="p">:</span>
661	<span class="k">for</span> <span class="n">data</span> <span class="ow">in</span> <span class="n">fdata</span><span class="p">:</span>
662	<span class="k">if</span> <span class="n">data</span> <span class="o">!=</span> <span class="s1">''</span> <span class="ow">and</span> <span class="n">data</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span> <span class="o">!=</span> <span class="s1">'!'</span><span class="p">:</span>
663	<span class="n">table128</span><span class="p">[</span><span class="n">data</span><span class="p">[</span><span class="mi">0</span><span class="p">:</span><span class="mi">3</span><span class="p">]]</span> <span class="o">=</span> <span class="n">data</span><span class="p">[</span><span class="mi">59</span><span class="p">:</span><span class="mi">65</span><span class="p">]</span><span class="o">.</span><span class="n">strip</span><span class="p">()</span>
664
665	<span class="k">return</span> <span class="n">table128</span></div>
666
667
668	<div class="viewcode-block" id="to_param_id"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.to_param_id">[docs]</a><span class="k">def</span> <span class="nf">to_param_id</span><span class="p">(</span><span class="n">pars</span><span class="p">,</span> <span class="n">table</span><span class="p">):</span>
669	<span class="sd">'''Transform parameter names to parameter ids with ECMWF grib table 128.</span>
670
671	<span class="sd"> Parameters</span>
672	<span class="sd"> ----------</span>
673	<span class="sd"> pars : str</span>
674	<span class="sd"> Addpar argument from CONTROL file in the format of</span>
675	<span class="sd"> parameter names instead of ids. The parameter short</span>
676	<span class="sd"> names are sepearted with "/" and they are passed as</span>
677	<span class="sd"> one single string.</span>
678
679	<span class="sd"> table : dict</span>
680	<span class="sd"> Contains the ECMWF grib table 128 information.</span>
681	<span class="sd"> The key is the parameter number and the value is the</span>
682	<span class="sd"> short name of the parameter.</span>
683
684	<span class="sd"> Return</span>
685	<span class="sd"> ------</span>
686	<span class="sd"> ipar : list of int</span>
687	<span class="sd"> List of addpar parameters from CONTROL file transformed to</span>
688	<span class="sd"> parameter ids in the format of integer.</span>
689	<span class="sd"> '''</span>
690	<span class="k">if</span> <span class="ow">not</span> <span class="n">pars</span><span class="p">:</span>
691	<span class="k">return</span> <span class="p">[]</span>
692	<span class="k">if</span> <span class="ow">not</span> <span class="nb">isinstance</span><span class="p">(</span><span class="n">pars</span><span class="p">,</span> <span class="nb">str</span><span class="p">):</span>
693	<span class="n">pars</span><span class="o">=</span><span class="nb">str</span><span class="p">(</span><span class="n">pars</span><span class="p">)</span>
694
695	<span class="n">cpar</span> <span class="o">=</span> <span class="n">pars</span><span class="o">.</span><span class="n">upper</span><span class="p">()</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="s1">'/'</span><span class="p">)</span>
696	<span class="n">ipar</span> <span class="o">=</span> <span class="p">[]</span>
697	<span class="k">for</span> <span class="n">par</span> <span class="ow">in</span> <span class="n">cpar</span><span class="p">:</span>
698	<span class="n">par</span> <span class="o">=</span> <span class="n">par</span><span class="o">.</span><span class="n">strip</span><span class="p">()</span>
699	<span class="k">for</span> <span class="n">k</span><span class="p">,</span> <span class="n">v</span> <span class="ow">in</span> <span class="n">table</span><span class="o">.</span><span class="n">items</span><span class="p">():</span>
700	<span class="k">if</span> <span class="n">par</span><span class="o">.</span><span class="n">isdigit</span><span class="p">():</span>
701	<span class="n">par</span> <span class="o">=</span> <span class="nb">str</span><span class="p">(</span><span class="nb">int</span><span class="p">(</span><span class="n">par</span><span class="p">))</span><span class="o">.</span><span class="n">zfill</span><span class="p">(</span><span class="mi">3</span><span class="p">)</span>
702	<span class="k">if</span> <span class="n">par</span> <span class="o">==</span> <span class="n">k</span> <span class="ow">or</span> <span class="n">par</span> <span class="o">==</span> <span class="n">v</span><span class="p">:</span>
703	<span class="n">ipar</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="nb">int</span><span class="p">(</span><span class="n">k</span><span class="p">))</span>
704	<span class="k">break</span>
705	<span class="k">else</span><span class="p">:</span>
706	<span class="nb">print</span><span class="p">(</span><span class="s1">'Warning: par '</span> <span class="o">+</span> <span class="n">par</span> <span class="o">+</span> <span class="s1">' not found in table 128'</span><span class="p">)</span>
707
708	<span class="k">return</span> <span class="n">ipar</span></div>
709
710	<div class="viewcode-block" id="to_param_id_with_tablenumber"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.to_param_id_with_tablenumber">[docs]</a><span class="k">def</span> <span class="nf">to_param_id_with_tablenumber</span><span class="p">(</span><span class="n">pars</span><span class="p">,</span> <span class="n">table</span><span class="p">):</span>
711	<span class="sd">'''Transform parameter names to parameter ids and add table id.</span>
712
713	<span class="sd"> Conversion with ECMWF grib table 128.</span>
714
715	<span class="sd"> Parameters</span>
716	<span class="sd"> ----------</span>
717	<span class="sd"> pars : str</span>
718	<span class="sd"> Addpar argument from CONTROL file in the format of</span>
719	<span class="sd"> parameter names instead of ids. The parameter short</span>
720	<span class="sd"> names are sepearted with "/" and they are passed as</span>
721	<span class="sd"> one single string.</span>
722
723	<span class="sd"> table : dict</span>
724	<span class="sd"> Contains the ECMWF grib table 128 information.</span>
725	<span class="sd"> The key is the parameter number and the value is the</span>
726	<span class="sd"> short name of the parameter.</span>
727
728	<span class="sd"> Return</span>
729	<span class="sd"> ------</span>
730	<span class="sd"> spar : str</span>
731	<span class="sd"> List of addpar parameters from CONTROL file transformed to</span>
732	<span class="sd"> parameter ids in the format of integer.</span>
733	<span class="sd"> '''</span>
734	<span class="k">if</span> <span class="ow">not</span> <span class="n">pars</span><span class="p">:</span>
735	<span class="k">return</span> <span class="p">[]</span>
736	<span class="k">if</span> <span class="ow">not</span> <span class="nb">isinstance</span><span class="p">(</span><span class="n">pars</span><span class="p">,</span> <span class="nb">str</span><span class="p">):</span>
737	<span class="n">pars</span><span class="o">=</span><span class="nb">str</span><span class="p">(</span><span class="n">pars</span><span class="p">)</span>
738
739	<span class="n">cpar</span> <span class="o">=</span> <span class="n">pars</span><span class="o">.</span><span class="n">upper</span><span class="p">()</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="s1">'/'</span><span class="p">)</span>
740	<span class="n">spar</span> <span class="o">=</span> <span class="p">[]</span>
741	<span class="k">for</span> <span class="n">par</span> <span class="ow">in</span> <span class="n">cpar</span><span class="p">:</span>
742	<span class="k">for</span> <span class="n">k</span><span class="p">,</span> <span class="n">v</span> <span class="ow">in</span> <span class="n">table</span><span class="o">.</span><span class="n">items</span><span class="p">():</span>
743	<span class="k">if</span> <span class="n">par</span><span class="o">.</span><span class="n">isdigit</span><span class="p">():</span>
744	<span class="n">par</span> <span class="o">=</span> <span class="nb">str</span><span class="p">(</span><span class="nb">int</span><span class="p">(</span><span class="n">par</span><span class="p">))</span><span class="o">.</span><span class="n">zfill</span><span class="p">(</span><span class="mi">3</span><span class="p">)</span>
745	<span class="k">if</span> <span class="n">par</span> <span class="o">==</span> <span class="n">k</span> <span class="ow">or</span> <span class="n">par</span> <span class="o">==</span> <span class="n">v</span><span class="p">:</span>
746	<span class="n">spar</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">k</span> <span class="o">+</span> <span class="s1">'.128'</span><span class="p">)</span>
747	<span class="k">break</span>
748	<span class="k">else</span><span class="p">:</span>
749	<span class="nb">print</span><span class="p">(</span><span class="s1">'</span><span class="se">\n\n\t\t</span><span class="s1">Warning: par '</span> <span class="o">+</span> <span class="n">par</span> <span class="o">+</span> <span class="s1">' not found in table 128</span><span class="se">\n\n</span><span class="s1">'</span><span class="p">)</span>
750
751	<span class="k">return</span> <span class="s1">'/'</span><span class="o">.</span><span class="n">join</span><span class="p">(</span><span class="n">spar</span><span class="p">)</span></div>
752
753	<div class="viewcode-block" id="get_list_as_string"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.get_list_as_string">[docs]</a><span class="k">def</span> <span class="nf">get_list_as_string</span><span class="p">(</span><span class="n">list_obj</span><span class="p">,</span> <span class="n">concatenate_sign</span><span class="o">=</span><span class="s1">', '</span><span class="p">):</span>
754	<span class="sd">'''Converts a list of arbitrary content into a single string.</span>
755
756	<span class="sd"> Parameters</span>
757	<span class="sd"> ----------</span>
758	<span class="sd"> list_obj : list of *</span>
759	<span class="sd"> A list with arbitrary content.</span>
760
761	<span class="sd"> concatenate_sign : str, optional</span>
762	<span class="sd"> A string which is used to concatenate the single</span>
763	<span class="sd"> list elements. Default value is ", ".</span>
764
765	<span class="sd"> Return</span>
766	<span class="sd"> ------</span>
767	<span class="sd"> str_of_list : str</span>
768	<span class="sd"> The content of the list as a single string.</span>
769	<span class="sd"> '''</span>
770
771	<span class="k">if</span> <span class="ow">not</span> <span class="nb">isinstance</span><span class="p">(</span><span class="n">list_obj</span><span class="p">,</span> <span class="nb">list</span><span class="p">):</span>
772	<span class="n">list_obj</span> <span class="o">=</span> <span class="nb">list</span><span class="p">(</span><span class="n">list_obj</span><span class="p">)</span>
773	<span class="n">str_of_list</span> <span class="o">=</span> <span class="n">concatenate_sign</span><span class="o">.</span><span class="n">join</span><span class="p">(</span><span class="nb">str</span><span class="p">(</span><span class="n">l</span><span class="p">)</span> <span class="k">for</span> <span class="n">l</span> <span class="ow">in</span> <span class="n">list_obj</span><span class="p">)</span>
774
775	<span class="k">return</span> <span class="n">str_of_list</span></div>
776
777	<div class="viewcode-block" id="make_dir"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.make_dir">[docs]</a><span class="k">def</span> <span class="nf">make_dir</span><span class="p">(</span><span class="n">directory</span><span class="p">):</span>
778	<span class="sd">'''Creates a directory.</span>
779
780	<span class="sd"> It gives a warning if the directory already exists and skips process.</span>
781	<span class="sd"> The program stops only if there is another problem.</span>
782
783	<span class="sd"> Parameters</span>
784	<span class="sd"> ----------</span>
785	<span class="sd"> directory : str</span>
786	<span class="sd"> The path to directory which should be created.</span>
787
788	<span class="sd"> Return</span>
789	<span class="sd"> ------</span>
790
791	<span class="sd"> '''</span>
792	<span class="k">try</span><span class="p">:</span>
793	<span class="n">os</span><span class="o">.</span><span class="n">makedirs</span><span class="p">(</span><span class="n">directory</span><span class="p">)</span>
794	<span class="k">except</span> <span class="ne">OSError</span> <span class="k">as</span> <span class="n">e</span><span class="p">:</span>
795	<span class="c1"># errno.EEXIST = directory already exists</span>
796	<span class="k">if</span> <span class="n">e</span><span class="o">.</span><span class="n">errno</span> <span class="o">==</span> <span class="n">errno</span><span class="o">.</span><span class="n">EEXIST</span><span class="p">:</span>
797	<span class="nb">print</span><span class="p">(</span><span class="s1">'WARNING: Directory </span><span class="si">{0}</span><span class="s1"> already exists!'</span><span class="o">.</span><span class="n">format</span><span class="p">(</span><span class="n">directory</span><span class="p">))</span>
798	<span class="k">else</span><span class="p">:</span>
799	<span class="k">raise</span> <span class="c1"># re-raise exception if a different error occured</span>
800
801	<span class="k">return</span></div>
802
803	<div class="viewcode-block" id="put_file_to_ecserver"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.put_file_to_ecserver">[docs]</a><span class="k">def</span> <span class="nf">put_file_to_ecserver</span><span class="p">(</span><span class="n">ecd</span><span class="p">,</span> <span class="n">filename</span><span class="p">,</span> <span class="n">target</span><span class="p">,</span> <span class="n">ecuid</span><span class="p">,</span> <span class="n">ecgid</span><span class="p">):</span>
804	<span class="sd">'''Uses the ecaccess-file-put command to send a file to the ECMWF servers.</span>
805
806	<span class="sd"> Note</span>
807	<span class="sd"> ----</span>
808	<span class="sd"> The return value is just for testing reasons. It does not have</span>
809	<span class="sd"> to be used from the calling function since the whole error handling</span>
810	<span class="sd"> is done in here.</span>
811
812	<span class="sd"> Parameters</span>
813	<span class="sd"> ----------</span>
814	<span class="sd"> ecd : str</span>
815	<span class="sd"> The path were the file is stored.</span>
816
817	<span class="sd"> filename : str</span>
818	<span class="sd"> The name of the file to send to the ECMWF server.</span>
819
820	<span class="sd"> target : str</span>
821	<span class="sd"> The target queue where the file should be sent to.</span>
822
823	<span class="sd"> ecuid : str</span>
824	<span class="sd"> The user id on ECMWF server.</span>
825
826	<span class="sd"> ecgid : str</span>
827	<span class="sd"> The group id on ECMWF server.</span>
828
829	<span class="sd"> Return</span>
830	<span class="sd"> ------</span>
831
832	<span class="sd"> '''</span>
833
834	<span class="k">try</span><span class="p">:</span>
835	<span class="n">subprocess</span><span class="o">.</span><span class="n">check_output</span><span class="p">([</span><span class="s1">'ecaccess-file-put'</span><span class="p">,</span>
836	<span class="n">ecd</span> <span class="o">+</span> <span class="s1">'/'</span> <span class="o">+</span> <span class="n">filename</span><span class="p">,</span>
837	<span class="n">target</span> <span class="o">+</span> <span class="s1">':/home/ms/'</span> <span class="o">+</span>
838	<span class="n">ecgid</span> <span class="o">+</span> <span class="s1">'/'</span> <span class="o">+</span> <span class="n">ecuid</span> <span class="o">+</span>
839	<span class="s1">'/'</span> <span class="o">+</span> <span class="n">filename</span><span class="p">],</span>
840	<span class="n">stderr</span><span class="o">=</span><span class="n">subprocess</span><span class="o">.</span><span class="n">STDOUT</span><span class="p">)</span>
841	<span class="k">except</span> <span class="n">subprocess</span><span class="o">.</span><span class="n">CalledProcessError</span> <span class="k">as</span> <span class="n">e</span><span class="p">:</span>
842	<span class="nb">print</span><span class="p">(</span><span class="s1">'... ERROR CODE: '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">e</span><span class="o">.</span><span class="n">returncode</span><span class="p">))</span>
843	<span class="nb">print</span><span class="p">(</span><span class="s1">'... ERROR MESSAGE:</span><span class="se">\n</span><span class="s1"> </span><span class="se">\t</span><span class="s1"> '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">e</span><span class="p">))</span>
844
845	<span class="nb">print</span><span class="p">(</span><span class="s1">'</span><span class="se">\n</span><span class="s1">... Do you have a valid ecaccess certification key?'</span><span class="p">)</span>
846	<span class="n">sys</span><span class="o">.</span><span class="n">exit</span><span class="p">(</span><span class="s1">'... ECACCESS-FILE-PUT FAILED!'</span><span class="p">)</span>
847	<span class="k">except</span> <span class="ne">OSError</span> <span class="k">as</span> <span class="n">e</span><span class="p">:</span>
848	<span class="nb">print</span><span class="p">(</span><span class="s1">'... ERROR CODE: '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">e</span><span class="o">.</span><span class="n">errno</span><span class="p">))</span>
849	<span class="nb">print</span><span class="p">(</span><span class="s1">'... ERROR MESSAGE:</span><span class="se">\n</span><span class="s1"> </span><span class="se">\t</span><span class="s1"> '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">e</span><span class="o">.</span><span class="n">strerror</span><span class="p">))</span>
850
851	<span class="nb">print</span><span class="p">(</span><span class="s1">'</span><span class="se">\n</span><span class="s1">... Most likely the ECACCESS library is not available!'</span><span class="p">)</span>
852	<span class="n">sys</span><span class="o">.</span><span class="n">exit</span><span class="p">(</span><span class="s1">'... ECACCESS-FILE-PUT FAILED!'</span><span class="p">)</span>
853
854	<span class="k">return</span></div>
855
856	<div class="viewcode-block" id="submit_job_to_ecserver"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.submit_job_to_ecserver">[docs]</a><span class="k">def</span> <span class="nf">submit_job_to_ecserver</span><span class="p">(</span><span class="n">target</span><span class="p">,</span> <span class="n">jobname</span><span class="p">):</span>
857	<span class="sd">'''Uses ecaccess-job-submit command to submit a job to the ECMWF server.</span>
858
859	<span class="sd"> Note</span>
860	<span class="sd"> ----</span>
861	<span class="sd"> The return value is just for testing reasons. It does not have</span>
862	<span class="sd"> to be used from the calling function since the whole error handling</span>
863	<span class="sd"> is done in here.</span>
864
865	<span class="sd"> Parameters</span>
866	<span class="sd"> ----------</span>
867	<span class="sd"> target : str</span>
868	<span class="sd"> The target where the file should be sent to, e.g. the queue.</span>
869
870	<span class="sd"> jobname : str</span>
871	<span class="sd"> The name of the jobfile to be submitted to the ECMWF server.</span>
872
873	<span class="sd"> Return</span>
874	<span class="sd"> ------</span>
875	<span class="sd"> job_id : int</span>
876	<span class="sd"> The id number of the job as a reference at the ecmwf server.</span>
877	<span class="sd"> '''</span>
878
879	<span class="k">try</span><span class="p">:</span>
880	<span class="n">job_id</span> <span class="o">=</span> <span class="n">subprocess</span><span class="o">.</span><span class="n">check_output</span><span class="p">([</span><span class="s1">'ecaccess-job-submit'</span><span class="p">,</span> <span class="s1">'-queueName'</span><span class="p">,</span>
881	<span class="n">target</span><span class="p">,</span> <span class="n">jobname</span><span class="p">])</span>
882
883	<span class="k">except</span> <span class="n">subprocess</span><span class="o">.</span><span class="n">CalledProcessError</span> <span class="k">as</span> <span class="n">e</span><span class="p">:</span>
884	<span class="nb">print</span><span class="p">(</span><span class="s1">'... ERROR CODE: '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">e</span><span class="o">.</span><span class="n">returncode</span><span class="p">))</span>
885	<span class="nb">print</span><span class="p">(</span><span class="s1">'... ERROR MESSAGE:</span><span class="se">\n</span><span class="s1"> </span><span class="se">\t</span><span class="s1"> '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">e</span><span class="p">))</span>
886
887	<span class="nb">print</span><span class="p">(</span><span class="s1">'</span><span class="se">\n</span><span class="s1">... Do you have a valid ecaccess certification key?'</span><span class="p">)</span>
888	<span class="n">sys</span><span class="o">.</span><span class="n">exit</span><span class="p">(</span><span class="s1">'... ECACCESS-JOB-SUBMIT FAILED!'</span><span class="p">)</span>
889	<span class="k">except</span> <span class="ne">OSError</span> <span class="k">as</span> <span class="n">e</span><span class="p">:</span>
890	<span class="nb">print</span><span class="p">(</span><span class="s1">'... ERROR CODE: '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">e</span><span class="o">.</span><span class="n">errno</span><span class="p">))</span>
891	<span class="nb">print</span><span class="p">(</span><span class="s1">'... ERROR MESSAGE:</span><span class="se">\n</span><span class="s1"> </span><span class="se">\t</span><span class="s1"> '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">e</span><span class="o">.</span><span class="n">strerror</span><span class="p">))</span>
892
893	<span class="nb">print</span><span class="p">(</span><span class="s1">'</span><span class="se">\n</span><span class="s1">... Most likely the ECACCESS library is not available!'</span><span class="p">)</span>
894	<span class="n">sys</span><span class="o">.</span><span class="n">exit</span><span class="p">(</span><span class="s1">'... ECACCESS-JOB-SUBMIT FAILED!'</span><span class="p">)</span>
895
896	<span class="k">return</span> <span class="n">job_id</span><span class="o">.</span><span class="n">decode</span><span class="p">()</span></div>
897
898
899	<div class="viewcode-block" id="get_informations"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.get_informations">[docs]</a><span class="k">def</span> <span class="nf">get_informations</span><span class="p">(</span><span class="n">filename</span><span class="p">):</span>
900	<span class="sd">'''Gets basic information from an example grib file.</span>
901
902	<span class="sd"> These information are important for later use and the</span>
903	<span class="sd"> initialization of numpy arrays for data storing.</span>
904
905	<span class="sd"> Parameters</span>
906	<span class="sd"> ----------</span>
907	<span class="sd"> filename : str</span>
908	<span class="sd"> Name of the file which will be opened to extract basic information.</span>
909
910	<span class="sd"> Return</span>
911	<span class="sd"> ------</span>
912	<span class="sd"> data : dict</span>
913	<span class="sd"> Contains basic informations of the ECMWF grib files, e.g.</span>
914	<span class="sd"> 'Ni', 'Nj', 'latitudeOfFirstGridPointInDegrees',</span>
915	<span class="sd"> 'longitudeOfFirstGridPointInDegrees', 'latitudeOfLastGridPointInDegrees',</span>
916	<span class="sd"> 'longitudeOfLastGridPointInDegrees', 'jDirectionIncrementInDegrees',</span>
917	<span class="sd"> 'iDirectionIncrementInDegrees', 'missingValue'</span>
918	<span class="sd"> '''</span>
919	<span class="kn">from</span> <span class="nn">eccodes</span> <span class="k">import</span> <span class="n">codes_grib_new_from_file</span><span class="p">,</span> <span class="n">codes_get</span><span class="p">,</span> <span class="n">codes_release</span>
920
921	<span class="n">data</span> <span class="o">=</span> <span class="p">{}</span>
922
923	<span class="c1"># --- open file ---</span>
924	<span class="nb">print</span><span class="p">(</span><span class="s2">"Opening file for getting information data --- </span><span class="si">%s</span><span class="s2">"</span> <span class="o">%</span> <span class="n">filename</span><span class="p">)</span>
925	<span class="k">with</span> <span class="nb">open</span><span class="p">(</span><span class="n">filename</span><span class="p">)</span> <span class="k">as</span> <span class="n">f</span><span class="p">:</span>
926	<span class="c1"># load first message from file</span>
927	<span class="n">gid</span> <span class="o">=</span> <span class="n">codes_grib_new_from_file</span><span class="p">(</span><span class="n">f</span><span class="p">)</span>
928
929	<span class="c1"># information needed from grib message</span>
930	<span class="n">keys</span> <span class="o">=</span> <span class="p">[</span>
931	<span class="s1">'Ni'</span><span class="p">,</span>
932	<span class="s1">'Nj'</span><span class="p">,</span>
933	<span class="s1">'latitudeOfFirstGridPointInDegrees'</span><span class="p">,</span>
934	<span class="s1">'longitudeOfFirstGridPointInDegrees'</span><span class="p">,</span>
935	<span class="s1">'latitudeOfLastGridPointInDegrees'</span><span class="p">,</span>
936	<span class="s1">'longitudeOfLastGridPointInDegrees'</span><span class="p">,</span>
937	<span class="s1">'jDirectionIncrementInDegrees'</span><span class="p">,</span>
938	<span class="s1">'iDirectionIncrementInDegrees'</span><span class="p">,</span>
939	<span class="s1">'missingValue'</span><span class="p">,</span>
940	<span class="p">]</span>
941
942	<span class="nb">print</span><span class="p">(</span><span class="s1">'</span><span class="se">\n</span><span class="s1">Informations are: '</span><span class="p">)</span>
943	<span class="k">for</span> <span class="n">key</span> <span class="ow">in</span> <span class="n">keys</span><span class="p">:</span>
944	<span class="c1"># Get the value of the key in a grib message.</span>
945	<span class="n">data</span><span class="p">[</span><span class="n">key</span><span class="p">]</span> <span class="o">=</span> <span class="n">codes_get</span><span class="p">(</span><span class="n">gid</span><span class="p">,</span><span class="n">key</span><span class="p">)</span>
946	<span class="nb">print</span><span class="p">(</span><span class="s2">"</span><span class="si">%s</span><span class="s2"> = </span><span class="si">%s</span><span class="s2">"</span> <span class="o">%</span> <span class="p">(</span><span class="n">key</span><span class="p">,</span><span class="n">data</span><span class="p">[</span><span class="n">key</span><span class="p">]))</span>
947
948	<span class="c1"># Free the memory for the message referred as gribid.</span>
949	<span class="n">codes_release</span><span class="p">(</span><span class="n">gid</span><span class="p">)</span>
950
951	<span class="k">return</span> <span class="n">data</span></div>
952
953
954	<div class="viewcode-block" id="get_dimensions"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.get_dimensions">[docs]</a><span class="k">def</span> <span class="nf">get_dimensions</span><span class="p">(</span><span class="n">info</span><span class="p">,</span> <span class="n">purefc</span><span class="p">,</span> <span class="n">dtime</span><span class="p">,</span> <span class="n">index_vals</span><span class="p">,</span> <span class="n">start_date</span><span class="p">,</span> <span class="n">end_date</span><span class="p">):</span>
955	<span class="sd">'''This function specifies the correct dimensions for x, y and t.</span>
956
957	<span class="sd"> Parameters</span>
958	<span class="sd"> ----------</span>
959	<span class="sd"> info : dict</span>
960	<span class="sd"> Contains basic informations of the ECMWF grib files, e.g.</span>
961	<span class="sd"> 'Ni', 'Nj', 'latitudeOfFirstGridPointInDegrees',</span>
962	<span class="sd"> 'longitudeOfFirstGridPointInDegrees', 'latitudeOfLastGridPointInDegrees',</span>
963	<span class="sd"> 'longitudeOfLastGridPointInDegrees', 'jDirectionIncrementInDegrees',</span>
964	<span class="sd"> 'iDirectionIncrementInDegrees', 'missingValue'</span>
965
966	<span class="sd"> purefc : int</span>
967	<span class="sd"> Switch for definition of pure forecast mode or not.</span>
968
969	<span class="sd"> dtime : str</span>
970	<span class="sd"> Time step in hours.</span>
971
972	<span class="sd"> index_vals : list of list of str</span>
973	<span class="sd"> Contains the values from the keys used for a distinct selection</span>
974	<span class="sd"> of grib messages in processing the grib files.</span>
975	<span class="sd"> Content looks like e.g.:</span>
976	<span class="sd"> index_vals[0]: ('20171106', '20171107', '20171108') ; date</span>
977	<span class="sd"> index_vals[1]: ('0', '1200', '1800', '600') ; time</span>
978	<span class="sd"> index_vals[2]: ('0', '12', '3', '6', '9') ; stepRange</span>
979
980	<span class="sd"> start_date : str</span>
981	<span class="sd"> The start date of the retrieval job.</span>
982
983	<span class="sd"> end_date : str</span>
984	<span class="sd"> The end date of the retrieval job.</span>
985
986	<span class="sd"> Return</span>
987	<span class="sd"> ------</span>
988	<span class="sd"> (ix, jy, it) : tuple of int</span>
989	<span class="sd"> Dimension in x-direction, y-direction and in time.</span>
990	<span class="sd"> '''</span>
991
992	<span class="n">ix</span> <span class="o">=</span> <span class="n">info</span><span class="p">[</span><span class="s1">'Ni'</span><span class="p">]</span>
993
994	<span class="n">jy</span> <span class="o">=</span> <span class="n">info</span><span class="p">[</span><span class="s1">'Nj'</span><span class="p">]</span>
995
996	<span class="k">if</span> <span class="ow">not</span> <span class="n">purefc</span><span class="p">:</span>
997	<span class="n">it</span> <span class="o">=</span> <span class="p">((</span><span class="n">end_date</span> <span class="o">-</span> <span class="n">start_date</span><span class="p">)</span><span class="o">.</span><span class="n">days</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)</span> <span class="o">*</span> <span class="mi">24</span> <span class="o">//</span> <span class="nb">int</span><span class="p">(</span><span class="n">dtime</span><span class="p">)</span>
998	<span class="k">else</span><span class="p">:</span>
999	<span class="c1"># #no of step * #no of times * #no of days</span>
1000	<span class="n">it</span> <span class="o">=</span> <span class="nb">len</span><span class="p">(</span><span class="n">index_vals</span><span class="p">[</span><span class="mi">2</span><span class="p">])</span> <span class="o"></span> <span class="nb">len</span><span class="p">(</span><span class="n">index_vals</span><span class="p">[</span><span class="mi">1</span><span class="p">])</span> <span class="o"></span> <span class="nb">len</span><span class="p">(</span><span class="n">index_vals</span><span class="p">[</span><span class="mi">0</span><span class="p">])</span>
1001
1002	<span class="k">return</span> <span class="p">(</span><span class="n">ix</span><span class="p">,</span> <span class="n">jy</span><span class="p">,</span> <span class="n">it</span><span class="p">)</span></div>
1003
1004
1005	<div class="viewcode-block" id="execute_subprocess"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.execute_subprocess">[docs]</a><span class="k">def</span> <span class="nf">execute_subprocess</span><span class="p">(</span><span class="n">cmd_list</span><span class="p">,</span> <span class="n">error_msg</span><span class="o">=</span><span class="s1">'SUBPROCESS FAILED!'</span><span class="p">):</span>
1006	<span class="sd">'''Executes a command line instruction via a subprocess.</span>
1007
1008	<span class="sd"> Error handling is done if an error occures.</span>
1009
1010	<span class="sd"> Parameters</span>
1011	<span class="sd"> ----------</span>
1012	<span class="sd"> cmd_list : list of str</span>
1013	<span class="sd"> A list of the components for the command line execution. Each</span>
1014	<span class="sd"> list entry is a single part of the command which is seperated from</span>
1015	<span class="sd"> the rest by a blank space.</span>
1016	<span class="sd"> E.g. ['mv', file1, file2]</span>
1017
1018	<span class="sd"> Return</span>
1019	<span class="sd"> ------</span>
1020	<span class="sd"> error_msg : str, optional</span>
1021	<span class="sd"> The possible error message if the subprocess failed.</span>
1022	<span class="sd"> By default it will just tell "SUBPROCESS FAILED!".</span>
1023	<span class="sd"> '''</span>
1024
1025	<span class="k">try</span><span class="p">:</span>
1026	<span class="n">subprocess</span><span class="o">.</span><span class="n">check_call</span><span class="p">(</span><span class="n">cmd_list</span><span class="p">)</span>
1027	<span class="k">except</span> <span class="n">subprocess</span><span class="o">.</span><span class="n">CalledProcessError</span> <span class="k">as</span> <span class="n">e</span><span class="p">:</span>
1028	<span class="nb">print</span><span class="p">(</span><span class="s1">'... ERROR CODE: '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">e</span><span class="o">.</span><span class="n">returncode</span><span class="p">))</span>
1029	<span class="nb">print</span><span class="p">(</span><span class="s1">'... ERROR MESSAGE:</span><span class="se">\n</span><span class="s1"> </span><span class="se">\t</span><span class="s1"> '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">e</span><span class="p">))</span>
1030
1031	<span class="n">sys</span><span class="o">.</span><span class="n">exit</span><span class="p">(</span><span class="s1">'... '</span> <span class="o">+</span> <span class="n">error_msg</span><span class="p">)</span>
1032	<span class="k">except</span> <span class="ne">OSError</span> <span class="k">as</span> <span class="n">e</span><span class="p">:</span>
1033	<span class="nb">print</span><span class="p">(</span><span class="s1">'... ERROR CODE: '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">e</span><span class="o">.</span><span class="n">errno</span><span class="p">))</span>
1034	<span class="nb">print</span><span class="p">(</span><span class="s1">'... ERROR MESSAGE:</span><span class="se">\n</span><span class="s1"> </span><span class="se">\t</span><span class="s1"> '</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">e</span><span class="o">.</span><span class="n">strerror</span><span class="p">))</span>
1035
1036	<span class="n">sys</span><span class="o">.</span><span class="n">exit</span><span class="p">(</span><span class="s1">'... '</span> <span class="o">+</span> <span class="n">error_msg</span><span class="p">)</span>
1037
1038	<span class="k">return</span></div>
1039
1040
1041	<div class="viewcode-block" id="generate_retrieval_period_boundary"><a class="viewcode-back" href="../Documentation/Api/api_python.html#tools.generate_retrieval_period_boundary">[docs]</a><span class="k">def</span> <span class="nf">generate_retrieval_period_boundary</span><span class="p">(</span><span class="n">c</span><span class="p">):</span>
1042	<span class="sd">'''Generates retrieval period boundary datetimes from CONTROL information.</span>
1043
1044	<span class="sd"> Parameters</span>
1045	<span class="sd"> ----------</span>
1046	<span class="sd"> c : ControlFile</span>
1047	<span class="sd"> Contains all the parameters of CONTROL file and</span>
1048	<span class="sd"> command line.</span>
1049
1050	<span class="sd"> Return</span>
1051	<span class="sd"> ------</span>
1052	<span class="sd"> start_period : datetime</span>
1053	<span class="sd"> The first timestamp of the actual retrieval period disregarding</span>
1054	<span class="sd"> the temporary times which were used for processing reasons.</span>
1055
1056	<span class="sd"> end_period : datetime</span>
1057	<span class="sd"> The last timestamp of the actual retrieval period disregarding</span>
1058	<span class="sd"> the temporary times which were used for processing reasons.</span>
1059	<span class="sd"> '''</span>
1060	<span class="c1"># generate start and end timestamp of the retrieval period</span>
1061	<span class="n">start_period</span> <span class="o">=</span> <span class="n">datetime</span><span class="o">.</span><span class="n">strptime</span><span class="p">(</span><span class="n">c</span><span class="o">.</span><span class="n">start_date</span> <span class="o">+</span> <span class="n">c</span><span class="o">.</span><span class="n">time</span><span class="p">[</span><span class="mi">0</span><span class="p">],</span> <span class="s1">'%Y%m</span><span class="si">%d</span><span class="s1">%H'</span><span class="p">)</span>
1062	<span class="n">start_period</span> <span class="o">=</span> <span class="n">start_period</span> <span class="o">+</span> <span class="n">timedelta</span><span class="p">(</span><span class="n">hours</span><span class="o">=</span><span class="nb">int</span><span class="p">(</span><span class="n">c</span><span class="o">.</span><span class="n">step</span><span class="p">[</span><span class="mi">0</span><span class="p">]))</span>
1063	<span class="n">end_period</span> <span class="o">=</span> <span class="n">datetime</span><span class="o">.</span><span class="n">strptime</span><span class="p">(</span><span class="n">c</span><span class="o">.</span><span class="n">end_date</span> <span class="o">+</span> <span class="n">c</span><span class="o">.</span><span class="n">time</span><span class="p">[</span><span class="o">-</span><span class="mi">1</span><span class="p">],</span> <span class="s1">'%Y%m</span><span class="si">%d</span><span class="s1">%H'</span><span class="p">)</span>
1064	<span class="n">end_period</span> <span class="o">=</span> <span class="n">end_period</span> <span class="o">+</span> <span class="n">timedelta</span><span class="p">(</span><span class="n">hours</span><span class="o">=</span><span class="nb">int</span><span class="p">(</span><span class="n">c</span><span class="o">.</span><span class="n">step</span><span class="p">[</span><span class="o">-</span><span class="mi">1</span><span class="p">]))</span>
1065
1066
1067	<span class="k">return</span> <span class="n">start_period</span><span class="p">,</span> <span class="n">end_period</span></div>
1068	</pre></div>
1069
1070	</div>
1071
1072	</div>
1073	<footer>
1074
1075
1076	<hr/>
1077
1078	<div role="contentinfo">
1079	<p>
1080	© Copyright 2019, Anne Philipp and Leopold Haimberger
1081
1082	</p>
1083	</div>
1084	Built with <a href="http://sphinx-doc.org/">Sphinx</a> using a <a href="https://github.com/rtfd/sphinx_rtd_theme">theme</a> provided by <a href="https://readthedocs.org">Read the Docs</a>.
1085
1086	</footer>
1087
1088	</div>
1089	</div>
1090
1091	</section>
1092
1093	</div>
1094
1095
1096
1097	<script type="text/javascript">
1098	jQuery(function () {
1099	SphinxRtdTheme.Navigation.enable(true);
1100	});
1101	</script>
1102
1103
1104
1105
1106
1107
1108	</body>
1109	</html>

Note: See TracBrowser for help on using the repository browser.

Download in other formats: