source: flex_extract.git/documentation/FORD/V5/Doc_original/sourcefile/grphreal.f.html @ ee06999

ctbtodev
Last change on this file since ee06999 was ee06999, checked in by pesei <petra seibert @…>, 5 years ago

Initial version of Fortran documentation

with FORD v5
  • Property mode set to 100644
File size: 54.0 KB
Line 
1<!DOCTYPE html>
2<html lang="en">
3  <head>
4    <meta charset="utf-8">
5    <meta http-equiv="X-UA-Compatible" content="IE=edge">
6    <meta name="viewport" content="width=device-width, initial-scale=1">
7         
8         <meta name="description" content="Calculation of vertical velocity for FLEXPART">
9   
10    <meta name="author" content="Leopold Haimberger" >
11    <link rel="icon" href="../favicon.png">
12
13    <title>grphreal.f &ndash; Flex_extract: Calculation of etadot</title>
14
15    <link href="../css/bootstrap.min.css" rel="stylesheet">
16    <link href="../css/pygments.css" rel="stylesheet">
17    <link href="../css/font-awesome.min.css" rel="stylesheet">
18    <link href="../css/local.css" rel="stylesheet">
19   
20
21    <!-- HTML5 shim and Respond.js for IE8 support of HTML5 elements and media queries -->
22    <!--[if lt IE 9]>
23      <script src="https://oss.maxcdn.com/html5shiv/3.7.2/html5shiv.min.js"></script>
24      <script src="https://oss.maxcdn.com/respond/1.4.2/respond.min.js"></script>
25    <![endif]-->
26   
27    <script src="../js/jquery-2.1.3.min.js"></script>
28    <script src="../js/svg-pan-zoom.min.js"></script>
29
30   
31    <script src="../tipuesearch/tipuesearch_content.js"></script>
32    <link  href="../tipuesearch/tipuesearch.css" rel="stylesheet">
33    <script src="../tipuesearch/tipuesearch_set.js"></script>
34    <script src="../tipuesearch/tipuesearch.js"></script>
35   
36
37  </head>
38
39  <body>
40
41    <!-- Fixed navbar -->
42    <nav class="navbar navbar-inverse navbar-fixed-top">
43      <div class="container">
44        <div class="navbar-header">
45          <button type="button" class="navbar-toggle collapsed" data-toggle="collapse" data-target="#navbar" aria-expanded="false" aria-controls="navbar">
46            <span class="sr-only">Toggle navigation</span>
47            <span class="icon-bar"></span>
48            <span class="icon-bar"></span>
49            <span class="icon-bar"></span>
50          </button>
51          <a class="navbar-brand" href="../index.html">Flex_extract: Calculation of etadot <small>7.1</small></a>
52        </div>
53        <div id="navbar" class="navbar-collapse collapse">
54          <ul class="nav navbar-nav">
55                               
56            <li class="dropdown hidden-xs visible-sm visible-md hidden-lg">
57              <a href="#" class="dropdown-toggle"
58              data-toggle="dropdown" role="button"
59              aria-haspopup="true"
60                 aria-expanded="false">Contents <span class="caret"></span></a>
61              <ul class="dropdown-menu">
62             
63            <li><a href="../lists/files.html">Source Files</a></li>
64                               
65                               
66            <li><a href="../lists/modules.html">Modules</a></li>
67                               
68           
69                               
70            <li><a href="../lists/procedures.html">Procedures</a></li>
71                               
72                                                               
73                               
74            <li><a href="../program/preconvert.html">Program</a></li>
75           
76            </ul>
77            </li>
78
79<li class="visible-xs hidden-sm visible-lg"><a href="../lists/files.html">Source Files</a></li>
80
81
82<li class="visible-xs hidden-sm visible-lg"><a href="../lists/modules.html">Modules</a></li>
83
84
85
86<li class="visible-xs hidden-sm visible-lg"><a href="../lists/procedures.html">Procedures</a></li>
87
88                             
89
90<li class="visible-xs hidden-sm visible-lg"><a href="../program/preconvert.html">Program</a></li>
91
92          </ul>
93       
94        <form action="../search.html" class="navbar-form navbar-right" role="search">
95        <div class="form-group">
96          <input type="text" class="form-control" placeholder="Search" name="q" id="tipue_search_input" autocomplete="off" required>
97        </div>
98<!--
99        <button type="submit" class="btn btn-default">Submit</button>
100-->
101        </form>
102       
103        </div><!--/.nav-collapse -->
104      </div>
105    </nav>
106
107    <div class="container">
108   
109 
110  <div class="row">
111    <h1>grphreal.f
112    <small>Source File</small>
113   
114    </h1>
115   
116<div class="row">
117  <div class="col-lg-12">
118<div class="well well-sm">
119  <ul class="list-inline" style="margin-bottom:0px;display:inline">
120     
121     
122     
123     
124   
125   
126     <li><i class="fa fa-list-ol"></i>
127       <a data-toggle="tooltip"
128          data-placement="bottom" data-html="true"
129          title="10.2% of total for source files.">136 statements</a>
130     </li> 
131     
132     
133    <li><i class="fa fa-code"></i><a href="../src/grphreal.f"> Source File</a></li>
134     
135  </ul>
136  <ol class="breadcrumb in-well text-right">
137 
138     <li class="active">grphreal.f</li>
139  </ol>
140</div>
141</div>
142</div>
143<script>
144  $(function () {
145  $('[data-toggle="tooltip"]').tooltip()
146  })
147</script>
148
149  </div>
150  <div class="row">
151    <div class="col-md-3 hidden-xs hidden-sm visible-md visible-lg">
152   
153<div id="sidebar">
154 
155
156
157
158
159
160
161<div class="panel panel-primary">
162  <div class="panel-heading text-left"><h3 class="panel-title"><a data-toggle="collapse" href="#mods-0">Modules</a></h3></div>
163  <div id="mods-0" class="panel-collapse collapse">
164    <div class="list-group">
165     
166      <a class="list-group-item" href="../module/grtoph.html">GRTOPH</a>
167     
168    </div>
169  </div>
170</div>
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186<div class="panel panel-primary">
187  <div class="panel-heading text-left"><h3 class="panel-title">Source Code</h3></div>
188  <div class="list-group">
189    <a class="list-group-item" href="../sourcefile/grphreal.f.html#src">grphreal.f</a>
190  </div>
191</div>
192
193
194  <hr>
195 
196
197<div class="panel panel-default">
198  <div class="panel-heading text-left"><h3 class="panel-title"><a data-toggle="collapse" href="#allfiles-0">All Source Files</a></h3></div>
199  <div id="allfiles-0" class="panel-collapse collapse">
200    <div class="list-group">
201     
202      <a class="list-group-item" href="../sourcefile/grphreal.f.html">grphreal.f</a>
203     
204      <a class="list-group-item" href="../sourcefile/phgrreal.f.html">phgrreal.f</a>
205     
206      <a class="list-group-item" href="../sourcefile/preconvert.f90.html">preconvert.f90</a>
207     
208      <a class="list-group-item" href="../sourcefile/rwgrib2.f90.html">rwGRIB2.f90</a>
209     
210    </div>
211  </div>
212</div>
213
214
215</div> 
216
217    </div>
218    <div class="col-md-9" id='text'>
219   
220   
221    <h3>This File Depends On</h3>
222   
223                <div class="depgraph"><?xml version="1.0" encoding="UTF-8" standalone="no"?>
224<!DOCTYPE svg PUBLIC "-//W3C//DTD SVG 1.1//EN"
225 "http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd">
226<!-- Generated by graphviz version 2.38.0 (20140413.2041)
227 -->
228<!-- Title: sourcefile~~grphreal.f~~EfferentGraph Pages: 1 -->
229<svg id="sourcefilegrphrealfEfferentGraph" width="160pt" height="32pt"
230 viewBox="0.00 0.00 160.00 32.00" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
231<g id="sourcefile~~grphreal.f~~EfferentGraph" class="graph" transform="scale(1 1) rotate(0) translate(4 28)">
232<title>sourcefile~~grphreal.f~~EfferentGraph</title>
233<polygon fill="white" stroke="none" points="-4,4 -4,-28 156,-28 156,4 -4,4"/>
234<!-- sourcefile~grphreal.f -->
235<g id="sourcefile~~grphreal.f~~EfferentGraph_node1" class="node"><title>sourcefile~grphreal.f</title>
236<polygon fill="none" stroke="black" points="152,-24 94,-24 94,-0 152,-0 152,-24"/>
237<text text-anchor="middle" x="123" y="-9.6" font-family="Helvetica,sans-Serif" font-size="10.50">grphreal.f</text>
238</g>
239<!-- sourcefile~phgrreal.f -->
240<g id="sourcefile~~grphreal.f~~EfferentGraph_node2" class="node"><title>sourcefile~phgrreal.f</title>
241<g id="a_sourcefile~~grphreal.f~~EfferentGraph_node2"><a xlink:href="../sourcefile/phgrreal.f.html" xlink:title="phgrreal.f">
242<polygon fill="#f0ad4e" stroke="#f0ad4e" points="58,-24 3.55271e-15,-24 3.55271e-15,-0 58,-0 58,-24"/>
243<text text-anchor="middle" x="29" y="-9.6" font-family="Helvetica,sans-Serif" font-size="10.50" fill="white">phgrreal.f</text>
244</a>
245</g>
246</g>
247<!-- sourcefile~phgrreal.f&#45;&gt;sourcefile~grphreal.f -->
248<g id="sourcefile~~grphreal.f~~EfferentGraph_edge1" class="edge"><title>sourcefile~phgrreal.f&#45;&gt;sourcefile~grphreal.f</title>
249<path fill="none" stroke="#000000" stroke-dasharray="5,2" d="M58.1029,-12C66.1425,-12 75.0461,-12 83.5904,-12"/>
250<polygon fill="#000000" stroke="#000000" points="83.7522,-15.5001 93.7521,-12 83.7521,-8.5001 83.7522,-15.5001"/>
251</g>
252</g>
253</svg>
254</div>
255                <div><a type="button" class="graph-help" data-toggle="modal" href="#graph-help-text">Help</a></div>
256                <div class="modal fade" id="graph-help-text" tabindex="-1" role="dialog">
257                  <div class="modal-dialog modal-lg" role="document">
258                    <div class="modal-content">
259                      <div class="modal-header">
260                        <button type="button" class="close" data-dismiss="modal" aria-label="Close"><span aria-hidden="true">&times;</span></button>
261                        <h4 class="modal-title" id="-graph-help-label">Graph Key</h4>
262                      </div>
263                      <div class="modal-body">
264                       
265    <p>Nodes of different colours represent the following: </p>
266    <?xml version="1.0" encoding="UTF-8" standalone="no"?>
267<!DOCTYPE svg PUBLIC "-//W3C//DTD SVG 1.1//EN"
268 "http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd">
269<!-- Generated by graphviz version 2.38.0 (20140413.2041)
270 -->
271<!-- Title: Graph Key Pages: 1 -->
272<svg width="190pt" height="32pt"
273 viewBox="0.00 0.00 190.00 32.00" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
274<g id="graph0" class="graph" transform="scale(1 1) rotate(0) translate(4 28)">
275<title>Graph Key</title>
276<polygon fill="white" stroke="none" points="-4,4 -4,-28 186,-28 186,4 -4,4"/>
277<!-- Source File -->
278<g id="node1" class="node"><title>Source File</title>
279<polygon fill="#f0ad4e" stroke="#f0ad4e" points="67,-24 0,-24 0,-0 67,-0 67,-24"/>
280<text text-anchor="middle" x="33.5" y="-9.6" font-family="Helvetica,sans-Serif" font-size="10.50" fill="white">Source File</text>
281</g>
282<!-- This Page&#39;s Entity -->
283<g id="node2" class="node"><title>This Page&#39;s Entity</title>
284<polygon fill="none" stroke="black" points="182,-24 85,-24 85,-0 182,-0 182,-24"/>
285<text text-anchor="middle" x="133.5" y="-9.6" font-family="Helvetica,sans-Serif" font-size="10.50">This Page&#39;s Entity</text>
286</g>
287</g>
288</svg>
289
290   
291    <p>Solid arrows point from a file to a file which depends upon it. A file
292    is dependent upon another if the latter must be compiled before the former
293    can be.
294    </p>
295   
296                      </div>
297                    </div>
298                  </div>
299                </div>
300               
301   
302     
303      <h3>Files Dependent On This One</h3>
304     
305                <div class="depgraph"><?xml version="1.0" encoding="UTF-8" standalone="no"?>
306<!DOCTYPE svg PUBLIC "-//W3C//DTD SVG 1.1//EN"
307 "http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd">
308<!-- Generated by graphviz version 2.38.0 (20140413.2041)
309 -->
310<!-- Title: sourcefile~~grphreal.f~~AfferentGraph Pages: 1 -->
311<svg id="sourcefilegrphrealfAfferentGraph" width="184pt" height="32pt"
312 viewBox="0.00 0.00 184.00 32.00" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
313<g id="sourcefile~~grphreal.f~~AfferentGraph" class="graph" transform="scale(1 1) rotate(0) translate(4 28)">
314<title>sourcefile~~grphreal.f~~AfferentGraph</title>
315<polygon fill="white" stroke="none" points="-4,4 -4,-28 180,-28 180,4 -4,4"/>
316<!-- sourcefile~grphreal.f -->
317<g id="sourcefile~~grphreal.f~~AfferentGraph_node1" class="node"><title>sourcefile~grphreal.f</title>
318<polygon fill="none" stroke="black" points="58,-24 3.55271e-15,-24 3.55271e-15,-0 58,-0 58,-24"/>
319<text text-anchor="middle" x="29" y="-9.6" font-family="Helvetica,sans-Serif" font-size="10.50">grphreal.f</text>
320</g>
321<!-- sourcefile~preconvert.f90 -->
322<g id="sourcefile~~grphreal.f~~AfferentGraph_node2" class="node"><title>sourcefile~preconvert.f90</title>
323<g id="a_sourcefile~~grphreal.f~~AfferentGraph_node2"><a xlink:href="../sourcefile/preconvert.f90.html" xlink:title="preconvert.f90">
324<polygon fill="#f0ad4e" stroke="#f0ad4e" points="176,-24 94,-24 94,-0 176,-0 176,-24"/>
325<text text-anchor="middle" x="135" y="-9.6" font-family="Helvetica,sans-Serif" font-size="10.50" fill="white">preconvert.f90</text>
326</a>
327</g>
328</g>
329<!-- sourcefile~grphreal.f&#45;&gt;sourcefile~preconvert.f90 -->
330<g id="sourcefile~~grphreal.f~~AfferentGraph_edge1" class="edge"><title>sourcefile~grphreal.f&#45;&gt;sourcefile~preconvert.f90</title>
331<path fill="none" stroke="#000000" stroke-dasharray="5,2" d="M58.1439,-12C66.0131,-12 74.8129,-12 83.5541,-12"/>
332<polygon fill="#000000" stroke="#000000" points="83.7293,-15.5001 93.7293,-12 83.7292,-8.5001 83.7293,-15.5001"/>
333</g>
334</g>
335</svg>
336</div>
337                <div><a type="button" class="graph-help" data-toggle="modal" href="#graph-help-text">Help</a></div>
338                <div class="modal fade" id="graph-help-text" tabindex="-1" role="dialog">
339                  <div class="modal-dialog modal-lg" role="document">
340                    <div class="modal-content">
341                      <div class="modal-header">
342                        <button type="button" class="close" data-dismiss="modal" aria-label="Close"><span aria-hidden="true">&times;</span></button>
343                        <h4 class="modal-title" id="-graph-help-label">Graph Key</h4>
344                      </div>
345                      <div class="modal-body">
346                       
347    <p>Nodes of different colours represent the following: </p>
348    <?xml version="1.0" encoding="UTF-8" standalone="no"?>
349<!DOCTYPE svg PUBLIC "-//W3C//DTD SVG 1.1//EN"
350 "http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd">
351<!-- Generated by graphviz version 2.38.0 (20140413.2041)
352 -->
353<!-- Title: Graph Key Pages: 1 -->
354<svg width="190pt" height="32pt"
355 viewBox="0.00 0.00 190.00 32.00" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
356<g id="graph0" class="graph" transform="scale(1 1) rotate(0) translate(4 28)">
357<title>Graph Key</title>
358<polygon fill="white" stroke="none" points="-4,4 -4,-28 186,-28 186,4 -4,4"/>
359<!-- Source File -->
360<g id="node1" class="node"><title>Source File</title>
361<polygon fill="#f0ad4e" stroke="#f0ad4e" points="67,-24 0,-24 0,-0 67,-0 67,-24"/>
362<text text-anchor="middle" x="33.5" y="-9.6" font-family="Helvetica,sans-Serif" font-size="10.50" fill="white">Source File</text>
363</g>
364<!-- This Page&#39;s Entity -->
365<g id="node2" class="node"><title>This Page&#39;s Entity</title>
366<polygon fill="none" stroke="black" points="182,-24 85,-24 85,-0 182,-0 182,-24"/>
367<text text-anchor="middle" x="133.5" y="-9.6" font-family="Helvetica,sans-Serif" font-size="10.50">This Page&#39;s Entity</text>
368</g>
369</g>
370</svg>
371
372   
373    <p>Solid arrows point from a file to a file which depends upon it. A file
374    is dependent upon another if the latter must be compiled before the former
375    can be.
376    </p>
377   
378                      </div>
379                    </div>
380                  </div>
381                </div>
382               
383     
384      <br>
385
386    <section class="visible-xs visible-sm hidden-md">
387     
388
389
390
391
392
393
394<div class="panel panel-primary">
395  <div class="panel-heading text-left"><h3 class="panel-title"><a data-toggle="collapse" href="#mods-1">Modules</a></h3></div>
396  <div id="mods-1" class="panel-collapse collapse">
397    <div class="list-group">
398     
399      <a class="list-group-item" href="../module/grtoph.html">GRTOPH</a>
400     
401    </div>
402  </div>
403</div>
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419<div class="panel panel-primary">
420  <div class="panel-heading text-left"><h3 class="panel-title">Source Code</h3></div>
421  <div class="list-group">
422    <a class="list-group-item" href="../sourcefile/grphreal.f.html#src">grphreal.f</a>
423  </div>
424</div>
425
426
427    </section>
428    <br class="visible-xs visible-sm hidden-md">
429
430    <section>
431      <h2><span class="anchor" id="src"></span>Source Code</h2>
432    <div class="hl"><pre><span></span><a name="ln-1"></a><span class="nl">     </span> <span class="k">MODULE </span><span class="n">GRTOPH</span>
433<a name="ln-2"></a>
434<a name="ln-3"></a><span class="nl">     </span> <span class="k">USE </span><span class="n">PHTOGR</span>
435<a name="ln-4"></a>
436<a name="ln-5"></a><span class="nl">     </span> <span class="k">CONTAINS</span>
437<a name="ln-6"></a><span class="c">C</span>
438<a name="ln-7"></a><span class="nl">     </span> <span class="k">SUBROUTINE </span><span class="n">GRPH213</span><span class="p">(</span><span class="n">CXMN</span><span class="p">,</span><span class="n">FELD</span><span class="p">,</span><span class="n">WSAVE</span><span class="p">,</span><span class="n">IFAX</span><span class="p">,</span><span class="n">Z</span><span class="p">,</span><span class="n">W</span><span class="p">,</span><span class="n">MLAT</span><span class="p">,</span>
439<a name="ln-8"></a><span class="nl">     </span><span class="gs">*</span><span class="n">MNAUF</span><span class="p">,</span><span class="n">MAXL</span><span class="p">,</span><span class="n">MAXB</span><span class="p">,</span><span class="n">MLEVEL</span><span class="p">)</span>
440<a name="ln-9"></a>
441<a name="ln-10"></a><span class="c">C     DIE ROUTINE F]HRT EINE TRANSFORMATION EINER</span>
442<a name="ln-11"></a><span class="c">C     FELDVARIABLEN VOM PHASENRAUM IN  DEN PHYSIKALISCHEN</span>
443<a name="ln-12"></a><span class="c">C     RAUM AUF KUGELKOORDINATEN DURCH</span>
444<a name="ln-13"></a><span class="c">C</span>
445<a name="ln-14"></a><span class="c">C     CXMN   = SPEKTRALKOEFFIZIENTEN IN DER REIHENFOLGE</span>
446<a name="ln-15"></a><span class="c">C              CX00,CX01,CX11,CX02,....CXMNAUFMNAUF</span>
447<a name="ln-16"></a><span class="c">C                   CXM              = FOURIERKOEFFIZIENTEN - nur ein Hilfsfeld</span>
448<a name="ln-17"></a><span class="c">C     FELD   = FELD DER METEOROLOGISCHEN VARIABLEN</span>
449<a name="ln-18"></a><span class="c">C                   WSAVE  = Working Array fuer Fouriertransformation</span>
450<a name="ln-19"></a><span class="c">C     Z      = LEGENDREFUNKTIONSWERTE</span>
451<a name="ln-20"></a><span class="c">C</span>
452<a name="ln-21"></a><span class="c">C     MNAUF    ANZAHL DER FOURIERKOEFFIZIENTEN</span>
453<a name="ln-22"></a><span class="c">C     MAXL     ANZAHL DER FUER DAS GITTER BENUTZTEN LAENGEN</span>
454<a name="ln-23"></a><span class="c">C     MAXB     ANZAHL DER FUER DAS GITTER BENOETIGTEN BREITEN</span>
455<a name="ln-24"></a><span class="c">C     MLEVEL   ANZAHL DER LEVELS, DIE TRANSFORMIERT WERDEN</span>
456<a name="ln-25"></a><span class="c">C</span>
457<a name="ln-26"></a><span class="nl">     </span> <span class="k">IMPLICIT </span><span class="kt">REAL</span> <span class="p">(</span><span class="n">A</span><span class="o">-</span><span class="n">H</span><span class="p">,</span><span class="n">O</span><span class="o">-</span><span class="n">Z</span><span class="p">)</span>
458<a name="ln-27"></a>
459<a name="ln-28"></a>
460<a name="ln-29"></a><span class="c">C                   Anzahl der Gitterpunkte pro Breitenkreis des reduzierten </span>
461<a name="ln-30"></a><span class="c">C                   Gauss&#39;schen Gitters</span>
462<a name="ln-31"></a><span class="nl">   INTE</span><span class="gs">G</span><span class="n">ER</span> <span class="n">MLAT</span><span class="p">(</span><span class="n">MAXB</span><span class="p">),</span><span class="n">ISIZE</span><span class="p">,</span><span class="n">IFAX</span><span class="p">(</span><span class="mi">10</span><span class="p">,</span><span class="n">MAXB</span><span class="p">)</span>
463<a name="ln-32"></a>                   
464<a name="ln-33"></a><span class="c">C    FELD DER LEGENDREPOLYNOME FUER EINE BREITE</span>
465<a name="ln-34"></a><span class="nl">     </span> <span class="kt">REAL</span><span class="o">*</span><span class="mi">8</span> <span class="n">Z</span><span class="p">(</span><span class="n">MAXB</span><span class="o">/</span><span class="mi">2</span><span class="p">,</span><span class="mi">0</span><span class="p">:((</span><span class="n">MNAUF</span><span class="o">+</span><span class="mi">3</span><span class="p">)</span><span class="o">*</span><span class="p">(</span><span class="n">MNAUF</span><span class="o">+</span><span class="mi">4</span><span class="p">))</span><span class="o">/</span><span class="mi">2</span><span class="p">)</span>
466<a name="ln-35"></a><span class="nl">     </span> 
467<a name="ln-36"></a><span class="c">C      LOGICAL*1 USED(((216*217)/2+1)*160)</span>
468<a name="ln-37"></a>
469<a name="ln-38"></a><span class="nl">     </span> <span class="k">DIMENSION </span><span class="n">CXMN</span><span class="p">(</span><span class="mi">0</span><span class="p">:(</span><span class="n">MNAUF</span><span class="o">+</span><span class="mi">1</span><span class="p">)</span><span class="o">*</span><span class="p">(</span><span class="n">MNAUF</span><span class="o">+</span><span class="mi">2</span><span class="p">)</span><span class="o">-</span><span class="mi">1</span><span class="p">,</span><span class="n">MLEVEL</span><span class="p">)</span>
470<a name="ln-39"></a><span class="nl">     </span> <span class="kt">REAL </span><span class="n">FELD</span><span class="p">(</span><span class="n">MAXL</span><span class="p">,</span><span class="n">MLEVEL</span><span class="p">)</span>
471<a name="ln-40"></a><span class="nl">     </span> <span class="k">DIMENSION </span><span class="n">WSAVE</span><span class="p">(</span><span class="mi">8</span><span class="o">*</span><span class="n">MAXB</span><span class="o">+</span><span class="mi">15</span><span class="p">,</span><span class="n">MAXB</span><span class="o">/</span><span class="mi">2</span><span class="p">)</span>
472<a name="ln-41"></a><span class="nl">     </span> <span class="kt">REAL</span><span class="o">*</span><span class="mi">8</span> <span class="n">W</span><span class="p">(</span><span class="n">MAXB</span><span class="p">)</span> 
473<a name="ln-42"></a><span class="nl">     </span> <span class="k">DIMENSION </span><span class="n">IND</span><span class="p">(</span><span class="n">MAXB</span><span class="p">)</span>
474<a name="ln-43"></a>                   
475<a name="ln-44"></a>                   
476<a name="ln-45"></a><span class="nl">     </span>       <span class="n">IND</span><span class="p">(</span><span class="mi">1</span><span class="p">)</span><span class="o">=</span><span class="mi">0</span>
477<a name="ln-46"></a><span class="nl">     </span> <span class="k">DO </span><span class="mi">6</span> <span class="n">J</span><span class="o">=</span><span class="mi">2</span><span class="p">,</span><span class="n">MAXB</span><span class="o">/</span><span class="mi">2</span>
478<a name="ln-47"></a><span class="nl">    IND</span><span class="gs">(</span><span class="n">j</span><span class="p">)</span><span class="o">=</span><span class="n">IND</span><span class="p">(</span><span class="n">J</span><span class="o">-</span><span class="mi">1</span><span class="p">)</span><span class="o">+</span><span class="n">MLAT</span><span class="p">(</span><span class="n">J</span><span class="o">-</span><span class="mi">1</span><span class="p">)</span>
479<a name="ln-48"></a><span class="nl"> 6   </span> <span class="k">CONTINUE</span>
480<a name="ln-49"></a><span class="c">!$OMP PARALLEL DO SCHEDULE(DYNAMIC)</span>
481<a name="ln-50"></a><span class="nl">     </span> <span class="k">DO </span><span class="mi">16</span> <span class="n">L</span><span class="o">=</span><span class="mi">1</span><span class="p">,</span><span class="n">MLEVEL</span>
482<a name="ln-51"></a><span class="nl">     </span>   <span class="k">CALL </span><span class="n">GRPHSUB</span><span class="p">(</span><span class="n">L</span><span class="p">,</span><span class="n">IND</span><span class="p">,</span><span class="n">CXMN</span><span class="p">,</span><span class="n">FELD</span><span class="p">,</span><span class="n">WSAVE</span><span class="p">,</span><span class="n">IFAX</span><span class="p">,</span><span class="n">Z</span><span class="p">,</span><span class="n">W</span><span class="p">,</span><span class="n">MLAT</span><span class="p">,</span>
483<a name="ln-52"></a><span class="nl">     </span><span class="gs">*</span><span class="n">MNAUF</span><span class="p">,</span><span class="n">MAXL</span><span class="p">,</span><span class="n">MAXB</span><span class="p">,</span><span class="n">MLEVEL</span><span class="p">)</span>
484<a name="ln-53"></a><span class="nl">16   </span> <span class="k">CONTINUE</span>
485<a name="ln-54"></a><span class="c">!$omp end parallel do</span>
486<a name="ln-55"></a>
487<a name="ln-56"></a>
488<a name="ln-57"></a><span class="nl">     </span> <span class="k">RETURN</span>
489<a name="ln-58"></a><span class="nl">     </span> <span class="k">END SUBROUTINE </span><span class="n">GRPH213</span>
490<a name="ln-59"></a><span class="c">C</span>
491<a name="ln-60"></a><span class="nl">     </span> <span class="k">SUBROUTINE </span><span class="n">GRPHSUB</span><span class="p">(</span><span class="n">L</span><span class="p">,</span><span class="n">IND</span><span class="p">,</span><span class="n">CXMN</span><span class="p">,</span><span class="n">FELD</span><span class="p">,</span><span class="n">WSAVE</span><span class="p">,</span><span class="n">IFAX</span><span class="p">,</span><span class="n">Z</span><span class="p">,</span><span class="n">W</span><span class="p">,</span><span class="n">MLAT</span><span class="p">,</span>
492<a name="ln-61"></a><span class="nl">     </span><span class="gs">*</span><span class="n">MNAUF</span><span class="p">,</span><span class="n">MAXL</span><span class="p">,</span><span class="n">MAXB</span><span class="p">,</span><span class="n">MLEVEL</span><span class="p">)</span>
493<a name="ln-62"></a>
494<a name="ln-63"></a><span class="c">C     DIE ROUTINE F]HRT EINE TRANSFORMATION EINER</span>
495<a name="ln-64"></a><span class="c">C     FELDVARIABLEN VOM PHASENRAUM IN  DEN PHYSIKALISCHEN</span>
496<a name="ln-65"></a><span class="c">C     RAUM AUF KUGELKOORDINATEN DURCH</span>
497<a name="ln-66"></a><span class="c">C</span>
498<a name="ln-67"></a><span class="c">C     CXMN   = SPEKTRALKOEFFIZIENTEN IN DER REIHENFOLGE</span>
499<a name="ln-68"></a><span class="c">C              CX00,CX01,CX11,CX02,....CXMNAUFMNAUF</span>
500<a name="ln-69"></a><span class="c">C                   CXM              = FOURIERKOEFFIZIENTEN - nur ein Hilfsfeld</span>
501<a name="ln-70"></a><span class="c">C     FELD   = FELD DER METEOROLOGISCHEN VARIABLEN</span>
502<a name="ln-71"></a><span class="c">C                   WSAVE  = Working Array fuer Fouriertransformation</span>
503<a name="ln-72"></a><span class="c">C     Z      = LEGENDREFUNKTIONSWERTE</span>
504<a name="ln-73"></a><span class="c">C</span>
505<a name="ln-74"></a><span class="c">C     MNAUF    ANZAHL DER FOURIERKOEFFIZIENTEN</span>
506<a name="ln-75"></a><span class="c">C     MAXL     ANZAHL DER FUER DAS GITTER BENUTZTEN LAENGEN</span>
507<a name="ln-76"></a><span class="c">C     MAXB     ANZAHL DER FUER DAS GITTER BENOETIGTEN BREITEN</span>
508<a name="ln-77"></a><span class="c">C     MLEVEL   ANZAHL DER LEVELS, DIE TRANSFORMIERT WERDEN</span>
509<a name="ln-78"></a><span class="c">C</span>
510<a name="ln-79"></a><span class="nl">     </span> <span class="k">IMPLICIT </span><span class="kt">REAL</span> <span class="p">(</span><span class="n">A</span><span class="o">-</span><span class="n">H</span><span class="p">,</span><span class="n">O</span><span class="o">-</span><span class="n">Z</span><span class="p">)</span>
511<a name="ln-80"></a>
512<a name="ln-81"></a><span class="c">C    FELD DER FOURIERKOEFFIZIENTEN</span>
513<a name="ln-82"></a><span class="nl">     </span> <span class="kt">REAL </span><span class="n">CXMS</span><span class="p">(</span><span class="mi">4</span><span class="o">*</span><span class="p">(</span><span class="n">MNAUF</span><span class="o">+</span><span class="mi">1</span><span class="p">))</span>
514<a name="ln-83"></a><span class="nl">     </span> <span class="kt">REAL </span><span class="n">CXMA</span><span class="p">(</span><span class="mi">4</span><span class="o">*</span><span class="p">(</span><span class="n">MNAUF</span><span class="o">+</span><span class="mi">1</span><span class="p">))</span>
515<a name="ln-84"></a><span class="nl">     </span> <span class="kt">REAL</span><span class="p">,</span><span class="k">ALLOCATABLE</span> <span class="kd">::</span> <span class="n">CXM</span><span class="p">(:,:)</span>
516<a name="ln-85"></a>
517<a name="ln-86"></a><span class="c">C                   Anzahl der Gitterpunkte pro Breitenkreis des reduzierten </span>
518<a name="ln-87"></a><span class="c">C                   Gauss&#39;schen Gitters</span>
519<a name="ln-88"></a><span class="nl">                   IN</span><span class="gs">T</span><span class="n">EGER</span> <span class="n">MLAT</span><span class="p">(</span><span class="n">MAXB</span><span class="p">),</span><span class="n">ISIZE</span>
520<a name="ln-89"></a>                   
521<a name="ln-90"></a><span class="c">C    FELD DER LEGENDREPOLYNOME FUER EINE BREITE</span>
522<a name="ln-91"></a><span class="nl">     </span> <span class="kt">REAL </span><span class="n">Z</span><span class="p">(</span><span class="n">MAXB</span><span class="o">/</span><span class="mi">2</span><span class="p">,</span><span class="mi">0</span><span class="p">:((</span><span class="n">MNAUF</span><span class="o">+</span><span class="mi">3</span><span class="p">)</span><span class="o">*</span><span class="p">(</span><span class="n">MNAUF</span><span class="o">+</span><span class="mi">4</span><span class="p">))</span><span class="o">/</span><span class="mi">2</span><span class="p">)</span>
523<a name="ln-92"></a><span class="nl">     </span> 
524<a name="ln-93"></a><span class="c">C      LOGICAL*1 USED(((216*217)/2+1)*160)</span>
525<a name="ln-94"></a>
526<a name="ln-95"></a><span class="nl">     </span> <span class="kt">REAL </span><span class="n">CXMN</span><span class="p">(</span><span class="mi">0</span><span class="p">:(</span><span class="n">MNAUF</span><span class="o">+</span><span class="mi">1</span><span class="p">)</span><span class="o">*</span><span class="p">(</span><span class="n">MNAUF</span><span class="o">+</span><span class="mi">2</span><span class="p">)</span><span class="o">-</span><span class="mi">1</span><span class="p">,</span><span class="n">MLEVEL</span><span class="p">)</span>
527<a name="ln-96"></a><span class="nl">     </span> <span class="kt">REAL </span><span class="n">FELD</span><span class="p">(</span><span class="n">MAXL</span><span class="p">,</span><span class="n">MLEVEL</span><span class="p">)</span>
528<a name="ln-97"></a><span class="nl">     </span> <span class="kt">REAL </span><span class="n">WSAVE</span><span class="p">(</span><span class="mi">8</span><span class="o">*</span><span class="n">MAXB</span><span class="o">+</span><span class="mi">15</span><span class="p">,</span><span class="n">MAXB</span><span class="o">/</span><span class="mi">2</span><span class="p">)</span>
529<a name="ln-98"></a><span class="nl">     </span> <span class="kt">INTEGER </span><span class="n">IFAX</span><span class="p">(</span><span class="mi">10</span><span class="p">,</span><span class="n">MAXB</span><span class="p">)</span>
530<a name="ln-99"></a><span class="nl">     </span> <span class="kt">REAL </span><span class="n">W</span><span class="p">(</span><span class="n">MAXB</span><span class="p">)</span>
531<a name="ln-100"></a><span class="nl">     </span> <span class="kt">INTEGER </span><span class="n">IND</span><span class="p">(</span><span class="n">MAXB</span><span class="p">)</span>
532<a name="ln-101"></a>                   
533<a name="ln-102"></a><span class="nl">     </span> <span class="k">ALLOCATE</span><span class="p">(</span><span class="n">CXM</span><span class="p">(</span> <span class="mi">4</span><span class="o">*</span><span class="n">MAXB</span><span class="p">,</span><span class="n">MAXB</span><span class="p">))</span>
534<a name="ln-103"></a><span class="nl">    DO</span> <span class="mi">5</span> <span class="n">J</span><span class="o">=</span><span class="mi">1</span><span class="p">,</span><span class="n">MAXB</span><span class="o">/</span><span class="mi">2</span>
535<a name="ln-104"></a><span class="nl">     </span>       <span class="n">CXMS</span><span class="p">(</span><span class="mi">1</span><span class="p">:</span><span class="n">MLAT</span><span class="p">(</span><span class="n">J</span><span class="p">))</span><span class="o">=</span><span class="n">FELD</span><span class="p">(</span><span class="n">IND</span><span class="p">(</span><span class="n">J</span><span class="p">)</span><span class="o">+</span><span class="mi">1</span><span class="p">:</span><span class="n">IND</span><span class="p">(</span><span class="n">J</span><span class="p">)</span><span class="o">+</span><span class="n">MLAT</span><span class="p">(</span><span class="n">J</span><span class="p">),</span><span class="n">L</span><span class="p">)</span>
536<a name="ln-105"></a><span class="nl">      </span><span class="gs">C</span><span class="k">ALL </span><span class="n">RFOUFTR</span><span class="p">(</span><span class="n">CXMS</span><span class="p">,</span><span class="n">WSAVE</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span><span class="n">J</span><span class="p">),</span><span class="n">IFAX</span><span class="p">(:,</span><span class="n">J</span><span class="p">),</span><span class="n">MNAUF</span><span class="p">,</span><span class="n">MLAT</span><span class="p">(</span><span class="n">J</span><span class="p">),</span><span class="mi">1</span><span class="p">)</span>
537<a name="ln-106"></a><span class="nl">     </span>       <span class="n">CXMA</span><span class="p">(</span><span class="mi">1</span><span class="p">:</span><span class="n">MLAT</span><span class="p">(</span><span class="n">J</span><span class="p">))</span><span class="o">=</span><span class="n">FELD</span><span class="p">(</span><span class="n">MAXL</span><span class="o">-</span><span class="n">IND</span><span class="p">(</span><span class="n">J</span><span class="p">)</span><span class="o">-</span><span class="n">MLAT</span><span class="p">(</span><span class="n">J</span><span class="p">)</span><span class="o">+</span><span class="mi">1</span><span class="p">:</span><span class="n">MAXL</span><span class="o">-</span><span class="n">IND</span><span class="p">(</span><span class="n">J</span><span class="p">),</span><span class="n">L</span><span class="p">)</span>
538<a name="ln-107"></a><span class="nl">      </span><span class="gs">C</span><span class="k">ALL </span><span class="n">RFOUFTR</span><span class="p">(</span><span class="n">CXMA</span><span class="p">,</span>
539<a name="ln-108"></a><span class="nl">     </span><span class="gs">*</span><span class="n">WSAVE</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span><span class="n">J</span><span class="p">),</span><span class="n">IFAX</span><span class="p">(:,</span><span class="n">J</span><span class="p">),</span><span class="n">MNAUF</span><span class="p">,</span><span class="n">MLAT</span><span class="p">(</span><span class="n">J</span><span class="p">),</span><span class="mi">1</span><span class="p">)</span>
540<a name="ln-109"></a><span class="nl">          </span><span class="gs">D</span><span class="n">O</span> <span class="mi">4</span> <span class="n">I</span><span class="o">=</span><span class="mi">1</span><span class="p">,</span><span class="mi">2</span><span class="o">*</span><span class="p">(</span><span class="n">MNAUF</span><span class="o">+</span><span class="mi">1</span><span class="p">)</span>
541<a name="ln-110"></a><span class="nl">      </span><span class="gs">    </span><span class="n">CXM</span><span class="p">(</span><span class="n">I</span><span class="p">,</span><span class="n">J</span><span class="p">)</span><span class="o">=</span><span class="n">CXMS</span><span class="p">(</span><span class="n">I</span><span class="p">)</span><span class="o">+</span><span class="n">CXMA</span><span class="p">(</span><span class="n">I</span><span class="p">)</span>
542<a name="ln-111"></a><span class="nl">      </span><span class="gs">    </span><span class="n">CXM</span><span class="p">(</span><span class="n">I</span><span class="p">,</span><span class="n">MAXB</span><span class="o">+</span><span class="mi">1</span><span class="o">-</span><span class="n">J</span><span class="p">)</span><span class="o">=</span><span class="n">CXMS</span><span class="p">(</span><span class="n">I</span><span class="p">)</span><span class="o">-</span><span class="n">CXMA</span><span class="p">(</span><span class="n">I</span><span class="p">)</span>
543<a name="ln-112"></a><span class="nl">4                 C</span><span class="gs">O</span><span class="n">NTINUE</span>
544<a name="ln-113"></a><span class="nl">    5</span>      <span class="k">CONTINUE</span>
545<a name="ln-114"></a><span class="nl">     </span>   <span class="k">CALL </span><span class="n">LGTR213</span><span class="p">(</span><span class="n">CXMN</span><span class="p">(</span><span class="mi">0</span><span class="p">,</span><span class="n">L</span><span class="p">),</span><span class="n">CXM</span><span class="p">,</span><span class="n">Z</span><span class="p">,</span><span class="n">W</span><span class="p">,</span><span class="n">MLAT</span><span class="p">,</span><span class="n">MNAUF</span><span class="p">,</span><span class="n">MAXB</span><span class="p">)</span>
546<a name="ln-115"></a>
547<a name="ln-116"></a><span class="nl">     </span> <span class="k">DEALLOCATE</span><span class="p">(</span><span class="n">CXM</span><span class="p">)</span>
548<a name="ln-117"></a>
549<a name="ln-118"></a><span class="nl">     </span> <span class="k">RETURN</span>
550<a name="ln-119"></a><span class="nl">     </span> <span class="k">END SUBROUTINE </span><span class="n">GRPHSUB</span>
551<a name="ln-120"></a><span class="c">C</span>
552<a name="ln-121"></a><span class="nl">     </span>   <span class="k">SUBROUTINE </span><span class="n">LGTR213</span><span class="p">(</span><span class="n">CXMN</span><span class="p">,</span><span class="n">CXM</span><span class="p">,</span><span class="n">Z</span><span class="p">,</span><span class="n">W</span><span class="p">,</span><span class="n">MLAT</span><span class="p">,</span><span class="n">MNAUF</span><span class="p">,</span><span class="n">MAXB</span><span class="p">)</span>
553<a name="ln-122"></a><span class="nl">     </span>   <span class="k">IMPLICIT </span><span class="kt">REAL</span> <span class="p">(</span><span class="n">A</span><span class="o">-</span><span class="n">H</span><span class="p">,</span><span class="n">O</span><span class="o">-</span><span class="n">Z</span><span class="p">)</span>
554<a name="ln-123"></a><span class="nl">     </span>   <span class="kt">INTEGER </span><span class="n">MLAT</span><span class="p">(</span><span class="n">MAXB</span><span class="p">)</span>
555<a name="ln-124"></a><span class="nl">     </span>   <span class="k">DIMENSION </span><span class="n">CXM</span><span class="p">(</span><span class="mi">0</span><span class="p">:</span><span class="mi">4</span><span class="o">*</span><span class="n">MAXB</span><span class="o">-</span><span class="mi">1</span><span class="p">,</span><span class="n">MAXB</span><span class="p">)</span>
556<a name="ln-125"></a><span class="nl">     </span>   <span class="k">DIMENSION </span><span class="n">CXMN</span><span class="p">(</span><span class="mi">0</span><span class="p">:</span><span class="mi">2</span><span class="o">*</span><span class="p">(((</span><span class="n">MNAUF</span><span class="o">+</span><span class="mi">1</span><span class="p">)</span><span class="o">*</span><span class="n">MNAUF</span><span class="p">)</span><span class="o">/</span><span class="mi">2</span><span class="o">+</span><span class="n">MNAUF</span><span class="p">)</span><span class="o">+</span><span class="mi">1</span><span class="p">)</span>
557<a name="ln-126"></a><span class="nl">     </span>   <span class="kt">REAL</span><span class="o">*</span><span class="mi">8</span> <span class="n">Z</span><span class="p">(</span><span class="n">MAXB</span><span class="o">/</span><span class="mi">2</span><span class="p">,</span><span class="mi">0</span><span class="p">:((</span><span class="n">MNAUF</span><span class="o">+</span><span class="mi">3</span><span class="p">)</span><span class="o">*</span><span class="p">(</span><span class="n">MNAUF</span><span class="o">+</span><span class="mi">4</span><span class="p">))</span><span class="o">/</span><span class="mi">2</span><span class="p">)</span>
558<a name="ln-127"></a><span class="nl">     </span>   <span class="kt">REAL</span><span class="o">*</span><span class="mi">8</span> <span class="n">W</span><span class="p">(</span><span class="n">MAXB</span><span class="p">),</span><span class="n">CR</span><span class="p">,</span><span class="n">CI</span><span class="p">,</span><span class="n">HILF</span>
559<a name="ln-128"></a><span class="nl">     </span>   <span class="kt">LOGICAL </span><span class="n">EVEN</span>
560<a name="ln-129"></a><span class="c">C</span>
561<a name="ln-130"></a><span class="c">C     DIESE ROUTINE BERECHNET DIE KFFKs CXMN</span>
562<a name="ln-131"></a><span class="c">C</span>
563<a name="ln-132"></a><span class="nl">                  LL</span><span class="gs">=</span><span class="mi">0</span>
564<a name="ln-133"></a><span class="nl">                  LL</span><span class="gs">P</span><span class="o">=</span><span class="mi">0</span>
565<a name="ln-134"></a><span class="nl">     </span> <span class="k">DO </span><span class="mi">1</span> <span class="n">I</span><span class="o">=</span><span class="mi">0</span><span class="p">,</span><span class="n">MNAUF</span>
566<a name="ln-135"></a><span class="nl">     </span>   <span class="n">KM</span><span class="o">=</span><span class="mi">0</span>
567<a name="ln-136"></a><span class="nl">    9</span>   <span class="n">KM</span><span class="o">=</span><span class="n">KM</span><span class="o">+</span><span class="mi">1</span>
568<a name="ln-137"></a><span class="nl">     </span>   <span class="k">IF</span><span class="p">(</span><span class="n">MLAT</span><span class="p">(</span><span class="n">KM</span><span class="p">).</span><span class="n">LE</span><span class="p">.</span><span class="mi">2</span><span class="o">*</span><span class="n">I</span><span class="p">)</span> <span class="k">THEN</span>
569<a name="ln-138"></a><span class="nl">     </span>      <span class="k">GOTO </span><span class="mi">9</span>
570<a name="ln-139"></a><span class="nl">     </span>   <span class="n">ENDIF</span>
571<a name="ln-140"></a><span class="nl">     </span>   <span class="k">DO </span><span class="mi">2</span> <span class="n">J</span><span class="o">=</span><span class="n">I</span><span class="p">,</span><span class="n">MNAUF</span>
572<a name="ln-141"></a><span class="nl">     </span>     <span class="n">CR</span><span class="o">=</span><span class="mi">0</span>
573<a name="ln-142"></a><span class="nl">     </span>     <span class="n">CI</span><span class="o">=</span><span class="mi">0</span>
574<a name="ln-143"></a><span class="nl">     </span>     <span class="n">EVEN</span><span class="o">=</span><span class="nb">MOD</span><span class="p">(</span><span class="n">I</span><span class="o">+</span><span class="n">J</span><span class="p">,</span><span class="mi">2</span><span class="p">).</span><span class="n">EQ</span><span class="p">.</span><span class="mi">0</span>
575<a name="ln-144"></a><span class="nl">     </span>     <span class="k">IF</span><span class="p">(</span><span class="n">EVEN</span><span class="p">)</span> <span class="k">THEN</span>
576<a name="ln-145"></a><span class="nl">     </span>       <span class="k">DO </span><span class="mi">3</span> <span class="n">K</span><span class="o">=</span><span class="n">KM</span><span class="p">,</span><span class="n">MAXB</span><span class="o">/</span><span class="mi">2</span>
577<a name="ln-146"></a><span class="nl">     </span>                        <span class="n">HILF</span><span class="o">=</span><span class="n">W</span><span class="p">(</span><span class="n">K</span><span class="p">)</span><span class="o">*</span><span class="n">Z</span><span class="p">(</span><span class="n">K</span><span class="p">,</span><span class="n">LLP</span><span class="p">)</span>
578<a name="ln-147"></a><span class="nl">     </span>            <span class="n">CR</span><span class="o">=</span><span class="n">CR</span><span class="o">+</span><span class="n">CXM</span><span class="p">(</span><span class="mi">2</span><span class="o">*</span><span class="n">I</span><span class="p">,</span><span class="n">K</span><span class="p">)</span><span class="o">*</span><span class="n">HILF</span>
579<a name="ln-148"></a><span class="nl">     </span>            <span class="n">CI</span><span class="o">=</span><span class="n">CI</span><span class="o">+</span><span class="n">CXM</span><span class="p">(</span><span class="mi">2</span><span class="o">*</span><span class="n">I</span><span class="o">+</span><span class="mi">1</span><span class="p">,</span><span class="n">K</span><span class="p">)</span><span class="o">*</span><span class="n">HILF</span>
580<a name="ln-149"></a><span class="nl">    3</span>      <span class="k">CONTINUE</span>
581<a name="ln-150"></a><span class="nl">     </span>     <span class="k">ELSE</span>
582<a name="ln-151"></a><span class="nl">     </span>       <span class="k">DO </span><span class="mi">4</span> <span class="n">K</span><span class="o">=</span><span class="n">KM</span><span class="p">,</span><span class="n">MAXB</span><span class="o">/</span><span class="mi">2</span>
583<a name="ln-152"></a><span class="nl">     </span>                        <span class="n">HILF</span><span class="o">=</span><span class="n">W</span><span class="p">(</span><span class="n">K</span><span class="p">)</span><span class="o">*</span><span class="n">Z</span><span class="p">(</span><span class="n">K</span><span class="p">,</span><span class="n">LLP</span><span class="p">)</span>
584<a name="ln-153"></a><span class="nl">     </span>            <span class="n">CR</span><span class="o">=</span><span class="n">CR</span><span class="o">+</span><span class="n">CXM</span><span class="p">(</span><span class="mi">2</span><span class="o">*</span><span class="n">I</span><span class="p">,</span><span class="n">MAXB</span><span class="o">+</span><span class="mi">1</span><span class="o">-</span><span class="n">K</span><span class="p">)</span><span class="o">*</span><span class="n">HILF</span>
585<a name="ln-154"></a><span class="nl">     </span>            <span class="n">CI</span><span class="o">=</span><span class="n">CI</span><span class="o">+</span><span class="n">CXM</span><span class="p">(</span><span class="mi">2</span><span class="o">*</span><span class="n">I</span><span class="o">+</span><span class="mi">1</span><span class="p">,</span><span class="n">MAXB</span><span class="o">+</span><span class="mi">1</span><span class="o">-</span><span class="n">K</span><span class="p">)</span><span class="o">*</span><span class="n">HILF</span>
586<a name="ln-155"></a><span class="nl">    4</span>      <span class="k">CONTINUE</span>
587<a name="ln-156"></a><span class="nl">     </span>     <span class="n">ENDIF</span>
588<a name="ln-157"></a><span class="nl">    5</span>     <span class="n">CXMN</span><span class="p">(</span><span class="mi">2</span><span class="o">*</span><span class="n">LL</span><span class="p">)</span><span class="o">=</span><span class="n">CR</span>
589<a name="ln-158"></a><span class="nl">     </span>     <span class="n">CXMN</span><span class="p">(</span><span class="mi">2</span><span class="o">*</span><span class="n">LL</span><span class="o">+</span><span class="mi">1</span><span class="p">)</span><span class="o">=</span><span class="n">CI</span>
590<a name="ln-159"></a><span class="nl">     </span>     <span class="n">LL</span><span class="o">=</span><span class="n">LL</span><span class="o">+</span><span class="mi">1</span>
591<a name="ln-160"></a><span class="nl">     </span>     <span class="n">LLP</span><span class="o">=</span><span class="n">LLP</span><span class="o">+</span><span class="mi">1</span>
592<a name="ln-161"></a><span class="nl">    2</span>   <span class="k">CONTINUE</span>
593<a name="ln-162"></a><span class="nl">     </span>   <span class="n">LLP</span><span class="o">=</span><span class="n">LLP</span><span class="o">+</span><span class="mi">2</span>
594<a name="ln-163"></a><span class="nl">    1</span> <span class="k">CONTINUE</span>
595<a name="ln-164"></a><span class="nl">     </span>   <span class="k">RETURN</span>
596<a name="ln-165"></a><span class="nl">     </span>   <span class="k">END SUBROUTINE </span><span class="n">LGTR213</span>
597<a name="ln-166"></a><span class="c">C</span>
598<a name="ln-167"></a>
599<a name="ln-168"></a><span class="c">C</span>
600<a name="ln-169"></a><span class="nl">     </span> <span class="k">SUBROUTINE </span><span class="n">RFOUFTR</span><span class="p">(</span><span class="n">CXM</span><span class="p">,</span><span class="n">TRIGS</span><span class="p">,</span><span class="n">IFAX</span><span class="p">,</span><span class="n">MNAUF</span><span class="p">,</span><span class="n">MAXL</span><span class="p">,</span><span class="nb">ISIGN</span><span class="p">)</span>
601<a name="ln-170"></a><span class="c">C     BERECHNET DIE FOURIERSUMME MIT EINEM FFT-ALGORITHMUS</span>
602<a name="ln-171"></a><span class="nl">     </span> <span class="k">IMPLICIT </span><span class="kt">REAL</span> <span class="p">(</span><span class="n">A</span><span class="o">-</span><span class="n">H</span><span class="p">,</span><span class="n">O</span><span class="o">-</span><span class="n">Z</span><span class="p">)</span>
603<a name="ln-172"></a><span class="nl">     </span> <span class="k">DIMENSION </span><span class="n">CXM</span><span class="p">(</span><span class="mi">0</span><span class="p">:</span><span class="mi">2</span><span class="o">*</span><span class="n">MAXL</span><span class="o">-</span><span class="mi">1</span><span class="p">)</span>
604<a name="ln-173"></a><span class="nl">     </span> <span class="k">DIMENSION </span><span class="n">FELD</span><span class="p">(</span><span class="n">MAXL</span><span class="p">),</span><span class="n">TRIGS</span><span class="p">(</span><span class="mi">2</span><span class="o">*</span><span class="n">MAXL</span><span class="p">)</span>
605<a name="ln-174"></a><span class="nl">     </span> <span class="k">DIMENSION </span><span class="n">WSAVE</span><span class="p">(</span><span class="n">MAXAUF</span><span class="p">)</span>
606<a name="ln-175"></a><span class="nl">     </span> <span class="kt">INTEGER </span><span class="n">IFAX</span><span class="p">(</span><span class="mi">10</span><span class="p">)</span>
607<a name="ln-176"></a>
608<a name="ln-177"></a>
609<a name="ln-178"></a><span class="c">C NORMIERUNG...</span>
610<a name="ln-179"></a><span class="nl">     </span> <span class="n">WSAVE</span><span class="p">(</span><span class="mi">1</span><span class="p">)</span><span class="o">=</span><span class="n">CXM</span><span class="p">(</span><span class="n">MAXL</span><span class="o">-</span><span class="mi">1</span><span class="p">)</span>
611<a name="ln-180"></a>
612<a name="ln-181"></a><span class="nl">     </span> <span class="n">CXM</span><span class="p">(</span><span class="mi">1</span><span class="p">:</span><span class="n">MAXL</span><span class="p">)</span><span class="o">=</span><span class="n">CXM</span><span class="p">(</span><span class="mi">0</span><span class="p">:</span><span class="n">MAXL</span><span class="o">-</span><span class="mi">1</span><span class="p">)</span><span class="o">/</span><span class="mi">2</span>
613<a name="ln-182"></a><span class="nl">     </span> <span class="n">CXM</span><span class="p">(</span><span class="mi">0</span><span class="p">)</span><span class="o">=</span><span class="n">WSAVE</span><span class="p">(</span><span class="mi">1</span><span class="p">)</span><span class="o">/</span><span class="mi">2</span>
614<a name="ln-183"></a><span class="c">!      CALL CFFTF(MAXL,CXM,WSAVE)</span>
615<a name="ln-184"></a><span class="nl">     </span> <span class="k">CALL </span><span class="n">FFT99</span><span class="p">(</span><span class="n">CXM</span><span class="p">,</span><span class="n">WSAVE</span><span class="p">,</span><span class="n">TRIGS</span><span class="p">,</span><span class="n">IFAX</span><span class="p">,</span><span class="mi">1</span><span class="p">,</span><span class="mi">1</span><span class="p">,</span><span class="n">MAXL</span><span class="p">,</span><span class="mi">1</span><span class="p">,</span><span class="o">-</span><span class="mi">1</span><span class="p">)</span>
616<a name="ln-185"></a><span class="nl">     </span> <span class="k">RETURN</span>
617<a name="ln-186"></a><span class="nl">     </span> <span class="k">END SUBROUTINE </span><span class="n">RFOUFTR</span>
618<a name="ln-187"></a>
619<a name="ln-188"></a><span class="nl">     </span> <span class="k">END MODULE </span><span class="n">GRTOPH</span>
620</pre></div>
621
622    </section>
623    </div>
624  </div>
625
626  <section class="visible-xs visible-sm hidden-md">
627    <hr>
628   
629
630<div class="panel panel-default">
631  <div class="panel-heading text-left"><h3 class="panel-title"><a data-toggle="collapse" href="#allfiles-1">All Source Files</a></h3></div>
632  <div id="allfiles-1" class="panel-collapse collapse">
633    <div class="list-group">
634     
635      <a class="list-group-item" href="../sourcefile/grphreal.f.html">grphreal.f</a>
636     
637      <a class="list-group-item" href="../sourcefile/phgrreal.f.html">phgrreal.f</a>
638     
639      <a class="list-group-item" href="../sourcefile/preconvert.f90.html">preconvert.f90</a>
640     
641      <a class="list-group-item" href="../sourcefile/rwgrib2.f90.html">rwGRIB2.f90</a>
642     
643    </div>
644  </div>
645</div>
646
647
648  </section>
649 
650    <hr>   
651    </div> <!-- /container -->
652    <footer>
653      <div class="container">
654      <div class="row">
655        <div class="col-xs-6 col-md-4"><p>&copy; 2019 <a rel="license" href="http://creativecommons.org/licenses/by/4.0/"><img alt="Creative Commons License" style="border-width:0" src="../80x15.png" /></a></p></div>
656        <div class="col-xs-6 col-md-4 col-md-push-4">
657          <p class="text-right">
658            Documentation generated by
659            <a href="https://github.com/cmacmackin/ford">FORD</a>
660           
661          </p>
662        </div>
663        <div class="col-xs-12 col-md-4 col-md-pull-4"><p class="text-center"> Flex_extract: Calculation of etadot was developed by Leopold Haimberger</p></div>
664      </div>
665      <br>
666      </div> <!-- /container -->   
667    </footer>
668
669    <!-- Bootstrap core JavaScript
670    ================================================== -->
671    <!-- Placed at the end of the document so the pages load faster -->
672<!--
673    <script src="https://ajax.googleapis.com/ajax/libs/jquery/1.11.1/jquery.min.js"></script>
674-->
675    <script src="../js/bootstrap.min.js"></script>
676    <!-- IE10 viewport hack for Surface/desktop Windows 8 bug -->
677    <script src="../js/ie10-viewport-bug-workaround.js"></script>
678
679    <!-- MathJax JavaScript
680    ================================================== -->
681    <!-- Placed at the end of the document so the pages load faster -->
682    <script type="text/x-mathjax-config">
683      MathJax.Hub.Config({
684        TeX: { extensions: ['AMSmath.js','AMSsymbols.js','noErrors.js','noUndefined.js'], equationNumbers: { autoNumber: 'AMS' } },
685        jax: ['input/TeX','input/MathML','output/HTML-CSS'],
686        extensions: ['tex2jax.js','mml2jax.js','MathMenu.js','MathZoom.js'],
687        'HTML-CSS': { 
688           styles: { '.MathJax_Display, .MathJax .mo, .MathJax .mi, .MathJax .mn': {color: '#000000 ! important'} }
689        }
690      });
691    </script>
692    <!-- <script src="https://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script> -->
693  </body>
694</html>
Note: See TracBrowser for help on using the repository browser.
hosted by ZAMG