Changeset 62 for palm/trunk/DOC/tec/technical_documentation.html
- Timestamp:
- Mar 13, 2007 2:52:40 AM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/DOC/tec/technical_documentation.html
r54 r62 1 1 <!DOCTYPE html PUBLIC "-//w3c//dtd html 4.0 transitional//en"> 2 2 <html><head> 3 4 5 <meta content="text/html; charset=ISO-8859-1" http-equiv="Content-Type"> 6 7 8 <meta content="Mozilla/4.73 [en] (X11; U; HP-UX B.10.20 9000/782) [Netscape]" name="GENERATOR"> 9 10 11 <meta content="Siegfried Raasch" name="Author"> 12 13 14 <meta content="parallel LES model" name="Keywords"><title>PALM technical/numerical documentation</title></head> 15 <body> 16 17 <font size="+4">PALM 3 <meta content="text/html; charset=ISO-8859-1" http-equiv="Content-Type"> <meta content="Mozilla/4.73 [en] (X11; U; HP-UX B.10.20 9000/782) [Netscape]" name="GENERATOR"> <meta content="Siegfried Raasch" name="Author"> <meta content="parallel LES model" name="Keywords"><title>PALM technical/numerical 4 documentation</title></head> 5 <body><font size="+4">PALM 18 6 </font>a 19 7 <b>pa</b>rallelized 20 8 <b>L</b>ES <b>m</b>odel 21 <br> 22 23 <hr width="100%"> 24 <br> 25 26 <font size="+2">technical/numerical documentation</font> 27 <br> 28 29 30 <br> 31 32 33 <br> 34 35 <table nosave="" cellpadding="0" cellspacing="0"> 36 37 <caption> <br> 38 39 </caption><tbody> 40 41 42 </tbody><tbody> 43 44 45 </tbody> <tbody> 46 47 <tr nosave=""> 48 49 <td nosave=""><b>Current model version:</b></td> 50 51 <td><span style="font-weight: bold;">3.1c</span></td> 52 53 </tr> 54 55 <tr nosave=""> 56 57 <td nosave=""><b>Last change of this document</b>: </td> 58 59 <td nosave=""><b>05/03/07</b></td> 60 61 </tr> 62 63 64 </tbody> 65 </table> 66 67 <hr width="100%"> 9 <br><hr width="100%"> 10 <br><font size="+2">technical/numerical documentation</font> 11 <br> 12 <br> 13 <br><table nosave="" cellpadding="0" cellspacing="0"> <caption> <br> </caption><tbody> 14 </tbody><tbody> </tbody> <tbody> <tr nosave=""> <td nosave=""><b>Current 15 model version:</b></td> <td><span style="font-weight: bold;">3.1c</span></td> </tr> 16 <tr nosave=""> <td nosave=""><b>Last 17 change of this document</b>: </td> <td nosave=""><b>$Id$</b></td> </tr> 18 </tbody> 19 </table><hr width="100%"> 68 20 <p><font size="+2">Table of contents</font> 69 <br> 70 71 72 <br> 73 74 75 </p> 76 77 <ul> 78 79 <a href="#Kapitel1.0">1.0</a> History of changes 80 81 <p><a href="#Kapitel2.0">2.0</a> How to change the model source 82 code </p> 83 84 85 <p><a href="#Kapitel3.0">3.0</a> Description of selected parts 21 <br> 22 <br> 23 </p><ul> <a href="#Kapitel1.0">1.0</a> 24 History of changes <p><a href="#Kapitel2.0">2.0</a> 25 How to change the model source 26 code </p> <p><a href="#Kapitel3.0">3.0</a> 27 Description of selected parts 86 28 of 87 29 the model source code <br> 88 89 30 </p> 90 91 </ul> 92 93 <b>Please note: Parts of this document are still in German. We 31 </ul><b>Please note: Parts of this document are still in 32 German. We 94 33 apologize 95 34 for any inconvenience.</b> 96 <br> 97 98 99 <h2><a name="Kapitel1.0"></a>1.0 History of changes</h2> 100 35 <br> 36 <h2><a name="Kapitel1.0"></a>1.0 37 History of changes</h2> 101 38 The following table documents all changes to the model since July, 7th 102 39 1998. Entries in the fourth column specify the type of the change: … … 107 44 respective 108 45 routines. 109 <br> 110 111 112 <br> 113 114 <table nosave="" style="width: 100%;" border="1"> 115 116 <caption> <br> 117 118 </caption><tbody> 119 120 121 </tbody><tbody> 122 123 124 </tbody> <tbody> 125 126 <tr> 127 128 <td><b>date</b></td> 129 130 <td><b>author</b></td> 131 132 <td><b>model</b> <br> 133 134 <b>version</b></td> 135 136 <td><b>type</b></td> 137 138 <td><b>description</b></td> 139 140 <td><b>affected routines</b></td> 141 142 </tr> 143 144 <tr nosave="" valign="top"> 145 146 <td>06/07/98</td> 147 148 <td>SR</td> 149 150 <td>1.0</td> 151 152 <td>C</td> 153 154 <td nosave="">Änderung von Default-Werten 46 <br> 47 <br><table nosave="" style="width: 100%;" border="1"> <caption> <br> </caption><tbody> 48 </tbody><tbody> </tbody> <tbody> <tr> <td><b>date</b></td> 49 <td><b>author</b></td> <td><b>model</b> 50 <br> <b>version</b></td> <td><b>type</b></td> 51 <td><b>description</b></td> <td><b>affected 52 routines</b></td> </tr> <tr nosave="" valign="top"> <td>06/07/98</td> <td>SR</td> 53 <td>1.0</td> <td>C</td> <td nosave="">Änderung 54 von Default-Werten 155 55 (initializing_actions, end_time, 156 56 fcl_factor, use_prior_plot1d_parameters), zusätzliche 157 57 Überprüfung 158 von Parametern auf Zulässigkeit (prandtl_layer, psolver, nx, ny, 58 von Parametern auf Zulässigkeit (prandtl_layer, psolver, nx, 59 ny, 159 60 nz, 160 fcl_factor), Höhe der Schnittebenen der Horizontalschnitte wird 61 fcl_factor), Höhe der Schnittebenen der Horizontalschnitte 62 wird 161 63 auf 162 Plots exakt ausgegeben (uv- bzw. w-Gitter).</td> 163 164 <td>modules, check_parameters, plot_2d</td> 165 166 </tr> 167 168 <tr nosave="" valign="top"> 169 170 <td>"</td> 171 172 <td>SR</td> 173 174 <td>1.0</td> 175 176 <td>N</td> 177 178 <td>Neues Modul test_variables wird in (fast) allen 64 Plots exakt ausgegeben (uv- bzw. w-Gitter).</td> <td>modules, 65 check_parameters, plot_2d</td> </tr> <tr nosave="" valign="top"> <td>"</td> <td>SR</td> 66 <td>1.0</td> <td>N</td> <td>Neues Modul 67 test_variables wird in (fast) allen 179 68 Programmteilen benutzt 180 69 (USE). <b>Wichtig</b>: dieses Modul steht in der Datei <tt>module_test.f90. 181 70 D</tt>iese Datei muß in der mrun- und add_library - 182 71 Konfigurationsdatei 183 bei der Variablen <tt>module_files</tt> aufgeführt werden.</td> 184 185 <td nosave="">(fast) alle <br> 186 187 Neu: module_test</td> 188 189 </tr> 190 191 <tr nosave="" valign="top"> 192 193 <td>16/07/98</td> 194 195 <td>SR</td> 196 197 <td>1.0</td> 198 199 <td>N</td> 200 201 <td nosave="">Anpassung an die T3E-Rechner in Juelich (cpp - 202 Direktiven)</td> 203 204 <td>advec_s_bc, check_open, check_parameters, cpu_zeitmessung, 72 bei der Variablen <tt>module_files</tt> 73 aufgeführt werden.</td> <td nosave="">(fast) 74 alle <br> 75 Neu: module_test</td> </tr> <tr nosave="" valign="top"> <td>16/07/98</td> <td>SR</td> 76 <td>1.0</td> <td>N</td> <td nosave="">Anpassung 77 an die T3E-Rechner in Juelich (cpp - 78 Direktiven)</td> <td>advec_s_bc, check_open, 79 check_parameters, cpu_zeitmessung, 205 80 init_pegrid, 206 81 local_getenv, local_system, local_tremain, local_tremain_ini, poisfft</td> 207 208 </tr> 209 210 <tr nosave="" valign="top"> 211 212 <td>24/07/98</td> 213 214 <td>MS</td> 215 216 <td>1.0</td> 217 218 <td>C</td> 219 220 <td nosave="">Anpassung des Druckloesealgorithmusses an 82 </tr> <tr nosave="" valign="top"> <td>24/07/98</td> 83 <td>MS</td> <td>1.0</td> <td>C</td> 84 <td nosave="">Anpassung des Druckloesealgorithmusses an 221 85 Vorgehensweise von 222 Schmidt, Schumann und Volkert (DFVLR-Mitteilung 84-15, S. 20-26) </td> 223 224 <td>poisfft</td> 225 226 </tr> 227 228 <tr nosave="" valign="top"> 229 230 <td>27/07/98</td> 231 232 <td>MS</td> 233 234 <td>1.0</td> 235 236 <td>C</td> 237 238 <td nosave="">Berechnung von l in maketri geschieht nun 86 Schmidt, Schumann und Volkert (DFVLR-Mitteilung 84-15, S. 87 20-26) </td> <td>poisfft</td> </tr> <tr nosave="" valign="top"> <td>27/07/98</td> 88 <td>MS</td> <td>1.0</td> <td>C</td> 89 <td nosave="">Berechnung von l in maketri geschieht nun 239 90 ausserhalb der k-Schleife, 240 da unabhaengig von k, dadurch Rechenzeitersparnis </td> 241 242 <td>poisfft</td> 243 244 </tr> 245 246 <tr nosave="" valign="top"> 247 248 <td>05/08/98</td> 249 250 <td>SR</td> 251 252 <td>1.0</td> 253 254 <td>C</td> 255 256 <td>Entfernung der Sonderbehandlung der Advektionsterme 91 da unabhaengig von k, dadurch Rechenzeitersparnis </td> <td>poisfft</td> 92 </tr> <tr nosave="" valign="top"> <td>05/08/98</td> 93 <td>SR</td> <td>1.0</td> <td>C</td> 94 <td>Entfernung der Sonderbehandlung der Advektionsterme 257 95 (Piascek-Williams) 258 bei k=1. Terme werden nun überall gleich berechnet. Evtl. leichte 96 bei k=1. Terme werden nun überall gleich berechnet. Evtl. 97 leichte 259 98 Erhöhung der Varianzen in Oberflächennähe.</td> 260 261 <td nosave="">advec_u_pw, advec_v_pw, advec_w_pw</td> 262 263 </tr> 264 265 <tr nosave="" valign="top"> 266 267 <td nosave="">"</td> 268 269 <td>SR</td> 270 271 <td>1.0</td> 272 273 <td>C</td> 274 275 <td>Mischungsweganpassung (s. adjust_mixing_length) erfolgt nur 99 <td nosave="">advec_u_pw, advec_v_pw, advec_w_pw</td> 100 </tr> <tr nosave="" valign="top"> <td nosave="">"</td> <td>SR</td> <td>1.0</td> 101 <td>C</td> <td>Mischungsweganpassung (s. 102 adjust_mixing_length) erfolgt nur 276 103 noch bei 277 k=1. Unrealistische sekundäre Maxima in den Km-Profilen oberhalb 104 k=1. Unrealistische sekundäre Maxima in den Km-Profilen 105 oberhalb 278 106 der 279 Prandtl-Schicht werden damit vermieden. Leichte Änderungen in den 280 oberflächennahen Turbulenzgrößen sind die Folge.</td> 281 282 <td>diffusivities</td> 283 284 </tr> 285 286 <tr nosave="" valign="top"> 287 288 <td>"</td> 289 290 <td>SR</td> 291 292 <td>1.0</td> 293 294 <td>B</td> 295 296 <td>Falsche Berechnung der benötigten CPU-Zeit pro Sekunde 107 Prandtl-Schicht werden damit vermieden. Leichte Änderungen in 108 den 109 oberflächennahen Turbulenzgrößen sind die 110 Folge.</td> <td>diffusivities</td> </tr> <tr nosave="" valign="top"> <td>"</td> <td>SR</td> 111 <td>1.0</td> <td>B</td> <td>Falsche 112 Berechnung der benötigten CPU-Zeit pro Sekunde 297 113 simulierter 298 Zeit bei Fortsetzungsläufen korrigiert. 299 300 <p>Berechnung von cross_uxmin(1) und cross_uxmax(1) aus den 114 Zeit bei Fortsetzungsläufen korrigiert. <p>Berechnung 115 von cross_uxmin(1) und cross_uxmax(1) aus den 301 116 Werten von 302 ug und vg entfernt. Lieferte falsche Wertebereiche für Plots von 117 ug und vg entfernt. Lieferte falsche Wertebereiche für Plots 118 von 303 119 Vertikalprofilen, 304 120 wenn im ersten Koordinatenkreuz nicht die Windgeschwindigkeiten 305 121 gezeichnet 306 wurden.</p> 307 308 </td> 309 310 <td nosave="">header, check_parameters</td> 311 312 </tr> 313 314 <tr nosave="" valign="top"> 315 316 <td nosave="">"</td> 317 318 <td>SR</td> 319 320 <td>1.0</td> 321 322 <td>N</td> 323 324 <td>Das erlaubte Maximum und Minimum der Rif-Zahl ist nun mittels 122 wurden.</p> </td> <td nosave="">header, 123 check_parameters</td> </tr> <tr nosave="" valign="top"> <td nosave="">"</td> <td>SR</td> 124 <td>1.0</td> <td>N</td> <td>Das 125 erlaubte Maximum und Minimum der Rif-Zahl ist nun mittels 325 126 der Initialparameter 326 127 rif_max und rif_min einstellbar (Werte waren bisher fest auf -2.0, +1.0 327 eingestellt).</td> 328 329 <td>check_parameters, header, init_1d_model, modules, parin, 128 eingestellt).</td> <td>check_parameters, header, 129 init_1d_model, modules, parin, 330 130 prandtl_fluxes, 331 read_var_list, write_var_list</td> 332 333 </tr> 334 335 <tr nosave="" valign="top"> 336 337 <td nosave="">23/09/98</td> 338 339 <td>SR</td> 340 341 <td>1.0</td> 342 343 <td>B</td> 344 345 <td>Initialisierung von Feldern im Fall von initializing_action = 131 read_var_list, write_var_list</td> </tr> <tr nosave="" valign="top"> <td nosave="">23/09/98</td> 132 <td>SR</td> <td>1.0</td> <td>B</td> 133 <td>Initialisierung von Feldern im Fall von initializing_action = 346 134 "<i>set_constant_profiles</i>" 347 und "<i>set_1d-model_profiles</i>" korrigiert. 348 349 <p>2D-Plotausgabe von u* und theta* korrigiert. Hier wurden 135 und "<i>set_1d-model_profiles</i>" korrigiert. <p>2D-Plotausgabe 136 von u* und theta* korrigiert. Hier wurden 350 137 bisher falsche 351 Felder herausgeschrieben. </p> 352 353 354 <p>Bei SOR-Aufrufen war bisher dp nicht initialisiert. 355 Entsprechend geändert.</p> 356 357 </td> 358 359 <td>init_3d_model, plot_2d, pres</td> 360 361 </tr> 362 363 <tr nosave="" valign="top"> 364 365 <td>"</td> 366 367 <td>SR</td> 368 369 <td>1.0</td> 370 371 <td>C</td> 372 373 <td nosave="">Ausgabeformate für Statistik-Ausgaben auf 138 Felder herausgeschrieben. </p> <p>Bei 139 SOR-Aufrufen war bisher dp nicht initialisiert. 140 Entsprechend geändert.</p> </td> <td>init_3d_model, 141 plot_2d, pres</td> </tr> <tr nosave="" valign="top"> <td>"</td> <td>SR</td> 142 <td>1.0</td> <td>C</td> <td nosave="">Ausgabeformate 143 für Statistik-Ausgaben auf 374 144 Datei RUN_CONTROL 375 um einige Kommastellen erweitert. 376 377 <p>Lieferung von ausführlicheren Informationen bei 145 um einige Kommastellen erweitert. <p>Lieferung von 146 ausführlicheren Informationen bei 378 147 Programmabbruch 379 durch Unterschreitung des minimal erlaubten Zeitschritts.</p> 380 381 </td> 382 383 <td>run_control, timestep</td> 384 385 </tr> 386 387 <tr nosave="" valign="top"> 388 389 <td nosave="">"</td> 390 391 <td>SR</td> 392 393 <td>1.0</td> 394 395 <td>N</td> 396 397 <td>Rechnungen mit geneigter Ebene zur Simulation von 148 durch Unterschreitung des minimal erlaubten Zeitschritts.</p> </td> 149 <td>run_control, timestep</td> </tr> <tr nosave="" valign="top"> <td nosave="">"</td> 150 <td>SR</td> <td>1.0</td> <td>N</td> 151 <td>Rechnungen mit geneigter Ebene zur Simulation von 398 152 Kaltluftabflüssen 399 153 sind erlaubt. Gestartet und gesteuert werden sie mit dem … … 401 155 Realisiert wird die geneigte Ebene durch entsprechende 402 156 Schrägstellung 403 des Schwerkraftvektors. Dadurch erhält auch die u-Komponente der 157 des Schwerkraftvektors. Dadurch erhält auch die u-Komponente 158 der 404 159 Windgeschwindigkeit 405 einen Anteil der Auftriebskraft. 406 407 <p>Diese Implementierung ist noch nicht abgeschlossen. Die 160 einen Anteil der Auftriebskraft. <p>Diese 161 Implementierung ist noch nicht abgeschlossen. Die 408 162 Simulation von 409 163 Kaltluftabflüssen erfordert bisher noch Eingriffe "von Hand".</p> 410 411 </td> 412 413 <td>buoyancy, check_parameters, header, leap_frog, modules, 164 </td> <td>buoyancy, check_parameters, header, leap_frog, 165 modules, 414 166 parin, prandtl_fluxes, 415 read_var_list, write_var_list</td> 416 417 </tr> 418 419 <tr nosave="" valign="top"> 420 421 <td>15/12/98</td> 422 423 <td>SR</td> 424 425 <td>1.0</td> 426 427 <td>N</td> 428 429 <td>Vereinbarung einer Cache-Layer zur schnelleren binären 167 read_var_list, write_var_list</td> </tr> <tr nosave="" valign="top"> <td>15/12/98</td> <td>SR</td> 168 <td>1.0</td> <td>N</td> <td>Vereinbarung 169 einer Cache-Layer zur schnelleren binären 430 170 I/O auf 431 T3E-Rechnern. 432 433 <p>Ermittelung wahrer I/O-Zeiten für binäre 171 T3E-Rechnern. <p>Ermittelung wahrer I/O-Zeiten 172 für binäre 434 173 Ein-/Ausgabe auf 435 T3E-Rechnern mit TIMEF.</p> 436 437 </td> 438 439 <td nosave="">check_open, <br> 440 441 cpu_zeitmessung, init_3d_model, write_3d_model</td> 442 443 </tr> 444 445 <tr nosave="" valign="top"> 446 447 <td nosave="">"</td> 448 449 <td>SR</td> 450 451 <td>1.0</td> 452 453 <td>B</td> 454 455 <td>Bekanntgabe von run_description_header an alle PE's. 456 457 <p>Diese Programmänderung erfordert auf T3E-Rechnern 174 T3E-Rechnern mit TIMEF.</p> </td> <td nosave="">check_open, 175 <br>cpu_zeitmessung, init_3d_model, write_3d_model</td> </tr> 176 <tr nosave="" valign="top"> <td nosave="">"</td> 177 <td>SR</td> <td>1.0</td> <td>B</td> 178 <td>Bekanntgabe von run_description_header an alle 179 PE's. <p>Diese Programmänderung erfordert auf 180 T3E-Rechnern 458 181 zwingend den 459 182 Einsatz der MPI-Version mpt.1.2.0.1. In dieser neuen Version tritt auch 460 der alte Fehler in MPI_ALLREDUCE nicht mehr auf, so daß ab jetzt 461 mpt.1.1.0.1 nicht mehr verwendet werden muß.</p> 462 463 </td> 464 465 <td>parles</td> 466 467 </tr> 468 469 <tr nosave="" valign="top"> 470 471 <td nosave="">05/02/99</td> 472 473 <td>SR</td> 474 475 <td>1.1</td> 476 477 <td>C/N</td> 478 479 <td>Implementierung des Upstream-Spline- Verfahrens zur Advektion 183 der alte Fehler in MPI_ALLREDUCE nicht mehr auf, so daß ab 184 jetzt 185 mpt.1.1.0.1 nicht mehr verwendet werden muß.</p> </td> 186 <td>parles</td> </tr> <tr nosave="" valign="top"> <td nosave="">05/02/99</td> 187 <td>SR</td> <td>1.1</td> <td>C/N</td> 188 <td>Implementierung des Upstream-Spline- Verfahrens zur Advektion 480 189 von Impuls 481 und skalaren Größen (neue UP's, s. rechts). Berechnung von 190 und skalaren Größen (neue UP's, s. rechts). 191 Berechnung von 482 192 Variablen 483 für spezielle Advektionsverfahren in eigenständiges UP 193 für spezielle Advektionsverfahren in eigenständiges 194 UP 484 195 ausgelagert 485 (<tt>init_advec</tt> + eigenes Modul <tt>advection</tt>). Reines 196 (<tt>init_advec</tt> + eigenes Modul <tt>advection</tt>). 197 Reines 486 198 Euler-/Upstream-Verfahren 487 kann gerechnet werden. 488 489 <p>Schnelle binäre I/O (s. 15/12/98) wieder abgeschaltet, da 199 kann gerechnet werden. <p>Schnelle binäre I/O 200 (s. 15/12/98) wieder abgeschaltet, da 490 201 Dateien 491 nicht wieder einlesbar sind. Druckstörung (p) wird zu Beginn auf 0 202 nicht wieder einlesbar sind. Druckstörung (p) wird zu Beginn 203 auf 0 492 204 gesetzt, damit sie keine Tendenzen im ersten Zeitschritt hervorruft.</p> 493 494 </td> 495 496 <td>advec_s_bc, check_open, check_parameters, flow_statistics, 205 </td> <td>advec_s_bc, check_open, check_parameters, 206 flow_statistics, 497 207 header, 498 208 init_3d_model, init_grid, leap_frog, modules, parin, plot_2d, plot_ts, 499 read_var_list, timestep, write_3d_binary, write_var_list 500 501 <p><b>Neu:</b> <br> 502 209 read_var_list, timestep, write_3d_binary, write_var_list <p><b>Neu:</b> 210 <br> 503 211 advec_s_ups, advec_u_ups, advec_v_ups, advec_w_ups, init_advec, 504 212 long_filter, 505 spline_x, spline_y, spline_z</p> 506 507 </td> 508 509 </tr> 510 511 <tr nosave="" valign="top"> 512 513 <td nosave="">24/02/99</td> 514 515 <td>SR</td> 516 517 <td>1.1a</td> 518 519 <td>C/N</td> 520 521 <td>Kleinere Ergänzungen / Änderungen für Upstream 213 spline_x, spline_y, spline_z</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">24/02/99</td> 214 <td>SR</td> <td>1.1a</td> <td>C/N</td> 215 <td>Kleinere Ergänzungen / Änderungen 216 für Upstream 522 217 - Spline, <br> 523 524 218 Diverse Verbesserungen bei Berechnung der Diffusionskoeffizienten 525 219 (exaktere … … 529 223 Wertebegrenzung 530 224 bei Spline - Berechnung zur Verhinderung von "Überschwingern"</td> 531 532 <td>check_parameters, diffusion_e, diffusion_pt, diffusivities, 225 <td>check_parameters, diffusion_e, diffusion_pt, diffusivities, 533 226 header, 534 227 init_advec, long_filter, modules, production_e, spline_x, spline_y, 535 spline_z</td> 536 537 </tr> 538 539 <tr nosave="" valign="top"> 540 541 <td>"</td> 542 543 <td nosave="">SR</td> 544 545 <td>1.1a</td> 546 547 <td>N/B</td> 548 549 <td>Erweiterung des des Bott - Chlond - Schemas auf Advektion der 228 spline_z</td> </tr> <tr nosave="" valign="top"> 229 <td>"</td> <td nosave="">SR</td> <td>1.1a</td> 230 <td>N/B</td> <td>Erweiterung des des Bott - Chlond - 231 Schemas auf Advektion der 550 232 TKE, 551 233 Korrektur der Dichteberechnung (war falsch bei Galilei-Transformation)</td> 552 553 <td>advec_s_bc, leap_frog</td> 554 555 </tr> 556 557 <tr nosave="" valign="top"> 558 559 <td nosave="">01/03/99</td> 560 561 <td>MS</td> 562 563 <td>1.1b</td> 564 565 <td>N</td> 566 567 <td>Kleinere Ergänzungen für Upstream - Spline: 568 569 <ul> 570 571 <li>Gradientenkontrolle für den nicht-parallelen Teil</li> 572 573 <li>statistische Auswertung über den prozentualen Anteil 234 <td>advec_s_bc, leap_frog</td> </tr> <tr nosave="" valign="top"> <td nosave="">01/03/99</td> 235 <td>MS</td> <td>1.1b</td> <td>N</td> 236 <td>Kleinere Ergänzungen für Upstream - 237 Spline: <ul> <li>Gradientenkontrolle 238 für den nicht-parallelen Teil</li> <li>statistische 239 Auswertung über den prozentualen Anteil 574 240 des 575 241 Upstream-Verfahrens 576 an der Gesamtadvektion fuer nicht-parallelen Teil</li> 577 578 579 </ul> 580 581 </td> 582 583 <td>modules, spline_x, spline_y, spline_z</td> 584 585 </tr> 586 587 <tr nosave="" valign="top"> 588 589 <td>"</td> 590 591 <td>SR</td> 592 593 <td>1.1b</td> 594 595 <td>N</td> 596 597 <td>Datenkompression für 3D-Plotausgaben</td> 598 599 <td nosave="">check_open, check_parameters, close_files, header, 242 an der Gesamtadvektion fuer nicht-parallelen Teil</li> </ul> 243 </td> <td>modules, spline_x, spline_y, spline_z</td> 244 </tr> <tr nosave="" valign="top"> <td>"</td> 245 <td>SR</td> <td>1.1b</td> <td>N</td> 246 <td>Datenkompression für 3D-Plotausgaben</td> <td nosave="">check_open, check_parameters, close_files, header, 600 247 modules, 601 parin, plot_3d 602 603 <p><b>Neu:</b> <br> 604 605 write_compressed</p> 606 607 </td> 608 609 </tr> 610 611 <tr nosave="" valign="top"> 612 613 <td nosave="">03/03/99</td> 614 615 <td>SR</td> 616 617 <td>1.1c</td> 618 619 <td>B/C</td> 620 621 <td>Schnelle binäre I/O auf T3E-Rechnern funktioniert jetzt. 248 parin, plot_3d <p><b>Neu:</b> <br> 249 write_compressed</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">03/03/99</td> 250 <td>SR</td> <td>1.1c</td> <td>B/C</td> 251 <td>Schnelle binäre I/O auf T3E-Rechnern funktioniert 252 jetzt. 622 253 Beim 623 254 binären Einlesen dürfen Records nicht mit dem Trick "<tt>READ 624 (..) idum</tt>" überlesen werden. Zeitmessungen auf T3E erfolgen 255 (..) idum</tt>" überlesen werden. Zeitmessungen auf T3E 256 erfolgen 625 257 jetzt 626 grundsätzlich mit <tt>TIMEF.</tt></td> 627 628 <td>check_open, cpu_zeitmessung, init_3d_model, modules, 629 write_3d_binary</td> 630 631 </tr> 632 633 <tr nosave="" valign="top"> 634 635 <td nosave="">25/03/99</td> 636 637 <td>SR</td> 638 639 <td>1.1d</td> 640 641 <td>N</td> 642 643 <td>Einbau einer Rayleigh-Dämpfungsschicht sowie flexiblere 258 grundsätzlich mit <tt>TIMEF.</tt></td> <td>check_open, 259 cpu_zeitmessung, init_3d_model, modules, 260 write_3d_binary</td> </tr> <tr nosave="" valign="top"> <td nosave="">25/03/99</td> 261 <td>SR</td> <td>1.1d</td> <td>N</td> 262 <td>Einbau einer Rayleigh-Dämpfungsschicht sowie 263 flexiblere 644 264 Wahl der 645 265 numerischen Begrenzer für das Upstream - Spline - Verfahren.</td> 646 647 <td>check_parameters, header, init_3d_model, leap_frog, modules, 266 <td>check_parameters, header, init_3d_model, leap_frog, modules, 648 267 parin, 649 read_var_list, spline_x, spline_y, spline_z, write_var_list</td> 650 651 </tr> 652 653 <tr nosave="" valign="top"> 654 655 <td nosave="">"</td> 656 657 <td>SR</td> 658 659 <td>1.1d</td> 660 661 <td>B</td> 662 663 <td>Korrektur des Zeitschrittkriteriums bei z=zp im 1D-Modell</td> 664 665 <td>init_1d_model</td> 666 667 </tr> 668 669 <tr nosave="" valign="top"> 670 671 <td>25/11/99</td> 672 673 <td>SR</td> 674 675 <td>1.1e</td> 676 677 <td>N</td> 678 679 <td nosave="">Partikelplots mit vtk-Graphiksoftware sind auf t3eh 268 read_var_list, spline_x, spline_y, spline_z, write_var_list</td> </tr> 269 <tr nosave="" valign="top"> <td nosave="">"</td> 270 <td>SR</td> <td>1.1d</td> <td>B</td> 271 <td>Korrektur des Zeitschrittkriteriums bei z=zp im 1D-Modell</td> 272 <td>init_1d_model</td> </tr> <tr nosave="" valign="top"> <td>25/11/99</td> <td>SR</td> 273 <td>1.1e</td> <td>N</td> <td nosave="">Partikelplots 274 mit vtk-Graphiksoftware sind auf t3eh 680 275 möglich. 681 276 Ausgabe erfolgt im dvr-Format. Bisher noch nicht nutzbar, da noch 682 277 Einstellmöglichkeiten 683 für die Partikelquelle(n) fehlen. <b>Wichtig:</b> Die Benutzung 278 für die Partikelquelle(n) fehlen. <b>Wichtig:</b> 279 Die Benutzung 684 280 der 685 281 vtk-Software erfordert die zusätzliche mrun-Option <tt>-g</tt>. 686 282 Das 687 Hauptprogramm <tt>parles</tt> wird in diesem Fall zum Unterprogramm 283 Hauptprogramm <tt>parles</tt> wird in diesem Fall zum 284 Unterprogramm 688 285 und 689 von einem c++-Programm aus aufgerufen.</td> 690 691 <td>header, init_3d_model, init_pegrid, leap_frog, modules, 286 von einem c++-Programm aus aufgerufen.</td> <td>header, 287 init_3d_model, init_pegrid, leap_frog, modules, 692 288 parin, parles, 693 read_var_list, write_var_list 694 695 <p><b>Neu:</b> <br> 696 697 advec_particles, init_particles, plot_particles</p> 698 699 </td> 700 701 </tr> 702 703 <tr nosave="" valign="top"> 704 705 <td>"</td> 706 707 <td>SR</td> 708 709 <td>1.1e</td> 710 711 <td nosave="">B/C</td> 712 713 <td>TIMEF-Zeitmessungen von Subroutine- auf Function-Aufruf 289 read_var_list, write_var_list <p><b>Neu:</b> 290 <br> 291 advec_particles, init_particles, plot_particles</p> </td> </tr> 292 <tr nosave="" valign="top"> <td>"</td> 293 <td>SR</td> <td>1.1e</td> <td nosave="">B/C</td> 294 <td>TIMEF-Zeitmessungen von Subroutine- auf Function-Aufruf 714 295 umgestellt. 715 Laufindexkorrektur in <tt>spline_y</tt>. Fehler trat bei nx /= ny auf.</td> 716 717 <td>cpu_zeitmessung, spline_y</td> 718 719 </tr> 720 721 <tr nosave="" valign="top"> 722 723 <td nosave="">28/12/99</td> 724 725 <td>SR</td> 726 727 <td>1.1f</td> 728 729 <td>N</td> 730 731 <td>Plot von Isooberflächen mit vtk-Graphiksoftware sind auf 296 Laufindexkorrektur in <tt>spline_y</tt>. Fehler trat bei 297 nx /= ny auf.</td> <td>cpu_zeitmessung, spline_y</td> 298 </tr> <tr nosave="" valign="top"> <td nosave="">28/12/99</td> <td>SR</td> <td>1.1f</td> 299 <td>N</td> <td>Plot von Isooberflächen mit 300 vtk-Graphiksoftware sind auf 732 301 t3eh 733 302 möglich. Ausgabe erfolgt im dvr-Format. Bisher noch nicht 734 303 allgemein 735 304 nutzbar, da noch diverse manuelle Eingriffe nötig sind.</td> 736 737 <td>header, leap_frog, modules, parin, plot_particles, 738 read_var_list, write_var_list 739 740 <p><b>Neu:</b> <br> 741 742 plot_isosurface</p> 743 744 </td> 745 746 </tr> 747 748 <tr nosave="" valign="top"> 749 750 <td nosave="">10/01/00</td> 751 752 <td>SR</td> 753 754 <td>1.2</td> 755 756 <td>C</td> 757 758 <td>Randbedingungen für w (=0) werden explizit gesetzt, und 305 <td>header, leap_frog, modules, parin, plot_particles, 306 read_var_list, write_var_list <p><b>Neu:</b> 307 <br> 308 plot_isosurface</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">10/01/00</td> 309 <td>SR</td> <td>1.2</td> <td>C</td> 310 <td>Randbedingungen für w (=0) werden explizit gesetzt, 311 und 759 312 zwar nicht 760 nur bei <tt>k = nzb, nzt+1</tt> sondern auch bei <tt>k = 313 nur bei <tt>k = nzb, nzt+1</tt> sondern auch 314 bei <tt>k = 761 315 nzt</tt>, 762 316 wo bislang der vom vorigen Zeitschritt durch den Drucklöser … … 764 318 Wert stand. Auch jetzt wird der Wert bei k = nzt anschließend 765 319 durch 766 den Drucklöser modifiziert. Soweit bis jetzt abzusehen, hat diese 767 Änderung nur minimale Auswirkungen auf den Prognoseverlauf. Evtl. 320 den Drucklöser modifiziert. Soweit bis jetzt abzusehen, hat 321 diese 322 Änderung nur minimale Auswirkungen auf den Prognoseverlauf. 323 Evtl. 768 324 gibt es stärkere Auswirkungen bei exzessivem Auftreten von 769 325 Schwerewellen 770 326 am oberen Rand (die dort aber eigentlich sowieso nicht 771 327 hingehören). <br> 772 773 Variablenübergabe an <tt>boundary_conds</tt> erfolgt <br> 774 775 über Modul und nicht mehr per Parameterliste. 776 777 <p>Bei Galilei-Transformation wird jetzt defaultmäßig 328 Variablenübergabe an <tt>boundary_conds</tt> 329 erfolgt <br> 330 über Modul und nicht mehr per Parameterliste. <p>Bei 331 Galilei-Transformation wird jetzt defaultmäßig 778 332 der geostrophische 779 333 Wind als Translationsgeschwindigkeit verwendet (neuer … … 782 336 höheren 783 337 geostrophischen Windgeschwindigkeiten unterdrückt. </p> 784 785 786 <p>3D-Feldvariablen (u, v, w, pt, e, km, kh) sowie einige 338 <p>3D-Feldvariablen (u, v, w, pt, e, km, kh) sowie einige 787 339 2D-Variablen 788 340 (usws, vsws, rif, shf) sind nun als Pointer deklariert, denen nach … … 792 344 3 Zeitebenen existieren die Targets ..._1 (z.B. u_1), ..._2, ..._3). 793 345 Durch 794 die jeweilige Neuzuordnung der Targets entfällt das Umspeichern 346 die jeweilige Neuzuordnung der Targets entfällt das 347 Umspeichern 795 348 der 796 Zeitebenen komplett (Zeitersparung: ca. 5% der Gesamtrechenzeit). 797 <br> 798 349 Zeitebenen komplett (Zeitersparung: ca. 5% der 350 Gesamtrechenzeit). <br> 799 351 Die Umordnung der Zeitebenen wurde ins neue UP <tt>swap_timelevel</tt> 800 ausgelagert. <br> 801 802 <b>Wichtig:</b> Die Übersetzungszeiten haben durch die 352 ausgelagert. <br> <b>Wichtig:</b> Die 353 Übersetzungszeiten haben durch die 803 354 Einführung 804 355 der Pointer t.w. sehr deutlich zugenommen (<tt>leap_frog</tt> 805 356 benötigt 806 jetzt z.B. ca. 60 Sekunden). Außerdem benötigt der Compiler 357 jetzt z.B. ca. 60 Sekunden). Außerdem benötigt der 358 Compiler 807 359 jetzt wesentlich mehr Speicher. Der "memory"-Parameter in den 808 360 Konfigurationsdateien 809 361 (.mrun.config, .al.config) muss nun mindestens den Wert 60 (MB) (t3eh, 810 t3ej2) bzw. 45 (t3eb) und 96 (vpp) haben. <br> 811 812 <b>Achtung:</b> bislang wurde diese Änderung nur für 362 t3ej2) bzw. 45 (t3eb) und 96 (vpp) haben. <br> <b>Achtung:</b> 363 bislang wurde diese Änderung nur für 813 364 das 814 reine Leapfrog-Verfahren getestet! </p> 815 816 817 <p>Neue Laderoption "<tt>-D preset=nan</tt>" empfohlen (s. 365 reine Leapfrog-Verfahren getestet! </p> <p>Neue 366 Laderoption "<tt>-D preset=nan</tt>" empfohlen (s. 818 367 Beispiel-Konfigurationsdatei). 819 Ebenso empfehlenswert für Debug-Läufe: Compileroption "<tt>-ei</tt>". 368 Ebenso empfehlenswert für Debug-Läufe: Compileroption 369 "<tt>-ei</tt>". 820 370 Diese bewirkt einen Programmabbruch, falls mit nicht initialisierten 821 371 Variablen 822 372 gearbeitet wird (ist aber zeitaufwendig in der Ausführung!).</p> 823 824 </td> 825 826 <td>boundary_conds, check_parameters, header, init_3d_model, 373 </td> <td>boundary_conds, check_parameters, header, 374 init_3d_model, 827 375 leap_frog, 828 modules, parin, read_var_list, timestep, write_var_list 829 830 <p><b>Neu:</b> <br> 831 832 swap_timelevel</p> 833 834 </td> 835 836 </tr> 837 838 <tr nosave="" valign="top"> 839 840 <td nosave="">10/01/00</td> 841 842 <td>SR</td> 843 844 <td>1.2</td> 845 846 <td>B</td> 847 848 <td>Der Asselin-Filter wird vor Umordnung der Zeitebenen 376 modules, parin, read_var_list, timestep, write_var_list <p><b>Neu:</b> 377 <br> 378 swap_timelevel</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">10/01/00</td> 379 <td>SR</td> <td>1.2</td> <td>B</td> 380 <td>Der Asselin-Filter wird vor Umordnung der Zeitebenen 849 381 aufgerufen. Dies 850 382 geschah bisher fälschlicherweise erst danach, sodass der 851 383 Zeitfilter 852 asymmetrisch arbeitete. <br> 853 854 <b>Achtung:</b> diese Änderung hat Auswirkungen auf den 384 asymmetrisch arbeitete. <br> <b>Achtung:</b> 385 diese Änderung hat Auswirkungen auf den 855 386 Prognoseverlauf. 856 387 Die Testergebnis-Datei wurde entsprechend geändert.</td> 857 858 <td>leap_frog</td> 859 860 </tr> 861 862 <tr nosave=""> 863 864 <td nosave="" valign="top">14/02/00</td> 865 866 <td nosave="" valign="top">MS</td> 867 868 <td nosave="" valign="top">1.2</td> 869 870 <td nosave="" valign="top">B/C</td> 871 872 <td><b>1D-Modell:</b> <br> 873 388 <td>leap_frog</td> </tr> <tr nosave=""> 389 <td nosave="" valign="top">14/02/00</td> <td nosave="" valign="top">MS</td> <td nosave="" valign="top">1.2</td> <td nosave="" valign="top">B/C</td> <td><b>1D-Modell:</b> 390 <br> 874 391 Der Asselin-Filter im 1D-Modell wird nun auch vor der Umordnung der 875 Zeitebenen aufgerufen (s.o.). 876 877 <p><b>3D-Modell:</b> <br> 878 879 TKE wurde bisher bei <tt>scalar_advec=bc-scheme</tt> einer zeitlichen 880 Filterung unterzogen. Dies wird nun unterbunden. </p> 881 882 883 <p>Die Divergenz der Impulsflüsse an der Prandtl - Schicht - 392 Zeitebenen aufgerufen (s.o.). <p><b>3D-Modell:</b> 393 <br> 394 TKE wurde bisher bei <tt>scalar_advec=bc-scheme</tt> einer 395 zeitlichen 396 Filterung unterzogen. Dies wird nun unterbunden. </p> <p>Die 397 Divergenz der Impulsflüsse an der Prandtl - Schicht - 884 398 Obergrenze 885 399 wird nun über eine ganze Gitterweite gebildet. Die bisherige … … 887 401 lieferte im Vergleich mit Ergebnissen anderer LES-Modelle zu 888 402 große 889 Impulsflüsse. </p> 890 891 892 <p>Bei der Berechnung der Diffusion der Temperatur werden die 403 Impulsflüsse. </p> <p>Bei der 404 Berechnung der Diffusion der Temperatur werden die 893 405 Diffusionskoeffizienten 894 406 nun wieder durch Mittelung für die entsprechenden … … 898 410 Upstream-Spline Advektionsverfahrens eingebaut wurde, bewirkte ein zu 899 411 rasches 900 Anwachsen der Grenzschicht mit der Zeit. </p> 901 902 903 <p><b>Achtung:</b> diese Änderungen haben Auswirkungen auf 412 Anwachsen der Grenzschicht mit der Zeit. </p> <p><b>Achtung:</b> 413 diese Änderungen haben Auswirkungen auf 904 414 den Prognoseverlauf. 905 Die Testergebnis-Datei wurde entsprechend geändert.</p> 906 907 </td> 908 909 <td nosave="" valign="top">init_1d_model, asselin_filter, 415 Die Testergebnis-Datei wurde entsprechend geändert.</p> </td> 416 <td nosave="" valign="top">init_1d_model, 417 asselin_filter, 910 418 diffusion_u, diffusion_v, <br> 911 912 diffusion_pt</td> 913 914 </tr> 915 916 <tr nosave="" valign="top"> 917 918 <td nosave="">17/02/00</td> 919 920 <td>MS</td> 921 922 <td>2.0</td> 923 924 <td>N</td> 925 926 <td><b>Allgemeine Beschreibung der Änderungen:</b> <br> 927 928 vollständige Implementierung des wolkenphysikalischen Moduls 929 930 <p>Zur Implementierung der Wolkenphysik waren umfangreiche 419 diffusion_pt</td> </tr> <tr nosave="" valign="top"> <td nosave="">17/02/00</td> 420 <td>MS</td> <td>2.0</td> <td>N</td> 421 <td><b>Allgemeine Beschreibung der Änderungen:</b> 422 <br> 423 vollständige Implementierung des wolkenphysikalischen 424 Moduls <p>Zur Implementierung der Wolkenphysik waren 425 umfangreiche 931 426 Änderungen 932 427 in zahlreichen Programmteilen notwendig (Beschreibung <a href="#Cloud_physics">siehe 933 428 unten</a> ). Im einzelnen wurden folgende Unprogramme 934 geändert: <br> 935 936 <b>modules</b>: <br> 937 938 1) Bereitstellung neuer 3D-, 2D- und 1D-Felder für die Prognose 429 geändert: <br> <b>modules</b>: 430 <br>1) Bereitstellung neuer 3D-, 2D- und 1D-Felder für 431 die Prognose 939 432 des Gesamtwassergehaltes <br> 940 941 433 2) "cloud_parameters" als neues Modul mit wolkenphysikalischen 942 Variablen <br> 943 944 <b>parin: </b>Initialisierungsparameterliste wurde erweitert 945 (s.u.) <br> 946 947 <font color="#000000"><b>read_var_list, write_var_list, 434 Variablen <br> <b>parin: </b>Initialisierungsparameterliste 435 wurde erweitert 436 (s.u.) <br> <font color="#000000"><b>read_var_list, 437 write_var_list, 948 438 write_3d_binary: </b>In/Output 949 439 der neuen Initialisierungsparamter und Felder für 950 Fortsetzungsläufe </font> <br> 951 952 <font color="#000000"><b>check_parameters:</b> Initialisierung 440 Fortsetzungsläufe </font> <br> <font color="#000000"><b>check_parameters:</b> 441 Initialisierung 953 442 des 954 Feuchteprofils und Ausgabe von neuen Vertikalprofilen: </font> <br> 955 956 <font color="#ff0000">vpt</font><font color="#000000">: virtuelle 443 Feuchteprofils und Ausgabe von neuen Vertikalprofilen: </font> 444 <br> <font color="#ff0000">vpt</font><font color="#000000">: virtuelle 957 445 poteniellen 958 Temperatur,</font> <br> 959 960 <font color="#ff0000">lpt</font><font color="#000000">: 446 Temperatur,</font> <br> <font color="#ff0000">lpt</font><font color="#000000">: 961 447 potentielle 962 Flüssigwassertemperatur,</font> <br> 963 964 <font color="#ff0000">q</font><font color="#000000">: 965 Gesamtwassergehalt,</font> <br> 966 967 <font color="#ff0000">qv</font><font color="#000000">: 448 Flüssigwassertemperatur,</font> <br> <font color="#ff0000">q</font><font color="#000000">: 449 Gesamtwassergehalt,</font> <br> <font color="#ff0000">qv</font><font color="#000000">: 968 450 spezifische 969 Feuchte,</font> <br> 970 971 <font color="#ff0000">ql</font><font color="#000000">: 972 Flüssigwassergehalt,</font> <br> 973 974 <font color="#ff0000">w"vpt", w*vpt*, wvpt</font><font color="#000000">: 451 Feuchte,</font> <br> <font color="#ff0000">ql</font><font color="#000000">: 452 Flüssigwassergehalt,</font> <br> <font color="#ff0000">w"vpt", w*vpt*, wvpt</font><font color="#000000">: 975 453 Auftriebsfluss (subskalig, aufgelöst, gesamt),</font> <br> 976 977 <font color="#ff0000">w"q", w*q*, wq</font><font color="#000000">: 978 Wasserfluss (subs., aufg., ges.),</font> <br> 979 980 <font color="#ff0000">w"qv", w*qv*, wqv</font><font color="#000000">: 981 Feuchtefluss (subs, aufg., ges.)</font> <br> 982 983 <b>init_3d_model</b>: Initialisierung aller 984 Feuchte/Wasserfelder <br> 985 986 <b>init_1d_model</b>: Anpassung an Rechnungen mit Feuchte <br> 987 988 <b>header</b>: Headerausgaben für Feuchtebehandlung und 454 <font color="#ff0000">w"q", w*q*, wq</font><font color="#000000">: 455 Wasserfluss (subs., aufg., ges.),</font> <br> <font color="#ff0000">w"qv", w*qv*, wqv</font><font color="#000000">: 456 Feuchtefluss (subs, aufg., ges.)</font> <br> <b>init_3d_model</b>: 457 Initialisierung aller 458 Feuchte/Wasserfelder <br> <b>init_1d_model</b>: 459 Anpassung an Rechnungen mit Feuchte <br> <b>header</b>: 460 Headerausgaben für Feuchtebehandlung und 989 461 Wolkenphysik 990 462 eingefügt, außerdem Formatnummern erweitert um Platz 991 463 für 992 weitere Ausgaben zu schaffen <br> 993 994 <b>leap_frog:</b> <br> 995 464 weitere Ausgaben zu schaffen <br> <b>leap_frog:</b> 465 <br> 996 466 1) alle prognostischen Gleichungen wurden in prognostic_equations 997 467 ausgelagert <br> 998 999 468 2) Aufrufe der neuen Unterprogramme (UP) compute_vpt, 1000 469 calc_liquid_water_content <br> 1001 1002 3) veraenderter Aufruf von diffusivities <br> 1003 1004 <b><tt>b</tt>oundary_conds, asselin_filter, swap_timelevel, 470 3) veraenderter Aufruf von diffusivities <br> <b><tt>b</tt>oundary_conds, 471 asselin_filter, swap_timelevel, 1005 472 flow_statistics: </b>Anpassung 1006 an Feuchtebehandlung <br> 1007 1008 <b>plot2d, plot3d</b>: Ausgabemöglichkeiten von q und ql <br> 1009 1010 <b>buoyancy</b>: Anpassung des Autriebstermes an 473 an Feuchtebehandlung <br> <b>plot2d, plot3d</b>: 474 Ausgabemöglichkeiten von q und ql <br> <b>buoyancy</b>: 475 Anpassung des Autriebstermes an 1011 476 Feuchtebehandlung 1012 477 durch neue Übergabeparameter: Temperatur + Nummer des 1013 478 entsprechenden 1014 mittleren Temperaturprofils <br> 1015 1016 <b>diffusion_e</b>: Anpassung der Mischungswegberechnung an 479 mittleren Temperaturprofils <br> <b>diffusion_e</b>: 480 Anpassung der Mischungswegberechnung an 1017 481 Feuchtebehandlung 1018 durch Variation der Übergabeparameter <br> 1019 1020 <b>production_e</b>: TKE-Produktion durch Auftriebskräfte an 482 durch Variation der Übergabeparameter <br> <b>production_e</b>: 483 TKE-Produktion durch Auftriebskräfte an 1021 484 Feuchtebehandlung 1022 angepasst <br> 1023 1024 <b>diffusivities</b>: siehe diffusion_e <br> 1025 1026 <b>prandtl-fluxes</b>: Prandtl-Schicht fuer den 485 angepasst <br> <b>diffusivities</b>: 486 siehe diffusion_e <br> <b>prandtl-fluxes</b>: 487 Prandtl-Schicht fuer den 1027 488 Gesamtwassergehalt 1028 eingefügt </p> 1029 1030 1031 <p>Zusätzlich wurden folgende Unterprogramme NEU 1032 eingeführt: <br> 1033 1034 <b>inti_cloud_physics</b>: Initialisierung wolkenphysikalischer 1035 Parameter <br> 1036 1037 <b>prognostic_equations</b>: Auslagerung der prognostischen 489 eingefügt </p> <p>Zusätzlich 490 wurden folgende Unterprogramme NEU 491 eingeführt: <br> <b>inti_cloud_physics</b>: 492 Initialisierung wolkenphysikalischer 493 Parameter <br> <b>prognostic_equations</b>: 494 Auslagerung der prognostischen 1038 495 Gleichungen 1039 aus leap_frog, zusätzliche porgnostische Gleichung für den1040 Gesamtwassergehalt <br> 1041 1042 <b>diffusion_s</b>:Berechnung der Diffusion skalarer496 aus leap_frog, zusätzliche porgnostische Gleichung 497 für den 498 Gesamtwassergehalt <br> <b>diffusion_s</b>: 499 Berechnung der Diffusion skalarer 1043 500 Größen 1044 501 (Temperatur und Gesamtwassergehalt), ersetzt diffusion_pt <br> 1045 1046 <b>calc_liquid_water_content</b>: Berechnung des 1047 Flüssigwassergehalts <br> 1048 1049 <b>calc_radiation</b>: Parametrisierung langwelliger 502 <b>calc_liquid_water_content</b>: Berechnung des 503 Flüssigwassergehalts <br> <b>calc_radiation</b>: 504 Parametrisierung langwelliger 1050 505 Strahlungsprozesse 1051 506 über des Schema der effektiven Emissivität <br> 1052 1053 <b>calc_precipitation</b>:Parametrisierung von 507 <b>calc_precipitation</b>:Parametrisierung von 1054 508 Niederschlagsprozessen 1055 mit Hilfe eines vereinfachten Kesslerschemas <br> 1056 1057 <b>impact_of_latent_heat</b>: Berücksichtigung der durch 509 mit Hilfe eines vereinfachten Kesslerschemas <br> <b>impact_of_latent_heat</b>: 510 Berücksichtigung der durch 1058 511 Niederschlagsprozesse 1059 512 zurückbleibenden latenten Waerme innerhalb eines 1060 Gittervolumens <br> 1061 1062 <b>compute_vpt</b>: Berechnung der virtuellen potentiellen 1063 Temperatur </p> 1064 1065 1066 <p><b>Liste der neuen Initialisierungsparameter:</b> <br> 1067 1068 1069 <table nosave="" border="1"> 1070 1071 <caption> <br> 1072 1073 </caption><tbody> 1074 1075 1076 </tbody><tbody> 1077 1078 1079 </tbody> <tbody> 1080 1081 <tr nosave=""> 1082 1083 <td nosave=""><b>Parameter</b></td> 1084 1085 <td nosave=""><b>Typ</b></td> 1086 1087 <td nosave="" width="4%"><b>Default</b></td> 1088 1089 <td nosave=""><b>Beschreibung</b></td> 1090 1091 </tr> 1092 1093 <tr nosave=""> 1094 1095 <td nosave=""><tt>moisture</tt></td> 1096 1097 <td>L</td> 1098 1099 <td>FALSE</td> 1100 1101 <td>Ab-/Zuschalten des Gesamtwassergerhaltes als neue 1102 progn. Variable</td> 1103 1104 </tr> 1105 1106 <tr> 1107 1108 <td><tt>cloud_physics</tt></td> 1109 1110 <td>L</td> 1111 1112 <td>FALSE</td> 1113 1114 <td>Ab-/Zuschalten des Kondensationsschemas</td> 1115 1116 </tr> 1117 1118 <tr nosave=""> 1119 1120 <td><tt>radiation</tt></td> 1121 1122 <td>L</td> 1123 1124 <td nosave="">FALSE</td> 1125 1126 <td>Ab-/Zuschalten des Strahlungsschemas</td> 1127 1128 </tr> 1129 1130 <tr> 1131 1132 <td><tt>precipitation</tt></td> 1133 1134 <td>L</td> 1135 1136 <td>FALSE</td> 1137 1138 <td>Ab-/Zuschalten der Niderschlagsparametrisierung</td> 1139 1140 </tr> 1141 1142 <tr> 1143 1144 <td><tt>bc_q_b</tt></td> 1145 1146 <td>C</td> 1147 1148 <td>'dirichlet'</td> 1149 1150 <td>untere Randbedingung für q (siehe auch bc_pt_b)</td> 1151 1152 </tr> 1153 1154 <tr> 1155 1156 <td><tt>bc_q_t</tt></td> 1157 1158 <td>C</td> 1159 1160 <td>'neumann'</td> 1161 1162 <td>obere RB für q (siehe auch bc_pt_t)</td> 1163 1164 </tr> 1165 1166 <tr> 1167 1168 <td><tt>q_surface</tt></td> 1169 1170 <td>R</td> 1171 1172 <td>0.0</td> 1173 1174 <td>Feuchtewert an der Erdoberfläche</td> 1175 1176 </tr> 1177 1178 <tr> 1179 1180 <td><tt>q_surface_initial_change</tt></td> 1181 1182 <td>R</td> 1183 1184 <td>0.0</td> 1185 1186 <td>vgl. pt_surface_initial_change</td> 1187 1188 </tr> 1189 1190 <tr> 1191 1192 <td><tt>q_vertical_gradient</tt></td> 1193 1194 <td>R(10)</td> 1195 1196 <td>0.0</td> 1197 1198 <td>vgl. pt_vertical_gradient</td> 1199 1200 </tr> 1201 1202 <tr nosave=""> 1203 1204 <td nosave=""><tt>q_vertical_gradient_level</tt></td> 1205 1206 <td>R(10)</td> 1207 1208 <td>10000.0</td> 1209 1210 <td>vgl. pt_vertical_gradient_level</td> 1211 1212 </tr> 1213 1214 <tr> 1215 1216 <td><tt>surface_pressure</tt></td> 1217 1218 <td>R</td> 1219 1220 <td>1013.25</td> 1221 1222 <td>Luftdruck an der Erdoberfläche</td> 1223 1224 </tr> 1225 1226 <tr> 1227 1228 <td><tt>surface_waterflux</tt></td> 1229 1230 <td>R</td> 1231 1232 <td>0.0</td> 1233 1234 <td>oberflächennaher Wasser/Feuchtefluss</td> 1235 1236 </tr> 1237 1238 1239 </tbody> 1240 1241 </table> 1242 1243 Fortsetzungsläufe mit Daten, die mit Modellversionen kleiner 2.0 513 Gittervolumens <br> <b>compute_vpt</b>: 514 Berechnung der virtuellen potentiellen 515 Temperatur </p> <p><b>Liste der neuen 516 Initialisierungsparameter:</b> <br> <table nosave="" border="1"> <caption> <br> </caption><tbody> 517 </tbody><tbody> </tbody> <tbody> <tr nosave=""> <td nosave=""><b>Parameter</b></td> 518 <td nosave=""><b>Typ</b></td> <td nosave="" width="4%"><b>Default</b></td> 519 <td nosave=""><b>Beschreibung</b></td> </tr> 520 <tr nosave=""> <td nosave=""><tt>moisture</tt></td> 521 <td>L</td> <td>FALSE</td> <td>Ab-/Zuschalten 522 des Gesamtwassergerhaltes als neue 523 progn. Variable</td> </tr> <tr> <td><tt>cloud_physics</tt></td> 524 <td>L</td> <td>FALSE</td> <td>Ab-/Zuschalten 525 des Kondensationsschemas</td> </tr> <tr nosave=""> 526 <td><tt>radiation</tt></td> <td>L</td> 527 <td nosave="">FALSE</td> <td>Ab-/Zuschalten 528 des Strahlungsschemas</td> </tr> <tr> <td><tt>precipitation</tt></td> 529 <td>L</td> <td>FALSE</td> <td>Ab-/Zuschalten 530 der Niderschlagsparametrisierung</td> </tr> <tr> <td><tt>bc_q_b</tt></td> 531 <td>C</td> <td>'dirichlet'</td> <td>untere 532 Randbedingung für q (siehe auch bc_pt_b)</td> </tr> 533 <tr> <td><tt>bc_q_t</tt></td> <td>C</td> 534 <td>'neumann'</td> <td>obere RB für q 535 (siehe auch bc_pt_t)</td> </tr> <tr> <td><tt>q_surface</tt></td> 536 <td>R</td> <td>0.0</td> <td>Feuchtewert 537 an der Erdoberfläche</td> </tr> <tr> <td><tt>q_surface_initial_change</tt></td> 538 <td>R</td> <td>0.0</td> <td>vgl. 539 pt_surface_initial_change</td> </tr> <tr> <td><tt>q_vertical_gradient</tt></td> 540 <td>R(10)</td> <td>0.0</td> <td>vgl. 541 pt_vertical_gradient</td> </tr> <tr nosave=""> 542 <td nosave=""><tt>q_vertical_gradient_level</tt></td> 543 <td>R(10)</td> <td>10000.0</td> <td>vgl. 544 pt_vertical_gradient_level</td> </tr> <tr> <td><tt>surface_pressure</tt></td> 545 <td>R</td> <td>1013.25</td> <td>Luftdruck 546 an der Erdoberfläche</td> </tr> <tr> <td><tt>surface_waterflux</tt></td> 547 <td>R</td> <td>0.0</td> <td>oberflächennaher 548 Wasser/Feuchtefluss</td> </tr> </tbody> </table> 549 Fortsetzungsläufe mit Daten, die mit Modellversionen kleiner 550 2.0 1244 551 erzeugt 1245 552 wurden, sind aufgrund der erweiterten Initialisierungsparameterliste … … 1247 554 zu den Einleseroutinen der Version 2.0. Ein Update auf Version 2.0 1248 555 sollte 1249 deshalb nach Abschluss eines kompletten Modelllaufs erfolgen. </p> 1250 1251 1252 <p><b>Achtung:</b> diese Änderungen haben Auswirkungen auf 556 deshalb nach Abschluss eines kompletten Modelllaufs erfolgen. 557 </p> <p><b>Achtung:</b> diese 558 Änderungen haben Auswirkungen auf 1253 559 den Prognoseverlauf. 1254 Die Testergebnis-Datei wurde entsprechend geändert.</p> 1255 1256 </td> 1257 1258 <td nosave="" valign="top">modules, parin, read_var_list, 560 Die Testergebnis-Datei wurde entsprechend geändert.</p> </td> 561 <td nosave="" valign="top">modules, parin, 562 read_var_list, 1259 563 check_parameters, 1260 564 init_3d_model, init_1d_model, header, leap_frog. buoyancy, diffusion_e, … … 1262 566 asselin_filter, 1263 567 swap_timelevels, flow_statistics, plot_2d, plot_3d, write_var_list, 1264 write_3d_binary 1265 1266 <p><b>Neu:</b> <br> 1267 568 write_3d_binary <p><b>Neu:</b> <br> 1268 569 init_cloud_physics, prognostic_equations, diffusion_s, 1269 570 calc_liquid_water_content, 1270 571 calc_radiation, calc_precipitation, impact_of_latent_heat, comute_vpt</p> 1271 1272 </td> 1273 1274 </tr> 1275 1276 <tr nosave=""> 1277 1278 <td nosave="" valign="top">18/02/00</td> 1279 1280 <td nosave="" valign="top">MS</td> 1281 1282 <td nosave="" valign="top">2.0</td> 1283 1284 <td nosave="" valign="top">B/N</td> 1285 1286 <td>Im Zuge der Umstellung auf Modellversion 2.0 wurden 572 </td> </tr> <tr nosave=""> <td nosave="" valign="top">18/02/00</td> <td nosave="" valign="top">MS</td> <td nosave="" valign="top">2.0</td> <td nosave="" valign="top">B/N</td> <td>Im Zuge der 573 Umstellung auf Modellversion 2.0 wurden 1287 574 noch kleine 1288 Änderungen an den Unterprogrammen diffusion_e, production_e und 575 Änderungen an den Unterprogrammen diffusion_e, production_e 576 und 1289 577 diffusivities 1290 578 vorgenommen. In allen genannten <br> 1291 1292 579 Unterprogrammen wurde bisher ein strengeres Kriterium zur Bestimmung 1293 der Schichtungsverhältnisse verwendet. Ein Modellvergleich zeigte 580 der Schichtungsverhältnisse verwendet. Ein Modellvergleich 581 zeigte 1294 582 jedoch, dass dadurch die Diffusion im Bereich der Inversion zu gross 1295 wird. 1296 1297 <p>Eine weitere Ergänzung betrifft die 583 wird. <p>Eine weitere Ergänzung betrifft die 1298 584 Advektionsverfahren: <br> 1299 1300 Das Bott-Chlond Advektionsverfahren ist nun auch für Rechnungen 1301 mit Feuchte/Wolkenphysik verfügbar. </p> 1302 1303 1304 <p><b>Achtung:</b> diese Änderungen haben Auswirkungen auf 585 Das Bott-Chlond Advektionsverfahren ist nun auch für 586 Rechnungen 587 mit Feuchte/Wolkenphysik verfügbar. </p> <p><b>Achtung:</b> 588 diese Änderungen haben Auswirkungen auf 1305 589 den Prognoseverlauf. 1306 Die Testergebnis-Datei wurde entsprechend geändert.</p> 1307 1308 </td> 1309 1310 <td nosave="" valign="top">diffusion_e, production_e, 1311 diffusivities, check_parameters</td> 1312 1313 </tr> 1314 1315 <tr nosave="" valign="top"> 1316 1317 <td nosave="">26/04/00</td> 1318 1319 <td>SR</td> 1320 1321 <td>2.0a</td> 1322 1323 <td>C</td> 1324 1325 <td>Vollständige Umstellung von vtk-Grafik-Software auf 590 Die Testergebnis-Datei wurde entsprechend geändert.</p> </td> 591 <td nosave="" valign="top">diffusion_e, 592 production_e, 593 diffusivities, check_parameters</td> </tr> <tr nosave="" valign="top"> <td nosave="">26/04/00</td> 594 <td>SR</td> <td>2.0a</td> <td>C</td> 595 <td>Vollständige Umstellung von vtk-Grafik-Software auf 1326 596 dvrp-Software. 1327 597 Grafik-Ausgaben im dvr-Format sind mittlerweile begrenzt nutzbar 1328 598 (Isooberflächen 1329 sind darstellbar, siehe neue Laufparameter <tt>dt_dvrp, threshold</tt>), 599 sind darstellbar, siehe neue Laufparameter <tt>dt_dvrp, 600 threshold</tt>), 1330 601 es wird aber in der Anwendung in naher Zukunft noch diverse 1331 602 Änderungen 1332 603 geben. Partikelausgabe vorerst nicht mehr möglich (Verlagerung 1333 schon).</td> 1334 1335 <td>header, init_particles, init_3d_model, leap_frog, modules, 604 schon).</td> <td>header, init_particles, init_3d_model, 605 leap_frog, modules, 1336 606 parin, parles, 1337 read_var_list, write_var_list <br> 1338 1339 <b>Gestrichen:</b> <br> 1340 1341 plot_isosurface, plot_particles 1342 1343 <p><b>Neu:</b> <br> 1344 1345 init_dvrp, plot_dvrp</p> 1346 1347 </td> 1348 1349 </tr> 1350 1351 <tr nosave="" valign="top"> 1352 1353 <td nosave="">26/04/00</td> 1354 1355 <td>SR</td> 1356 1357 <td>2.0a</td> 1358 1359 <td>C</td> 1360 1361 <td>Durch zusätzlichen zyklischen Rand im Bott-Chlond-Schema 607 read_var_list, write_var_list <br> <b>Gestrichen:</b> 608 <br> 609 plot_isosurface, plot_particles <p><b>Neu:</b> 610 <br> 611 init_dvrp, plot_dvrp</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">26/04/00</td> 612 <td>SR</td> <td>2.0a</td> <td>C</td> 613 <td>Durch zusätzlichen zyklischen Rand im 614 Bott-Chlond-Schema 1362 615 sehr 1363 616 viele SENDRECV-Aufrufe eingespart. <br> 1364 1365 Änderung des Namens einer eingelesenen Environment-Variablen in <tt>check_open</tt> 1366 von <tt>remote_addres</tt> nach <tt>return_addres</tt>. <br> 1367 1368 Prozessor-Topologie kann durch Benutzer vorgegeben werden. Dazu neue 1369 Initialisierungsparameter <tt>npex</tt> und <tt>npey</tt>. Durch 617 Änderung des Namens einer eingelesenen Environment-Variablen 618 in <tt>check_open</tt> 619 von <tt>remote_addres</tt> nach <tt>return_addres</tt>. 620 <br>Prozessor-Topologie kann durch Benutzer vorgegeben werden. 621 Dazu neue 622 Initialisierungsparameter <tt>npex</tt> und <tt>npey</tt>. 623 Durch 1370 624 Änderungen 1371 625 am Header-Format wurde Testergebnis-Datei geändert.</td> 1372 1373 <td>advec_s_bc, check_open, header, init_pegrid, modules, parin, 626 <td>advec_s_bc, check_open, header, init_pegrid, modules, parin, 1374 627 read_var_list, 1375 write_var_list</td> 1376 1377 </tr> 1378 1379 <tr nosave="" valign="top"> 1380 1381 <td nosave="">26/04/00</td> 1382 1383 <td>SR</td> 1384 1385 <td>2.0a</td> 1386 1387 <td>N</td> 1388 1389 <td>Rechnungen mit geneigter Oberfläche erlauben jetzt die 628 write_var_list</td> </tr> <tr nosave="" valign="top"> <td nosave="">26/04/00</td> 629 <td>SR</td> <td>2.0a</td> <td>N</td> 630 <td>Rechnungen mit geneigter Oberfläche erlauben jetzt 631 die 1390 632 Vorgabe 1391 633 eines stabil geschichteten Temperaturprofils mit konstantem Gradienten. 1392 634 Hinzuschalten von Feuchte funktioniert für Hangrechnungen noch 1393 nicht.</td> 1394 1395 <td>advec_s_bc, buoyancy, init_3d_model, modules, 1396 prognostic_equations 1397 1398 <p><b>Neu:</b> <br> 1399 1400 init_slope</p> 1401 1402 </td> 1403 1404 </tr> 1405 1406 <tr nosave="" valign="top"> 1407 1408 <td>13/06/00</td> 1409 1410 <td>MS</td> 1411 1412 <td>2.0a</td> 1413 1414 <td nosave="">B</td> 1415 1416 <td>Es werden nun auch die seitlichen Ränder für ql 635 nicht.</td> <td>advec_s_bc, buoyancy, init_3d_model, 636 modules, 637 prognostic_equations <p><b>Neu:</b> <br> 638 init_slope</p> </td> </tr> <tr nosave="" valign="top"> <td>13/06/00</td> <td>MS</td> 639 <td>2.0a</td> <td nosave="">B</td> <td>Es 640 werden nun auch die seitlichen Ränder für ql 1417 641 gesetzt. 1418 642 Ein Nichtsetzen der Randbedingungen führte zu Fehlern bei der 1419 643 Ausgabe 1420 von pt</td> 1421 1422 <td nosave="">calc_liquid_water_content</td> 1423 1424 </tr> 1425 1426 <tr nosave="" valign="top"> 1427 1428 <td nosave="">03/07/00</td> 1429 1430 <td>SR</td> 1431 1432 <td>2.0b</td> 1433 1434 <td>C</td> 1435 1436 <td>Für eine Reihe von Unterprogrammen, die innerhalb 644 von pt</td> <td nosave="">calc_liquid_water_content</td> 645 </tr> <tr nosave="" valign="top"> <td nosave="">03/07/00</td> <td>SR</td> <td>2.0b</td> 646 <td>C</td> <td>Für eine Reihe von 647 Unterprogrammen, die innerhalb 1437 648 von <tt>prognostic_equation</tt> 1438 aufgerufen werden und an die Pointer-Variablen übergeben werden, 649 aufgerufen werden und an die Pointer-Variablen übergeben 650 werden, 1439 651 sind 1440 652 jetzt explizite Interfaces vereinbart (neues Modul <tt>pointer_interfaces</tt>). 1441 653 Auch im Vereinbarungsteil der entsprechenden Unterprogramme sind diese 1442 Übergabeparameter jetzt aus Konsistenzgründen (und weil es 654 Übergabeparameter jetzt aus Konsistenzgründen (und 655 weil es 1443 656 sonst 1444 657 Laufzeitfehler gibt) als Pointer-Variablen deklariert. Mit dieser 1445 658 Änderung 1446 ist der große Speicher- und CPU-Zeit-Bedarf beim Übersetzen 1447 von prognostic_equations wieder auf ein erträgliches Maß 659 ist der große Speicher- und CPU-Zeit-Bedarf beim 660 Übersetzen 661 von prognostic_equations wieder auf ein erträgliches 662 Maß 1448 663 reduziert. <br> 1449 1450 664 Im Rahmen dieser Änderungen wurde auch an diffusion_e ein 1451 665 Hilfsfeld 1452 666 weniger übergeben. <br> 1453 1454 Aus nicht vollständig gekärten Gründen haben sich die 1455 Zahlen in der Testergebnis-Datei minimal geändert (und zwar die 667 Aus nicht vollständig gekärten Gründen haben 668 sich die 669 Zahlen in der Testergebnis-Datei minimal geändert (und zwar 670 die 1456 671 Spalte 1457 672 mit der Gesamtdivergenz nach Aufruf des Drucklösers).</td> 1458 1459 <td>advec_s_bc, buoyancy, diffusion_e, diffusion_s, diffusion_u, 673 <td>advec_s_bc, buoyancy, diffusion_e, diffusion_s, diffusion_u, 1460 674 diffusion_v, 1461 675 diffusion_w, disturb_field, modules, production_e, prognostic_equations</td> 1462 1463 </tr> 1464 1465 <tr nosave="" valign="top"> 1466 1467 <td nosave="">04/07/00</td> 1468 1469 <td>SR</td> 1470 1471 <td>2.0b</td> 1472 1473 <td>B</td> 1474 1475 <td>Diriclet-Randbedingungen für Temperatur und 676 </tr> <tr nosave="" valign="top"> <td nosave="">04/07/00</td> <td>SR</td> <td>2.0b</td> 677 <td>B</td> <td>Diriclet-Randbedingungen für 678 Temperatur und 1476 679 Flüssigwassergehalt 1477 werden gesetzt. Dies wäre bereits nach Einführung der Pointer 680 werden gesetzt. Dies wäre bereits nach Einführung der 681 Pointer 1478 682 zwingend notwendig gewesen. Rechnungen mit vorgegebener 1479 683 Oberflächentemperatur 1480 schlugen deshalb bislang fehl.</td> 1481 1482 <td>boundary_conds</td> 1483 1484 </tr> 1485 1486 <tr nosave="" valign="top"> 1487 1488 <td nosave="">07/09/00</td> 1489 1490 <td>MS</td> 1491 1492 <td>2.0b</td> 1493 1494 <td>B</td> 1495 1496 <td>Die virtuelle potenielle Temperatur wird nun auch für 684 schlugen deshalb bislang fehl.</td> <td>boundary_conds</td> 685 </tr> <tr nosave="" valign="top"> <td nosave="">07/09/00</td> <td>MS</td> <td>2.0b</td> 686 <td>B</td> <td>Die virtuelle potenielle Temperatur 687 wird nun auch für 1497 688 k=nzb und 1498 689 k=nzt+1 berechnet, um die Randwerte auch für Ausgabezwecke 1499 bereichtzustellen. </td> 1500 1501 <td>compute_vpt</td> 1502 1503 </tr> 1504 1505 <tr nosave="" valign="top"> 1506 1507 <td nosave="">28/12/00</td> 1508 1509 <td>SR</td> 1510 1511 <td>2.1</td> 1512 1513 <td>C/N</td> 1514 1515 <td>PALM erlaubt jetzt den Einsatz optionaler Software-Pakete, 690 bereichtzustellen. </td> <td>compute_vpt</td> 691 </tr> <tr nosave="" valign="top"> <td nosave="">28/12/00</td> <td>SR</td> <td>2.1</td> 692 <td>C/N</td> <td>PALM erlaubt jetzt den Einsatz 693 optionaler Software-Pakete, 1516 694 die zwar 1517 zum Modell gehören, standardmäßig aber nicht mit 695 zum Modell gehören, standardmäßig aber 696 nicht mit 1518 697 übersetzt 1519 698 werden, um so unter anderem Kompilationszeit einzusparen. Siehe neues 1520 699 Kapitel 1521 3.7 in der Modelldokumentation. 1522 1523 <p>Einsatz der dvrp-Software wurde vollständig 700 3.7 in der Modelldokumentation. <p>Einsatz der 701 dvrp-Software wurde vollständig 1524 702 überarbeitet. 1525 703 Sie ist nun als Software-Paket optional im Modell einsetzbar. Siehe 1526 704 überarbeitetes 1527 Kapitel 4.5.6 in der Modelldokumentation. </p> 1528 1529 1530 <p>Die Prognose der Partikelverlagerung ist ebenfalls in ein 705 Kapitel 4.5.6 in der Modelldokumentation. </p> <p>Die 706 Prognose der Partikelverlagerung ist ebenfalls in ein 1531 707 optionales 1532 Software-Paket ausgelagert.</p> 1533 1534 </td> 1535 1536 <td>advec_particles, header, init_3d_model, init_dvrp, 708 Software-Paket ausgelagert.</p> </td> <td>advec_particles, 709 header, init_3d_model, init_dvrp, 1537 710 init_particles, 1538 711 init_pegrid, modules, parin, parles, plot_dvrp, prognostic_equations, 1539 712 read_var_list, 1540 write_var_list 1541 1542 <p><b>Neu:</b> <br> 1543 1544 package_parin</p> 1545 1546 </td> 1547 1548 </tr> 1549 1550 <tr nosave="" valign="top"> 1551 1552 <td nosave="">28/12/00</td> 1553 1554 <td>SR</td> 1555 1556 <td>2.1</td> 1557 1558 <td>B</td> 1559 1560 <td>MPI_FINALIZE ans Ende des Hauptprogramms geschoben (stand 713 write_var_list <p><b>Neu:</b> <br> 714 package_parin</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">28/12/00</td> 715 <td>SR</td> <td>2.1</td> <td>B</td> 716 <td>MPI_FINALIZE ans Ende des Hauptprogramms geschoben (stand 1561 717 vorher vor 1562 718 cpu_auswertung, wo noch ein Barrier-Aufruf steht). comm2d war zu Beginn 1563 719 undefiniert und wird nun erst einmal gleich MPI_COMM_WORLD gesetzt.</td> 1564 1565 <td>parles</td> 1566 1567 </tr> 1568 1569 <tr nosave="" valign="top"> 1570 1571 <td nosave="">02/01/01</td> 1572 1573 <td>SR</td> 1574 1575 <td>2.1a</td> 1576 1577 <td>C</td> 1578 1579 <td>Offene Dateien werden sobald möglich geschlossen. 720 <td>parles</td> </tr> <tr nosave="" valign="top"> <td nosave="">02/01/01</td> 721 <td>SR</td> <td>2.1a</td> <td>C</td> 722 <td>Offene Dateien werden sobald möglich geschlossen. 1580 723 Wiedereröffnung 1581 724 bestimmter Dateien mit POSITION='APPEND' möglich. In diesem … … 1583 726 wurde das Unterprogramm close_files in close_file umbenannt und hat nun 1584 727 1 Argument. <br> 1585 1586 728 Anpassung der Diffusionsparametrisierung im 1D-Modell an das 3D-Modell.</td> 1587 1588 <td>advec_particles, check_cpu_time, check_open, cpu_auswertung, 729 <td>advec_particles, check_cpu_time, check_open, cpu_auswertung, 1589 730 init_1d_model, 1590 731 init_3d_model, init_particles, modules.f90, parin, parles, plot_2d, 1591 write_3d_binary 1592 1593 <p><b>Neu:</b> <br> 1594 1595 close_file </p> 1596 1597 1598 <p><b>Eliminiert:</b> <br> 1599 1600 close_files</p> 1601 1602 </td> 1603 1604 </tr> 1605 1606 <tr nosave="" valign="top"> 1607 1608 <td nosave="">02/01/01</td> 1609 1610 <td>SR</td> 1611 1612 <td>2.1a</td> 1613 1614 <td>B</td> 1615 1616 <td>Beseitigung kleiner Fehler.</td> 1617 1618 <td>check_cpu_time, diffusion_e, diffusivities, long_filter, 1619 production_e</td> 1620 1621 </tr> 1622 1623 <tr nosave="" valign="top"> 1624 1625 <td>05/01/01</td> 1626 1627 <td>SR</td> 1628 1629 <td>2.1b</td> 1630 1631 <td>N</td> 1632 1633 <td nosave="">Neues Software-Paket zur Berechnung von Spektren im 1634 Ortsraum</td> 1635 1636 <td>check_open, header, init_3d_model, leap_frog, modules, 732 write_3d_binary <p><b>Neu:</b> <br> 733 close_file </p> <p><b>Eliminiert:</b> 734 <br> 735 close_files</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">02/01/01</td> 736 <td>SR</td> <td>2.1a</td> <td>B</td> 737 <td>Beseitigung kleiner Fehler.</td> <td>check_cpu_time, 738 diffusion_e, diffusivities, long_filter, 739 production_e</td> </tr> <tr nosave="" valign="top"> <td>05/01/01</td> <td>SR</td> 740 <td>2.1b</td> <td>N</td> <td nosave="">Neues 741 Software-Paket zur Berechnung von Spektren im 742 Ortsraum</td> <td>check_open, header, init_3d_model, 743 leap_frog, modules, 1637 744 package_parin, 1638 parles, read_var_list, write_3d_binary, write_var_list 1639 1640 <p><b>Neu:</b> <br> 1641 1642 calc_spectra, plot_spectra</p> 1643 1644 </td> 1645 1646 </tr> 1647 1648 <tr valign="top"> 1649 1650 <td>25/01/01</td> 1651 1652 <td>SR</td> 1653 1654 <td>2.1c</td> 1655 1656 <td>C/N</td> 1657 1658 <td>Drucklöser kann jetzt die Singleton-FFT benutzen. Modul <tt>singleton</tt> 745 parles, read_var_list, write_3d_binary, write_var_list <p><b>Neu:</b> 746 <br> 747 calc_spectra, plot_spectra</p> </td> </tr> <tr valign="top"> <td>25/01/01</td> <td>SR</td> 748 <td>2.1c</td> <td>C/N</td> <td>Drucklöser 749 kann jetzt die Singleton-FFT benutzen. Modul <tt>singleton</tt> 1659 750 wurde deshalb von <tt>calc_spectra.f90</tt> nach <tt>modules.f90</tt> 1660 751 geschoben 1661 (neuer Initialisierungsparameter <tt>fft_method</tt>). <br>1662 1663 Flüsse können jetzt bei k=1unabhängig von der752 (neuer Initialisierungsparameter <tt>fft_method</tt>). 753 <br>Flüsse können jetzt bei k=1 754 unabhängig von der 1664 755 Verwendung 1665 756 einer Prandtl-Schicht vorgegeben werden (neuer 1666 Initialisierungsparameter <tt>use_surface_fluxes</tt>). <br>1667 1668 Modul <tt>test_variables</tt> wurde aus allenProgrammteilen entfernt.1669 Hauptprogramm wurde von parles nach <tt>palm</tt> umbenannt. <br>1670 757 Initialisierungsparameter <tt>use_surface_fluxes</tt>). 758 <br>Modul <tt>test_variables</tt> wurde aus allen 759 Programmteilen entfernt. 760 Hauptprogramm wurde von parles nach <tt>palm</tt> 761 umbenannt. <br> 1671 762 Modell kann auf DEC-Workstations eingesetzt werden (-D <tt>dec</tt>). 1672 <br> 1673 1674 Zusätzliche Zeitmessungen mit dvrp-Software. Schreiben von 763 <br>Zusätzliche Zeitmessungen mit dvrp-Software. 764 Schreiben von 1675 765 Partikel-Informationen 1676 766 ist optional (neuer Parameter <tt>write_particle_informations</tt>).</td> 1677 1678 <td>advec_particles, calc_spectra, check_parameters, 767 <td>advec_particles, calc_spectra, check_parameters, 1679 768 cpu_zeitmessung, diffusion_s, 1680 769 header, init_3d_model, init_dvrp, init_particles, leap_frog, modules, 1681 770 package_parin, 1682 771 parin, poisfft, read_var_list, swap_timelevel, write_var_list 1683 1684 <p><b>Neu:</b> <br> 1685 1686 palm </p> 1687 1688 1689 <p><b>Eliminiert:</b> <br> 1690 1691 parles, module_test</p> 1692 1693 </td> 1694 1695 </tr> 1696 1697 <tr valign="top"> 1698 1699 <td>25/01/01</td> 1700 1701 <td>SR</td> 1702 1703 <td>2.1c</td> 1704 1705 <td>B</td> 1706 1707 <td>Fehler beim öffnen von Unit 23 beseitigt. Modul <tt>interface</tt> 772 <p><b>Neu:</b> <br> 773 palm </p> <p><b>Eliminiert:</b> <br> 774 parles, module_test</p> </td> </tr> <tr valign="top"> <td>25/01/01</td> <td>SR</td> 775 <td>2.1c</td> <td>B</td> <td>Fehler 776 beim öffnen von Unit 23 beseitigt. Modul <tt>interface</tt> 1708 777 fehlte in disturb_field. <br> 1709 1710 778 Wertebereich von theta* in Prandtl-Fluxes eingeschränkt, weil 1711 sonst auf einigen Rechnern overflow droht, wenn u*=0.</td> 1712 1713 <td>check_open, disturb_field, prandtl_fluxes</td> 1714 1715 </tr> 1716 1717 <tr valign="top"> 1718 1719 <td>30/01/01</td> 1720 1721 <td>SR</td> 1722 1723 <td>2.1d</td> 1724 1725 <td>N/C</td> 1726 1727 <td>Prognose eines passiven Skalars ist möglich (anstatt 779 sonst auf einigen Rechnern overflow droht, wenn u*=0.</td> <td>check_open, 780 disturb_field, prandtl_fluxes</td> </tr> <tr valign="top"> <td>30/01/01</td> <td>SR</td> 781 <td>2.1d</td> <td>N/C</td> <td>Prognose 782 eines passiven Skalars ist möglich (anstatt 1728 783 Feuchte), 1729 784 dazu neuer Initialisierungsparameter <tt>passive_scalar</tt>. 1730 1731 <p>Falls kein Pfad für ftpcopy existiert, werden in <tt>check_open</tt> 1732 die Dateinamen für AVS-Koordinaten- und Datendatei auf "unknown" 1733 gesetzt.</p> 1734 1735 </td> 1736 1737 <td>asselin_filter, boundary_conds, check_open, check_parameters, 785 <p>Falls kein Pfad für ftpcopy existiert, werden in <tt>check_open</tt> 786 die Dateinamen für AVS-Koordinaten- und Datendatei auf 787 "unknown" 788 gesetzt.</p> </td> <td>asselin_filter, 789 boundary_conds, check_open, check_parameters, 1738 790 flow_statistics, 1739 791 header, init_1d_model, init_3d_model, modules, parin, plot_2d, … … 1741 793 prandtl_fluxes, prognostic_equations, read_var_list, swap_timelevel, 1742 794 write_3d_binary, 1743 write_var_list</td> 1744 1745 </tr> 1746 1747 <tr valign="top"> 1748 1749 <td>30/01/01</td> 1750 1751 <td>SR</td> 1752 1753 <td>2.1d</td> 1754 1755 <td>B</td> 1756 1757 <td>String-Ausgabe für use_surface_fluxes korrigiert. 795 write_var_list</td> </tr> <tr valign="top"> <td>30/01/01</td> 796 <td>SR</td> <td>2.1d</td> <td>B</td> 797 <td>String-Ausgabe für use_surface_fluxes korrigiert. 1758 798 Version vom 1759 799 25/01/01 brach deswegen bei Fortsetzungsläufen ab. <br> 1760 1761 xy-Plotausgabe von q und ql korrigiert.</td> 1762 1763 <td>plot_2d, write_var_list</td> 1764 1765 </tr> 1766 1767 <tr nosave="" valign="top"> 1768 1769 <td nosave="">30/03/01</td> 1770 1771 <td>SR</td> 1772 1773 <td>2.2</td> 1774 1775 <td>N/C</td> 1776 1777 <td>Sämtliche 3D-Felder der Zeitebene t+dt wurden entfernt. 800 xy-Plotausgabe von q und ql korrigiert.</td> <td>plot_2d, 801 write_var_list</td> </tr> <tr nosave="" valign="top"> <td nosave="">30/03/01</td> 802 <td>SR</td> <td>2.2</td> <td>N/C</td> 803 <td>Sämtliche 3D-Felder der Zeitebene t+dt wurden 804 entfernt. 1778 805 Ebenso 1779 806 die Arbeitsfelder work und work1. Hilfsfelder wurden soweit … … 1783 810 eine Reihe von Übergabeparametern gestrichen werden (z.B. beim 1784 811 Aufruf 1785 von <tt>pres, advec_s_bc, poisfft, production_e, sor </tt>sowie den 812 von <tt>pres, advec_s_bc, poisfft, production_e, sor </tt>sowie 813 den 1786 814 Spline-Unterprogrammen). 1787 815 Der Long-Filter wurde in das Upstream-Spline-Verfahren integriert. Da 1788 816 der 1789 817 Zeitfilter jetzt direkt in der prognostischen Gleichung angewendet wird 1790 (bisher geschah dies erst nach Aufruf des Drucklösers), haben sich 818 (bisher geschah dies erst nach Aufruf des Drucklösers), haben 819 sich 1791 820 die Zahlen in der Testergebnis-Datei geändert. Die mittleren 1792 821 Profile 1793 sind aber unbeeinflusst. 1794 1795 <p>Eine weitere kleinere Änderung der Testergebnis-Daten 822 sind aber unbeeinflusst. <p>Eine weitere kleinere 823 Änderung der Testergebnis-Daten 1796 824 wird dadurch 1797 825 verursacht, dass jetzt der Parameter <tt>adjust_mixing_length</tt> 1798 826 defaultmäßig <tt>.FALSE.</tt> 1799 ist und der Mischungsweg grundsätzlich zusätzlich auf <tt>0.7 827 ist und der Mischungsweg grundsätzlich zusätzlich 828 auf <tt>0.7 1800 829 * zu</tt> begrenzt wird. Falls <tt>adjust_mixing_length=T</tt>, 1801 wird jetzt der Mischungsweg in <tt>diffusivities</tt> analog zu <tt>diffusion_e</tt> 830 wird jetzt der Mischungsweg in <tt>diffusivities</tt> 831 analog zu <tt>diffusion_e</tt> 1802 832 an allen Gitterpunkten modifiziert. <b>Achtung: </b>Diese 1803 833 Änderungen 1804 können bei anderen Simulationen t.w. erhebliche Auswirkungen auf 834 können bei anderen Simulationen t.w. erhebliche Auswirkungen 835 auf 1805 836 die 1806 oberflächennahen Ergebnisse haben. </p> 1807 1808 1809 <p><b>Kleinere Änderungen:</b> In <tt>init_cloud_physics</tt> 837 oberflächennahen Ergebnisse haben. </p> <p><b>Kleinere 838 Änderungen:</b> In <tt>init_cloud_physics</tt> 1810 839 wird 1811 der Wert der Variablen <tt>surface_pressure</tt> in hPa belassen 840 der Wert der Variablen <tt>surface_pressure</tt> in hPa 841 belassen 1812 842 (bisher 1813 843 Umwandlung in Pa). <br> 1814 1815 844 Die Namelist-Namen der Software-Pakete sind t.w. 1816 845 geändert. <br> 1817 1818 Unit 14 (binäre Ausgabe fürFortsetzungsläufe) wird846 Unit 14 (binäre Ausgabe für 847 Fortsetzungsläufe) wird 1819 848 jetzt im Hauptprogramm geschlossen, damit auf diese Datei noch 1820 849 benutzergesteuerte 1821 850 Ausgaben erfolgen können. <br> 1822 1823 851 Laufparameter werden in keinem Fall mehr auf Unit 14 ausgegeben bzw. 1824 852 von Unit 13 gelesen, d.h. sie gelten jetzt tatsächlich nur 1825 853 für 1826 854 den jeweils aktuellen Lauf. <br> 1827 1828 855 Partikeladvektion funktioniert jetzt auch zusammen mit der 1829 Galilei-Transformation. </p> 1830 1831 1832 <p>Restliche deutsche Variablen-, Unterprogramm- und Modulnamen 856 Galilei-Transformation. </p> <p>Restliche 857 deutsche Variablen-, Unterprogramm- und Modulnamen 1833 858 wurden 1834 859 ins Englische übersetzt. Davon sind fast alle Programmteile 1835 860 betroffen 1836 (nicht in rechter Spalte aufgeführt).</p> 1837 1838 </td> 1839 1840 <td>advec_particles, advec_s_bc, advec_s_ups, advec_u_ups, 861 (nicht in rechter Spalte aufgeführt).</p> </td> <td>advec_particles, 862 advec_s_bc, advec_s_ups, advec_u_ups, 1841 863 advec_v_ups, 1842 864 advec_w_ups, boundary_conds, calc_spectra, check_parameters, … … 1847 869 prognostic_equations, 1848 870 read_var_list, sor, swap_timelevel, transpose_*, write_3d_binary, 1849 write_var_list 1850 1851 <p><b>Neu:</b> <br> 1852 1853 cpu_log, cpu_statistics </p> 1854 1855 1856 <p><b>Eliminiert:</b> <br> 1857 1858 asselin_filter, cpu_auswertung, cpu_zeitmessung, long_filter</p> 1859 1860 </td> 1861 1862 </tr> 1863 1864 <tr nosave="" valign="top"> 1865 1866 <td nosave="">12/07/01</td> 1867 1868 <td>SR</td> 1869 1870 <td>2.2a</td> 1871 1872 <td>N/C</td> 1873 1874 <td>Defaultwert der unteren Randbedingung für die TKE ist ab 871 write_var_list <p><b>Neu:</b> <br> 872 cpu_log, cpu_statistics </p> <p><b>Eliminiert:</b> 873 <br> 874 asselin_filter, cpu_auswertung, cpu_zeitmessung, long_filter</p> </td> 875 </tr> <tr nosave="" valign="top"> <td nosave="">12/07/01</td> <td>SR</td> <td>2.2a</td> 876 <td>N/C</td> <td>Defaultwert der unteren 877 Randbedingung für die TKE ist ab 1875 878 sofort <tt>bc_e_b 1876 879 = 'neumann'</tt>. <br> 1877 1878 880 Partikelquelle ist über Paketparameter steuerbar. 1879 881 Partikeleigenschaften … … 1881 883 ist 1882 884 jetzt ebenfalls über Paketparameter steuerbar. <br> 1883 1884 885 Allen Modulen (bis auf singleton) wurde eine SAVE-Anweisung 1885 886 hinzugefügt, 1886 damit die durch sie vereinbarten Variablen während der Rechnung in 1887 keinem Fall undefiniert werden können (diese Gefahr besteht z.B. 887 damit die durch sie vereinbarten Variablen während der 888 Rechnung in 889 keinem Fall undefiniert werden können (diese Gefahr besteht 890 z.B. 1888 891 auf 1889 SGI-Origin-Maschinen).</td> 1890 1891 <td>advec_particles, header, init_dvrp, init_particles, 892 SGI-Origin-Maschinen).</td> <td>advec_particles, header, 893 init_dvrp, init_particles, 1892 894 modules.f90, package_parin, 1893 plot_dvrp, user_interface</td> 1894 1895 </tr> 1896 1897 <tr> 1898 1899 <td>12/07/01</td> 1900 1901 <td>SR</td> 1902 1903 <td>2.2a</td> 1904 1905 <td>B</td> 1906 1907 <td>Verschiebung der OPEN-Anweisung für Unit 33 1908 (AVS-FLD-Datei).</td> 1909 1910 <td>check_open</td> 1911 1912 </tr> 1913 1914 <tr nosave="" valign="top"> 1915 1916 <td nosave="">20/07/01</td> 1917 1918 <td>SR</td> 1919 1920 <td>2.3</td> 1921 1922 <td>N</td> 1923 1924 <td>Einbau des Mehrgitterverfahrens zur Lösung der 1925 Poisson-Gleichung.</td> 1926 1927 <td>check_parameters, exchange_horiz, header, init_grid, 895 plot_dvrp, user_interface</td> </tr> <tr> <td>12/07/01</td> 896 <td>SR</td> <td>2.2a</td> <td>B</td> 897 <td>Verschiebung der OPEN-Anweisung für Unit 33 898 (AVS-FLD-Datei).</td> <td>check_open</td> </tr> 899 <tr nosave="" valign="top"> <td nosave="">20/07/01</td> 900 <td>SR</td> <td>2.3</td> <td>N</td> 901 <td>Einbau des Mehrgitterverfahrens zur Lösung der 902 Poisson-Gleichung.</td> <td>check_parameters, 903 exchange_horiz, header, init_grid, 1928 904 init_pegrid, modules, 1929 parin, pres, run_control 1930 1931 <p><b>Neu:</b> <br> 1932 1933 poismg</p> 1934 1935 </td> 1936 1937 </tr> 1938 1939 <tr nosave="" valign="top"> 1940 1941 <td nosave="">21/08/01</td> 1942 1943 <td>SR</td> 1944 1945 <td>2.3a</td> 1946 1947 <td>N/C</td> 1948 1949 <td>Ergänzung der dvrp-Software: Partikel können mit 905 parin, pres, run_control <p><b>Neu:</b> <br> 906 poismg</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">21/08/01</td> 907 <td>SR</td> <td>2.3a</td> <td>N/C</td> 908 <td>Ergänzung der dvrp-Software: Partikel 909 können mit 1950 910 Schwänzen 1951 911 versehen werden. Die Ausgabe der dvrp-Plotdaten durch einen separaten 1952 912 PE 1953 ist möglich (mrun-Option "<tt>-p dvrp_graphics+1PE</tt>"). Der 913 ist möglich (mrun-Option "<tt>-p dvrp_graphics+1PE</tt>"). 914 Der 1954 915 Kommunikator <tt>MPI_COMM_WORLD</tt> 1955 wurde deshalb durch <tt>comm_palm</tt> ersetzt. 1956 1957 <p>Reduzierung des Mischungsweges imWandbereich kann wahlweise916 wurde deshalb durch <tt>comm_palm</tt> 917 ersetzt. <p>Reduzierung des Mischungsweges im 918 Wandbereich kann wahlweise 1958 919 abgeschaltet 1959 920 werden. Für km und kh wird jetzt eine Neumann-Randbedingung am 1960 921 unteren 1961 922 Rand verwendet. Die TKE-Energieproduktionsberechnung durch Scherung des 1962 Grundstroms am unteren Rand wurde verbessert. </p> 1963 1964 1965 <p>Generelle Verwendung der bodennahen vertikalen 923 Grundstroms am unteren Rand wurde verbessert. </p> <p>Generelle 924 Verwendung der bodennahen vertikalen 1966 925 Impulsflüsse in 1967 926 den Diffusionstermen von u und v, immer wenn diese als Randbedingung 1968 927 vorgegeben 1969 sind (gemäß <tt>use_surface_fluxes</tt>, - bisher wurde 928 sind (gemäß <tt>use_surface_fluxes</tt>, 929 - bisher wurde 1970 930 dies 1971 931 nur bei eingeschalteter Prandtl-Schicht gemacht). Schubspannungen … … 1973 933 jetzt als Randbedingungen vorgegeben werden (allerdings unter 1974 934 zuhilfenahme 1975 benutzereigener Software). </p> 1976 1977 1978 <p>Zusätzliche Schnittstellen für benutzereigene 935 benutzereigener Software). </p> <p>Zusätzliche 936 Schnittstellen für benutzereigene 1979 937 Software in 1980 allen Tendenztermen, in <tt>flow_statistics,</tt> und für 938 allen Tendenztermen, in <tt>flow_statistics,</tt> und 939 für 1981 940 Partikel- 1982 941 bzw. dvrp-Programmteile (Bestimmung von Partikeleigenschaften, 1983 942 Festlegung 1984 von Farbtabellen). </p> 1985 1986 1987 <p><b>Achtung: </b>Diese Änderungen können bei 943 von Farbtabellen). </p> <p><b>Achtung: </b>Diese 944 Änderungen können bei 1988 945 Simulationen mit 1989 mittlerem Wind durch die Änderungen im bereich des unteren Randes 1990 t.w. erhebliche Auswirkungen auf die oberflächennahen Ergebnisse 1991 haben.</p> 1992 1993 </td> 1994 1995 <td>advec_particles, diffusion_e, diffusion_u, diffusion_v, 946 mittlerem Wind durch die Änderungen im bereich des unteren 947 Randes 948 t.w. erhebliche Auswirkungen auf die oberflächennahen 949 Ergebnisse 950 haben.</p> </td> <td>advec_particles, diffusion_e, 951 diffusion_u, diffusion_v, 1996 952 diffusivities, 1997 953 flow_statistics, header, init_dvrp, init_particles, init_pegrid, … … 2000 956 prognostic_equations, 2001 957 read_var_list, swap_timelevel, user_interface, write_var_list</td> 2002 2003 </tr> 2004 2005 <tr nosave="" valign="top"> 2006 2007 <td nosave="">21/08/01</td> 2008 2009 <td>SR</td> 2010 2011 <td>2.3a</td> 2012 2013 <td>B</td> 2014 2015 <td>Falsche Positionierung von MPI_ALLREDUCE in check_cpu_time 2016 korrigiert.</td> 2017 2018 <td>check_cpu_time</td> 2019 2020 </tr> 2021 2022 <tr nosave="" valign="top"> 2023 2024 <td nosave="">04/09/01</td> 2025 2026 <td>SR</td> 2027 2028 <td>2.3b</td> 2029 2030 <td>N</td> 2031 2032 <td>Zusätzliche Profilausgaben für 958 </tr> <tr nosave="" valign="top"> <td nosave="">21/08/01</td> <td>SR</td> <td>2.3a</td> 959 <td>B</td> <td>Falsche Positionierung von 960 MPI_ALLREDUCE in check_cpu_time 961 korrigiert.</td> <td>check_cpu_time</td> </tr> 962 <tr nosave="" valign="top"> <td nosave="">04/09/01</td> 963 <td>SR</td> <td>2.3b</td> <td>N</td> 964 <td>Zusätzliche Profilausgaben für 2033 965 Energieproduktionsterme möglich. 2034 Stördruck wird bei der FFT-Methode nicht mehr durch Aufsummieren 966 Stördruck wird bei der FFT-Methode nicht mehr durch 967 Aufsummieren 2035 968 der 2036 Werte zu den einzelnen Zeitschritten ermittelt.</td> 2037 2038 <td>check_parameters, flow_statistics, modules, pres</td> 2039 2040 </tr> 2041 2042 <tr nosave="" valign="top"> 2043 2044 <td nosave="">04/09/01</td> 2045 2046 <td>SR</td> 2047 2048 <td>2.3b</td> 2049 2050 <td>B</td> 2051 2052 <td>Austausch der Geisterränder für die zeitgefilterten 969 Werte zu den einzelnen Zeitschritten ermittelt.</td> <td>check_parameters, 970 flow_statistics, modules, pres</td> </tr> <tr nosave="" valign="top"> <td nosave="">04/09/01</td> 971 <td>SR</td> <td>2.3b</td> <td>B</td> 972 <td>Austausch der Geisterränder für die 973 zeitgefilterten 2053 974 Felder 2054 975 direkt nach Durchführung der Filterung (war fehlerhaft seit … … 2056 977 der Zeitebene t+dt (Version 2.2) und führte zu sehr kleinen 2057 978 Störungen 2058 an den seitlichen Rändern der Teilgebiete). 2059 2060 <p><b>Achtung:</b> Der Inhalt der Testergebnis-Datei ändert 979 an den seitlichen Rändern der Teilgebiete). <p><b>Achtung:</b> 980 Der Inhalt der Testergebnis-Datei ändert 2061 981 sich durch 2062 diese Fehlerkorrektur.</p> 2063 2064 </td> 2065 2066 <td>prognostic_equations</td> 2067 2068 </tr> 2069 2070 <tr nosave="" valign="top"> 2071 2072 <td nosave="">09/11/01</td> 2073 2074 <td>SR</td> 2075 2076 <td>2.3c</td> 2077 2078 <td>N</td> 2079 2080 <td>Farbe entlang der Partikelschwänze kann sich 982 diese Fehlerkorrektur.</p> </td> <td>prognostic_equations</td> 983 </tr> <tr nosave="" valign="top"> <td nosave="">09/11/01</td> <td>SR</td> <td>2.3c</td> 984 <td>N</td> <td>Farbe entlang der 985 Partikelschwänze kann sich 2081 986 ändern. Einschalten 2082 987 der Partikeladvektion bei Fortsetzungsläufen möglich.</td> 2083 2084 <td>advec_particles, init_particles, modules, package_parin, 2085 plot_dvrp</td> 2086 2087 </tr> 2088 2089 <tr nosave="" valign="top"> 2090 2091 <td nosave="">16/04/02</td> 2092 2093 <td>SR</td> 2094 2095 <td>2.3d</td> 2096 2097 <td>N</td> 2098 2099 <td>Vorgabe von Randbedingungen für Partikeladvektion 988 <td>advec_particles, init_particles, modules, package_parin, 989 plot_dvrp</td> </tr> <tr nosave="" valign="top"> 990 <td nosave="">16/04/02</td> <td>SR</td> 991 <td>2.3d</td> <td>N</td> <td>Vorgabe 992 von Randbedingungen für Partikeladvektion 2100 993 möglich. 2101 Partikeldaten können für spätere Analysen auf Datei 994 Partikeldaten können für spätere Analysen 995 auf Datei 2102 996 geschrieben 2103 werden (Unit 85). <tt>PARTICLE</tt>-Datentyp enthält 997 werden (Unit 85). <tt>PARTICLE</tt>-Datentyp 998 enthält 2104 999 Informationen 2105 über Partikelgeschwindigkeitskomponenten und Startposition. 2106 2107 <p>Skalartransport sowie entsprechendeDatenausgaben können1000 über Partikelgeschwindigkeitskomponenten und 1001 Startposition. <p>Skalartransport sowie entsprechende 1002 Datenausgaben können 2108 1003 durch 2109 1004 eigene Parameter gesteuert werden (nicht mehr über die 2110 Feuchteparameter). </p> 2111 2112 2113 <p>Im Fall von pdims(1)=1 (eindimensionales virtuelles 1005 Feuchteparameter). </p> <p>Im Fall von 1006 pdims(1)=1 (eindimensionales virtuelles 2114 1007 Prozessornetz in 2115 1008 y-Richtung) werden die zyklischen Randbedingungen in x-Richtung durch 2116 1009 direktes 2117 Umspeichern anstatt durch <tt>sendrecv </tt>realisiert. Ebenso 1010 Umspeichern anstatt durch <tt>sendrecv </tt>realisiert. 1011 Ebenso 2118 1012 werden 2119 die Transponierungen xz, yz, zx und zy eingespart (es wird aber 1013 die Transponierungen xz, yz, zx und zy eingespart (es wird 1014 aber 2120 1015 weiter 2121 umsortiert). </p> 2122 2123 2124 <p>Im Fall von 3D-Plotausgaben Aufruf von ftpcopy-Script durch 1016 umsortiert). </p> <p>Im Fall von 1017 3D-Plotausgaben Aufruf von ftpcopy-Script durch 2125 1018 batch_scp-Script 2126 1019 ersetzt (dient der Ermittelung von Dateizyklusnummern als Information 2127 1020 für 2128 die AVS-fld-Datei). </p> 2129 2130 2131 <p>Bei fehlerhafter Eröffnung einer CPU_MEASURES-Datei 1021 die AVS-fld-Datei). </p> <p>Bei fehlerhafter 1022 Eröffnung einer CPU_MEASURES-Datei 2132 1023 werden erneute 2133 Eröffnungsversuche durchgeführt.</p> 2134 2135 </td> 2136 2137 <td>advec_particles, check_open, check_parameters, 1024 Eröffnungsversuche durchgeführt.</p> </td> 1025 <td>advec_particles, check_open, check_parameters, 2138 1026 exchange_horiz, flow_statistics, 2139 1027 header, init_particles, modules, package_parin, parin, plot_2d, … … 2141 1029 prognostic_equations, read_var_list, transpose_xz, transpose_yz, 2142 1030 transpose_zx, 2143 transpose_zy, write_var_list</td> 2144 2145 </tr> 2146 2147 <tr nosave="" valign="top"> 2148 2149 <td nosave="">02/05/02</td> 2150 2151 <td>SR</td> 2152 2153 <td>2.3e</td> 2154 2155 <td>B</td> 2156 2157 <td>Wiedereinführung der 3D-Felder für die Zeitebene 1031 transpose_zy, write_var_list</td> </tr> <tr nosave="" valign="top"> <td nosave="">02/05/02</td> 1032 <td>SR</td> <td>2.3e</td> <td>B</td> 1033 <td>Wiedereinführung der 3D-Felder für die 1034 Zeitebene 2158 1035 t+dt, die 2159 1036 notwendig ist, damit in den Diffusionstermen bei Leapfrog-Zeitschritten 2160 1037 mit der korrekten Zeitebene t-dt gearbeitet wird. Seit Version 2.2 2161 1038 wurde 2162 fälschlicherweise die Zeitebene t+dt verwendet. Asselin-Filter ist 1039 fälschlicherweise die Zeitebene t+dt verwendet. Asselin-Filter 1040 ist 2163 1041 jetzt wieder eigenständiges Unterprogramm. <br> 2164 2165 1042 Die Wiedereinführung dieser dritten Zeitebene ist auch 2166 1043 Voraussetzung 2167 1044 für die in der nächsten Version geplante skalare 2168 Optimierung. 2169 2170 <p>Fehler bei Berechnung von Ausgabezeitpunkten eliminiert (trat 1045 Optimierung. <p>Fehler bei Berechnung von 1046 Ausgabezeitpunkten eliminiert (trat 2171 1047 auf, wenn 2172 bei Fortsetzungsläufen von Ausgabeabständen 0.0 auf von Null 2173 verschiedene Werte gewechselt werden sollte). </p> 2174 2175 2176 <p><b>Achtung:</b> Der Inhalt der Testergebnis-Datei ändert 1048 bei Fortsetzungsläufen von Ausgabeabständen 0.0 auf 1049 von Null 1050 verschiedene Werte gewechselt werden sollte). </p> <p><b>Achtung:</b> 1051 Der Inhalt der Testergebnis-Datei ändert 2177 1052 sich durch 2178 diese Fehlerkorrektur.</p> 2179 2180 </td> 2181 2182 <td>init_3d_model, leap_frog, modules, prognostic_equations, 2183 swap_timelevel 2184 2185 <p><b>Neu:</b> <br> 2186 2187 asselin_filter</p> 2188 2189 </td> 2190 2191 </tr> 2192 2193 <tr nosave="" valign="top"> 2194 2195 <td>02/05/02</td> 2196 2197 <td>SR</td> 2198 2199 <td>2.3e</td> 2200 2201 <td>N/C</td> 2202 2203 <td>Kleinere Anpassungen an IBM-Regatta-Systeme in 2204 check_parameters, cpu_log. 2205 2206 <p>Modul singleton in eigenständige Datei 2207 überführt. </p> 2208 2209 2210 <p>Global_min_max arbeitet mit REAL*4, um Komplikationen auf 1053 diese Fehlerkorrektur.</p> </td> <td>init_3d_model, 1054 leap_frog, modules, prognostic_equations, 1055 swap_timelevel <p><b>Neu:</b> <br> 1056 asselin_filter</p> </td> </tr> <tr nosave="" valign="top"> <td>02/05/02</td> <td>SR</td> 1057 <td>2.3e</td> <td>N/C</td> <td>Kleinere 1058 Anpassungen an IBM-Regatta-Systeme in 1059 check_parameters, cpu_log. <p>Modul singleton in 1060 eigenständige Datei 1061 überführt. </p> <p>Global_min_max 1062 arbeitet mit REAL*4, um Komplikationen auf 2211 1063 32-bit-Rechnern 2212 1064 zu vermeiden (Datentyp MPI_2REAL müsste dort sonst auf 2213 1065 MPI_2DOUBLE_PRECISION 2214 1066 geändert werden). <br> 2215 2216 1067 Horizontale Geschwindigkeitskomponenten für Partikeladvektion 2217 1068 werden exakt zwischen den jeweils benachbarten horizontalen 2218 1069 Gitterflächen 2219 interpoliert.</p> 2220 2221 </td> 2222 2223 <td nosave="">advec_particles, check_parameters, cpu_log, 1070 interpoliert.</p> </td> <td nosave="">advec_particles, 1071 check_parameters, cpu_log, 2224 1072 global_min_max, 2225 modules 2226 2227 <p><b>Neu:</b> <br> 2228 2229 singleton</p> 2230 2231 </td> 2232 2233 </tr> 2234 2235 <tr nosave="" valign="top"> 2236 2237 <td nosave="">11/06/02</td> 2238 2239 <td>SR</td> 2240 2241 <td>2.4</td> 2242 2243 <td>N/C</td> 2244 2245 <td>Optimierung und Anpassungen für einzelne Knoten der 2246 IBM-Regatta-Systeme. 2247 2248 <p><b>Skalare (Cache) Optimierung:</b> Tendenzterme werden 1073 modules <p><b>Neu:</b> <br> 1074 singleton</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">11/06/02</td> 1075 <td>SR</td> <td>2.4</td> <td>N/C</td> 1076 <td>Optimierung und Anpassungen für einzelne Knoten der 1077 IBM-Regatta-Systeme. <p><b>Skalare (Cache) Optimierung:</b> 1078 Tendenzterme werden 2249 1079 innerhalb einer 2250 1080 großen (i,j)-Schleife berechnet, die nun jeweils die gesamte … … 2252 1082 Gleichung umfasst. Schleifenindices i und j werden an die 2253 1083 Tendenzunterprogramme 2254 als Argumente übergeben. Die Tendenzunterprogramme sind als Module 2255 geschrieben, die mittels Technik der überladenen Funktionen auch 1084 als Argumente übergeben. Die Tendenzunterprogramme sind als 1085 Module 1086 geschrieben, die mittels Technik der überladenen Funktionen 1087 auch 2256 1088 ohne 2257 diese Indices i und j aufgerufen werden können, und dann wie in 1089 diese Indices i und j aufgerufen werden können, und dann wie 1090 in 2258 1091 früheren 2259 Versionen funktionieren (d.h. in ihnen laufen die Schleifen über 1092 Versionen funktionieren (d.h. in ihnen laufen die Schleifen 1093 über 2260 1094 alle 2261 1095 3 Dimensionen). Solche Teile der prognostischen Gleichungen, die … … 2265 1099 mit Piascek-Williams-Advektion gearbeitet wird, werden alle 2266 1100 prognostischen 2267 Gleichungen in einer einzigen (i,j)-Schleife gerechnet. </p> 2268 2269 2270 <p><b>Kommunikationsoptimierung:</b> Zusätzliche 1101 Gleichungen in einer einzigen (i,j)-Schleife gerechnet. </p> <p><b>Kommunikationsoptimierung:</b> 1102 Zusätzliche 2271 1103 Implementierung einer 2272 1104 1D-Gebietszerlegung in x-Richtung. Mit dem neuen direkten … … 2279 1111 dieser 2280 1112 1D-Zerlegung durch einfache zyklische Randbedingungen ersetzt. </p> 2281 2282 2283 <p><b>Kleinere Anpassungen:</b> String-Vergleiche erfordern beim 1113 <p><b>Kleinere Anpassungen:</b> String-Vergleiche 1114 erfordern beim 2284 1115 IBM-Compiler 2285 in der Regel den Einsatz der TRIM-Funktion, um überschüssige 1116 in der Regel den Einsatz der TRIM-Funktion, um 1117 überschüssige 2286 1118 Blanks am Ende des Strings abzuschneiden (sonst liefert der Vergleich 2287 1119 .F.). <br> 2288 2289 1120 Lokal allokierte Felder t.w. in sogenannte automatische Felder 2290 1121 umgewandelt 2291 (günstiger, um Gefahr von eventuellen Speicherlecks zu vermeiden). 1122 (günstiger, um Gefahr von eventuellen Speicherlecks zu 1123 vermeiden). 2292 1124 Einige Feldoperationen wurden wegen schlechter Performance auf der IBM 2293 entfernt (s. pres). </p> 2294 2295 2296 <p><b>Achtung: </b>Die Übersetzung des Modells erfordert 1125 entfernt (s. pres). </p> <p><b>Achtung: </b>Die 1126 Übersetzung des Modells erfordert 2297 1127 aufgrund 2298 1128 der Vielzahl nunmehr eingesetzter Module die Verwendung des 2299 1129 make-Mechanismus, 2300 um die korrekte Reihenfolge bei der Übersetzung der Programmteile 2301 zu gewährleisten. Die Beschreibung dieser Abhängigkeiten 1130 um die korrekte Reihenfolge bei der Übersetzung der 1131 Programmteile 1132 zu gewährleisten. Die Beschreibung dieser 1133 Abhängigkeiten 2302 1134 erfolgt 2303 in der Datei<tt> Makefile</tt>, die zusammen mit den Programmdateien 1135 in der Datei<tt> Makefile</tt>, die zusammen mit den 1136 Programmdateien 2304 1137 abgelegt 2305 ist und ebenfalls unter RCS-Verwaltung steht. </p> 2306 2307 2308 <p>Im User-Interface ist das Unterprogramm<tt> user_actions </tt>nun 1138 ist und ebenfalls unter RCS-Verwaltung steht. </p> <p>Im 1139 User-Interface ist das Unterprogramm<tt> user_actions </tt>nun 2309 1140 ebenfalls 2310 als Modul geschrieben und bedient sich der Methode der überladenen 1141 als Modul geschrieben und bedient sich der Methode der 1142 überladenen 2311 1143 Funktionen. Existierende benutzereigene Software muss entsprechend 2312 1144 angepasst 2313 werden.</p> 2314 2315 </td> 2316 2317 <td>advec_s_pw, advec_s_up, advec_u_pw, advec_u_up, advec_v_pw, 1145 werden.</p> </td> <td>advec_s_pw, advec_s_up, 1146 advec_u_pw, advec_u_up, advec_v_pw, 2318 1147 advec_v_up, 2319 1148 advec_w_pw, advec_w_up, buoyancy, calc_precipitation, calc_radiation, … … 2327 1156 plot_3d, plot_spectra, poisfft, pres, production_e, 2328 1157 prognostic_equations, 2329 user_interface, 2330 2331 <p><b>Neu:</b> <br> 2332 2333 fft_for_1d_decomp, Makefile, poisfft_hybrid, </p> 2334 2335 </td> 2336 2337 </tr> 2338 2339 <tr nosave="" valign="top"> 2340 2341 <td nosave="">11/06/02</td> 2342 2343 <td>SR</td> 2344 2345 <td>2.4</td> 2346 2347 <td>B</td> 2348 2349 <td>Fehlende Variablenvereinbarung in <tt>advec_particles</tt> 2350 ergänzt.</td> 2351 2352 <td>advec_particles</td> 2353 2354 </tr> 2355 2356 <tr> 2357 2358 <td valign="top">12/09/02</td> 2359 2360 <td valign="top">SR</td> 2361 2362 <td valign="top">2.4a</td> 2363 2364 <td valign="top">N</td> 2365 2366 <td valign="top">Partikel können mit Trägheit versehen 1158 user_interface, <p><b>Neu:</b> <br> 1159 fft_for_1d_decomp, Makefile, poisfft_hybrid, </p> </td> 1160 </tr> <tr nosave="" valign="top"> <td nosave="">11/06/02</td> <td>SR</td> <td>2.4</td> 1161 <td>B</td> <td>Fehlende Variablenvereinbarung in <tt>advec_particles</tt> 1162 ergänzt.</td> <td>advec_particles</td> </tr> 1163 <tr> <td valign="top">12/09/02</td> <td valign="top">SR</td> <td valign="top">2.4a</td> 1164 <td valign="top">N</td> <td valign="top">Partikel 1165 können mit Trägheit versehen 2367 1166 werden und 2368 1167 einen Dichteunterschied zum umgebenden Fluid haben. Zur Steuerung dient 2369 1168 der neue Paketparameter <tt>density_ratio</tt>.</td> 2370 2371 <td valign="top">advec_particles, header, init_particles, 2372 modules, package_parin</td> 2373 2374 </tr> 2375 2376 <tr> 2377 2378 <td valign="top">12/09/02</td> 2379 2380 <td valign="top">SR</td> 2381 2382 <td valign="top">2.4a</td> 2383 2384 <td valign="top">B</td> 2385 2386 <td valign="top">Fehler in Berechnung der Anfangsprofile von pt 1169 <td valign="top">advec_particles, header, init_particles, 1170 modules, package_parin</td> </tr> <tr> <td valign="top">12/09/02</td> <td valign="top">SR</td> 1171 <td valign="top">2.4a</td> <td valign="top">B</td> 1172 <td valign="top">Fehler in Berechnung der Anfangsprofile 1173 von pt 2387 1174 und q entfernt 2388 1175 (obere Feldgrenze von 10 konnte in <tt>pt_vertical_gradient_level_ind</tt> 2389 überschritten werden). Fehler führte unter Umständen zu 2390 "segmentation fault" bei Verwendung sehr großer Gitterpunktzahlen 2391 in z-Richtung. 2392 2393 <p>Fehler in Berechnung der v-Komponente am unteren Rand 1176 überschritten werden). Fehler führte unter 1177 Umständen zu 1178 "segmentation fault" bei Verwendung sehr großer 1179 Gitterpunktzahlen 1180 in z-Richtung. <p>Fehler in Berechnung der v-Komponente am 1181 unteren Rand 2394 1182 entfernt, die 2395 für den Scherungsproduktionsterm der TKE verwendet wird (Variable 2396 v_0). Fälschlicherweise wurde zur Berechnung u(k=1) statt v(k=1) 2397 verwendet.</p> 2398 2399 </td> 2400 2401 <td valign="top">check_parameters, modules, production_e</td> 2402 2403 </tr> 2404 2405 <tr> 2406 2407 <td valign="top">19/12/02</td> 2408 2409 <td valign="top">SR</td> 2410 2411 <td valign="top">2.5</td> 2412 2413 <td valign="top">N</td> 2414 2415 <td valign="top">Restart times can be set by user with new 1183 für den Scherungsproduktionsterm der TKE verwendet wird 1184 (Variable 1185 v_0). Fälschlicherweise wurde zur Berechnung u(k=1) statt 1186 v(k=1) 1187 verwendet.</p> </td> <td valign="top">check_parameters, 1188 modules, production_e</td> </tr> <tr> <td valign="top">19/12/02</td> <td valign="top">SR</td> 1189 <td valign="top">2.5</td> <td valign="top">N</td> 1190 <td valign="top">Restart times can be set by user with new 2416 1191 runtime (d3par) 2417 parameters <tt>restart_time</tt> and <tt>dt_restart</tt>. Run 1192 parameters <tt>restart_time</tt> and <tt>dt_restart</tt>. 1193 Run 2418 1194 description 2419 1195 header is written on file CONTINUE_RUN. Output of cpu statistics … … 2423 1199 Output of warnings in the job protocol in case of negative measured 2424 1200 cpu-times. 2425 Remaining cpu-time is also evaluated on IBM-Regatta.</td> 2426 2427 <td valign="top">cpu_log, cpu_statistics, 1201 Remaining cpu-time is also evaluated on IBM-Regatta.</td> <td valign="top">cpu_log, cpu_statistics, 2428 1202 check_parameters,header, leap_frog, 2429 local_tremain, local_tremain_ini, modules, palm, 2430 2431 <p><b>new:</b> <br> 2432 2433 check_for_restart </p> 2434 2435 2436 <p><b>deleted:</b> <br> 2437 2438 check_cpu_time</p> 2439 2440 </td> 2441 2442 </tr> 2443 2444 <tr> 2445 2446 <td valign="top"> <br> 2447 2448 </td> 2449 2450 <td valign="top"> <br> 2451 2452 </td> 2453 2454 <td valign="top"> <br> 2455 2456 </td> 2457 2458 <td valign="top">C</td> 2459 2460 <td valign="top">Unit 14 (BINOUT) must be opened using a special 1203 local_tremain, local_tremain_ini, modules, palm, <p><b>new:</b> 1204 <br> 1205 check_for_restart </p> <p><b>deleted:</b> <br> 1206 check_cpu_time</p> </td> </tr> <tr> <td valign="top"> <br> </td> <td valign="top"> <br> </td> <td valign="top"> <br> </td> <td valign="top">C</td> <td valign="top">Unit 1207 14 (BINOUT) must be opened using a special 2461 1208 process 2462 1209 id string, because on IBM, the PE rank differs between communicators … … 2466 1213 written on file. As the consequence, these files could not be read by 2467 1214 restart-jobs. 2468 Reading of array <tt>hom</tt> moved from <tt>init_3d_model </tt>to <tt>read_var_list</tt> 2469 (binary version number incremented to 2.2). 2470 2471 <p><b>Further optimization for IBM-Regatta-systems:</b> <br> 2472 1215 Reading of array <tt>hom</tt> moved from <tt>init_3d_model 1216 </tt>to <tt>read_var_list</tt> 1217 (binary version number incremented to 2.2). <p><b>Further 1218 optimization for IBM-Regatta-systems:</b> <br> 2473 1219 Additional optimization of the hybrid-solver for multinode usage 2474 1220 (overlapping 2475 of communication and computation). </p> 2476 2477 2478 <p>Further cache optimization by using strides and joining loops (<tt>hybrid_solver, 2479 pres</tt> and <tt>timestep</tt>). Joining of MPI_ALLREDUCE calls in 2480 timestep. </p> 2481 2482 2483 <p>In case of 1d-decomposition along x only a part of the 1221 of communication and computation). </p> <p>Further cache 1222 optimization by using strides and joining loops (<tt>hybrid_solver, 1223 pres</tt> and <tt>timestep</tt>). Joining of 1224 MPI_ALLREDUCE calls in 1225 timestep. </p> <p>In case of 1d-decomposition along x 1226 only a part of the 2484 1227 integral divisor 2485 1228 conditions is checked, on IBM hosts a 1d-decomposition along x is the … … 2489 1232 1d-decomposition along x, the hybrid-solver does not force a 2490 1233 1d-decomposition 2491 any more. </p> 2492 2493 2494 <p>Array notation changed to do-loop constructs due to better 2495 performance. </p> 2496 2497 2498 <p>Cyclic boundary conditions along <tt>y</tt> used instead of 1234 any more. </p> <p>Array notation changed to do-loop 1235 constructs due to better 1236 performance. </p> <p>Cyclic boundary conditions along <tt>y</tt> 1237 used instead of 2499 1238 sendrecv 2500 in case of a 1d-decomposition along <tt>x</tt>. SENDRECV replaced by 1239 in case of a 1d-decomposition along <tt>x</tt>. SENDRECV 1240 replaced by 2501 1241 nonblocking 2502 1242 routines ISEND and IRECV in <tt>exchange_horiz.f90</tt>. </p> 2503 2504 2505 <p>Speed optimization by removing MINVAL/MAXVAL calls and by 1243 <p>Speed optimization by removing MINVAL/MAXVAL calls and by 2506 1244 handling <br> 2507 2508 1245 the "abs" case in a different way than the min/max cases (routine <tt>global_min_max</tt>).</p> 2509 2510 </td> 2511 2512 <td valign="top">asselin_filter, check_open, exchange_horiz, 1246 </td> <td valign="top">asselin_filter, check_open, 1247 exchange_horiz, 2513 1248 global_min_max, 2514 1249 init_3d_model, init_pegrid, parin, poisfft_hybrid, read_var_list, 2515 1250 timestep, 2516 write_3d_binary, write_var_list</td> 2517 2518 </tr> 2519 2520 <tr> 2521 2522 <td valign="top"><br> 2523 2524 </td> 2525 2526 <td valign="top"><br> 2527 2528 </td> 2529 2530 <td valign="top"><br> 2531 2532 </td> 2533 2534 <td valign="top">B</td> 2535 2536 <td valign="top">Correction of mixing length term (l(k)/ll(k)). 1251 write_3d_binary, write_var_list</td> </tr> <tr> <td valign="top"><br> </td> <td valign="top"><br> 1252 </td> <td valign="top"><br> </td> <td valign="top">B</td> <td valign="top">Correction 1253 of mixing length term (l(k)/ll(k)). 2537 1254 The condition 2538 1255 kh=3*km in the unstable case is now also exactly met in the wall 2539 1256 adjustment 2540 1257 region. Factor 0.7 in wall adjustment part replaced by variable <tt>wall_adjustment_factor</tt>, 2541 which is set to 1.8 in <tt>modules.f90</tt>. The factor 0.7 was the 1258 which is set to 1.8 in <tt>modules.f90</tt>. The factor 1259 0.7 was the 2542 1260 possible 2543 reason for 2-delta-x-waves, which were observed since version 2.1d. 2544 2545 <p>Calculation of deformation tensor re-designed (<tt>production_e</tt>). 2546 </p> 2547 2548 2549 <p>STOP statements replaced by call of new subroutine <tt>local_stop, 2550 </tt>where 1261 reason for 2-delta-x-waves, which were observed since version 2.1d. <p>Calculation 1262 of deformation tensor re-designed (<tt>production_e</tt>). </p> 1263 <p>STOP statements replaced by call of new subroutine <tt>local_stop, 1264 </tt>where 2551 1265 MPI_FINALIZE is called before STOP in case of a parallel environment. </p> 2552 2553 2554 <p>ISO2D parameter dp set to true for ibm hosts. </p> 2555 2556 2557 <p>tend=p added in <tt>pres</tt> after calling sor method.</p> 2558 2559 </td> 2560 2561 <td valign="top">advec_s_bc, buoyancy, check_open, 1266 <p>ISO2D parameter dp set to true for ibm hosts. </p> <p>tend=p 1267 added in <tt>pres</tt> after calling sor method.</p> 1268 </td> <td valign="top">advec_s_bc, buoyancy, 1269 check_open, 2562 1270 check_parameters, close_file, 2563 1271 coriolis, diffusion_e, diffusivities, fft_for_1d_decomp, … … 2565 1273 init_3d_model, init_dvrp, init_grid, init_particles, init_pegrid, 2566 1274 parin, 2567 plot_2d, poisfft, pres, read_var_list, user_interface, 2568 2569 <p><b>new:</b> <br> 2570 2571 local_stop, production_e</p> 2572 2573 </td> 2574 2575 </tr> 2576 2577 <tr> 2578 2579 <td valign="top">03/03/03</td> 2580 2581 <td valign="top">SR</td> 2582 2583 <td valign="top">2.5a</td> 2584 2585 <td valign="top">N</td> 2586 2587 <td valign="top">Particle velocities are also stored in array 1275 plot_2d, poisfft, pres, read_var_list, user_interface, <p><b>new:</b> 1276 <br> 1277 local_stop, production_e</p> </td> </tr> <tr> 1278 <td valign="top">03/03/03</td> <td valign="top">SR</td> 1279 <td valign="top">2.5a</td> <td valign="top">N</td> 1280 <td valign="top">Particle velocities are also stored in 1281 array 2588 1282 particles in 2589 1283 case of zero density ratio. Steering of variables by dvrp browser 2590 included.</td> 2591 2592 <td valign="top">advec_particles, init_dvrp, modules,</td> 2593 2594 </tr> 2595 2596 <tr> 2597 2598 <td valign="top"><br> 2599 2600 </td> 2601 2602 <td valign="top"><br> 2603 2604 </td> 2605 2606 <td valign="top"><br> 2607 2608 </td> 2609 2610 <td valign="top">C</td> 2611 2612 <td valign="top">AVS data format changed from float to xdr_float 1284 included.</td> <td valign="top">advec_particles, 1285 init_dvrp, modules,</td> </tr> <tr> <td valign="top"><br> </td> <td valign="top"><br> 1286 </td> <td valign="top"><br> </td> <td valign="top">C</td> <td valign="top">AVS 1287 data format changed from float to xdr_float 2613 1288 (needed 2614 1289 on linux machines <br> 2615 2616 1290 due to the little/big endian problem). Updates for new version of dvrp 2617 software (e.g. using module dvrp is now mandatory).</td> 2618 2619 <td valign="top">close_file, init_dvrp, leap_frog, plot_dvrp</td> 2620 2621 </tr> 2622 2623 <tr> 2624 2625 <td valign="top"><br> 2626 2627 </td> 2628 2629 <td valign="top"><br> 2630 2631 </td> 2632 2633 <td valign="top"><br> 2634 2635 </td> 2636 2637 <td valign="top">B</td> 2638 2639 <td valign="top">Error in particle inertia part removed (exp_arg 1291 software (e.g. using module dvrp is now mandatory).</td> <td valign="top">close_file, init_dvrp, leap_frog, plot_dvrp</td> 1292 </tr> <tr> <td valign="top"><br> </td> 1293 <td valign="top"><br> </td> <td valign="top"><br> </td> <td valign="top">B</td> 1294 <td valign="top">Error in particle inertia part removed 1295 (exp_arg 2640 1296 must not 2641 1297 contain the timestep). <br> 2642 2643 1298 Error in calculation of the vertical flux of resolved scale energy 2644 1299 (profile 57) removed. Displacement for integers in mpi_particle_type 2645 1300 reduced 2646 from 16 to 8 on ibm.</td> 2647 2648 <td valign="top">advec_particles, flow_statistics, init_particles</td> 2649 2650 </tr> 2651 2652 <tr> 2653 2654 <td valign="top">12/03/03</td> 2655 2656 <td valign="top">SR</td> 2657 2658 <td valign="top">2.6</td> 2659 2660 <td valign="top">N</td> 2661 2662 <td valign="top"><b>Version optimized for NEC-SX6 parallel-vector 2663 machines. </b> 2664 2665 <p>There are two main changes. A new vectorizable routine <tt>prognostic_equations_vec</tt> 1301 from 16 to 8 on ibm.</td> <td valign="top">advec_particles, 1302 flow_statistics, init_particles</td> </tr> <tr> <td valign="top">12/03/03</td> <td valign="top">SR</td> 1303 <td valign="top">2.6</td> <td valign="top">N</td> 1304 <td valign="top"><b>Version optimized for NEC-SX6 1305 parallel-vector 1306 machines. </b> <p>There are two main changes. A 1307 new vectorizable routine <tt>prognostic_equations_vec</tt> 2666 1308 is added, where cache optimizations are undone. The tendency 2667 1309 subroutines 2668 called by <tt>prognostic_equations_vec</tt> contain the full 3d-loops 1310 called by <tt>prognostic_equations_vec</tt> contain the 1311 full 3d-loops 2669 1312 (compared 2670 1313 with the cache-optimized versions where only loops over k are carried 2671 1314 out). 2672 Additionally, a new pressure solver <tt>poisfft_hybrid_vec</tt> is 1315 Additionally, a new pressure solver <tt>poisfft_hybrid_vec</tt> 1316 is 2673 1317 added, 2674 1318 which requires a 1d-domain-decomposition (like poisfft_hybrid). In this … … 2676 1320 optionally 2677 1321 available by choosing <tt>fft_method </tt>= <i>'temperton-algorithm'</i>. 2678 </p> 2679 2680 2681 <p>Additional changes in <tt>flow_statistics</tt> to allow 2682 better vectorization. </p> 2683 2684 2685 <p>Small changes in routines <tt>diffusivities</tt> and <tt>pres</tt> 1322 </p> <p>Additional changes in <tt>flow_statistics</tt> 1323 to allow 1324 better vectorization. </p> <p>Small changes in routines <tt>diffusivities</tt> 1325 and <tt>pres</tt> 2686 1326 which 2687 1327 caused run time errors on IBM and NEC due to compiler problems. </p> 2688 2689 2690 <p>Reading of environment variable <tt>tasks_per_node</tt> moved 1328 <p>Reading of environment variable <tt>tasks_per_node</tt> 1329 moved 2691 1330 from routine <tt>poisfft_hybrid</tt> 2692 to routine <tt>parin</tt>. </p> 2693 2694 2695 <p><b>Changes in makefile:</b> modules and user_interface now 1331 to routine <tt>parin</tt>. </p> <p><b>Changes 1332 in makefile:</b> modules and user_interface now 2696 1333 depend on 2697 1334 the f90 files, dependency of singleton added, LDFLAGS moved to the end 2698 of PROG rule, temperton_fft added.</p> 2699 2700 </td> 2701 2702 <td valign="top">advec_s_pw, advec_s_up, advec_u_pw, advec_u_up, 1335 of PROG rule, temperton_fft added.</p> </td> <td valign="top">advec_s_pw, advec_s_up, advec_u_pw, advec_u_up, 2703 1336 advec_v_pw, 2704 1337 advec_v_up, advec_w_pw, advec_w_up, buoyancy, calc_precipitation, … … 2710 1343 init_pegrid, leap_frog, local_tremain, modules, parin, poisfft, 2711 1344 poisfft_hybrid, 2712 pres, production_e, prognostic_equations, user_interface 2713 2714 <p><b>new:</b> <br> 2715 2716 temperton_fft</p> 2717 2718 </td> 2719 2720 </tr> 2721 2722 <tr> 2723 2724 <td valign="top"><br> 2725 2726 </td> 2727 2728 <td valign="top"><br> 2729 2730 </td> 2731 2732 <td valign="top"><br> 2733 2734 </td> 2735 2736 <td valign="top">B</td> 2737 2738 <td valign="top">tend=p added in routine <tt>pres</tt> after 1345 pres, production_e, prognostic_equations, user_interface <p><b>new:</b> 1346 <br> 1347 temperton_fft</p> </td> </tr> <tr> <td valign="top"><br> </td> <td valign="top"><br> 1348 </td> <td valign="top"><br> </td> <td valign="top">B</td> <td valign="top">tend=p 1349 added in routine <tt>pres</tt> after 2739 1350 calling sor 2740 method.</td> 2741 2742 <td valign="top">pres</td> 2743 2744 </tr> 2745 2746 <tr> 2747 2748 <td valign="top">14/03/03</td> 2749 2750 <td valign="top">SR</td> 2751 2752 <td valign="top">2.6a</td> 2753 2754 <td valign="top">N</td> 2755 2756 <td valign="top">Additional vector optimization for NEC-SX6. 2757 2758 <p>Optional system-specific random number generator available 1351 method.</td> <td valign="top">pres</td> </tr> 1352 <tr> <td valign="top">14/03/03</td> <td valign="top">SR</td> <td valign="top">2.6a</td> 1353 <td valign="top">N</td> <td valign="top">Additional 1354 vector optimization for NEC-SX6. <p>Optional system-specific 1355 random number generator available 2759 1356 (for equal 2760 distributed numbers). </p> 2761 2762 2763 <p>Define strings in all relevant subroutines changed from </p> 2764 2765 2766 <p><tt> #if defined( define_string )</tt> 2767 </p> 2768 2769 2770 <p>to </p> 2771 2772 2773 <p><tt> #if defined( __define_string )</tt> 2774 </p> 2775 2776 2777 <p><b>Caution: </b>The relevant subroutines are NOT listed in 1357 distributed numbers). </p> <p>Define strings in all 1358 relevant subroutines changed from </p> <p><tt> 1359 #if defined( define_string )</tt> </p> <p>to </p> 1360 <p><tt> #if 1361 defined( __define_string )</tt> </p> <p><b>Caution: 1362 </b>The relevant subroutines are NOT listed in 2778 1363 the right 2779 column!</p> 2780 2781 </td> 2782 2783 <td valign="top">asselin_filter, check_parameters, diffusion_e, 1364 column!</p> </td> <td valign="top">asselin_filter, 1365 check_parameters, diffusion_e, 2784 1366 diffusivities, 2785 1367 disturb_field, global_min_max, header, local_tremain, modules, parin, 2786 1368 poisfft_hybrid, 2787 pres, read_var_list, write_var_list</td> 2788 2789 </tr> 2790 2791 <tr> 2792 2793 <td valign="top"><br> 2794 2795 </td> 2796 2797 <td valign="top"><br> 2798 2799 </td> 2800 2801 <td valign="top"><br> 2802 2803 </td> 2804 2805 <td valign="top">B</td> 2806 2807 <td valign="top">Error in particle boundary condition removed 1369 pres, read_var_list, write_var_list</td> </tr> <tr> 1370 <td valign="top"><br> </td> <td valign="top"><br> </td> <td valign="top"><br> 1371 </td> <td valign="top">B</td> <td valign="top">Error in particle boundary condition removed 2808 1372 (velocity must 2809 be inverted in case of reflection)</td> 2810 2811 <td valign="top">advec_particles</td> 2812 2813 </tr> 2814 2815 <tr> 2816 2817 <td valign="top">16/04/03</td> 2818 2819 <td valign="top">SR</td> 2820 2821 <td valign="top">2.6b</td> 2822 2823 <td valign="top">N</td> 2824 2825 <td valign="top">Time series output of Monin Obukhov length. 2826 2827 <p>Temperton fft can now be used for all hosts and every domain 1373 be inverted in case of reflection)</td> <td valign="top">advec_particles</td> 1374 </tr> <tr> <td valign="top">16/04/03</td> 1375 <td valign="top">SR</td> <td valign="top">2.6b</td> 1376 <td valign="top">N</td> <td valign="top">Time 1377 series output of Monin Obukhov length. <p>Temperton fft can now 1378 be used for all hosts and every domain 2828 1379 decomposition. 2829 1380 Abort in case of Temperton fft, if number of gridpoints along x and/or 2830 y contain illegal factors. </p> 2831 2832 2833 <p>Index values for the extrema found in global_min_max are 1381 y contain illegal factors. </p> <p>Index values 1382 for the extrema found in global_min_max are 2834 1383 limited to 2835 1384 the range 0..nx, 0..ny (on IBM machines -1 and nx+1/ny+1 occured which 2836 1385 produced different RUN_CONTROL output compared to other machines). </p> 2837 2838 2839 <p>Output format of iteration count in routine run_control 2840 enlarged.</p> 2841 2842 </td> 2843 2844 <td valign="top">check_parameters, fft_for_1d_decomp, 1386 <p>Output format of iteration count in routine run_control 1387 enlarged.</p> </td> <td valign="top">check_parameters, 1388 fft_for_1d_decomp, 2845 1389 global_min_max, init_pegrid, 2846 modules, plot_ts, poisfft, run_control, temperton_fft</td> 2847 2848 </tr> 2849 2850 <tr> 2851 2852 <td valign="top"><br> 2853 2854 </td> 2855 2856 <td valign="top"><br> 2857 2858 </td> 2859 2860 <td valign="top"><br> 2861 2862 </td> 2863 2864 <td valign="top">B</td> 2865 2866 <td valign="top">Additional checks of variables <tt>hybrid_solver</tt> 1390 modules, plot_ts, poisfft, run_control, temperton_fft</td> </tr> 1391 <tr> <td valign="top"><br> </td> <td valign="top"><br> </td> <td valign="top"><br> 1392 </td> <td valign="top">B</td> <td valign="top">Additional checks of variables <tt>hybrid_solver</tt> 2867 1393 and <tt>host</tt> 2868 (<tt>hybrid_solver</tt> must not be used for a 2d-domain-decomposition,2869 1394 (<tt>hybrid_solver</tt> must not be used for a 1395 2d-domain-decomposition, <tt>host</tt> 2870 1396 should be set by setting environment variable <tt>localhost</tt>). 2871 2872 <p>Header output for mixing length limitations revised.</p> 2873 2874 </td> 2875 2876 <td valign="top">check_parameters, header</td> 2877 2878 </tr> 2879 2880 <tr> 2881 2882 <td valign="top">09/05/03</td> 2883 2884 <td valign="top">SR</td> 2885 2886 <td valign="top">2.7</td> 2887 2888 <td valign="top">N</td> 2889 2890 <td valign="top"><b>Version running on Linux Clusters using MPICH 1397 <p>Header output for mixing length limitations revised.</p> 1398 </td> <td valign="top">check_parameters, header</td> 1399 </tr> <tr> <td valign="top">09/05/03</td> 1400 <td valign="top">SR</td> <td valign="top">2.7</td> 1401 <td valign="top">N</td> <td valign="top"><b>Version 1402 running on Linux Clusters using MPICH 2891 1403 and Intel 2892 FORTRAN compiler (ifc)</b> 2893 2894 <p>So far, only absolutely neccessary changes have been done 1404 FORTRAN compiler (ifc)</b> <p>So far, only absolutely 1405 neccessary changes have been done 2895 1406 (possible 2896 optimizations will follow in a later version): </p> 2897 2898 2899 <p>New time measurements for Linux (ifc) environment added. </p> 2900 2901 2902 <p>Some parameters in MPI calls had to be modified (arrays had to 1407 optimizations will follow in a later version): </p> <p>New 1408 time measurements for Linux (ifc) environment added. </p> <p>Some 1409 parameters in MPI calls had to be modified (arrays had to 2903 1410 be replaced 2904 1411 by the first element of the regarding array) in order to fulfill f90 … … 2906 1413 and rank requirements, since on the MUK-cluster a FORTRAN90-version of 2907 1414 MPI is used. This also required to replace "mpif.h" by using a special 2908 module (named <tt>mpi</tt>). I did not find out how character 1415 module (named <tt>mpi</tt>). I did not find out 1416 how character 2909 1417 strings 2910 1418 can be send with MPI. Therefore, these strings are transformed to … … 2912 1420 before they are send, and transformed back to characters on the 2913 1421 receiving 2914 PE. </p> 2915 2916 2917 <p>On the MUK-cluster, only PE0 is able to read the values of 1422 PE. </p> <p>On the MUK-cluster, only PE0 is able to read 1423 the values of 2918 1424 environment 2919 1425 variables. Therefore, these values are communicated via broadcast to 2920 1426 the 2921 other PEs. </p> 2922 2923 2924 <p>1d-decomposition is set as the default on Linux Clusters. 1427 other PEs. </p> <p>1d-decomposition is set as the default 1428 on Linux Clusters. 2925 1429 Character 2926 strings have to be transformed to integer </p> 2927 2928 2929 <p>Smaller changes: </p> 2930 2931 2932 <p>Batch_scp paths for IBM, NEC and Linux added</p> 2933 2934 </td> 2935 2936 <td valign="top">check_open, cpu_log, flow_statistics, 1430 strings have to be transformed to integer </p> <p>Smaller 1431 changes: </p> <p>Batch_scp paths for IBM, NEC and Linux 1432 added</p> </td> <td valign="top">check_open, 1433 cpu_log, flow_statistics, 2937 1434 init_pegrid, local_getenv, 2938 modules, palm, poisfft_hybrid</td> 2939 2940 </tr> 2941 2942 <tr> 2943 2944 <td valign="top"><br> 2945 2946 </td> 2947 2948 <td valign="top"><br> 2949 2950 </td> 2951 2952 <td valign="top"><br> 2953 2954 </td> 2955 2956 <td valign="top">B</td> 2957 2958 <td valign="top">Measurements on IBM are now using function <tt>irtc</tt>, 1435 modules, palm, poisfft_hybrid</td> </tr> <tr> <td valign="top"><br> </td> <td valign="top"><br> 1436 </td> <td valign="top"><br> </td> <td valign="top">B</td> <td valign="top">Measurements 1437 on IBM are now using function <tt>irtc</tt>, 2959 1438 which allows correct measurements for jobs running over the 24:00 2960 timeline.</td> 2961 2962 <td valign="top">cpu_log, local_tremain, local_tremain_ini</td> 2963 2964 </tr> 2965 2966 <tr> 2967 2968 <td valign="top">01/08/03</td> 2969 2970 <td valign="top">SR</td> 2971 2972 <td valign="top">2.7a</td> 2973 2974 <td valign="top">B</td> 2975 2976 <td valign="top">Check that the number processors is also an 1439 timeline.</td> <td valign="top">cpu_log, 1440 local_tremain, local_tremain_ini</td> </tr> <tr> <td valign="top">01/08/03</td> <td valign="top">SR</td> 1441 <td valign="top">2.7a</td> <td valign="top">B</td> 1442 <td valign="top">Check that the number processors is also 1443 an 2977 1444 integral divisor 2978 of the number of gridpoints along y in case of a 1d-decomposition2979 along x. 2980 2981 <p>Error concerning the multinode-version ofpoisfft_hybrid1445 of the number of gridpoints along y in case of a 1446 1d-decomposition 1447 along x. <p>Error concerning the multinode-version of 1448 poisfft_hybrid 2982 1449 removed. It 2983 1450 was caused by the Linux changes of version 2.7 and led to program abort 2984 due to MPI errors in <tt>MPI_ALLTOALL</tt>. </p> 2985 2986 2987 <p>No abort on t3e-systems if system specific routines are used. </p> 2988 2989 2990 <p>Array<tt> sums </tt>is initialized in<tt> init_3d_model </tt>before 2991 the first call of<tt> pres</tt>. Call of<tt> init_cloud_physics </tt>moved 2992 before<tt> init_particles</tt>.</p> 2993 2994 </td> 2995 2996 <td valign="top">fft_1dd_init, init_3d_model, init_pegrid, 2997 poisfft_hybrid</td> 2998 2999 </tr> 3000 3001 <tr> 3002 3003 <td valign="top"><br> 3004 3005 </td> 3006 3007 <td valign="top"><br> 3008 3009 </td> 3010 3011 <td valign="top"><br> 3012 3013 </td> 3014 3015 <td valign="top">N</td> 3016 3017 <td valign="top">Interactive steering of dvrp-graphic features by 1451 due to MPI errors in <tt>MPI_ALLTOALL</tt>. </p> <p>No 1452 abort on t3e-systems if system specific routines are used. </p> <p>Array<tt> 1453 sums </tt>is initialized in<tt> init_3d_model </tt>before 1454 the first call of<tt> pres</tt>. Call of<tt> 1455 init_cloud_physics </tt>moved 1456 before<tt> init_particles</tt>.</p> </td> <td valign="top">fft_1dd_init, init_3d_model, init_pegrid, 1457 poisfft_hybrid</td> </tr> <tr> <td valign="top"><br> 1458 </td> <td valign="top"><br> </td> <td valign="top"><br> </td> <td valign="top">N</td> 1459 <td valign="top">Interactive steering of dvrp-graphic 1460 features by 3018 1461 dvrp-browser 3019 1462 plugin extended (position of slicers, threshold values of isosurfaces) 3020 1463 . New dvrp_graphics package parameter<tt> slicer_range_limits_dvrp</tt>.</td> 3021 3022 <td valign="top">init_dvrp, modules, package_parin, plot_dvrp</td> 3023 3024 </tr> 3025 3026 <tr> 3027 3028 <td valign="top">29/10/03</td> 3029 3030 <td valign="top">SR</td> 3031 3032 <td valign="top">2.7b</td> 3033 3034 <td valign="top">N</td> 3035 3036 <td valign="top">In the multigrid method, on a defined level, 1464 <td valign="top">init_dvrp, modules, package_parin, 1465 plot_dvrp</td> </tr> <tr> <td valign="top">29/10/03</td> 1466 <td valign="top">SR</td> <td valign="top">2.7b</td> 1467 <td valign="top">N</td> <td valign="top">In 1468 the multigrid method, on a defined level, 3037 1469 data are gathered 3038 1470 on PE0 and further calculations are carried out only on this PE. New<tt> 3039 1471 d3par </tt>parameter <tt>mg_switch_to_pe0_level</tt>. 3040 3041 <p>Particle groups implemented. New data type<tt> 1472 <p>Particle groups implemented. New data type<tt> 3042 1473 particle_groups_type </tt>and data type<tt> particle_type </tt>modified. 3043 1474 Version 3044 1475 numbers are 3045 1476 output on the particle files. New<tt> particles_par </tt>parameter<tt> 3046 diameter</tt>. Parameter name<tt> uniform_psize</tt> changed to<tt> 1477 diameter</tt>. Parameter name<tt> uniform_psize</tt> 1478 changed to<tt> 3047 1479 dvrp_psize</tt>. 3048 1480 Parameters<tt> density_ratio </tt>is now an array.</p> 3049 3050 </td> 3051 3052 <td valign="top">advec_particles, check_open, exchange_horiz, 1481 </td> <td valign="top">advec_particles, check_open, 1482 exchange_horiz, 3053 1483 header, init_particles, 3054 init_pegrid, modules, package_parin, parin, poismg</td> 3055 3056 </tr> 3057 3058 <tr> 3059 3060 <td valign="top"><br> 3061 3062 </td> 3063 3064 <td valign="top"><br> 3065 3066 </td> 3067 3068 <td valign="top"><br> 3069 3070 </td> 3071 3072 <td valign="top">C</td> 3073 3074 <td valign="top">Random number generator from numerical recipes 1484 init_pegrid, modules, package_parin, parin, poismg</td> </tr> 1485 <tr> <td valign="top"><br> </td> <td valign="top"><br> </td> <td valign="top"><br> 1486 </td> <td valign="top">C</td> <td valign="top">Random number generator from numerical recipes 3075 1487 is now rewritten 3076 1488 as a module. In restart runs the seed of this generator is set to the 3077 1489 last 3078 1490 values of the previous run in order to keep the sequence of the random 3079 numbers. New binary version 2.2. 3080 3081 <p>Buffer is flushed for file RUN_CONTROL immediately after every 1491 numbers. New binary version 2.2. <p>Buffer is flushed for file 1492 RUN_CONTROL immediately after every 3082 1493 output 3083 on IBM and Linux cluster.</p> 3084 3085 </td> 3086 3087 <td valign="top">advec_particles, disturb_field, init_3d_model, 1494 on IBM and Linux cluster.</p> </td> <td valign="top">advec_particles, 1495 disturb_field, init_3d_model, 3088 1496 init_particles, 3089 random_function, random_gauss, run_control, write_3d_binary</td> 3090 3091 </tr> 3092 3093 <tr> 3094 3095 <td valign="top"><br> 3096 3097 </td> 3098 3099 <td valign="top"><br> 3100 3101 </td> 3102 3103 <td valign="top"><br> 3104 3105 </td> 3106 3107 <td valign="top">E</td> 3108 3109 <td valign="top">In routine buoyance, horizontal mean temperature 1497 random_function, random_gauss, run_control, write_3d_binary</td> </tr> 1498 <tr> <td valign="top"><br> </td> <td valign="top"><br> </td> <td valign="top"><br> 1499 </td> <td valign="top">E</td> <td valign="top">In routine buoyance, horizontal mean 1500 temperature 3110 1501 is now 3111 1502 taken from array hom instead of array sums (otherwise inconsistence in 3112 case of using more than one statistical region). 3113 3114 <p>Information about the time of the last timestep change added 1503 case of using more than one statistical region). <p>Information 1504 about the time of the last timestep change added 3115 1505 to the 3116 restart file. </p> 3117 3118 3119 <p>System call of<tt> batch_scp </tt>on IBM corrected.</p> 3120 3121 </td> 3122 3123 <td valign="top">buoyancy, check_open, read_var_list, 3124 write_var_list</td> 3125 3126 </tr> 3127 3128 <tr> 3129 3130 <td style="vertical-align: top;">28/01/04</td> 3131 3132 <td style="vertical-align: top;">SR</td> 3133 3134 <td style="vertical-align: top;">2.8</td> 3135 3136 <td style="vertical-align: top;">N</td> 3137 3138 <td style="vertical-align: top;">Runge-Kutta schemes (2nd and 3rd 1506 restart file. </p> <p>System call of<tt> batch_scp </tt>on 1507 IBM corrected.</p> </td> <td valign="top">buoyancy, 1508 check_open, read_var_list, 1509 write_var_list</td> </tr> <tr> <td style="vertical-align: top;">28/01/04</td> <td style="vertical-align: top;">SR</td> <td style="vertical-align: top;">2.8</td> <td style="vertical-align: top;">N</td> <td style="vertical-align: top;">Runge-Kutta schemes (2nd and 1510 3rd 3139 1511 order) 3140 1512 for time integration implemented. The third order scheme is the new 3141 1513 default. 3142 1514 It allows the timestep to be 0.9*CFL, which is much larger than for the 3143 former default leapfrog scheme. The intermediate steps, which are 1515 former default leapfrog scheme. The intermediate steps, which 1516 are 3144 1517 part of these schemes, are realized by an additional loop within routine<span style="font-family: monospace;"> 3145 1518 time_integration </span>(currently<span style="font-family: monospace;"> 3146 leap_frog</span>, but this routine will be renamed soon). Steering 1519 leap_frog</span>, but this routine will be renamed soon). 1520 Steering 3147 1521 variables<span style="font-family: monospace;"> 3148 at </span>and<span style="font-family: monospace;"> bt </span>have 3149 been replaced by array<span style="font-family: monospace;"> sct</span>. 1522 at </span>and<span style="font-family: monospace;"> 1523 bt </span>have 1524 been replaced by array<span style="font-family: monospace;"> 1525 sct</span>. 3150 1526 , which is particularly used in the prognostic equations. Values of 3151 1527 this 3152 1528 steering array are calculated within the new routine<span style="font-family: monospace;"> 3153 timestep_scheme_steering</span>. 3154 3155 <p>When using Runge-Kutta schemes, the timestep increment 1529 timestep_scheme_steering</span>. <p>When using Runge-Kutta 1530 schemes, the timestep increment 3156 1531 is freely 3157 1532 allowed to adjust after each timestep (the older schemes are using some 3158 restrictions, see routine<span style="font-family: monospace;"> timestep</span>). 3159 Also, routine<span style="font-family: monospace;"> run_control </span>is 1533 restrictions, see routine<span style="font-family: monospace;"> 1534 timestep</span>). 1535 Also, routine<span style="font-family: monospace;"> 1536 run_control </span>is 3160 1537 not automatically called in case of timestep changes, when Runge-Kutta 3161 schemes are switched on. </p> 3162 3163 3164 <p>The old leapfrog scheme is still implemented and should 1538 schemes are switched on. </p> <p>The old leapfrog scheme 1539 is still implemented and should 3165 1540 produce the 3166 same results as in the previous version(s)! </p> 3167 3168 3169 <p>Depending on the user experience, further adjustments to the 1541 same results as in the previous version(s)! </p> <p>Depending 1542 on the user experience, further adjustments to the 3170 1543 Runge-Kutta 3171 1544 schemes may be necessary within the next minor versions. </p> 3172 3173 3174 <p><span style="font-weight: bold;">Attention:</span> <br> 3175 1545 <p><span style="font-weight: bold;">Attention:</span> 1546 <br> 3176 1547 The content of the testresult-file (run control output) has 3177 1548 significantly 3178 changed due to the new default timestep scheme!</p> 3179 3180 </td> 3181 3182 <td style="vertical-align: top;">check_parameters, header, 1549 changed due to the new default timestep scheme!</p> </td> <td style="vertical-align: top;">check_parameters, header, 3183 1550 init_3d_model, 3184 1551 init_rankine, leap_frog, modules, prognostic_equations, read_var_list, 3185 swap_timelevel, timestep, write_var_list 3186 3187 <p><span style="font-weight: bold;">new:</span> <br> 3188 3189 timestep_scheme_steering</p> 3190 3191 </td> 3192 3193 </tr> 3194 3195 <tr> 3196 3197 <td style="vertical-align: top;"><br> 3198 3199 </td> 3200 3201 <td style="vertical-align: top;"><br> 3202 3203 </td> 3204 3205 <td style="vertical-align: top;"><br> 3206 3207 </td> 3208 3209 <td style="vertical-align: top;">C</td> 3210 3211 <td style="vertical-align: top;">Output of particle infos in 1552 swap_timelevel, timestep, write_var_list <p><span style="font-weight: bold;">new:</span> <br> 1553 timestep_scheme_steering</p> </td> </tr> <tr> 1554 <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">C</td> <td style="vertical-align: top;">Output of particle infos in 3212 1555 subroutine<span style="font-family: monospace;"> 3213 allocate_prt_memory </span>on demand only. 3214 3215 <p>Type<span style="font-family: monospace;"> log </span>(used 1556 allocate_prt_memory </span>on demand only. <p>Type<span style="font-family: monospace;"> log </span>(used 3216 1557 for 3217 1558 cpu time measurements) changed to<span style="font-family: monospace;"> 3218 logpoint </span>due to name conflict with FORTRAN intrinsic<span style="font-family: monospace;"> 3219 log.</span></p> 3220 3221 </td> 3222 3223 <td style="vertical-align: top;">advec_particles, cpu_log, </td> 3224 3225 </tr> 3226 3227 <tr> 3228 3229 <td style="vertical-align: top;"><br> 3230 3231 </td> 3232 3233 <td style="vertical-align: top;"><br> 3234 3235 </td> 3236 3237 <td style="vertical-align: top;"><br> 3238 3239 </td> 3240 3241 <td style="vertical-align: top;">E</td> 3242 3243 <td style="vertical-align: top;">In case of opening unit 80, a 1559 logpoint </span>due to name conflict with FORTRAN 1560 intrinsic<span style="font-family: monospace;"> 1561 log.</span></p> </td> <td style="vertical-align: top;">advec_particles, 1562 cpu_log, </td> </tr> <tr> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">E</td> <td style="vertical-align: top;">In case of opening unit 80, a 3244 1563 barrier 3245 1564 is set only for the first call of<span style="font-family: monospace;"> 3246 1565 check_open </span>(from routine<span style="font-family: monospace;"> 3247 init_particles</span>), in order to avoid the possibility of hanging 1566 init_particles</span>), in order to avoid the possibility of 1567 hanging 3248 1568 jobs, 3249 1569 which may occur if unit 80 is opened within routine<span style="font-family: monospace;"> 3250 allocate_prt_memory</span>.</td> 3251 3252 <td style="vertical-align: top;">check_open</td> 3253 3254 </tr> 3255 3256 <tr> 3257 3258 <td style="vertical-align: top;">30/01/04</td> 3259 3260 <td style="vertical-align: top;">SR</td> 3261 3262 <td style="vertical-align: top;">2.8a</td> 3263 3264 <td style="vertical-align: top;">N</td> 3265 3266 <td style="vertical-align: top;">In order to prepare the code for 1570 allocate_prt_memory</span>.</td> <td style="vertical-align: top;">check_open</td> </tr> 1571 <tr> <td style="vertical-align: top;">30/01/04</td> 1572 <td style="vertical-align: top;">SR</td> <td style="vertical-align: top;">2.8a</td> <td style="vertical-align: top;">N</td> <td style="vertical-align: top;">In order to prepare the code 1573 for 3267 1574 the simulation 3268 1575 of flow around buildings (to be realized in one of the next major … … 3270 1577 the lower<span style="font-family: monospace;"> k </span>index<span style="font-family: monospace;"> 3271 1578 nzb </span>has been replaced by a two-dimensional array<span style="font-family: monospace;"> 3272 nzb_2d </span>in many of the three-dimensional loops. So far, all 1579 nzb_2d </span>in many of the three-dimensional loops. So 1580 far, all 3273 1581 elements of this array are set to<span style="font-family: monospace;"> 3274 nzb</span>. 3275 3276 <p>This change should not effect simulation results and cpu time 1582 nzb</span>. <p>This change should not effect simulation 1583 results and cpu time 3277 1584 in any 3278 way.</p> 3279 3280 </td> 3281 3282 <td style="vertical-align: top;">advec_s_pw, advec_s_up, 1585 way.</p> </td> <td style="vertical-align: top;">advec_s_pw, 1586 advec_s_up, 3283 1587 advec_u_pw, advec_u_up, 3284 1588 advec_v_pw, advec_v_up, advec_w_pw, advec_w_up, asselin_filter, … … 3289 1593 diffusivities, 3290 1594 impact_of_latent_heat, init_3d_model, modules, production_e, 3291 prognostic_equations, </td> 3292 3293 </tr> 3294 3295 <tr> 3296 3297 <td style="vertical-align: top;"><br> 3298 3299 </td> 3300 3301 <td style="vertical-align: top;"><br> 3302 3303 </td> 3304 3305 <td style="vertical-align: top;"><br> 3306 3307 </td> 3308 3309 <td style="vertical-align: top;">C</td> 3310 3311 <td style="vertical-align: top;">Euler step informations are 1595 prognostic_equations, </td> </tr> <tr> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">C</td> <td style="vertical-align: top;">Euler step informations are 3312 1596 removed in 3313 1597 the run control output in case of Runge-Kutta schemes. Also, the 3314 1598 timestep 3315 increment limitation is removed for the Runge-Kutta schemes. 3316 3317 <p><span style="font-weight: bold;">Attention:</span> <br> 3318 1599 increment limitation is removed for the Runge-Kutta schemes. <p><span style="font-weight: bold;">Attention:</span> <br> 3319 1600 The content of the testresult-file (run control output) has changed 3320 due to this modification.</p> 3321 3322 </td> 3323 3324 <td style="vertical-align: top;">run_control, timestep</td> 3325 3326 </tr> 3327 3328 <tr> 3329 3330 <td style="vertical-align: top;"><br> 3331 3332 </td> 3333 3334 <td style="vertical-align: top;"><br> 3335 3336 </td> 3337 3338 <td style="vertical-align: top;"><br> 3339 3340 </td> 3341 3342 <td style="vertical-align: top;">E</td> 3343 3344 <td style="vertical-align: top;">Velocity gradients at the 1601 due to this modification.</p> </td> <td style="vertical-align: top;">run_control, timestep</td> 1602 </tr> <tr> <td style="vertical-align: top;"><br> 1603 </td> <td style="vertical-align: top;"><br> </td> 1604 <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">E</td> <td style="vertical-align: top;">Velocity gradients at the 3345 1605 surface calculated 3346 1606 in routine production_e are now limited. In case of inhomogeneous 3347 1607 surface 3348 1608 heating small diffusivities (<span style="font-family: monospace;">km</span>) 3349 sometimes caused very small timesteps due to this problem.</td> 3350 3351 <td style="vertical-align: top;">production_e</td> 3352 3353 </tr> 3354 3355 <tr> 3356 3357 <td style="vertical-align: top;">30/04/04<br> 3358 3359 </td> 3360 3361 <td style="vertical-align: top;">SR<br> 3362 3363 </td> 3364 3365 <td style="vertical-align: top;">2.8b<br> 3366 3367 </td> 3368 3369 <td style="vertical-align: top;">N<br> 3370 3371 </td> 3372 3373 <td style="vertical-align: top;">The number of processors along 1609 sometimes caused very small timesteps due to this problem.</td> <td style="vertical-align: top;">production_e</td> </tr> 1610 <tr> <td style="vertical-align: top;">30/04/04<br> 1611 </td> <td style="vertical-align: top;">SR<br> 1612 </td> <td style="vertical-align: top;">2.8b<br> 1613 </td> <td style="vertical-align: top;">N<br> </td> 1614 <td style="vertical-align: top;">The number of processors 1615 along 3374 1616 the respective directions of the virtual processor grid must not be 3375 1617 integral divisors of the number of gridpoints along x, y or z any more. … … 3379 1621 Although this causes a load imbalance, only runs with very small 3380 1622 numbers of processors will feel a significant decrease of performance. 3381 In order to be able to carry out the transpositions, array<span style="font-family: monospace;"> d </span>(and also working array<span style="font-family: monospace;"> tend</span>) have the same size on 1623 In order to be able to carry out the transpositions, array<span style="font-family: monospace;"> d </span>(and also 1624 working array<span style="font-family: monospace;"> tend</span>) 1625 have the same size on 3382 1626 all processors. This new feature is switched on by setting the 3383 inipar-parameter<span style="font-family: monospace;"> grid_matching </span>= 3384 <span style="font-style: italic;">'match'</span>, which is the 3385 default now.<br> 3386 3387 <br> 3388 3389 Routine<span style="font-family: monospace;"> poisfft </span>is now 1627 inipar-parameter<span style="font-family: monospace;"> 1628 grid_matching </span>= <span style="font-style: italic;">'match'</span>, 1629 which is the 1630 default now.<br> <br> 1631 Routine<span style="font-family: monospace;"> poisfft </span>is 1632 now 3390 1633 the 3391 1634 default for solving the Poisson-equation for the case of 2d- as well as 3392 1635 1d-domain-decompositions and has been changed from a subroutine to a 3393 1636 module. This routine now includes most of the optimization from the 3394 (Ketelsen) routine<span style="font-family: monospace;"> poisfft_hybrid 3395 </span>and additionally allows 1d-decompositions along y.<br> 3396 3397 <span style="font-weight: bold;">ATTENTION:</span> Routine<span style="font-family: monospace;"> poisfft_hybrid </span>might still be 1637 (Ketelsen) routine<span style="font-family: monospace;"> 1638 poisfft_hybrid </span>and additionally allows 1d-decompositions 1639 along y.<br> <span style="font-weight: bold;">ATTENTION:</span> 1640 Routine<span style="font-family: monospace;"> 1641 poisfft_hybrid </span>might still be 3398 1642 faster when running PALM on more than one node of the IBM-Regatta, 3399 although the difference is not very big. The <span style="font-family: monospace;">poisfft</span>-code for a 1643 although the difference is not very big. The <span style="font-family: monospace;">poisfft</span>-code 1644 for a 3400 1645 decomposition along y has still a poor performance and needs further 3401 optimization.<br> 3402 3403 <br> 3404 1646 optimization.<br> <br> 3405 1647 The spectra-package is now able to use all the available FFT-methods 3406 1648 (selected by inipar-parameter<span style="font-family: monospace;"> … … 3408 1650 Performance of the transpositions needed for calculating the spectra 3409 1651 has been improved by using new direct transpositions from z to y and 3410 from y to x instead of two subsequent transpositions.<br> 3411 3412 <br> 3413 3414 Runge-Kutta scheme implemented in the 1d-version (<span style="font-family: monospace;">init_1d_model</span>). Diffusion 1652 from y to x instead of two subsequent transpositions.<br> <br> 1653 Runge-Kutta scheme implemented in the 1d-version (<span style="font-family: monospace;">init_1d_model</span>). 1654 Diffusion 3415 1655 quantity arrays are allocated for previous timelevels only in case of 3416 leapfrog scheme.<br> 3417 3418 <br> 3419 1656 leapfrog scheme.<br> <br> 3420 1657 This version is adapted for use on COMPAQ-DECALPHA systems (main 3421 1658 changes: calling of SHAPE function within subroutine arguments (e.g.<span style="font-family: monospace;"> singleton</span>) 3422 1659 causes compiler errors, particle data type modified for better 3423 1660 alignment, DATA statements from modules removed because they caused 3424 internal compiler errors).<br> 3425 3426 </td> 3427 3428 <td style="vertical-align: top;">advec_s_ups, advec_u_ups, 1661 internal compiler errors).<br> </td> <td style="vertical-align: top;">advec_s_ups, advec_u_ups, 3429 1662 advec_v_ups, advec_w_ups, calc_spectra, check_parameters, fft_xy, 3430 1663 header, init_1d_model, init_3d_model, init_particles, init_pegrid, 3431 1664 modules, parin, plot_2d, poisfft, prandtl_fluxes, pres, read_var_list, 3432 1665 singleton, spline_x, spline_y, spline_z, time_integration, transpose, 3433 write_3d_binary, write_var_list<br> 3434 3435 </td> 3436 3437 </tr> 3438 3439 <tr> 3440 3441 <td style="vertical-align: top;"><br> 3442 3443 </td> 3444 3445 <td style="vertical-align: top;"><br> 3446 3447 </td> 3448 3449 <td style="vertical-align: top;"><br> 3450 3451 </td> 3452 3453 <td style="vertical-align: top;">C<br> 3454 3455 </td> 3456 3457 <td style="vertical-align: top;">FFT routines renamed. Module<span style="font-family: monospace;"> fft_for_1d_decomp </span>renamed<span style="font-family: monospace;"> fft_xy</span>. Subroutine<span style="font-family: monospace;"> leap_frog </span>renamed<span style="font-family: monospace;"> time_integration</span>.<br> 3458 3459 <br> 3460 1666 write_3d_binary, write_var_list<br> </td> </tr> <tr> 1667 <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">C<br> </td> <td style="vertical-align: top;">FFT routines renamed. Module<span style="font-family: monospace;"> fft_for_1d_decomp </span>renamed<span style="font-family: monospace;"> fft_xy</span>. 1668 Subroutine<span style="font-family: monospace;"> leap_frog 1669 </span>renamed<span style="font-family: monospace;"> 1670 time_integration</span>.<br> <br> 3461 1671 For all transpositions, MPI_ALLTOALL is used instead of MPI_ALLTOALLV. 3462 The names of many transposition indices have been changed.<br> 3463 3464 <br> 3465 3466 Inipar-parameter<span style="font-family: monospace;"> impulse_advec </span>renamed<span style="font-family: monospace;"> momentum_advec</span>.<br> 3467 3468 <br> 3469 3470 Reading of restart data is moved from<span style="font-family: monospace;"> init_3d_model </span>to new routine <span style="font-family: monospace;">read_3d_binary</span>.<br> 3471 3472 <br> 3473 3474 All former transpose routines collected into one file.<br> 3475 3476 </td> 3477 3478 <td style="vertical-align: top;">check_parameters, 1672 The names of many transposition indices have been changed.<br> <br> 1673 Inipar-parameter<span style="font-family: monospace;"> 1674 impulse_advec </span>renamed<span style="font-family: monospace;"> momentum_advec</span>.<br> 1675 <br> 1676 Reading of restart data is moved from<span style="font-family: monospace;"> init_3d_model </span>to 1677 new routine <span style="font-family: monospace;">read_3d_binary</span>.<br> 1678 <br> 1679 All former transpose routines collected into one file.<br> </td> 1680 <td style="vertical-align: top;">check_parameters, 3479 1681 flow_statistics, header, init_1d_model, init_3d_model, init_advec, 3480 1682 init_pegrid, palm, parin, poisfft, poisfft_hybrid, production_e, 3481 prognostic_equations, read_var_list, sline_x, spline_y, write_var_list<span style="font-weight: bold;"><br> 3482 3483 <br> 3484 1683 prognostic_equations, read_var_list, sline_x, spline_y, write_var_list<span style="font-weight: bold;"><br> <br> 3485 1684 new:</span><br> 3486 3487 fft_xy, read_3d_binary, time_integration, transpose<br> 3488 3489 <br> 3490 3491 <span style="font-weight: bold;">deleted:</span><br> 3492 1685 fft_xy, read_3d_binary, time_integration, transpose<br> <br> 1686 <span style="font-weight: bold;">deleted:</span><br> 3493 1687 fft_for_1d_decomp, leap_frog, transpose_xy, transpose_xz, transpose_yx, 3494 transpose_yz, transpose_zx, transpose_zy<br> 3495 3496 </td> 3497 3498 </tr> 3499 3500 <tr> 3501 3502 <td style="vertical-align: top;"><br> 3503 3504 </td> 3505 3506 <td style="vertical-align: top;"><br> 3507 3508 </td> 3509 3510 <td style="vertical-align: top;"><br> 3511 3512 </td> 3513 3514 <td style="vertical-align: top;">E<br> 3515 3516 </td> 3517 3518 <td style="vertical-align: top;">Small errors in the 1688 transpose_yz, transpose_zx, transpose_zy<br> </td> </tr> 1689 <tr> <td style="vertical-align: top;"><br> </td> 1690 <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">E<br> </td> <td style="vertical-align: top;">Small errors in the 3519 1691 particle-package removed, which sometimes caused aborts of restart runs 3520 on IBM.<br> 3521 3522 <br> 3523 1692 on IBM.<br> <br> 3524 1693 Small error in calculating the time series profiles removed (array 3525 1694 bounds of some arrays needed for profil-output had to be increased).<br> 3526 3527 <br> 3528 1695 <br> 3529 1696 Function fft (singleton method) replaced by subroutine fftn due to 3530 problems with 64-bit mode on IBM<br> 3531 3532 </td> 3533 3534 <td style="vertical-align: top;">advec_particles, 3535 check_parameters, <br> 3536 3537 </td> 3538 3539 </tr> 3540 3541 <tr> 3542 3543 <td style="vertical-align: top;">26/03/05<br> 3544 3545 </td> 3546 3547 <td style="vertical-align: top;">SR<br> 3548 3549 </td> 3550 3551 <td style="vertical-align: top;">2.9<br> 3552 3553 </td> 3554 3555 <td style="vertical-align: top;">N<br> 3556 3557 </td> 3558 3559 <td style="vertical-align: top;">Non-cyclic lateral boundary 1697 problems with 64-bit mode on IBM<br> </td> <td style="vertical-align: top;">advec_particles, 1698 check_parameters, <br> </td> </tr> <tr> <td style="vertical-align: top;">26/03/05<br> </td> 1699 <td style="vertical-align: top;">SR<br> </td> 1700 <td style="vertical-align: top;">2.9<br> </td> 1701 <td style="vertical-align: top;">N<br> </td> <td style="vertical-align: top;">Non-cyclic lateral boundary 3560 1702 conditions implemented.<br> 3561 3562 1703 Dirichlet/Neumann conditions at the inflow/outflow can be used either 3563 1704 along x or along y. Cyclic conditions have to be used along the … … 3565 1706 upper bound nxl (nyn) of array u (v) is enlarged by one gridpoint on 3566 1707 the rightmost (northest) processor. These array index informations are 3567 passed to routine<span style="font-family: monospace;"> exchange_horiz </span>and<span style="font-family: monospace;"> disturb_field </span>by two 1708 passed to routine<span style="font-family: monospace;"> 1709 exchange_horiz </span>and<span style="font-family: monospace;"> 1710 disturb_field </span>by two 3568 1711 additional arguments. Boundary conditions are set with new 3569 inipar-parameters<span style="font-family: monospace;"> bc_lr </span>and<span style="font-family: monospace;"> bc_ns</span>. A damping layer is 3570 needed at the outflow in order to damp reflections (see new parameters<span style="font-family: monospace;"> km_damp_max </span>and<span style="font-family: monospace;"> outflow_damping_width</span>) and an 3571 additional disturbance is imposed near the inflow (see new parameters<span style="font-family: monospace;"> inflow_disturbance_begin </span>and<span style="font-family: monospace;"> inflow_disturbance_end</span>). 1712 inipar-parameters<span style="font-family: monospace;"> 1713 bc_lr </span>and<span style="font-family: monospace;"> 1714 bc_ns</span>. A damping layer is 1715 needed at the outflow in order to damp reflections (see new parameters<span style="font-family: monospace;"> km_damp_max </span>and<span style="font-family: monospace;"> outflow_damping_width</span>) 1716 and an 1717 additional disturbance is imposed near the inflow (see new parameters<span style="font-family: monospace;"> 1718 inflow_disturbance_begin </span>and<span style="font-family: monospace;"> inflow_disturbance_end</span>). 3572 1719 Non-cyclic conditions require the use of the multigrid-method for 3573 solving the Poisson equation. Argument<span style="font-family: monospace;"> range </span>added to routine<span style="font-family: monospace;"> boundary_conds </span>(where the 3574 non-cyclic conditions are set). <br> 3575 3576 <br> 3577 1720 solving the Poisson equation. Argument<span style="font-family: monospace;"> range </span>added 1721 to routine<span style="font-family: monospace;"> 1722 boundary_conds </span>(where the 1723 non-cyclic conditions are set). <br> <br> 3578 1724 Minor changes:<br> 3579 3580 Output of parameter<span style="font-family: monospace;"> datform </span>on 3581 &GLOBAL-Namelist (iso2d). This is set to <span style="font-style: italic;">'big_endian'</span> on IBM- and 3582 NEC-machines.<br> 3583 3584 <br> 3585 3586 Number of particle groups are additionally calculated in<span style="font-family: monospace;"> header </span>(formerly unknown).<br> 3587 3588 </td> 3589 3590 <td style="vertical-align: top;">advec_u_pw, 1725 Output of parameter<span style="font-family: monospace;"> 1726 datform </span>on 1727 &GLOBAL-Namelist (iso2d). This is set to <span style="font-style: italic;">'big_endian'</span> on 1728 IBM- and 1729 NEC-machines.<br> <br> 1730 Number of particle groups are additionally calculated in<span style="font-family: monospace;"> header </span>(formerly 1731 unknown).<br> </td> <td style="vertical-align: top;">advec_u_pw, 3591 1732 advec_u_up, advec_v_pw, advec_v_up, boundary_conds, 3592 1733 calc_liquid_water_content, check_parameters, close_file, coriolis, … … 3595 1736 init_rankine, modules, parin, poismg, pres, prognostic_equations, 3596 1737 read_var_list, sor, time_integration, timestep, write_var_list<br> 3597 3598 </td> 3599 3600 </tr> 3601 3602 <tr> 3603 3604 <td style="vertical-align: top;"><br> 3605 3606 </td> 3607 3608 <td style="vertical-align: top;"><br> 3609 3610 </td> 3611 3612 <td style="vertical-align: top;"><br> 3613 3614 </td> 3615 3616 <td style="vertical-align: top;">C<br> 3617 3618 </td> 3619 3620 <td style="vertical-align: top;"><span style="font-family: monospace;">pres </span>is called by default only 1738 </td> </tr> <tr> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">C<br> </td> <td style="vertical-align: top;"><span style="font-family: monospace;">pres </span>is called 1739 by default only 3621 1740 at the last Runge-Kutta-substep, which spares a big amount of CPU-time 3622 1741 (see also new parameter<span style="font-family: monospace;"> 3623 call_psolver_at_all_substeps</span>). No pressure term in the momentum 3624 equations in case of Runge-Kutta-schemes.<br> 3625 3626 <span style="font-weight: bold;">Attention:</span> <br> 3627 1742 call_psolver_at_all_substeps</span>). No pressure term in the 1743 momentum 1744 equations in case of Runge-Kutta-schemes.<br> <span style="font-weight: bold;">Attention:</span> <br> 3628 1745 The content of the testresult-file (run control output) has changed 3629 due to this modification.<br> 3630 3631 <br> 3632 3633 Program speedup by removing the exchange of ghost points in routines <span style="font-family: monospace;">calc_liquid_water_content</span>, <span style="font-family: monospace;">diffusivities</span>, .<br> 3634 3635 <br> 3636 3637 Default value for residual limit increased from 1E-6 to 1E-4. Routine<span style="font-family: monospace;"> prolong </span>in<span style="font-family: monospace;"> poismg </span>simplified (one call of<span style="font-family: monospace;"> exchange_horiz </span>spared).<br> 3638 3639 </td> 3640 3641 <td style="vertical-align: top;">calc_liquid_water_content, 1746 due to this modification.<br> <br> 1747 Program speedup by removing the exchange of ghost points in routines <span style="font-family: monospace;">calc_liquid_water_content</span>, 1748 <span style="font-family: monospace;">diffusivities</span>, 1749 .<br> <br> 1750 Default value for residual limit increased from 1E-6 to 1E-4. Routine<span style="font-family: monospace;"> prolong </span>in<span style="font-family: monospace;"> poismg </span>simplified 1751 (one call of<span style="font-family: monospace;"> 1752 exchange_horiz </span>spared).<br> </td> <td style="vertical-align: top;">calc_liquid_water_content, 3642 1753 diffusivities, modules, poismg, pres, time_integration, 3643 timestep_scheme_steering<br> 3644 3645 </td> 3646 3647 </tr> 3648 3649 <tr> 3650 3651 <td style="vertical-align: top;"><br> 3652 3653 </td> 3654 3655 <td style="vertical-align: top;"><br> 3656 3657 </td> 3658 3659 <td style="vertical-align: top;"><br> 3660 3661 </td> 3662 3663 <td style="vertical-align: top;">E<br> 3664 3665 </td> 3666 3667 <td style="vertical-align: top;">Calculation of vertical particle 3668 velocity (with inertia) corrected, exp_arg had a wrong sign.<br> 3669 3670 <br> 3671 1754 timestep_scheme_steering<br> </td> </tr> <tr> 1755 <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">E<br> </td> <td style="vertical-align: top;">Calculation of vertical 1756 particle 1757 velocity (with inertia) corrected, exp_arg had a wrong sign.<br> <br> 3672 1758 (e)**1.5 replaced by e*SQRT(e) in init_1d_model because of wrong 3673 1759 results on NEC machines (as a side effect, the new calculation is much 3674 faster).<br> 3675 3676 <br> 3677 1760 faster).<br> <br> 3678 1761 Two errors concerning switch to PE0 (multigrid method) removed from<span style="font-family: monospace;"> init_pegrid</span>.<br> 3679 3680 <br> 3681 3682 Two errors in<span style="font-family: monospace;"> plot_2d </span>removed: 1762 <br> 1763 Two errors in<span style="font-family: monospace;"> 1764 plot_2d </span>removed: 3683 1765 local_2d is allocated with upper bound<span style="font-family: monospace;"> nzt</span> (former<span style="font-family: monospace;"> nzt+1 </span>was 3684 1766 wrong), additional barrier in case that PE0 gathers data from the other 3685 1767 PEs, this barrier had to be set at the end of the gathering because 3686 otherwise PE0 sometimes received wrong data on tag 0.<br> 3687 3688 <br> 3689 3690 Use of module indices removed from routine<span style="font-family: monospace;"> split_1dd </span>in<span style="font-family: monospace;"> poisfft </span>since it caused 3691 errors in case of<span style="font-family: monospace;"> nx /= ny </span>and 3692 a 1D-decomposition along y.<br> 3693 3694 <br> 3695 1768 otherwise PE0 sometimes received wrong data on tag 0.<br> <br> 1769 Use of module indices removed from routine<span style="font-family: monospace;"> split_1dd </span>in<span style="font-family: monospace;"> poisfft </span>since 1770 it caused 1771 errors in case of<span style="font-family: monospace;"> nx 1772 /= ny </span>and 1773 a 1D-decomposition along y.<br> <br> 3696 1774 Wrong re-definition of dx2 and dy2 in<span style="font-family: monospace;"> sor </span>removed.<br> 3697 3698 </td> 3699 3700 <td style="vertical-align: top;">advec_particles, init_1d_model, 3701 init_pegrid, plot_2d, poisfft, sor<br> 3702 3703 </td> 3704 3705 </tr> 3706 3707 <tr> 3708 3709 <td style="vertical-align: top;">23/04/05<br> 3710 3711 </td> 3712 3713 <td style="vertical-align: top;">SR<br> 3714 3715 </td> 3716 3717 <td style="vertical-align: top;">2.9a<br> 3718 3719 </td> 3720 3721 <td style="vertical-align: top;">C<br> 3722 3723 </td> 3724 3725 <td style="vertical-align: top;">Additional check for 1775 </td> <td style="vertical-align: top;">advec_particles, 1776 init_1d_model, 1777 init_pegrid, plot_2d, poisfft, sor<br> </td> </tr> <tr> 1778 <td style="vertical-align: top;">23/04/05<br> </td> 1779 <td style="vertical-align: top;">SR<br> </td> 1780 <td style="vertical-align: top;">2.9a<br> </td> 1781 <td style="vertical-align: top;">C<br> </td> <td style="vertical-align: top;">Additional check for 3726 1782 incompatibilities between non-cyclic lateral boundary conditions and 3727 other schemes.<br> 3728 3729 <br> 3730 3731 fcl_factor renamed cfl_factor.<br> 3732 3733 <br> 3734 1783 other schemes.<br> <br> 1784 fcl_factor renamed cfl_factor.<br> <br> 3735 1785 New local array sums_ll declared in flow_statistics instead of 3736 temporarily using sums_l.<br> 3737 3738 <br> 3739 1786 temporarily using sums_l.<br> <br> 3740 1787 crmax (maximum number of crosses allowed in a plot of vertical 3741 profiles) increased from 20 to 100.<br> 3742 3743 </td> 3744 3745 <td style="vertical-align: top;">calc_spectra, check_parameters, 3746 flow_statistics, header, modules, parin, timestep<br> 3747 3748 </td> 3749 3750 </tr> 3751 3752 <tr> 3753 3754 <td style="vertical-align: top;"><br> 3755 3756 </td> 3757 3758 <td style="vertical-align: top;"><br> 3759 3760 </td> 3761 3762 <td style="vertical-align: top;"><br> 3763 3764 </td> 3765 3766 <td style="vertical-align: top;">E<br> 3767 3768 </td> 3769 3770 <td style="vertical-align: top;">Revised calculation of output 1788 profiles) increased from 20 to 100.<br> </td> <td style="vertical-align: top;">calc_spectra, check_parameters, 1789 flow_statistics, header, modules, parin, timestep<br> </td> 1790 </tr> <tr> <td style="vertical-align: top;"><br> 1791 </td> <td style="vertical-align: top;"><br> </td> 1792 <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">E<br> </td> <td style="vertical-align: top;">Revised calculation of output 3771 1793 time counters regarding a possible decrease of the output time interval 3772 in case of restart runs.<br> 3773 3774 <br> 3775 1794 in case of restart runs.<br> <br> 3776 1795 Error removed in Dirichlet bottom boundary conditions for pt and q in 3777 case of Runge-Kutta schemes.<br> 3778 3779 <br> 3780 3781 <span style="font-weight: bold;">batch_scp</span> calls (needed 1796 case of Runge-Kutta schemes.<br> <br> <span style="font-weight: bold;">batch_scp</span> calls 1797 (needed 3782 1798 for determining correct filenames for AVS-fld-files) is given the 3783 remote username as an additional argument.<br> 3784 3785 <br> 3786 3787 Default setting of outflow_damping_width corrected.<br> 3788 3789 <br> 3790 1799 remote username as an additional argument.<br> <br> 1800 Default setting of outflow_damping_width corrected.<br> <br> 3791 1801 Initial horizontal velocities at the lowest grid levels in the 1d-model 3792 1802 are set to a very small value in order to avoid wrong results and the 3793 resulting too small timesteps.<br> 3794 3795 <br> 3796 1803 resulting too small timesteps.<br> <br> 3797 1804 Implicit counters i renamed i9 in modules due to declaration conflicts 3798 with i in other subroutines (reported by Portland compiler).<br> 3799 3800 </td> 3801 3802 <td style="vertical-align: top;">advec_particles, boundary_conds, 1805 with i in other subroutines (reported by Portland compiler).<br> </td> 1806 <td style="vertical-align: top;">advec_particles, 1807 boundary_conds, 3803 1808 check_open, check_parameters, init_1d_model, modules, time_integration<br> 3804 3805 </td> 3806 3807 </tr> 3808 3809 <tr> 3810 3811 <td style="vertical-align: top;">18/05/05<br> 3812 3813 </td> 3814 3815 <td style="vertical-align: top;">SR<br> 3816 3817 </td> 3818 3819 <td style="vertical-align: top;">2.10<br> 3820 3821 </td> 3822 3823 <td style="vertical-align: top;">N<br> 3824 3825 </td> 3826 3827 <td style="vertical-align: top;">NetCDF support implemented.<br> 3828 1809 </td> </tr> <tr> <td style="vertical-align: top;">18/05/05<br> </td> 1810 <td style="vertical-align: top;">SR<br> </td> 1811 <td style="vertical-align: top;">2.10<br> </td> 1812 <td style="vertical-align: top;">N<br> </td> <td style="vertical-align: top;">NetCDF support implemented.<br> 3829 1813 Profiles, time series, spectra, cross-sections, 3d-volume-data and 3830 1814 particle data can now be output in NetCDF format (file handles 3831 1815 101-108). This output can be switched on with the new d3par-parameter <span style="font-family: monospace;">data_output_format</span>. 3832 1816 So far, this support is available on IBM systems at HLRN and on the NEC 3833 system at DKRZ. Output in format suitable for graphic software <span style="font-weight: bold;">profil</span>, <span style="font-weight: bold;">iso2d</span> and <span style="font-weight: bold;">avs</span> is still possible.<br> 3834 3835 <br> 3836 1817 system at DKRZ. Output in format suitable for graphic software <span style="font-weight: bold;">profil</span>, <span style="font-weight: bold;">iso2d</span> and <span style="font-weight: bold;">avs</span> is still 1818 possible.<br> <br> 3837 1819 Output of cross sections for qv, vpt and lwp (liquid water path).<br> 3838 3839 </td> 3840 3841 <td style="vertical-align: top;">advec_particles, check_open, 1820 </td> <td style="vertical-align: top;">advec_particles, 1821 check_open, 3842 1822 check_parameters, close_file, header, modules, parin, plot_1d, plot_2d, 3843 1823 plot_3d, plot_spectra, plot_ts, read_var_list, write_var_list<br> 3844 3845 <br> 3846 3847 <span style="font-weight: bold;">new:</span><br> 3848 3849 netcdf<br> 3850 3851 </td> 3852 3853 </tr> 3854 3855 <tr> 3856 3857 <td style="vertical-align: top;"><br> 3858 3859 </td> 3860 3861 <td style="vertical-align: top;"><br> 3862 3863 </td> 3864 3865 <td style="vertical-align: top;"><br> 3866 3867 </td> 3868 3869 <td style="vertical-align: top;">C<br> 3870 3871 </td> 3872 3873 <td style="vertical-align: top;">In <span style="font-family: monospace;">calc_spectra</span>, n is assigned the 1824 <br> <span style="font-weight: bold;">new:</span><br> 1825 netcdf<br> </td> </tr> <tr> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">C<br> </td> <td style="vertical-align: top;">In <span style="font-family: monospace;">calc_spectra</span>, n 1826 is assigned the 3874 1827 number of spectra (formerly, the number of spectra was given by n-1)<br> 3875 3876 <br> 3877 1828 <br> 3878 1829 Abort if <span style="font-family: monospace;">poisfft_hybrid</span> 3879 is called in a non-parallel environment.<br> 3880 3881 <br> 3882 3883 Default values of spectra package parameters <span style="font-family: monospace;">pl_spectra</span> and <span style="font-family: monospace;">spectra_direction</span> changed.<br> 3884 3885 </td> 3886 3887 <td style="vertical-align: top;">calc_spectra, check_parameters, 3888 modules, pres, read_3d_binary, write_3d_binary<br> 3889 3890 </td> 3891 3892 </tr> 3893 3894 <tr> 3895 3896 <td style="vertical-align: top;"><br> 3897 3898 </td> 3899 3900 <td style="vertical-align: top;"><br> 3901 3902 </td> 3903 3904 <td style="vertical-align: top;"><br> 3905 3906 </td> 3907 3908 <td style="vertical-align: top;">E<br> 3909 3910 </td> 3911 3912 <td style="vertical-align: top;">Missing argument in ffty 3913 (non-parallel case) added.<br> 3914 3915 <br> 3916 3917 Error in output of particle inertia information removed.<br> 3918 3919 </td> 3920 3921 <td style="vertical-align: top;">poisfft, header<br> 3922 3923 </td> 3924 3925 </tr> 3926 3927 <tr> 3928 3929 <td style="vertical-align: top;">05/07/05<br> 3930 3931 </td> 3932 3933 <td style="vertical-align: top;">SR<br> 3934 3935 </td> 3936 3937 <td style="vertical-align: top;">2.11<br> 3938 3939 </td> 3940 3941 <td style="vertical-align: top;">N/C<br> 3942 3943 </td> 3944 3945 <td style="vertical-align: top;">New cloud physics code 1830 is called in a non-parallel environment.<br> <br> 1831 Default values of spectra package parameters <span style="font-family: monospace;">pl_spectra</span> and <span style="font-family: monospace;">spectra_direction</span> 1832 changed.<br> </td> <td style="vertical-align: top;">calc_spectra, 1833 check_parameters, 1834 modules, pres, read_3d_binary, write_3d_binary<br> </td> </tr> 1835 <tr> <td style="vertical-align: top;"><br> </td> 1836 <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">E<br> </td> <td style="vertical-align: top;">Missing argument in ffty 1837 (non-parallel case) added.<br> <br> 1838 Error in output of particle inertia information removed.<br> </td> 1839 <td style="vertical-align: top;">poisfft, header<br> 1840 </td> </tr> <tr> <td style="vertical-align: top;">05/07/05<br> </td> 1841 <td style="vertical-align: top;">SR<br> </td> 1842 <td style="vertical-align: top;">2.11<br> </td> 1843 <td style="vertical-align: top;">N/C<br> </td> 1844 <td style="vertical-align: top;">New cloud physics code 3946 1845 implemented, which explicitly simulates single cloud droplets, 3947 1846 including droplet growth by condensation and collision. Using this code 3948 requires to set the new inipar-parameter <span style="font-family: monospace;">cloud_droplets</span> = .TRUE. and to use the particle package (by using the <span style="font-weight: bold;">mrun</span>-option <span style="font-family: monospace;">-p particles</span>). Coupling between droplets and the thermodynamic quantities (potential 3949 temperature and specific humidity) is regarded in the new subroutine <span style="font-family: monospace;">interaction_droplets_ptq</span>. The real number of droplets in a grid cell can steered by the new inipar-parameter <span style="font-family: monospace;">initial_weighting_factor</span>.<br> 3950 3951 3952 <br> 3953 3954 3955 <span style="font-weight: bold;">This code will be further extended and probably changed in the near future, so the current status is that of a test version!</span><br> 3956 3957 3958 <br> 3959 3960 3961 In connection with this new code implementation, several changes has been made to the existing code:<br> 3962 3963 3964 <br> 3965 3966 3967 The particle (droplet) size is now given as radius (previously diameter). E.g., the name of the respective package parameter <span style="font-family: monospace;">diameter</span> has changed to <span style="font-family: monospace;">radius</span>.<br> 3968 3969 3970 <br> 3971 3972 3973 The current weighting_factor and radius are stored as additional particle attributes in the type structure <span style="font-family: monospace;">particle_type</span>. The MPI data type <span style="font-family: monospace;">mpi_particle_type</span>, needed for exchanging particles between subdomains, has been changed correspondingly.<br> 3974 3975 3976 <br> 3977 3978 3979 <span style="font-family: monospace;">gas_constant</span> is renamed <span style="font-family: monospace;">r_d</span>, <span style="font-family: monospace;">latent_heat</span> is renamed <span style="font-family: monospace;">r_v</span>, the allowed string length of <span style="font-family: monospace;">pl2d</span> and <span style="font-family: monospace;">pl3d</span> is increased from 6 to 10.<br> 3980 3981 3982 <br> 3983 3984 3985 New 3d arrays <span style="font-family: monospace;">ql_1</span>, <span style="font-family: monospace;">ql_2</span>, <span style="font-family: monospace;">ql_v</span> and <span style="font-family: monospace;">ql_vp</span>, <span style="font-family: monospace;">ql</span> is now a pointer, new pointer <span style="font-family: monospace;">ql_c</span>.<br> 3986 3987 3988 <br> 3989 3990 3991 Data format for unit 85 (particle data output) is changed.</td> 3992 3993 <td style="vertical-align: top;">advec_particles.f90 1847 requires to set the new inipar-parameter <span style="font-family: monospace;">cloud_droplets</span> 1848 = .TRUE. and to use the particle package (by using the <span style="font-weight: bold;">mrun</span>-option <span style="font-family: monospace;">-p particles</span>). 1849 Coupling between droplets and the thermodynamic quantities (potential 1850 temperature and specific humidity) is regarded in the new subroutine <span style="font-family: monospace;">interaction_droplets_ptq</span>. 1851 The real number of droplets in a grid cell can steered by the new 1852 inipar-parameter <span style="font-family: monospace;">initial_weighting_factor</span>.<br> 1853 <br> <span style="font-weight: bold;">This code 1854 will be further extended and probably changed in the near future, so 1855 the current status is that of a test version!</span><br> <br> 1856 In connection with this new code implementation, several changes has 1857 been made to the existing code:<br> <br> 1858 The particle (droplet) size is now given as radius (previously 1859 diameter). E.g., the name of the respective package parameter <span style="font-family: monospace;">diameter</span> has 1860 changed to <span style="font-family: monospace;">radius</span>.<br> 1861 <br> 1862 The current weighting_factor and radius are stored as additional 1863 particle attributes in the type structure <span style="font-family: monospace;">particle_type</span>. 1864 The MPI data type <span style="font-family: monospace;">mpi_particle_type</span>, 1865 needed for exchanging particles between subdomains, has been changed 1866 correspondingly.<br> <br> <span style="font-family: monospace;">gas_constant</span> is 1867 renamed <span style="font-family: monospace;">r_d</span>, 1868 <span style="font-family: monospace;">latent_heat</span> 1869 is renamed <span style="font-family: monospace;">r_v</span>, 1870 the allowed string length of <span style="font-family: monospace;">pl2d</span> 1871 and <span style="font-family: monospace;">pl3d</span> 1872 is increased from 6 to 10.<br> <br> 1873 New 3d arrays <span style="font-family: monospace;">ql_1</span>, 1874 <span style="font-family: monospace;">ql_2</span>, <span style="font-family: monospace;">ql_v</span> and <span style="font-family: monospace;">ql_vp</span>, <span style="font-family: monospace;">ql</span> is now a 1875 pointer, new pointer <span style="font-family: monospace;">ql_c</span>.<br> 1876 <br> 1877 Data format for unit 85 (particle data output) is changed.</td> <td style="vertical-align: top;">advec_particles.f90 3994 1878 check_open.f90, header.f90, init_3d_model.f90, init_cloud_physics.f90, 3995 1879 init_particles.f90, modules.f90, package_parin.f90, parin.f90, 3996 1880 plot_2d.f90, plot_3d.f90, read_var_list.f90, time_integration.f90, 3997 write_var_list.f90<br> 3998 3999 4000 <br> 4001 4002 4003 <span style="font-weight: bold;">new:</span><br> 4004 4005 4006 interaction_droplets_ptq.f90</td> 4007 4008 </tr> 4009 4010 <tr> 4011 4012 <td style="vertical-align: top;"><br> 4013 4014 </td> 4015 4016 <td style="vertical-align: top;"><br> 4017 4018 </td> 4019 4020 <td style="vertical-align: top;"><br> 4021 4022 </td> 4023 4024 <td style="vertical-align: top;">C<br> 4025 4026 </td> 4027 4028 <td style="vertical-align: top;">Default value of <span style="font-family: monospace;">call_psolver_at_all_substeps</span> is <span style="font-style: italic;">.TRUE.</span>, because small scale waves occured with <span style="font-style: italic;">.FALSE.</span></td> 4029 4030 <td style="vertical-align: top;">modules.f90</td> 4031 4032 </tr> 4033 4034 <tr> 4035 4036 <td style="vertical-align: top;"><br> 4037 4038 </td> 4039 4040 <td style="vertical-align: top;"><br> 4041 4042 </td> 4043 4044 <td style="vertical-align: top;"><br> 4045 4046 </td> 4047 4048 <td style="vertical-align: top;">E<br> 4049 4050 </td> 4051 4052 <td style="vertical-align: top;">2d- and 3d-NetCDF-files are not opened if there is no output for them (<span style="font-family: monospace;">pl2d</span> = <span style="font-family: monospace;">pl3d</span> = <span style="font-style: italic;">' '</span>).<br> 4053 4054 4055 <br> 4056 4057 1881 write_var_list.f90<br> <br> <span style="font-weight: bold;">new:</span><br> 1882 interaction_droplets_ptq.f90</td> </tr> <tr> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">C<br> </td> <td style="vertical-align: top;">Default value of <span style="font-family: monospace;">call_psolver_at_all_substeps</span> 1883 is <span style="font-style: italic;">.TRUE.</span>, 1884 because small scale waves occured with <span style="font-style: italic;">.FALSE.</span></td> <td style="vertical-align: top;">modules.f90</td> </tr> 1885 <tr> <td style="vertical-align: top;"><br> </td> 1886 <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">E<br> </td> <td style="vertical-align: top;">2d- and 3d-NetCDF-files are not 1887 opened if there is no output for them (<span style="font-family: monospace;">pl2d</span> = <span style="font-family: monospace;">pl3d</span> = <span style="font-style: italic;">' '</span>).<br> <br> 4058 1888 No opening of units 101 - 108 in case of missing NetCDF support.</td> 4059 4060 <td style="vertical-align: top;">check_open.f90, netcdf.f90<br> 4061 4062 </td> 4063 4064 </tr> 4065 4066 <tr> 4067 4068 <td style="vertical-align: top;">30/06/05<br> 4069 4070 </td> 4071 4072 <td style="vertical-align: top;">GS<br> 4073 4074 </td> 4075 4076 <td style="vertical-align: top;">2.11a<br> 4077 4078 </td> 4079 4080 <td style="vertical-align: top;">C<br> 4081 4082 </td> 4083 4084 <td style="vertical-align: top;">Scalars ug and vg have been 1889 <td style="vertical-align: top;">check_open.f90, netcdf.f90<br> 1890 </td> </tr> <tr> <td style="vertical-align: top;">30/06/05<br> </td> 1891 <td style="vertical-align: top;">GS<br> </td> 1892 <td style="vertical-align: top;">2.11a<br> </td> 1893 <td style="vertical-align: top;">C<br> </td> <td style="vertical-align: top;">Scalars ug and vg have been 4085 1894 changed into 4086 1895 arrays in order to allow the specification of a geostrophic wind that 4087 1896 depends on height (baroclinicity). The initial profiles of the the 4088 1897 u- and v-component of the geostrophic wind are initialized by 4089 specifying the new initialization parameters <span style="font-family: monospace;">ug_surface</span>, 4090 <span style="font-family: monospace;">ug_vertical_gradient</span>, <span style="font-family: monospace;">ug_vertical_gradient_level</span> and <span style="font-family: monospace;">vg_surface</span>, 4091 <span style="font-family: monospace;">vg_vertical_gradient</span>, <span style="font-family: monospace;">vg_vertical_gradient_level</span>, respectively. <br> 4092 4093 <br> 4094 4095 The former initialization parameters <span style="font-family: monospace;">ug</span> and <span style="font-family: monospace;">vg</span> have been removed. <br> 4096 4097 </td> 4098 4099 <td style="vertical-align: top;">boundary_conds, 1898 specifying the new initialization parameters <span style="font-family: monospace;">ug_surface</span>, <span style="font-family: monospace;">ug_vertical_gradient</span>, 1899 <span style="font-family: monospace;">ug_vertical_gradient_level</span> 1900 and <span style="font-family: monospace;">vg_surface</span>, 1901 <span style="font-family: monospace;">vg_vertical_gradient</span>, 1902 <span style="font-family: monospace;">vg_vertical_gradient_level</span>, 1903 respectively. <br> <br> 1904 The former initialization parameters <span style="font-family: monospace;">ug</span> and <span style="font-family: monospace;">vg</span> have been 1905 removed. <br> </td> <td style="vertical-align: top;">boundary_conds, 4100 1906 check_parameters, coriolis, header, init_1d_model, modules, parin, 4101 1907 prognostic_equations, read_var_list, spline_z, write_var_list<br> 4102 4103 </td> 4104 4105 </tr> 4106 4107 <tr> 4108 4109 <td style="vertical-align: top;"><br> 4110 4111 </td> 4112 4113 <td style="vertical-align: top;"><br> 4114 4115 </td> 4116 4117 <td style="vertical-align: top;"><br> 4118 4119 </td> 4120 4121 <td style="vertical-align: top;">E<br> 4122 4123 </td> 4124 4125 <td style="vertical-align: top;">Error in the calculation of u_0 1908 </td> </tr> <tr> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">E<br> </td> <td style="vertical-align: top;">Error in the calculation of u_0 4126 1909 and 4127 v_0 in production_e.f90 removed.<br> 4128 4129 <br> 4130 4131 <span style="font-weight: bold;">Attention:</span> <br> 4132 1910 v_0 in production_e.f90 removed.<br> <br> <span style="font-weight: bold;">Attention:</span> <br> 4133 1911 The content of the testresult-file (run control output) has changed 4134 due to this modification.<br> 4135 4136 </td> 4137 4138 <td style="vertical-align: top;">production_e<br> 4139 4140 </td> 4141 4142 </tr> 4143 <tr> 4144 <td style="text-align: left; vertical-align: top;">20/10/05</td> 4145 <td style="vertical-align: top;">SR</td> 4146 <td style="vertical-align: top;">2.11b</td> 4147 <td style="vertical-align: top;">N/C</td> 4148 <td style="vertical-align: top;">Collision process for droplet growth completed.<br> 4149 <br> 4150 <span style="font-weight: bold;">The droplet code is still under development and requires further speed optimization!<br> 4151 <br> 4152 </span>Number of particles really used is additionally output on the netcdf particle data file.<span style="font-weight: bold;"></span></td> 4153 <td style="vertical-align: top;">advec_particles, modules, netcdf</td> 4154 </tr> 4155 <tr> 4156 <td style="vertical-align: top;"> </td> 4157 <td style="vertical-align: top;"> </td> 4158 <td style="vertical-align: top;"> </td> 4159 <td style="vertical-align: top;">E</td> 4160 <td style="vertical-align: top;">The 1912 due to this modification.<br> </td> <td style="vertical-align: top;">production_e<br> </td> 1913 </tr> <tr> <td style="text-align: left; vertical-align: top;">20/10/05</td> 1914 <td style="vertical-align: top;">SR</td> <td style="vertical-align: top;">2.11b</td> <td style="vertical-align: top;">N/C</td> <td style="vertical-align: top;">Collision process for droplet 1915 growth completed.<br> <br> <span style="font-weight: bold;">The droplet code is still under 1916 development and requires further speed optimization!<br> <br> 1917 </span>Number of particles really used is additionally output on 1918 the netcdf particle data file.<span style="font-weight: bold;"></span></td> 1919 <td style="vertical-align: top;">advec_particles, modules, 1920 netcdf</td> </tr> <tr> <td style="vertical-align: top;"> </td> <td style="vertical-align: top;"> </td> <td style="vertical-align: top;"> </td> <td style="vertical-align: top;">E</td> <td style="vertical-align: top;">The 4161 1921 last PE in a row is not allowed to have more grid points than the other 4162 1922 PES (only less). Jobs crashed in these cases that the last PE has more 4163 1923 grid points. Therefore, the number of gridpoints along x (and y) must 4164 now meet a special condition which is checked in case that <span style="font-family: monospace;">grid_matching</span> = <span style="font-style: italic;">'match'</span> is used.<br> 4165 <br> 4166 Error removed in calculating y axis data for yz cross sections. Error in output of netcdf yz-cross-sections removed.<br> 4167 <br> 4168 2*r replaced by r in the exponential term of the particle momentum equation.<br> 4169 Error in output of yz-slice information removed.</td> 4170 <td style="vertical-align: top;">advec_particles, header, init_pegrid, netcdf, plot_2d</td> 4171 </tr> 4172 <tr> 4173 <td style="vertical-align: top;">06/12/05</td> 4174 <td style="vertical-align: top;">SR</td> 4175 <td style="vertical-align: top;">2.11c</td> 4176 <td style="vertical-align: top;">N</td> 4177 <td style="vertical-align: top;">Output of ql profile is allowed in case of using cloud droplets.</td> 4178 <td style="vertical-align: top;">check_parameters, flow_statistics, modules</td> 4179 </tr> 4180 <tr> 4181 <td style="vertical-align: top;"> </td> 4182 <td style="vertical-align: top;"> </td> 4183 <td style="vertical-align: top;"> </td> 4184 <td style="vertical-align: top;">E</td> 4185 <td style="vertical-align: top;">The horizontally averaged pt profile is calculated within<span style="font-family: monospace;"> buoyancy </span>only 1924 now meet a special condition which is checked in case that <span style="font-family: monospace;">grid_matching</span> = 1925 <span style="font-style: italic;">'match'</span> is 1926 used.<br> <br> 1927 Error removed in calculating y axis data for yz cross sections. Error 1928 in output of netcdf yz-cross-sections removed.<br> <br> 1929 2*r replaced by r in the exponential term of the particle momentum 1930 equation.<br> 1931 Error in output of yz-slice information removed.</td> <td style="vertical-align: top;">advec_particles, header, 1932 init_pegrid, netcdf, plot_2d</td> </tr> <tr> <td style="vertical-align: top;">06/12/05</td> <td style="vertical-align: top;">SR</td> <td style="vertical-align: top;">2.11c</td> <td style="vertical-align: top;">N</td> <td style="vertical-align: top;">Output of ql profile is allowed 1933 in case of using cloud droplets.</td> <td style="vertical-align: top;">check_parameters, 1934 flow_statistics, modules</td> </tr> <tr> <td style="vertical-align: top;"> </td> <td style="vertical-align: top;"> </td> <td style="vertical-align: top;"> </td> <td style="vertical-align: top;">E</td> <td style="vertical-align: top;">The horizontally averaged pt 1935 profile is calculated within<span style="font-family: monospace;"> 1936 buoyancy </span>only 4186 1937 in case of the first respective intermediate timestep. This is done in 4187 1938 order to spare communication time and to produce identical model 4188 results with jobs which are calling<span style="font-family: monospace;"> flow_statistics </span>at different time intervals.<br> 4189 <br> 4190 Error in netcdf variable declaration for<span style="font-family: monospace;"> statistic_regions </span>/= <span style="font-style: italic;">0</span> removed.</td> 4191 <td style="vertical-align: top;">buoyancy, netcdf</td> 4192 </tr><tr><td style="vertical-align: top;">23/02/06</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.0</td><td style="vertical-align: top;">N</td><td><span style="font-weight: bold;">Version for simulating flows over topography (buildings, mountains, etc.)</span><br style="font-weight: bold;"><br>Besides 1939 results with jobs which are calling<span style="font-family: monospace;"> flow_statistics </span>at 1940 different time intervals.<br> <br> 1941 Error in netcdf variable declaration for<span style="font-family: monospace;"> statistic_regions </span>/= 1942 <span style="font-style: italic;">0</span> removed.</td> 1943 <td style="vertical-align: top;">buoyancy, netcdf</td> 1944 </tr><tr><td style="vertical-align: top;">23/02/06</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.0</td><td style="vertical-align: top;">N</td><td><span style="font-weight: bold;">Version for simulating flows over 1945 topography (buildings, mountains, etc.)</span><br style="font-weight: bold;"><br>Besides 4193 1946 the implementation of topography, this version contains several other 4194 improvements. The most important ones are listed here:<br><ul><li>The cloud droplet code has been optimized (it now runs about 100 times faster than before).</li><li>Output of time-averaged 2d/3d data (sections / volume data) is possible.</li><li>Output of 2d-section data averaged along the direction normal to this section can be done.</li><li>The user-interface has been extended in order to allow 2d-/3d-output of user-defined quantities.</li></ul>The kind of topography to be used is controlled by the new initialization parameter <span style="font-family: monospace;">topography</span>. Allowed values are <span style="font-style: italic;">'flat'</span>, <span style="font-style: italic;">'single_building'</span>, and <span style="font-style: italic;">'read_from_file'</span>. The user can define his own special topography setting within the new user-interface routine <span style="font-family: monospace;">user_init_grid</span>. For the single building case, the building (rectangular) size and position can be controlled with parameters <span style="font-family: monospace;">building_height</span>, <span style="font-family: monospace;">building_length_x</span>, <span style="font-family: monospace;">building_length_y</span>, <span style="font-family: monospace;">building_wall_left</span>, and <span style="font-family: monospace;">building_wall_south</span>.<br> 4195 <br> 4196 The topography code is already optimized (vectorized) for use on the NEC-SX6.<br><br>Output of time-averaged 2d-/3d-data can be enabled by appending the string <span style="font-style: italic;">'_av'</span> to the respective output quantities given with the new runtime-parameter <span style="font-family: monospace;">data_output</span> 1947 improvements. The most important ones are listed here:<br><ul><li>The 1948 cloud droplet code has been optimized (it now runs about 100 times 1949 faster than before).</li><li>Output of time-averaged 2d/3d 1950 data (sections / volume data) is possible.</li><li>Output 1951 of 2d-section data averaged along the direction normal to this section 1952 can be done.</li><li>The user-interface has been extended 1953 in order to allow 2d-/3d-output of user-defined quantities.</li></ul>The 1954 kind of topography to be used is controlled by the new initialization 1955 parameter <span style="font-family: monospace;">topography</span>. 1956 Allowed values are <span style="font-style: italic;">'flat'</span>, 1957 <span style="font-style: italic;">'single_building'</span>, 1958 and <span style="font-style: italic;">'read_from_file'</span>. 1959 The user can define his own special topography setting within the new 1960 user-interface routine <span style="font-family: monospace;">user_init_grid</span>. 1961 For the single building case, the building (rectangular) size 1962 and position can be controlled with parameters <span style="font-family: monospace;">building_height</span>, 1963 <span style="font-family: monospace;">building_length_x</span>, 1964 <span style="font-family: monospace;">building_length_y</span>, 1965 <span style="font-family: monospace;">building_wall_left</span>, 1966 and <span style="font-family: monospace;">building_wall_south</span>.<br> 1967 <br> 1968 The topography code is already optimized (vectorized) for use on the 1969 NEC-SX6.<br><br>Output of time-averaged 2d-/3d-data can be 1970 enabled by appending the string <span style="font-style: italic;">'_av'</span> 1971 to the respective output quantities given with the new 1972 runtime-parameter <span style="font-family: monospace;">data_output</span> 4197 1973 (see further below for the (name)changes in the parameters for data 4198 1974 output). Example: if a time-average of the vertical cross-section of … … 4202 1978 the averaging interval and the temporal distance of the single 4203 1979 instantaneous sections (volumes) entering into this averaging is 4204 controlled by the new runtime-parameters <span style="font-family: monospace;">averaging_interval</span>and <span style="font-family: monospace;">dt_averaging_input</span>. <br><br>An average normal to the direction of a cross-section can be output by setting the respective cross section index to <span style="font-style: italic;">-1</span> (example: <span style="font-family: monospace;">section_xz</span> = <span style="font-style: italic;">-1</span>).<br><br>New user-interface routines allowing the output of user-defined quantities are <span style="font-family: monospace;">user_check_data_output</span>, <span style="font-family: monospace;">user_define_netcdf_grid</span>, <span style="font-family: monospace;">user_data_output_2d</span>, <span style="font-family: monospace;">user_data_output_3d</span>, and <span style="font-family: monospace;">user_3d_data_averaging</span>. These quantities can be calculated using the new <span style="font-family: monospace;">CASE</span>-Entry <span style="font-style: italic;">'after integration'</span> in the user-interface routine <span style="font-family: monospace;">user_actions</span>.<br><br>Optimization 1980 controlled by the new runtime-parameters <span style="font-family: monospace;">averaging_interval</span>and 1981 <span style="font-family: monospace;">dt_averaging_input</span>. 1982 <br><br>An average normal to the direction of a 1983 cross-section can be output by setting the respective cross section 1984 index to <span style="font-style: italic;">-1</span> 1985 (example: <span style="font-family: monospace;">section_xz</span> 1986 = <span style="font-style: italic;">-1</span>).<br><br>New 1987 user-interface routines allowing the output of user-defined quantities 1988 are <span style="font-family: monospace;">user_check_data_output</span>, 1989 <span style="font-family: monospace;">user_define_netcdf_grid</span>, 1990 <span style="font-family: monospace;">user_data_output_2d</span>, 1991 <span style="font-family: monospace;">user_data_output_3d</span>, 1992 and <span style="font-family: monospace;">user_3d_data_averaging</span>. 1993 These quantities can be calculated using the new <span style="font-family: monospace;">CASE</span>-Entry <span style="font-style: italic;">'after integration'</span> 1994 in the user-interface routine <span style="font-family: monospace;">user_actions</span>.<br><br>Optimization 4205 1995 of the cloud droplet code is mainly done by re-sorting all particles 4206 1996 after every timestep in a way that now all particles which belong to 4207 1997 the same grid-box (k,j,i) are stored consecutively (wthout any strides) 4208 in memory.<br> 4209 <br> 4210 Code adaptation for the IBM at the department of atmospheric sciences, Yonsei university (ibmy). Affected subroutines: <span style="font-family: monospace;">fft_xy</span>, <span style="font-family: monospace;">netcdf</span>, <br><br><span style="font-weight: bold;">Further new features:</span><br><ul><li>Use of particle tails now requires setting of the new package-parameter <span style="font-family: monospace;">use_particle_tails</span>. The new parameter skip_particles_for_tail allows only every <span style="font-family: monospace;">skip_particle_for_tail</span>'th parameter having a tail.</li><li>Initialization-parameter <span style="font-family: monospace;">e_min</span> can be used to set a lower limit for the subgrid-scale TKE (i.e. to guarantee a lower limit for the flow's Reynolds-number).</li><li>A 1998 in memory.<br> <br> 1999 Code adaptation for the IBM at the department of atmospheric sciences, 2000 Yonsei university (ibmy). Affected subroutines: <span style="font-family: monospace;">fft_xy</span>, <span style="font-family: monospace;">netcdf</span>, <br><br><span style="font-weight: bold;">Further new features:</span><br><ul><li>Use 2001 of particle tails now requires setting of the new package-parameter <span style="font-family: monospace;">use_particle_tails</span>. 2002 The new parameter skip_particles_for_tail allows only every <span style="font-family: monospace;">skip_particle_for_tail</span>'th 2003 parameter having a tail.</li><li>Initialization-parameter <span style="font-family: monospace;">e_min</span> can be 2004 used to set a lower limit for the subgrid-scale TKE (i.e. to guarantee 2005 a lower limit for the flow's Reynolds-number).</li><li>A 4211 2006 conservation of the volume flow (through the complete xz- and 4212 2007 yz-sections of the total domain) can be enforced by the new 4213 initialization-parameter <span style="font-family: monospace;">conserve_volume_flow</span>.</li><li>Data output can be skipped for a given time interval from simulation start (t=0) using new runtime-parameters <span style="font-family: monospace;">skip_time_dopr</span>, <span style="font-family: monospace;">skip_time_dosp</span>, <span style="font-family: monospace;">skip_time_do2d_xy</span>, <span style="font-family: monospace;">skip_time_do2d_xz</span>, <span style="font-family: monospace;">skip_time_do2d_yz</span>, and <span style="font-family: monospace;">skip_time_do3d</span>.<br> 4214 </li><li>By default, NetCDF output is now using 64-bit offset format 2008 initialization-parameter <span style="font-family: monospace;">conserve_volume_flow</span>.</li><li>Data 2009 output can be skipped for a given time interval from simulation start 2010 (t=0) using new runtime-parameters <span style="font-family: monospace;">skip_time_dopr</span>, 2011 <span style="font-family: monospace;">skip_time_dosp</span>, 2012 <span style="font-family: monospace;">skip_time_do2d_xy</span>, 2013 <span style="font-family: monospace;">skip_time_do2d_xz</span>, 2014 <span style="font-family: monospace;">skip_time_do2d_yz</span>, 2015 and <span style="font-family: monospace;">skip_time_do3d</span>.<br> 2016 </li><li>By default, NetCDF output is now using 64-bit 2017 offset format 4215 2018 (large file support) on most machines. The user can switch to the 4216 32-bit offset format with the new runtime-parameter <span style="font-family: monospace;">netcdf_64bit</span>. Units 111:113 and 116 are opened for NetCDF output of time-averaged 2d-3d data.<br>4217 </li> 4218 2019 32-bit offset format with the new runtime-parameter <span style="font-family: monospace;">netcdf_64bit</span>. 2020 Units 111:113 and 116 are opened for NetCDF output of 2021 time-averaged 2d-3d data.<br> </li> <li>Data 4219 2022 logging routines are added (see file data_log.f90), which can be used 4220 2023 for debugging purposes. Output is done on unit 20 (local file 4221 name/directory name DATA_LOG).</li> 4222 <li>New routine <span style="font-family: monospace;">exchange_horiz_2d_int</span> for ghostpoint exchange of 2d-integer arrays. Routine <span style="font-family: monospace;">exchange_horiz_2d</span> is extended for non-cyclic boundary conditions in the non-parallel case.</li> 4223 <li>Two different methods for calculating the mixing length and 2024 name/directory name DATA_LOG).</li> <li>New routine <span style="font-family: monospace;">exchange_horiz_2d_int</span> 2025 for ghostpoint exchange of 2d-integer arrays. Routine <span style="font-family: monospace;">exchange_horiz_2d</span> 2026 is extended for non-cyclic boundary conditions in the non-parallel case.</li> 2027 <li>Two different methods for calculating the mixing length and 4224 2028 the dissipation can be used in the 1d-model. These are steered by the 4225 new initializing-parameters <span style="font-family: monospace;">dissipation_1d</span> and <span style="font-family: monospace;">mixing_length_1d</span>.<br> 4226 </li> 4227 </ul></td><td style="vertical-align: top;">advec_particles, advec_s_pw, 2029 new initializing-parameters <span style="font-family: monospace;">dissipation_1d</span> 2030 and <span style="font-family: monospace;">mixing_length_1d</span>.<br> 2031 </li> 2032 </ul></td><td style="vertical-align: top;">advec_particles, 2033 advec_s_pw, 4228 2034 advec_s_up, advec_u_pw, advec_u_up, advec_v_pw, advec_v_up, advec_w_pw, 4229 2035 advec_w_up, boundary_conds, buoyancy, check_open, check_parameters, … … 4235 2041 prandtl_fluxes, pres, production_e, prognostic_equations, 4236 2042 read_3d_binary, read_var_list, run_control, time_integration, 4237 user_interface, write_3d_binary, write_var_list<br> 4238 <br> 4239 <span style="font-weight: bold;">new:</span><br> 4240 average_3d_data, data_log, sum_up_3d_data<br> 4241 <br> 4242 <span style="font-weight: bold;">renamed:<br> 4243 </span>plot_dvrp -> data_output_dvrp<br> 2043 user_interface, write_3d_binary, write_var_list<br> <br> <span style="font-weight: bold;">new:</span><br> 2044 average_3d_data, data_log, sum_up_3d_data<br> <br> <span style="font-weight: bold;">renamed:<br> </span>plot_dvrp 2045 -> data_output_dvrp<br> 4244 2046 plot_spectra -> data_output_spectra<br> 4245 2047 plot_ts -> data_output_tseries<br> 4246 plot_1d -> data_output_profiles<br> 4247 <span style="font-weight: bold;"></span>plot_2d ->data_output_2d<br>2048 plot_1d -> data_output_profiles<br> <span style="font-weight: bold;"></span>plot_2d -> 2049 data_output_2d<br> 4248 2050 plot_3d -> data_output_3d<br> 4249 </td></tr><tr><td style="vertical-align: top;"> </td><td> </td><td> </td><td style="vertical-align: top;">C</td><td>In parallel mode, all PEs are opening the same copy of the NAMELIST-parameter file <span style="font-family: monospace;">PARIN</span>, i.e. from now on, the file attribute "<span style="font-family: monospace;">npe</span>" must be removed from the corresponding file connection statements in the <span style="font-weight: bold;">mrun</span> configuration file (.mrun.config).<br><br>A 2051 </td></tr><tr><td style="vertical-align: top;"> </td><td> </td><td> </td><td style="vertical-align: top;">C</td><td>In 2052 parallel mode, all PEs are opening the same copy of the 2053 NAMELIST-parameter file <span style="font-family: monospace;">PARIN</span>, 2054 i.e. from now on, the file attribute "<span style="font-family: monospace;">npe</span>" must be 2055 removed from the corresponding file connection statements in the <span style="font-weight: bold;">mrun</span> configuration 2056 file (.mrun.config).<br><br>A 4250 2057 large number of parameter names and local filenames have been changed 4251 2058 (see list below). Most of them are parameters concerning data … … 4253 2060 to "data_output" and "do", respectively. In addition to these changes, 4254 2061 the names of many internal PALM parameters and variables have also 4255 changed.<br><br>Restart runs now require setting of the initialization-parameter <span style="font-family: monospace;">initializing_actions</span> = <span style="font-style: italic;">'read_restart_data'</span>.<br> 4256 <br> 2062 changed.<br><br>Restart runs now require setting of the 2063 initialization-parameter <span style="font-family: monospace;">initializing_actions</span> 2064 = <span style="font-style: italic;">'read_restart_data'</span>.<br> 2065 <br> 4257 2066 Concerning particles, the parameters <span style="font-family: monospace;">psl</span>, <span style="font-family: monospace;">psr</span>, <span style="font-family: monospace;">pdx</span>, 4258 2067 etc., which are controlling the position of the particle source and the 4259 2068 number of particles within this source, are now 1d-arrays. This allows 4260 2069 to define different particle sources for different particle groups.<br> 4261 <br> 4262 The number of particle groups has to be set by the user (no automatic calculation from the value of parameter <span style="font-family: monospace;">density_ratio</span> any more)<span style="font-style: italic;"></span>. Variable <span style="font-family: monospace;">idum</span> in type <span style="font-family: monospace;">particle_type</span> has been renamed <span style="font-family: monospace;">tail_id</span>.<br><br>Output for units 15 and 17 flushed in <span style="font-family: monospace;">init_1d_model</span>.<br><br>List of changed parameter names:<br><table style="text-align: left; width: 559px; height: 831px;" border="1" cellpadding="2" cellspacing="2"><tbody><tr><td style="font-weight: bold;" align="undefined" valign="undefined">Old name</td><td align="undefined" valign="undefined"><span style="font-weight: bold;">New name</span></td><td align="undefined" valign="undefined"><span style="font-weight: bold;">Comments</span></td></tr><tr><td align="undefined" valign="undefined">average_period_pl1d</td><td align="undefined" valign="undefined">averaging_interval_pr</td><td align="undefined" valign="undefined"><br> 4263 </td></tr><tr><td style="vertical-align: top;">average_period_pr1d</td><td style="vertical-align: top;">---</td><td style="vertical-align: top;">averaging of printed profiles is not possible any more</td></tr><tr><td align="undefined" valign="undefined">average_period_sp</td><td align="undefined" valign="undefined">averaging_interval_sp</td><td align="undefined" valign="undefined"><br> 2070 <br> 2071 The number of particle groups has to be set by the user (no automatic 2072 calculation from the value of parameter <span style="font-family: monospace;">density_ratio</span> 2073 any more)<span style="font-style: italic;"></span>. 2074 Variable <span style="font-family: monospace;">idum</span> 2075 in type <span style="font-family: monospace;">particle_type</span> 2076 has been renamed <span style="font-family: monospace;">tail_id</span>.<br><br>Output 2077 for units 15 and 17 flushed in <span style="font-family: monospace;">init_1d_model</span>.<br><br>List 2078 of changed parameter names:<br><table style="text-align: left; width: 559px; height: 831px;" border="1" cellpadding="2" cellspacing="2"><tbody><tr><td style="font-weight: bold;" align="undefined" valign="undefined">Old name</td><td align="undefined" valign="undefined"><span style="font-weight: bold;">New name</span></td><td align="undefined" valign="undefined"><span style="font-weight: bold;">Comments</span></td></tr><tr><td align="undefined" valign="undefined">average_period_pl1d</td><td align="undefined" valign="undefined">averaging_interval_pr</td><td align="undefined" valign="undefined"><br> 2079 </td></tr><tr><td style="vertical-align: top;">average_period_pr1d</td><td style="vertical-align: top;">---</td><td style="vertical-align: top;">averaging of printed profiles 2080 is not possible any more</td></tr><tr><td align="undefined" valign="undefined">average_period_sp</td><td align="undefined" valign="undefined">averaging_interval_sp</td><td align="undefined" valign="undefined"><br> 4264 2081 </td></tr><tr><td align="undefined" valign="undefined">dt_average</td><td align="undefined" valign="undefined">dt_averaging_input_pr</td><td align="undefined" valign="undefined"><br> 4265 2082 </td></tr><tr><td align="undefined" valign="undefined">dt_pl1d</td><td align="undefined" valign="undefined">dt_dopr</td><td align="undefined" valign="undefined"><br> … … 4275 2092 </td></tr><tr><td align="undefined" valign="undefined">plts</td><td align="undefined" valign="undefined">data_output_ts</td><td align="undefined" valign="undefined"><br> 4276 2093 </td></tr><tr><td align="undefined" valign="undefined">pl1d</td><td align="undefined" valign="undefined">data_output_pr</td><td align="undefined" valign="undefined"><br> 4277 </td></tr><tr><td style="vertical-align: top;">pl2d</td><td style="vertical-align: top;">data_output</td><td style="vertical-align: top;">old parameters pl2d and pl3d are joined to new parameter data_output</td></tr><tr><td align="undefined" valign="undefined">pl3d</td><td align="undefined" valign="undefined">data_output</td><td align="undefined" valign="undefined"><br> 2094 </td></tr><tr><td style="vertical-align: top;">pl2d</td><td style="vertical-align: top;">data_output</td><td style="vertical-align: top;">old parameters pl2d and pl3d 2095 are joined to new parameter data_output</td></tr><tr><td align="undefined" valign="undefined">pl3d</td><td align="undefined" valign="undefined">data_output</td><td align="undefined" valign="undefined"><br> 4278 2096 </td></tr><tr><td align="undefined" valign="undefined">pl2d_at_begin</td><td align="undefined" valign="undefined">do2d_at_begin</td><td align="undefined" valign="undefined"><br> 4279 2097 </td></tr><tr><td align="undefined" valign="undefined">pl3d_at_begin</td><td align="undefined" valign="undefined">do3d_at_begin</td><td align="undefined" valign="undefined"><br> … … 4292 2110 run_control, time_integration, timestep, write_compressed, 4293 2111 write_var_list<br> 4294 </td></tr><tr><td> </td><td> </td><td> </td><td style="vertical-align: top;">E</td><td>Error in the particle release (defaults of <span style="font-family: monospace;">psl</span>, <span style="font-family: monospace;">psr</span>, <span style="font-family: monospace;">psn</span>, <span style="font-family: monospace;">pss</span>) at the PE boundary has been fixed.<br> 4295 <br> 4296 Error removed getting the variable ids for <span style="font-family: monospace;">_ext</span> for profiles and time series.<br><br>NetCDF calls in <span style="font-family: monospace;">data_output_tseries</span> only if NetCDF output is switched on.<br><br>Output of the initial geostrophic wind profile corrected.<br><br>Maximum possibel timestep for Runge-Kutta-schemes reduced, error concerning initialization of <span style="font-family: monospace;">l_black</span> removed (both in <span style="font-family: monospace;">init_1d_model</span>).<br> 4297 </td><td style="vertical-align: top;">advec_particles, data_output_tseries, header, init_particles, init_1d_model, netcdf<br> 4298 </td></tr><tr><td align="left" valign="top">26/04/06</td><td align="left" valign="top">SR</td><td align="left" valign="top">3.0a</td><td align="left" valign="top">N</td><td align="left" valign="top">OpenMP parallelization for the cache-optimized PALM code.<br><br>OpenMP 2112 </td></tr><tr><td> </td><td> </td><td> </td><td style="vertical-align: top;">E</td><td>Error in 2113 the particle release (defaults of <span style="font-family: monospace;">psl</span>, <span style="font-family: monospace;">psr</span>, <span style="font-family: monospace;">psn</span>, <span style="font-family: monospace;">pss</span>) at the PE 2114 boundary has been fixed.<br> <br> 2115 Error removed getting the variable ids for <span style="font-family: monospace;">_ext</span> for 2116 profiles and time series.<br><br>NetCDF calls in <span style="font-family: monospace;">data_output_tseries</span> 2117 only if NetCDF output is switched on.<br><br>Output of the 2118 initial geostrophic wind profile corrected.<br><br>Maximum 2119 possibel timestep for Runge-Kutta-schemes reduced, error concerning 2120 initialization of <span style="font-family: monospace;">l_black</span> 2121 removed (both in <span style="font-family: monospace;">init_1d_model</span>).<br> 2122 </td><td style="vertical-align: top;">advec_particles, 2123 data_output_tseries, header, init_particles, init_1d_model, netcdf<br> 2124 </td></tr><tr><td align="left" valign="top">26/04/06</td><td align="left" valign="top">SR</td><td align="left" valign="top">3.0a</td><td align="left" valign="top">N</td><td align="left" valign="top">OpenMP parallelization for 2125 the cache-optimized PALM code.<br><br>OpenMP 4299 2126 directives have been added where necessary. Still further tuning will 4300 2127 be necessary to get optimum performance. Other parts of the code than 4301 the cache-optimized code still need to be (OpenMP-) parallelized. <br><br>So far, threads are activated by setting the <span style="font-weight: bold;">mrun</span>-options <span style="font-family: monospace;">-T</span> <span style="font-style: italic;"><number_of_threads></span> and <span style="font-family: monospace;">-O</span>, i.e. when the Option <span style="font-family: monospace;">-O</span> is given, the <span style="font-family: monospace;">-T</span> argument is interpreted as the number of threads to be used. When <span style="font-family: monospace;">-O</span> is given, the number of MPI tasks per node is automatically set to <span style="font-style: italic;">1</span>.<br><br>Beside 2128 the cache-optimized code still need to be (OpenMP-) parallelized. <br><br>So 2129 far, threads are activated by setting the <span style="font-weight: bold;">mrun</span>-options <span style="font-family: monospace;">-T</span> <span style="font-style: italic;"><number_of_threads></span> 2130 and <span style="font-family: monospace;">-O</span>, 2131 i.e. when the Option <span style="font-family: monospace;">-O</span> 2132 is given, the <span style="font-family: monospace;">-T</span> 2133 argument is interpreted as the number of threads to be used. When <span style="font-family: monospace;">-O</span> is given, 2134 the number of MPI tasks per node is automatically set to <span style="font-style: italic;">1</span>.<br><br>Beside 4302 2135 the additional compiler directives, the main code change required by 4303 2136 the OpenMP parallelization is that global sums are now calculated by 4304 2137 first calculating local sums on each thread. As the next step, local 4305 2138 sums are calculated from these thread sums for each MPI task before the 4306 global sum is computed via MPI_ALLREDUCE. For storing the thread 2139 global sum is computed via MPI_ALLREDUCE. For storing the 2140 thread 4307 2141 sums, arrays sums_l and sums_l_l now have three dimensions where the 4308 2142 number of elements of the third dimension is equal to the number of … … 4310 2144 cpu_statistics, diffusivities, flow_statistics, header, init_pegrid, 4311 2145 init_3d_model, modules, prandtl_fluxes, pres, production_e</td></tr><tr><td align="left" valign="top"> </td><td align="left" valign="top">ML</td><td align="left" valign="top"> </td><td align="left" valign="top">C</td><td align="left" valign="top">Provisional 4312 correction for Piacsek & Williams advection scheme: keep u and v 4313 zero one layer below the topography (in case of set_1d-model_profiles).</td><td align="left" valign="top">init_3d_model</td></tr><tr><td align="left" valign="top"> </td><td align="left" valign="top">SR/ML</td><td align="left" valign="top"> </td><td align="left" valign="top">E</td><td align="left" valign="top">In <span style="font-family: monospace;">flow_statistics.f90</span> nzb_s_inner(j,i) replaced by nzb in determination of z_i.<br><br>Errors removed in the computation of the diabatic mixing length (<span style="font-family: monospace;">init_1d_model.f90</span>).<br><br>Error removed in extend mode when checking whether the selected cross sections match those in the already existing NetCDF file (<span style="font-family: monospace;">netcdf.f90</span>).</td><td align="left" valign="top">flow_statistics, init_1d_model, netcdf</td></tr><tr><td align="left" valign="top">02/06/06</td><td align="left" valign="top">SR</td><td align="left" valign="top">3.0b</td><td align="left" valign="top">N</td><td align="left" valign="top">This version is adapted to the Sun Fire X4600 system at TIT (to be used by setting the <span style="font-weight: bold;">mrun</span> option <span style="font-family: monospace;">-h lctit</span>).</td><td align="left" valign="top">cpu_log, init_1d_model, local_tremain, locaö_tremain_ini, run_control</td></tr><tr><td align="left" valign="top"> </td><td align="left" valign="top"> </td><td align="left" valign="top"> </td><td align="left" valign="top">C</td><td align="left" valign="top">NetCDF 2d- 2146 correction for Piacsek & Williams advection scheme: keep u and 2147 v 2148 zero one layer below the topography (in case of set_1d-model_profiles).</td><td align="left" valign="top">init_3d_model</td></tr><tr><td align="left" valign="top"> </td><td align="left" valign="top">SR/ML</td><td align="left" valign="top"> </td><td align="left" valign="top">E</td><td align="left" valign="top">In <span style="font-family: monospace;">flow_statistics.f90</span> 2149 nzb_s_inner(j,i) replaced by nzb in determination of z_i.<br><br>Errors 2150 removed in the computation of the diabatic mixing length (<span style="font-family: monospace;">init_1d_model.f90</span>).<br><br>Error 2151 removed in extend mode when checking whether the selected cross 2152 sections match those in the already existing NetCDF file (<span style="font-family: monospace;">netcdf.f90</span>).</td><td align="left" valign="top">flow_statistics, 2153 init_1d_model, netcdf</td></tr><tr><td align="left" valign="top">02/06/06</td><td align="left" valign="top">SR</td><td align="left" valign="top">3.0b</td><td align="left" valign="top">N</td><td align="left" valign="top">This version is adapted to the Sun Fire X4600 2154 system at TIT (to be used by setting the <span style="font-weight: bold;">mrun</span> option <span style="font-family: monospace;">-h lctit</span>).</td><td align="left" valign="top">cpu_log, init_1d_model, 2155 local_tremain, locaö_tremain_ini, run_control</td></tr><tr><td align="left" valign="top"> </td><td align="left" valign="top"> </td><td align="left" valign="top"> </td><td align="left" valign="top">C</td><td align="left" valign="top">NetCDF 2d- 4314 2156 and 3d-datasets now contain the exact coordinates of the variables 4315 2157 along x and y with respect to the positions where they are defined on 4316 the staggered grid.<br><br>NetCDF datasets of cross sections (e.g. <span style="font-family: monospace;">DATA_2D_XY_NETCDF</span>) are only opened if the respective cross sections are really requested by the settings of parameter <span style="font-family: monospace;">data_output</span>.</td><td align="left" valign="top">check_parameters, data_output_2d, data_output_3d, modules, netcdf, user_interface</td></tr><tr><td align="left" valign="top"> </td><td align="left" valign="top"> </td><td align="left" valign="top"> </td><td align="left" valign="top">E</td><td align="left" valign="top">p is assigned to <span style="font-family: monospace;">to_be_resorted</span> instead of w</td><td align="left" valign="top">data_output_2d</td></tr><tr><td align="left" valign="top">04/08/06</td><td align="left" valign="top">SR</td><td align="left" valign="top">3.1</td><td align="left" valign="top">N</td><td align="left" valign="top">Subgrid 2158 the staggered grid.<br><br>NetCDF datasets of cross 2159 sections (e.g. <span style="font-family: monospace;">DATA_2D_XY_NETCDF</span>) 2160 are only opened if the respective cross sections are really requested 2161 by the settings of parameter <span style="font-family: monospace;">data_output</span>.</td><td align="left" valign="top">check_parameters, 2162 data_output_2d, data_output_3d, modules, netcdf, user_interface</td></tr><tr><td align="left" valign="top"> </td><td align="left" valign="top"> </td><td align="left" valign="top"> </td><td align="left" valign="top">E</td><td align="left" valign="top">p is assigned to <span style="font-family: monospace;">to_be_resorted</span> 2163 instead of w</td><td align="left" valign="top">data_output_2d</td></tr><tr><td align="left" valign="top">04/08/06</td><td align="left" valign="top">SR</td><td align="left" valign="top">3.1</td><td align="left" valign="top">N</td><td align="left" valign="top">Subgrid 4317 2164 scale velocities can (optionally) be included for calculating the 4318 2165 particle advection, using the method of Weil et al. (2004, JAS, 61, 4319 2166 2877-2887). This method is switched on by the new particle package 4320 parameter <span style="font-family: monospace;">use_sgs_for_particles</span>. This also forces the Euler/upstream method to be used for time advancement of the TKE (see new parameter <span style="font-family: monospace;">use_upstream_for_tke</span>). The minimum timestep during the sub-timesteps is controlled by parameter <span style="font-family: monospace;">dt_min_part</span>. The data type <span style="font-family: monospace;">particle_type</span> has additional new attributes <span style="font-family: monospace;">e_m</span>, <span style="font-family: monospace;">dt_sum</span>, and <span style="font-family: monospace;">speed_x/y/z_sgs</span>.<br><br>Other new particle features: output of particle quantities as timeseries in NetCDF format (on local file <span style="font-family: monospace;">DATA_1D_PTS_NETCDF</span>). The output time interval is controlled by parameter <span style="font-family: monospace;">dt_dopts</span>. Particle advection can be switched off after some time using parameter <span style="font-family: monospace;">end_time_prel</span>. More than one particle per point can be started with parameter <span style="font-family: monospace;">particles_per_point</span>.<br><br>Additional parameter in function <span style="font-family: monospace;">random_gauss</span> which limits the range of the created random numbers to five times the standard deviation (=1). Seeds <span style="font-family: monospace;">iran</span> and <span style="font-family: monospace;">iran_part</span> are stored for restart runs.</td><td align="left" valign="top">advec_particles, 2167 parameter <span style="font-family: monospace;">use_sgs_for_particles</span>. 2168 This also forces the Euler/upstream method to be used for time 2169 advancement of the TKE (see new parameter <span style="font-family: monospace;">use_upstream_for_tke</span>). 2170 The minimum timestep during the sub-timesteps is controlled by 2171 parameter <span style="font-family: monospace;">dt_min_part</span>. 2172 The data type <span style="font-family: monospace;">particle_type</span> 2173 has additional new attributes <span style="font-family: monospace;">e_m</span>, 2174 <span style="font-family: monospace;">dt_sum</span>, 2175 and <span style="font-family: monospace;">speed_x/y/z_sgs</span>.<br><br>Other 2176 new particle features: output of particle quantities as timeseries in 2177 NetCDF format (on local file <span style="font-family: monospace;">DATA_1D_PTS_NETCDF</span>). 2178 The output time interval is controlled by parameter <span style="font-family: monospace;">dt_dopts</span>. 2179 Particle advection can be switched off after some time using parameter <span style="font-family: monospace;">end_time_prel</span>. 2180 More than one particle per point can be started with parameter <span style="font-family: monospace;">particles_per_point</span>.<br><br>Additional 2181 parameter in function <span style="font-family: monospace;">random_gauss</span> 2182 which limits the range of the created random numbers to five times the 2183 standard deviation (=1). Seeds <span style="font-family: monospace;">iran</span> 2184 and <span style="font-family: monospace;">iran_part</span> 2185 are stored for restart runs.</td><td align="left" valign="top">advec_particles, 4321 2186 check_open, check_parameters, close_file, diffusion_e, 4322 2187 disturb_heatflux, header, init_3d_model, init_particles, modules, 4323 2188 netcdf, package_parin, parin, prognostic_equations, random_gauss, 4324 2189 read_3d_binary, read_var_list, time_integration, write_3d_binary, 4325 write_var_list<br><br><span style="font-weight: bold;">new:</span><br>data_output_ptseries</td></tr><tr><td align="left" valign="top"> </td><td align="left" valign="top"> </td><td align="left" valign="top"> </td><td align="left" valign="top">C</td><td align="left" valign="top">Call of subroutine <span style="font-family: monospace;">header</span> is moved after call of subroutine <span style="font-family: monospace;">init_3d_model</span>. Generation of <span style="font-family: monospace;">run_description_header</span> is moved from routines <span style="font-family: monospace;">palm</span> and <span style="font-family: monospace;">header</span> to <span style="font-family: monospace;">check_parameters</span>. Determination of the number of particle groups is moved from <span style="font-family: monospace;">header</span> to <span style="font-family: monospace;">init_particles</span>.<br><br><span style="font-family: monospace;">izuf</span> renamed <span style="font-family: monospace;">iran</span>.</td><td align="left" valign="top">check_parameters, disturb_field, disturb_heatflux, header, init_particles, modules, palm</td></tr><tr><td align="left" valign="top"> </td><td align="left" valign="top"> </td><td align="left" valign="top"> </td><td align="left" valign="top">E</td><td align="left" valign="top">Variables 2190 write_var_list<br><br><span style="font-weight: bold;">new:</span><br>data_output_ptseries</td></tr><tr><td align="left" valign="top"> </td><td align="left" valign="top"> </td><td align="left" valign="top"> </td><td align="left" valign="top">C</td><td align="left" valign="top">Call of subroutine <span style="font-family: monospace;">header</span> is moved 2191 after call of subroutine <span style="font-family: monospace;">init_3d_model</span>. 2192 Generation of <span style="font-family: monospace;">run_description_header</span> 2193 is moved from routines <span style="font-family: monospace;">palm</span> 2194 and <span style="font-family: monospace;">header</span> 2195 to <span style="font-family: monospace;">check_parameters</span>. 2196 Determination of the number of particle groups is moved from <span style="font-family: monospace;">header</span> to <span style="font-family: monospace;">init_particles</span>.<br><br><span style="font-family: monospace;">izuf</span> renamed <span style="font-family: monospace;">iran</span>.</td><td align="left" valign="top">check_parameters, 2197 disturb_field, disturb_heatflux, header, init_particles, modules, palm</td></tr><tr><td align="left" valign="top"> </td><td align="left" valign="top"> </td><td align="left" valign="top"> </td><td align="left" valign="top">E</td><td align="left" valign="top">Variables 4326 2198 do2d_unit and do3d_unit now defined as 2d-arrays. Before, in some 4327 2199 cases, the units of variables in the 2d- and 3d-NetCDF-files have been 4328 wrong.<br><br>In routine <span style="font-family: monospace;">poisfft</span> default setting of the thread number tn in case of not using OpenMP.<br><br>Error removed in the non-parallel part of routine <span style="font-family: monospace;">flow_statistics</span> (number of arguments of array <span style="font-family: monospace;">sums_l</span>).<br><br>Error removed which appeared if the user had decreased the value of <span style="font-family: monospace;">dt_dopr</span> within a job chain and if simulataneous time averaging of profiles was switched on.</td><td align="left" valign="top">check_parameters, data_output_profiles, flow_statistics, modules, netcdf, poisfft</td></tr><tr><td style="vertical-align: top;">22/08/06</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.1a</td><td style="vertical-align: top;">N</td><td style="vertical-align: top;">In case of more than one particle group, seperate output of timeseries for each of the groups.<br><br>New initial parameter <span style="font-family: monospace;">dz_max</span>, which limits the vertical gridspacing in case of a vertically stretched grid.</td><td style="vertical-align: top;">data_output_ptseries, header, init_grid, modules, netcdf, parin, read_var_list, write_var_list</td></tr><tr><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Disturbances are imposed only for the last Runge-Kutta-substep. <br><span style="font-weight: bold;">Attention:</span> <br> 4329 2200 wrong.<br><br>In routine <span style="font-family: monospace;">poisfft</span> default 2201 setting of the thread number tn in case of not using OpenMP.<br><br>Error 2202 removed in the non-parallel part of routine <span style="font-family: monospace;">flow_statistics</span> 2203 (number of arguments of array <span style="font-family: monospace;">sums_l</span>).<br><br>Error 2204 removed which appeared if the user had decreased the value of <span style="font-family: monospace;">dt_dopr</span> within 2205 a job chain and if simulataneous time averaging of profiles was 2206 switched on.</td><td align="left" valign="top">check_parameters, 2207 data_output_profiles, flow_statistics, modules, netcdf, poisfft</td></tr><tr><td style="vertical-align: top;">22/08/06</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.1a</td><td style="vertical-align: top;">N</td><td style="vertical-align: top;">In case of more than one 2208 particle group, seperate output of timeseries for each of the groups.<br><br>New 2209 initial parameter <span style="font-family: monospace;">dz_max</span>, 2210 which limits the vertical gridspacing in case of a vertically stretched 2211 grid.</td><td style="vertical-align: top;">data_output_ptseries, 2212 header, init_grid, modules, netcdf, parin, read_var_list, write_var_list</td></tr><tr><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Disturbances are imposed only 2213 for the last Runge-Kutta-substep. <br><span style="font-weight: bold;">Attention:</span> <br> 4330 2214 The content of the testresult-file (run control output) has changed 4331 due to this modification.<br><br>Output of xz and yz cross sections now up to gridpoint nzt+1.<br><br>Default settings of particle start positions changed.</td><td style="vertical-align: top;">check_open, close_file, data_output_2d, init_particles, netcdf, time_integration</td></tr><tr><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;">E</td><td style="vertical-align: top;">Initialisation of all tendency arrays (t.._m) needed for the Runge-Kutta schemes.<br><br>Bugfix in index array initialization for line- or point-like topography structures.<br><br>Bugfix: yv coordinates are now used for yz cross sections, where neccessary.</td><td style="vertical-align: top;">init_3d_model, init_grid, netcdf</td></tr><tr><td style="vertical-align: top;">12/02/07</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.1b</td><td style="vertical-align: top;"> N</td><td style="vertical-align: top;"><span style="font-weight: bold;">Attention:</span><br>This is the last revision before switching to the subversion revison control system!<br><br>Informative output to the job protocoll in case of restarts.</td><td style="vertical-align: top;">check_for_restart</td></tr><tr><td> </td><td> </td><td> </td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Values of environment variables are now read from local file <span style="font-family: monospace;">ENVPAR</span> instead of getting them by a system call.<br><br>For interpolation in <span style="font-family: monospace;">advec_particles</span>, allways level k is used.</td><td style="vertical-align: top;">advec_particles, local_tremain_ini, parin</td></tr><tr><td> </td><td> </td><td> </td><td style="vertical-align: top;">E</td><td style="vertical-align: top;">Several bugfixes in <span style="font-family: monospace;">advec_particles</span>: new particles are released only if <span style="font-family: monospace;">end_time_prel</span> > <span style="font-family: monospace;">simulated_time</span>, transfer 2215 due to this modification.<br><br>Output of xz and yz cross 2216 sections now up to gridpoint nzt+1.<br><br>Default settings 2217 of particle start positions changed.</td><td style="vertical-align: top;">check_open, close_file, 2218 data_output_2d, init_particles, netcdf, time_integration</td></tr><tr><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;">E</td><td style="vertical-align: top;">Initialisation of all tendency 2219 arrays (t.._m) needed for the Runge-Kutta schemes.<br><br>Bugfix 2220 in index array initialization for line- or point-like topography 2221 structures.<br><br>Bugfix: yv coordinates are now used for 2222 yz cross sections, where neccessary.</td><td style="vertical-align: top;">init_3d_model, init_grid, netcdf</td></tr><tr><td style="vertical-align: top;">12/02/07</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.1b</td><td style="vertical-align: top;"> N</td><td style="vertical-align: top;"><span style="font-weight: bold;">Attention:</span><br>This 2223 is the last revision before switching to the subversion revison control 2224 system!<br><br>Informative output to the job protocoll in 2225 case of restarts.</td><td style="vertical-align: top;">check_for_restart</td></tr><tr><td> </td><td> </td><td> </td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Values of environment variables 2226 are now read from local file <span style="font-family: monospace;">ENVPAR</span> 2227 instead of getting them by a system call.<br><br>For 2228 interpolation in <span style="font-family: monospace;">advec_particles</span>, 2229 allways level k is used.</td><td style="vertical-align: top;">advec_particles, 2230 local_tremain_ini, parin</td></tr><tr><td> </td><td> </td><td> </td><td style="vertical-align: top;">E</td><td style="vertical-align: top;">Several bugfixes in <span style="font-family: monospace;">advec_particles</span>: 2231 new particles are released only if <span style="font-family: monospace;">end_time_prel</span> 2232 > <span style="font-family: monospace;">simulated_time</span>, transfer 4332 2233 of particles when x < -0.5*dx (0.0 before), etc., index i,j used 4333 2234 instead of cartesian (x,y) coordinate to check for transfer because 4334 2235 this failed under very rare conditions, calculation of number of 4335 particles with same radius as the current particle (cloud droplet code).<br><br>Allocation of <span style="font-family: monospace;">tail_mask</span> and <span style="font-family: monospace;">new_tail_id</span> in case of restart-runs. "__" added in a cpp-directive. (both <span style="font-family: monospace;">init_particles</span>)</td><td style="vertical-align: top;">advec_particles, init_particles</td></tr><tr><td style="vertical-align: top;">02/03/07</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.1c</td><td style="vertical-align: top;">N</td><td style="vertical-align: top;">A heatflux can be prescribed at the top with new inipar parameters <span style="font-family: monospace;">top_heatflux</span> and <span style="font-family: monospace;">use_top_fluxes</span>. New 2d-arrays <span style="font-family: monospace;">qswst</span>, <span style="font-family: monospace;">qswst_m</span>, <span style="font-family: monospace;">tswst</span> and <span style="font-family: monospace;">tswst_m</span> are used to store this flux. Use of fluxes are controlled with new index variable <span style="font-family: monospace;">nzt_diff</span>. A Neumann boundary condition for temperature can be applied under these conditions.<br><br>Additionally, a Dirichlet condition for temperature can be used at the top.<br><br><span style="font-weight: bold;">Attention:</span> <br> 4336 2236 particles with same radius as the current particle (cloud droplet code).<br><br>Allocation 2237 of <span style="font-family: monospace;">tail_mask</span> 2238 and <span style="font-family: monospace;">new_tail_id</span> 2239 in case of restart-runs. "__" added in a cpp-directive. (both <span style="font-family: monospace;">init_particles</span>)</td><td style="vertical-align: top;">advec_particles, init_particles</td></tr><tr><td style="vertical-align: top;">02/03/07</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.1c</td><td style="vertical-align: top;">N</td><td style="vertical-align: top;">A heatflux can be prescribed at 2240 the top with new inipar parameters <span style="font-family: monospace;">top_heatflux</span> 2241 and <span style="font-family: monospace;">use_top_fluxes</span>. 2242 New 2d-arrays <span style="font-family: monospace;">qswst</span>, 2243 <span style="font-family: monospace;">qswst_m</span>, 2244 <span style="font-family: monospace;">tswst</span> 2245 and <span style="font-family: monospace;">tswst_m</span> 2246 are used to store this flux. Use of fluxes are controlled with new 2247 index variable <span style="font-family: monospace;">nzt_diff</span>. 2248 A Neumann boundary condition for temperature can be applied under these 2249 conditions.<br><br>Additionally, a Dirichlet condition for 2250 temperature can be used at the top.<br><br><span style="font-weight: bold;">Attention:</span> <br> 4337 2251 The content of the testresult-file (run control output) has changed 4338 2252 due to this modification.</td><td style="vertical-align: top;">check_parameters, 4339 2253 diffusion_s, flow_statistics, header, init_grid, init_3d_model, 4340 2254 modules, parin, production_e, prognostic_equations, read_var_list, 4341 read_3d_binary, swap_timelevel, write_var_list, write_3d_binary</td></tr><tr><td> </td><td> </td><td> </td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Prognostic equations for all scalars are now solved up to gridpoint <span style="font-family: monospace;">nzt</span> (formerly <span style="font-family: monospace;">nzt-1</span>). Boundary conditions for scalars at top adjusted respectively (now applied only at <span style="font-family: monospace;">nzt+1</span>).<br><br>The default top boundary condition for temperature has been renamed to <span style="font-style: italic;">'initial_gradient'</span>.<br><br>Calls of <span style="font-family: monospace;">dvrp_output_local</span>, which were commented out for a long time, are now activated for all streams.</td><td style="vertical-align: top;">advec_s_pw, 2255 read_3d_binary, swap_timelevel, write_var_list, write_3d_binary</td></tr><tr><td> </td><td> </td><td> </td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Prognostic equations for all 2256 scalars are now solved up to gridpoint <span style="font-family: monospace;">nzt</span> (formerly <span style="font-family: monospace;">nzt-1</span>). 2257 Boundary conditions for scalars at top adjusted respectively (now 2258 applied only at <span style="font-family: monospace;">nzt+1</span>).<br><br>The 2259 default top boundary condition for temperature has been renamed to <span style="font-style: italic;">'initial_gradient'</span>.<br><br>Calls 2260 of <span style="font-family: monospace;">dvrp_output_local</span>, 2261 which were commented out for a long time, are now activated for all 2262 streams.</td><td style="vertical-align: top;">advec_s_pw, 4342 2263 boundary_conds, calc_precipitation, check_parameters, diffusion_e, 4343 2264 diffusion_s, impact_of_latent_heat, init_dvrp, init_pt_anomaly, 4344 modules, production_e, prognostic_equations, spline_z</td></tr><tr><td> </td><td> </td><td> </td><td style="vertical-align: top;">E</td><td style="vertical-align: top;">Bugfix: 3d-array <span style="font-family: monospace;">p</span> is not a pointer any more.<br><br>Bugfix in <span style="font-family: monospace;">init_particles</span>: <span style="font-family: monospace;">MPI_REAL</span> argument in <span style="font-family: monospace;">MPI_ALLREDUCE</span> replaced by<br><span style="font-family: monospace;">MPI_INTEGER</span> (caused error on NEC only).<br><br>Bugfix: <span style="font-family: monospace;">ddzw</span> now dimensioned <span style="font-family: monospace;">1:nzt"+1".</span></td><td style="vertical-align: top;">diffusion_e, diffusion_s, diffusion_u, diffusion_v, diffusion_w, init_particles, modules</td></tr> 4345 4346 4347 </tbody> 4348 </table> 4349 4350 <b><blink>Attention:</blink></b> If<tt> make </tt>is not used 2265 modules, production_e, prognostic_equations, spline_z</td></tr><tr><td> </td><td> </td><td> </td><td style="vertical-align: top;">E</td><td style="vertical-align: top;">Bugfix: 3d-array <span style="font-family: monospace;">p</span> is not a 2266 pointer any more.<br><br>Bugfix in <span style="font-family: monospace;">init_particles</span>: 2267 <span style="font-family: monospace;">MPI_REAL</span> 2268 argument in <span style="font-family: monospace;">MPI_ALLREDUCE</span> 2269 replaced by<br><span style="font-family: monospace;">MPI_INTEGER</span> 2270 (caused error on NEC only).<br><br>Bugfix: <span style="font-family: monospace;">ddzw</span> now 2271 dimensioned <span style="font-family: monospace;">1:nzt"+1".</span></td><td style="vertical-align: top;">diffusion_e, diffusion_s, 2272 diffusion_u, diffusion_v, diffusion_w, init_particles, modules</td></tr> 2273 </tbody> 2274 </table> <b><blink>Attention:</blink></b> 2275 If<tt> make </tt>is not used 4351 2276 for 4352 2277 compiling, 4353 <b>all routines</b> have to be re-compiled after every change 2278 <b>all routines</b> have to be re-compiled after every 2279 change 4354 2280 to <b>modules.f90!</b> 4355 <br> 4356 4357 4358 <br> 4359 4360 4361 <h2><a name="Kapitel2.0"></a>2.0 How to change the model source 2281 <br> 2282 <br> 2283 <h2><a name="Kapitel2.0"></a>2.0 How 2284 to change the model source 4362 2285 code</h2> 4363 4364 2286 Änderungen, Erweiterungen oder Fehlerkorrekturen am Modellcode 4365 2287 (gemeint 4366 2288 sind alle Modellteile, die nicht zur benutzereigenen Software 4367 2289 gehören) 4368 erfordern vom jeweiligen Benutzer Gruppenzugehörigkeit zur Gruppe 2290 erfordern vom jeweiligen Benutzer Gruppenzugehörigkeit zur 2291 Gruppe 4369 2292 <tt>"palm"</tt>. 4370 2293 Zu dieser Gruppe gehören alle Benutzer, die aktiv an der 4371 2294 Modellentwicklung 4372 2295 mitwirken (im weiteren "Modellarbeitsgruppe" genannt). 4373 <p>Die einzelnen Quelltextdateien des Modells werden mit RCS (Revision 2296 <p>Die einzelnen Quelltextdateien des Modells werden mit RCS 2297 (Revision 4374 2298 Control System) verwaltet. Grundlegende Kenntnisse über dieses 4375 2299 System 4376 2300 werden im weiteren vorausgesetzt. 4377 </p> 4378 4379 <p>Änderungen am Modellcode erfordern die Durchführung 2301 </p><p>Änderungen am Modellcode erfordern die 2302 Durchführung 4380 2303 bestimmter 4381 2304 Aktionen in einer festgelegten Reihenfolge, die nun beschrieben werden 4382 2305 sollen. 4383 </p> 4384 4385 <ol> 4386 4387 <li>Konzeptplanung.</li> 4388 4389 <li>Diskussion der geplanten Änderungen in der 4390 Modellarbeitsgruppe.</li> 4391 4392 <li>Auschecken der benötigten Quelltextdateien, 2306 </p><ol> <li>Konzeptplanung.</li> <li>Diskussion 2307 der geplanten Änderungen in der 2308 Modellarbeitsgruppe.</li> <li>Auschecken der 2309 benötigten Quelltextdateien, 4393 2310 Durchführung der 4394 2311 Änderungen und Test des Modells. Die Modul-Datei <tt>modules.f90</tt> 4395 2312 sollte nur dann ausgecheckt werden, wenn dort aufgeführte 4396 2313 Variablen 4397 verändert werden. Neu geschaffene Variablen sollten während 2314 verändert werden. Neu geschaffene Variablen sollten 2315 während 4398 2316 der 4399 2317 Testphase in ein ebenfalls temporär neu zu schaffendes Modul … … 4403 2321 entsprechender USE-Anweisungen). Die in den einzelnen Dateien 4404 2322 durchgeführten 4405 Arbeiten müssen in den Kopfzeilen (unter "Aktuelle Aenderungen") 2323 Arbeiten müssen in den Kopfzeilen (unter "Aktuelle 2324 Aenderungen") 4406 2325 der 4407 Dateien vermerkt werden. Dies tut man am besten noch <b>bevor</b> man 2326 Dateien vermerkt werden. Dies tut man am besten noch <b>bevor</b> 2327 man 4408 2328 die 4409 2329 eigentlichen Änderungen durchführt (sonst droht das 4410 Vergessen...).</li> 4411 4412 <li>Zum Abschluß der Tests sollte mit dem Modell eine 2330 Vergessen...).</li> <li>Zum Abschluß der Tests 2331 sollte mit dem Modell eine 4413 2332 Kontrollrechnung 4414 2333 durchgeführt und die Ergebnisse mit denen des sogenannten … … 4418 2337 durchzuführen, die auch zum Testen der Modellinstallation 4419 2338 verwendet 4420 werden (s. <a href="Anwendung/Kapitel_5.0.html">Kapitel 5.0</a> der 2339 werden (s. <a href="Anwendung/Kapitel_5.0.html">Kapitel 2340 5.0</a> der 4421 2341 Anwenderdokumentation). 4422 2342 Die Ergebnisdatei mit dem lokalen Namen <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM-1/Dokumentationen/Anwendung/Kapitel_3.4.html#RUN_CONTROL">RUN_CONTROL</a> … … 4430 2350 so 4431 2351 ist unbedingt eine Diskussion in der Modellarbeitsgruppe notwendig.</li> 4432 4433 <li>Nach Durchführung der Tests müssen beiBedarf die2352 <li>Nach Durchführung der Tests müssen bei 2353 Bedarf die 4434 2354 vorläufig 4435 im temporären Modul eingetragenen Variablen in die eigentlichen 2355 im temporären Modul eingetragenen Variablen in die 2356 eigentlichen 4436 2357 Module 4437 2358 des Modells übertragen werden. Dazu ist die Modul-Datei <tt>modules.f90</tt> 4438 2359 auszuchecken. Nach Eintrag der Variablen sollte das Modell noch einmal 4439 testweise übersetzt werden, um Sytaxfehler auszuschließen 2360 testweise übersetzt werden, um Sytaxfehler 2361 auszuschließen 4440 2362 (evtl. 4441 2363 auch noch einmal einen Probelauf starten). Die Verwendung des 4442 2364 temporären 4443 2365 Moduls ist aus allen entsprechenden Programmteilen zu entfernen.</li> 4444 4445 <li>Vergabe einer neuen Versionsnummer (Variable <tt>version</tt>im2366 <li>Vergabe einer neuen Versionsnummer (Variable <tt>version</tt> 2367 im 4446 2368 Modul <tt>steuer</tt>), 4447 Bsp.: 2.3c. Kleinere Änderungen werden durch angehängten 2369 Bsp.: 2.3c. Kleinere Änderungen werden durch 2370 angehängten 4448 2371 kleinen 4449 Buchstaben kenntlich gemacht, größere Änderungen durch 2372 Buchstaben kenntlich gemacht, größere 2373 Änderungen durch 4450 2374 die Zahl hinter dem Punkt. Eine Inkrementierung der Zahl vor dem Punkt 4451 geschieht nur bei ganz grundsätzlichen Revisionen des Modells, die4452 weite Teile des Modellcodes betreffen.</li> 4453 4454 <li>Nunkönnen die ausgecheckten Dateien wieder eingecheckt2375 geschieht nur bei ganz grundsätzlichen Revisionen des Modells, 2376 die 2377 weite Teile des Modellcodes betreffen.</li> <li>Nun 2378 können die ausgecheckten Dateien wieder eingecheckt 4455 2379 werden. 4456 2380 Vor … … 4462 2386 der vom RCS-System vergebenen Versionsnummer (dafür sorgt die 4463 2387 $Log: 4464 ... $ - Zeile, die in jeder Quelltextdatei steht).</li> 4465 4466 <li>Eine Zusammenfassung der durchgeführten Arbeiten ist im <a href="#Kapitel1.0">Kapitel 4467 1.0</a> einzutragen.</li> 4468 4469 <li>Die Modellarbeitsgruppe ist über den Abschluß der 2388 ... $ - Zeile, die in jeder Quelltextdatei steht).</li> <li>Eine 2389 Zusammenfassung der durchgeführten Arbeiten ist im <a href="#Kapitel1.0">Kapitel 2390 1.0</a> einzutragen.</li> <li>Die 2391 Modellarbeitsgruppe ist über den Abschluß der 4470 2392 durchgeführten 4471 2393 Arbeiten zu informieren.</li> 4472 4473 </ol> 4474 4475 <h2> 4476 <a name="Kapitel3.0"></a>3.0 Description of selected parts of the 2394 </ol><h2> 2395 <a name="Kapitel3.0"></a>3.0 Description 2396 of selected parts of the 4477 2397 model source code</h2> 4478 4479 2398 <b>This chapter is still under construction (Feb. 04).</b> 4480 <p>Dieser Abschnitt verweist auf vorhandene Beschreibungen der genauen 4481 Funktionsweise einzelner, spezieller Modellteile. Dazu gehört z.B. 2399 <p>Dieser Abschnitt verweist auf vorhandene Beschreibungen der 2400 genauen 2401 Funktionsweise einzelner, spezieller Modellteile. Dazu gehört 2402 z.B. 4482 2403 die Beschreibung bestimmter numerischer Verfahren (auch ihr 4483 2404 theoretischer … … 4488 2409 Aufgaben (z.B. der Berechnung horizontal gemittelter 4489 2410 Größen). 4490 <br> 4491 4492 Die Beschreibungen liegen in jeweils getrennten Dokumenten, in der 2411 <br>Die Beschreibungen liegen in jeweils getrennten Dokumenten, 2412 in der 4493 2413 Regel als TeX- und/oder Postscript-Datei vor und sind über die 4494 2414 folgenden 4495 2415 Links zugänglich. 4496 <br> 4497 2416 <br> 2417 </p><ul> <li> <a name="UPS"></a>upstream-spline 2418 advection scheme ( <a href="Upstream_Spline/ups.ps">.ps</a> 2419 )</li> <li> <a name="Cloud_physics"></a>cloud 2420 physics module ( <a href="Cloud_physics/wolken.pdf">.pdf</a> 2421 )</li> <li> <a name="Bitkompression"></a>Datenkompression 2422 mit dem Verfahren 2423 der <a href="Bitkompression/kompression.html">Bitverschiebung</a></li> 2424 <li> <a name="Runge-kutta"></a>Runge-Kutta 2425 time integration scheme ( <a href="numerik.heiko/zeitschrittverfahren.pdf">.pdf</a> 2426 )</li> 2427 </ul><br> 4498 2428 4499 </p> 4500 4501 <ul> 4502 4503 <li> <a name="UPS"></a>upstream-spline advection scheme ( <a href="Upstream_Spline/ups.ps">.ps</a> 4504 )</li> 4505 4506 <li> <a name="Cloud_physics"></a>cloud physics module ( <a href="Cloud_physics/wolken.pdf">.pdf</a> 4507 )</li> 4508 4509 <li> <a name="Bitkompression"></a>Datenkompression mit dem Verfahren 4510 der <a href="Bitkompression/kompression.html">Bitverschiebung</a></li> 4511 4512 <li> <a name="Runge-kutta"></a>Runge-Kutta time integration scheme ( 4513 <a href="numerik.heiko/zeitschrittverfahren.pdf">.pdf</a> 4514 )</li> 4515 4516 </ul> 4517 4518 <br> 4519 4520 4521 <br> 4522 4523 4524 <br> 4525 4526 4527 <br> 4528 4529 4530 <br> 4531 4532 4533 <br> 4534 4535 4536 <br> 4537 4538 4539 <br> 4540 4541 2429 <br> 2430 <br> 2431 <br> 2432 <br> 2433 <br> 2434 <br> 2435 <br> 4542 2436 </body></html>
Note: See TracChangeset
for help on using the changeset viewer.