Changeset 392 for palm/trunk/DOC/tec
- Timestamp:
- Sep 24, 2009 10:39:14 AM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/DOC/tec/technical_documentation.html
r238 r392 1 1 <!DOCTYPE html PUBLIC "-//w3c//dtd html 4.0 transitional//en"> 2 2 <html><head> 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 documentation</title></head> 4 < body><font size="+4">PALM 3 4 <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 documentation</title></head><body><font size="+4">PALM 5 5 </font>a 6 6 <b>pa</b>rallelized … … 12 12 <br><table nosave="" cellpadding="0" cellspacing="0"> <caption> <br> </caption><tbody> 13 13 </tbody><tbody> </tbody> <tbody> <tr nosave=""> <td nosave=""><b>Current 14 model version:</b></td> <td><span style="font-weight: bold;">3. 6 (rev 235)</span></td> </tr>14 model version:</b></td> <td><span style="font-weight: bold;">3.7 (rev 392)</span></td> </tr> 15 15 <tr nosave=""> <td nosave=""><b>Last 16 16 change of this document</b>: </td> <td nosave=""><b>$Id$</b></td> </tr> … … 28 28 the model source code <br> 29 29 </p> 30 </ul><b>Please note: Parts of this document are still in 31 German. We 32 apologize 33 for any inconvenience.</b> 34 <br> 30 </ul><b></b><br> 31 35 32 <h2><a name="Kapitel1.0"></a>1.0 36 33 History of changes</h2> … … 50 47 <td><b>description</b></td> <td><b>affected 51 48 routines</b></td> </tr> <tr nosave="" valign="top"> <td>06/07/98</td> <td>SR</td> 52 <td>1.0</td> <td>C</td> <td nosave=""> Änderung49 <td>1.0</td> <td>C</td> <td nosave="">Änderung 53 50 von Default-Werten 54 51 (initializing_actions, end_time, 55 fcl_factor, use_prior_plot1d_parameters), zus ätzliche56 Überprüfung57 von Parametern auf Zul ässigkeit (prandtl_layer, psolver, nx,52 fcl_factor, use_prior_plot1d_parameters), zusätzliche 53 Überprüfung 54 von Parametern auf Zulässigkeit (prandtl_layer, psolver, nx, 58 55 ny, 59 56 nz, 60 fcl_factor), H öhe der Schnittebenen der Horizontalschnitte57 fcl_factor), Höhe der Schnittebenen der Horizontalschnitte 61 58 wird 62 59 auf … … 67 64 Programmteilen benutzt 68 65 (USE). <b>Wichtig</b>: dieses Modul steht in der Datei <tt>module_test.f90. 69 D</tt>iese Datei mu ßin der mrun- und add_library -66 D</tt>iese Datei muß in der mrun- und add_library - 70 67 Konfigurationsdatei 71 68 bei der Variablen <tt>module_files</tt> 72 aufgef ührt werden.</td> <td nosave="">(fast)69 aufgeführt werden.</td> <td nosave="">(fast) 73 70 alle <br> 74 71 Neu: module_test</td> </tr> <tr nosave="" valign="top"> <td>16/07/98</td> <td>SR</td> … … 93 90 <td>Entfernung der Sonderbehandlung der Advektionsterme 94 91 (Piascek-Williams) 95 bei k=1. Terme werden nun überall gleich berechnet. Evtl.92 bei k=1. Terme werden nun überall gleich berechnet. Evtl. 96 93 leichte 97 Erh öhung der Varianzen in Oberflächennähe.</td>94 Erhöhung der Varianzen in Oberflächennähe.</td> 98 95 <td nosave="">advec_u_pw, advec_v_pw, advec_w_pw</td> 99 96 </tr> <tr nosave="" valign="top"> <td nosave="">"</td> <td>SR</td> <td>1.0</td> … … 101 98 adjust_mixing_length) erfolgt nur 102 99 noch bei 103 k=1. Unrealistische sekund äre Maxima in den Km-Profilen100 k=1. Unrealistische sekundäre Maxima in den Km-Profilen 104 101 oberhalb 105 102 der 106 Prandtl-Schicht werden damit vermieden. Leichte Änderungen in103 Prandtl-Schicht werden damit vermieden. Leichte Änderungen in 107 104 den 108 oberfl ächennahen Turbulenzgrößen sind die105 oberflächennahen Turbulenzgrößen sind die 109 106 Folge.</td> <td>diffusivities</td> </tr> <tr nosave="" valign="top"> <td>"</td> <td>SR</td> 110 107 <td>1.0</td> <td>B</td> <td>Falsche 111 Berechnung der ben ötigten CPU-Zeit pro Sekunde108 Berechnung der benötigten CPU-Zeit pro Sekunde 112 109 simulierter 113 Zeit bei Fortsetzungsl äufen korrigiert. <p>Berechnung110 Zeit bei Fortsetzungsläufen korrigiert. <p>Berechnung 114 111 von cross_uxmin(1) und cross_uxmax(1) aus den 115 112 Werten von 116 ug und vg entfernt. Lieferte falsche Wertebereiche f ür Plots113 ug und vg entfernt. Lieferte falsche Wertebereiche für Plots 117 114 von 118 115 Vertikalprofilen, … … 137 134 Felder herausgeschrieben. </p> <p>Bei 138 135 SOR-Aufrufen war bisher dp nicht initialisiert. 139 Entsprechend ge ändert.</p> </td> <td>init_3d_model,136 Entsprechend geändert.</p> </td> <td>init_3d_model, 140 137 plot_2d, pres</td> </tr> <tr nosave="" valign="top"> <td>"</td> <td>SR</td> 141 138 <td>1.0</td> <td>C</td> <td nosave="">Ausgabeformate 142 f ür Statistik-Ausgaben auf139 für Statistik-Ausgaben auf 143 140 Datei RUN_CONTROL 144 141 um einige Kommastellen erweitert. <p>Lieferung von 145 ausf ührlicheren Informationen bei142 ausführlicheren Informationen bei 146 143 Programmabbruch 147 144 durch Unterschreitung des minimal erlaubten Zeitschritts.</p> </td> … … 149 146 <td>SR</td> <td>1.0</td> <td>N</td> 150 147 <td>Rechnungen mit geneigter Ebene zur Simulation von 151 Kaltluftabfl üssen148 Kaltluftabflüssen 152 149 sind erlaubt. Gestartet und gesteuert werden sie mit dem 153 150 Intialisierungsparameter <tt>alpha_surface</tt>. 154 151 Realisiert wird die geneigte Ebene durch entsprechende 155 Schr ägstellung156 des Schwerkraftvektors. Dadurch erh ält auch die u-Komponente152 Schrägstellung 153 des Schwerkraftvektors. Dadurch erhält auch die u-Komponente 157 154 der 158 155 Windgeschwindigkeit … … 160 157 Implementierung ist noch nicht abgeschlossen. Die 161 158 Simulation von 162 Kaltluftabfl üssen erfordert bisher noch Eingriffe "von Hand".</p>159 Kaltluftabflüssen erfordert bisher noch Eingriffe "von Hand".</p> 163 160 </td> <td>buoyancy, check_parameters, header, leap_frog, 164 161 modules, … … 166 163 read_var_list, write_var_list</td> </tr> <tr nosave="" valign="top"> <td>15/12/98</td> <td>SR</td> 167 164 <td>1.0</td> <td>N</td> <td>Vereinbarung 168 einer Cache-Layer zur schnelleren bin ären165 einer Cache-Layer zur schnelleren binären 169 166 I/O auf 170 167 T3E-Rechnern. <p>Ermittelung wahrer I/O-Zeiten 171 f ür binäre168 für binäre 172 169 Ein-/Ausgabe auf 173 170 T3E-Rechnern mit TIMEF.</p> </td> <td nosave="">check_open, … … 176 173 <td>SR</td> <td>1.0</td> <td>B</td> 177 174 <td>Bekanntgabe von run_description_header an alle 178 PE's. <p>Diese Programm änderung erfordert auf175 PE's. <p>Diese Programmänderung erfordert auf 179 176 T3E-Rechnern 180 177 zwingend den 181 178 Einsatz der MPI-Version mpt.1.2.0.1. In dieser neuen Version tritt auch 182 der alte Fehler in MPI_ALLREDUCE nicht mehr auf, so da ßab179 der alte Fehler in MPI_ALLREDUCE nicht mehr auf, so daß ab 183 180 jetzt 184 mpt.1.1.0.1 nicht mehr verwendet werden mu ß.</p> </td>181 mpt.1.1.0.1 nicht mehr verwendet werden muß.</p> </td> 185 182 <td>parles</td> </tr> <tr nosave="" valign="top"> <td nosave="">05/02/99</td> 186 183 <td>SR</td> <td>1.1</td> <td>C/N</td> 187 184 <td>Implementierung des Upstream-Spline- Verfahrens zur Advektion 188 185 von Impuls 189 und skalaren Gr ößen (neue UP's, s. rechts).186 und skalaren Größen (neue UP's, s. rechts). 190 187 Berechnung von 191 188 Variablen 192 f ür spezielle Advektionsverfahren in eigenständiges189 für spezielle Advektionsverfahren in eigenständiges 193 190 UP 194 191 ausgelagert … … 196 193 Reines 197 194 Euler-/Upstream-Verfahren 198 kann gerechnet werden. <p>Schnelle bin äre I/O195 kann gerechnet werden. <p>Schnelle binäre I/O 199 196 (s. 15/12/98) wieder abgeschaltet, da 200 197 Dateien 201 nicht wieder einlesbar sind. Druckst örung (p) wird zu Beginn198 nicht wieder einlesbar sind. Druckstörung (p) wird zu Beginn 202 199 auf 0 203 200 gesetzt, damit sie keine Tendenzen im ersten Zeitschritt hervorruft.</p> … … 212 209 spline_x, spline_y, spline_z</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">24/02/99</td> 213 210 <td>SR</td> <td>1.1a</td> <td>C/N</td> 214 <td>Kleinere Erg änzungen / Änderungen215 f ür Upstream211 <td>Kleinere Ergänzungen / Änderungen 212 für Upstream 216 213 - Spline, <br> 217 214 Diverse Verbesserungen bei Berechnung der Diffusionskoeffizienten … … 221 218 auch unterhalb des entsprechenden Gitterpunktes stabil ist), 222 219 Wertebegrenzung 223 bei Spline - Berechnung zur Verhinderung von " Überschwingern"</td>220 bei Spline - Berechnung zur Verhinderung von "Überschwingern"</td> 224 221 <td>check_parameters, diffusion_e, diffusion_pt, diffusivities, 225 222 header, … … 233 230 <td>advec_s_bc, leap_frog</td> </tr> <tr nosave="" valign="top"> <td nosave="">01/03/99</td> 234 231 <td>MS</td> <td>1.1b</td> <td>N</td> 235 <td>Kleinere Erg änzungen für Upstream -232 <td>Kleinere Ergänzungen für Upstream - 236 233 Spline: <ul> <li>Gradientenkontrolle 237 f ür den nicht-parallelen Teil</li> <li>statistische238 Auswertung über den prozentualen Anteil234 für den nicht-parallelen Teil</li> <li>statistische 235 Auswertung über den prozentualen Anteil 239 236 des 240 237 Upstream-Verfahrens … … 243 240 </tr> <tr nosave="" valign="top"> <td>"</td> 244 241 <td>SR</td> <td>1.1b</td> <td>N</td> 245 <td>Datenkompression f ür 3D-Plotausgaben</td> <td nosave="">check_open, check_parameters, close_files, header,242 <td>Datenkompression für 3D-Plotausgaben</td> <td nosave="">check_open, check_parameters, close_files, header, 246 243 modules, 247 244 parin, plot_3d <p><b>Neu:</b> <br> 248 245 write_compressed</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">03/03/99</td> 249 246 <td>SR</td> <td>1.1c</td> <td>B/C</td> 250 <td>Schnelle bin äre I/O auf T3E-Rechnern funktioniert247 <td>Schnelle binäre I/O auf T3E-Rechnern funktioniert 251 248 jetzt. 252 249 Beim 253 bin ären Einlesen dürfen Records nicht mit dem Trick "<tt>READ254 (..) idum</tt>" überlesen werden. Zeitmessungen auf T3E250 binären Einlesen dürfen Records nicht mit dem Trick "<tt>READ 251 (..) idum</tt>" überlesen werden. Zeitmessungen auf T3E 255 252 erfolgen 256 253 jetzt 257 grunds ätzlich mit <tt>TIMEF.</tt></td> <td>check_open,254 grundsätzlich mit <tt>TIMEF.</tt></td> <td>check_open, 258 255 cpu_zeitmessung, init_3d_model, modules, 259 256 write_3d_binary</td> </tr> <tr nosave="" valign="top"> <td nosave="">25/03/99</td> 260 257 <td>SR</td> <td>1.1d</td> <td>N</td> 261 <td>Einbau einer Rayleigh-D ämpfungsschicht sowie258 <td>Einbau einer Rayleigh-Dämpfungsschicht sowie 262 259 flexiblere 263 260 Wahl der 264 numerischen Begrenzer f ür das Upstream - Spline - Verfahren.</td>261 numerischen Begrenzer für das Upstream - Spline - Verfahren.</td> 265 262 <td>check_parameters, header, init_3d_model, leap_frog, modules, 266 263 parin, … … 272 269 <td>1.1e</td> <td>N</td> <td nosave="">Partikelplots 273 270 mit vtk-Graphiksoftware sind auf t3eh 274 m öglich.271 möglich. 275 272 Ausgabe erfolgt im dvr-Format. Bisher noch nicht nutzbar, da noch 276 Einstellm öglichkeiten277 f ür die Partikelquelle(n) fehlen. <b>Wichtig:</b>273 Einstellmöglichkeiten 274 für die Partikelquelle(n) fehlen. <b>Wichtig:</b> 278 275 Die Benutzung 279 276 der 280 vtk-Software erfordert die zus ätzliche mrun-Option <tt>-g</tt>.277 vtk-Software erfordert die zusätzliche mrun-Option <tt>-g</tt>. 281 278 Das 282 279 Hauptprogramm <tt>parles</tt> wird in diesem Fall zum … … 296 293 nx /= ny auf.</td> <td>cpu_zeitmessung, spline_y</td> 297 294 </tr> <tr nosave="" valign="top"> <td nosave="">28/12/99</td> <td>SR</td> <td>1.1f</td> 298 <td>N</td> <td>Plot von Isooberfl ächen mit295 <td>N</td> <td>Plot von Isooberflächen mit 299 296 vtk-Graphiksoftware sind auf 300 297 t3eh 301 m öglich. Ausgabe erfolgt im dvr-Format. Bisher noch nicht298 möglich. Ausgabe erfolgt im dvr-Format. Bisher noch nicht 302 299 allgemein 303 nutzbar, da noch diverse manuelle Eingriffe n ötig sind.</td>300 nutzbar, da noch diverse manuelle Eingriffe nötig sind.</td> 304 301 <td>header, leap_frog, modules, parin, plot_particles, 305 302 read_var_list, write_var_list <p><b>Neu:</b> … … 307 304 plot_isosurface</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">10/01/00</td> 308 305 <td>SR</td> <td>1.2</td> <td>C</td> 309 <td>Randbedingungen f ür w (=0) werden explizit gesetzt,306 <td>Randbedingungen für w (=0) werden explizit gesetzt, 310 307 und 311 308 zwar nicht … … 313 310 bei <tt>k = 314 311 nzt</tt>, 315 wo bislang der vom vorigen Zeitschritt durch den Druckl öser312 wo bislang der vom vorigen Zeitschritt durch den Drucklöser 316 313 ermittelte 317 Wert stand. Auch jetzt wird der Wert bei k = nzt anschlie ßend314 Wert stand. Auch jetzt wird der Wert bei k = nzt anschließend 318 315 durch 319 den Druckl öser modifiziert. Soweit bis jetzt abzusehen, hat316 den Drucklöser modifiziert. Soweit bis jetzt abzusehen, hat 320 317 diese 321 Änderung nur minimale Auswirkungen auf den Prognoseverlauf.318 Änderung nur minimale Auswirkungen auf den Prognoseverlauf. 322 319 Evtl. 323 gibt es st ärkere Auswirkungen bei exzessivem Auftreten von320 gibt es stärkere Auswirkungen bei exzessivem Auftreten von 324 321 Schwerewellen 325 322 am oberen Rand (die dort aber eigentlich sowieso nicht 326 hingeh ören). <br>327 Variablen übergabe an <tt>boundary_conds</tt>323 hingehören). <br> 324 Variablenübergabe an <tt>boundary_conds</tt> 328 325 erfolgt <br> 329 über Modul und nicht mehr per Parameterliste. <p>Bei330 Galilei-Transformation wird jetzt defaultm äßig326 über Modul und nicht mehr per Parameterliste. <p>Bei 327 Galilei-Transformation wird jetzt defaultmäßig 331 328 der geostrophische 332 329 Wind als Translationsgeschwindigkeit verwendet (neuer 333 330 Initialisierungsparameter <tt>use_ug_for_galilei_tr</tt>). 334 Dadurch wird das verst ärkte Auftreten von Schwerewellen bei335 h öheren336 geostrophischen Windgeschwindigkeiten unterdr ückt. </p>331 Dadurch wird das verstärkte Auftreten von Schwerewellen bei 332 höheren 333 geostrophischen Windgeschwindigkeiten unterdrückt. </p> 337 334 <p>3D-Feldvariablen (u, v, w, pt, e, km, kh) sowie einige 338 335 2D-Variablen 339 336 (usws, vsws, rif, shf) sind nun als Pointer deklariert, denen nach 340 337 jedem 341 Zeitschritt neue Targets zugeordnet werden (f ür die342 ben ötigten 338 Zeitschritt neue Targets zugeordnet werden (für die 339 benötigten 343 340 3 Zeitebenen existieren die Targets ..._1 (z.B. u_1), ..._2, ..._3). 344 341 Durch 345 die jeweilige Neuzuordnung der Targets entf ällt das342 die jeweilige Neuzuordnung der Targets entfällt das 346 343 Umspeichern 347 344 der … … 350 347 Die Umordnung der Zeitebenen wurde ins neue UP <tt>swap_timelevel</tt> 351 348 ausgelagert. <br> <b>Wichtig:</b> Die 352 Übersetzungszeiten haben durch die353 Einf ührung349 Übersetzungszeiten haben durch die 350 Einführung 354 351 der Pointer t.w. sehr deutlich zugenommen (<tt>leap_frog</tt> 355 ben ötigt356 jetzt z.B. ca. 60 Sekunden). Au ßerdem benötigt der352 benötigt 353 jetzt z.B. ca. 60 Sekunden). Außerdem benötigt der 357 354 Compiler 358 355 jetzt wesentlich mehr Speicher. Der "memory"-Parameter in den … … 360 357 (.mrun.config, .al.config) muss nun mindestens den Wert 60 (MB) (t3eh, 361 358 t3ej2) bzw. 45 (t3eb) und 96 (vpp) haben. <br> <b>Achtung:</b> 362 bislang wurde diese Änderung nur für359 bislang wurde diese Änderung nur für 363 360 das 364 361 reine Leapfrog-Verfahren getestet! </p> <p>Neue 365 362 Laderoption "<tt>-D preset=nan</tt>" empfohlen (s. 366 363 Beispiel-Konfigurationsdatei). 367 Ebenso empfehlenswert f ür Debug-Läufe: Compileroption364 Ebenso empfehlenswert für Debug-Läufe: Compileroption 368 365 "<tt>-ei</tt>". 369 366 Diese bewirkt einen Programmabbruch, falls mit nicht initialisierten 370 367 Variablen 371 gearbeitet wird (ist aber zeitaufwendig in der Ausf ührung!).</p>368 gearbeitet wird (ist aber zeitaufwendig in der Ausführung!).</p> 372 369 </td> <td>boundary_conds, check_parameters, header, 373 370 init_3d_model, … … 379 376 <td>Der Asselin-Filter wird vor Umordnung der Zeitebenen 380 377 aufgerufen. Dies 381 geschah bisher f älschlicherweise erst danach, sodass der378 geschah bisher fälschlicherweise erst danach, sodass der 382 379 Zeitfilter 383 380 asymmetrisch arbeitete. <br> <b>Achtung:</b> 384 diese Änderung hat Auswirkungen auf den381 diese Änderung hat Auswirkungen auf den 385 382 Prognoseverlauf. 386 Die Testergebnis-Datei wurde entsprechend ge ändert.</td>383 Die Testergebnis-Datei wurde entsprechend geändert.</td> 387 384 <td>leap_frog</td> </tr> <tr nosave=""> 388 385 <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> … … 394 391 zeitlichen 395 392 Filterung unterzogen. Dies wird nun unterbunden. </p> <p>Die 396 Divergenz der Impulsfl üsse an der Prandtl - Schicht -393 Divergenz der Impulsflüsse an der Prandtl - Schicht - 397 394 Obergrenze 398 wird nun über eine ganze Gitterweite gebildet. Die bisherige395 wird nun über eine ganze Gitterweite gebildet. Die bisherige 399 396 Methode 400 397 lieferte im Vergleich mit Ergebnissen anderer LES-Modelle zu 401 gro ße402 Impulsfl üsse. </p> <p>Bei der398 große 399 Impulsflüsse. </p> <p>Bei der 403 400 Berechnung der Diffusion der Temperatur werden die 404 401 Diffusionskoeffizienten 405 nun wieder durch Mittelung f ür die entsprechenden 402 nun wieder durch Mittelung für die entsprechenden 406 403 Gitterpunkte 407 404 bereitgestllt. Das Maximumkriterium, das im Zuge der Implementierung … … 410 407 rasches 411 408 Anwachsen der Grenzschicht mit der Zeit. </p> <p><b>Achtung:</b> 412 diese Änderungen haben Auswirkungen auf409 diese Änderungen haben Auswirkungen auf 413 410 den Prognoseverlauf. 414 Die Testergebnis-Datei wurde entsprechend ge ändert.</p> </td>411 Die Testergebnis-Datei wurde entsprechend geändert.</p> </td> 415 412 <td nosave="" valign="top">init_1d_model, 416 413 asselin_filter, … … 418 415 diffusion_pt</td> </tr> <tr nosave="" valign="top"> <td nosave="">17/02/00</td> 419 416 <td>MS</td> <td>2.0</td> <td>N</td> 420 <td><b>Allgemeine Beschreibung der Änderungen:</b>417 <td><b>Allgemeine Beschreibung der Änderungen:</b> 421 418 <br> 422 vollst ändige Implementierung des wolkenphysikalischen419 vollständige Implementierung des wolkenphysikalischen 423 420 Moduls <p>Zur Implementierung der Wolkenphysik waren 424 421 umfangreiche 425 Änderungen422 Änderungen 426 423 in zahlreichen Programmteilen notwendig (Beschreibung <a href="#Cloud_physics">siehe 427 424 unten</a> ). Im einzelnen wurden folgende Unprogramme 428 ge ändert: <br> <b>modules</b>: 429 <br>1) Bereitstellung neuer 3D-, 2D- und 1D-Felder f ür425 geändert: <br> <b>modules</b>: 426 <br>1) Bereitstellung neuer 3D-, 2D- und 1D-Felder für 430 427 die Prognose 431 428 des Gesamtwassergehaltes <br> … … 436 433 write_var_list, 437 434 write_3d_binary: </b>In/Output 438 der neuen Initialisierungsparamter und Felder f ür439 Fortsetzungsl äufe </font> <br> <font color="#000000"><b>check_parameters:</b>435 der neuen Initialisierungsparamter und Felder für 436 Fortsetzungsläufe </font> <br> <font color="#000000"><b>check_parameters:</b> 440 437 Initialisierung 441 438 des … … 445 442 Temperatur,</font> <br> <font color="#ff0000">lpt</font><font color="#000000">: 446 443 potentielle 447 Fl üssigwassertemperatur,</font> <br> <font color="#ff0000">q</font><font color="#000000">:444 Flüssigwassertemperatur,</font> <br> <font color="#ff0000">q</font><font color="#000000">: 448 445 Gesamtwassergehalt,</font> <br> <font color="#ff0000">qv</font><font color="#000000">: 449 446 spezifische 450 447 Feuchte,</font> <br> <font color="#ff0000">ql</font><font color="#000000">: 451 Fl üssigwassergehalt,</font> <br> <font color="#ff0000">w"vpt", w*vpt*, wvpt</font><font color="#000000">:452 Auftriebsfluss (subskalig, aufgel öst, gesamt),</font> <br>448 Flüssigwassergehalt,</font> <br> <font color="#ff0000">w"vpt", w*vpt*, wvpt</font><font color="#000000">: 449 Auftriebsfluss (subskalig, aufgelöst, gesamt),</font> <br> 453 450 <font color="#ff0000">w"q", w*q*, wq</font><font color="#000000">: 454 451 Wasserfluss (subs., aufg., ges.),</font> <br> <font color="#ff0000">w"qv", w*qv*, wqv</font><font color="#000000">: … … 457 454 Feuchte/Wasserfelder <br> <b>init_1d_model</b>: 458 455 Anpassung an Rechnungen mit Feuchte <br> <b>header</b>: 459 Headerausgaben f ür Feuchtebehandlung und456 Headerausgaben für Feuchtebehandlung und 460 457 Wolkenphysik 461 eingef ügt, außerdem Formatnummern erweitert um Platz462 f ür458 eingefügt, außerdem Formatnummern erweitert um Platz 459 für 463 460 weitere Ausgaben zu schaffen <br> <b>leap_frog:</b> 464 461 <br> … … 471 468 flow_statistics: </b>Anpassung 472 469 an Feuchtebehandlung <br> <b>plot2d, plot3d</b>: 473 Ausgabem öglichkeiten von q und ql <br> <b>buoyancy</b>:470 Ausgabemöglichkeiten von q und ql <br> <b>buoyancy</b>: 474 471 Anpassung des Autriebstermes an 475 472 Feuchtebehandlung 476 durch neue Übergabeparameter: Temperatur + Nummer des473 durch neue Übergabeparameter: Temperatur + Nummer des 477 474 entsprechenden 478 475 mittleren Temperaturprofils <br> <b>diffusion_e</b>: 479 476 Anpassung der Mischungswegberechnung an 480 477 Feuchtebehandlung 481 durch Variation der Übergabeparameter <br> <b>production_e</b>:482 TKE-Produktion durch Auftriebskr äfte an478 durch Variation der Übergabeparameter <br> <b>production_e</b>: 479 TKE-Produktion durch Auftriebskräfte an 483 480 Feuchtebehandlung 484 481 angepasst <br> <b>diffusivities</b>: … … 486 483 Prandtl-Schicht fuer den 487 484 Gesamtwassergehalt 488 eingef ügt </p> <p>Zusätzlich485 eingefügt </p> <p>Zusätzlich 489 486 wurden folgende Unterprogramme NEU 490 eingef ührt: <br> <b>inti_cloud_physics</b>:487 eingeführt: <br> <b>inti_cloud_physics</b>: 491 488 Initialisierung wolkenphysikalischer 492 489 Parameter <br> <b>prognostic_equations</b>: 493 490 Auslagerung der prognostischen 494 491 Gleichungen 495 aus leap_frog, zus ätzliche porgnostische Gleichung496 f ür den492 aus leap_frog, zusätzliche porgnostische Gleichung 493 für den 497 494 Gesamtwassergehalt <br> <b>diffusion_s</b>: 498 495 Berechnung der Diffusion skalarer 499 Gr ößen496 Größen 500 497 (Temperatur und Gesamtwassergehalt), ersetzt diffusion_pt <br> 501 498 <b>calc_liquid_water_content</b>: Berechnung des 502 Fl üssigwassergehalts <br> <b>calc_radiation</b>:499 Flüssigwassergehalts <br> <b>calc_radiation</b>: 503 500 Parametrisierung langwelliger 504 501 Strahlungsprozesse 505 über des Schema der effektiven Emissivität <br>502 über des Schema der effektiven Emissivität <br> 506 503 <b>calc_precipitation</b>:Parametrisierung von 507 504 Niederschlagsprozessen 508 505 mit Hilfe eines vereinfachten Kesslerschemas <br> <b>impact_of_latent_heat</b>: 509 Ber ücksichtigung der durch506 Berücksichtigung der durch 510 507 Niederschlagsprozesse 511 zur ückbleibenden latenten Waerme innerhalb eines508 zurückbleibenden latenten Waerme innerhalb eines 512 509 Gittervolumens <br> <b>compute_vpt</b>: 513 510 Berechnung der virtuellen potentiellen … … 529 526 der Niderschlagsparametrisierung</td> </tr> <tr> <td><tt>bc_q_b</tt></td> 530 527 <td>C</td> <td>'dirichlet'</td> <td>untere 531 Randbedingung f ür q (siehe auch bc_pt_b)</td> </tr>528 Randbedingung für q (siehe auch bc_pt_b)</td> </tr> 532 529 <tr> <td><tt>bc_q_t</tt></td> <td>C</td> 533 <td>'neumann'</td> <td>obere RB f ür q530 <td>'neumann'</td> <td>obere RB für q 534 531 (siehe auch bc_pt_t)</td> </tr> <tr> <td><tt>q_surface</tt></td> 535 532 <td>R</td> <td>0.0</td> <td>Feuchtewert 536 an der Erdoberfl äche</td> </tr> <tr> <td><tt>q_surface_initial_change</tt></td>533 an der Erdoberfläche</td> </tr> <tr> <td><tt>q_surface_initial_change</tt></td> 537 534 <td>R</td> <td>0.0</td> <td>vgl. 538 535 pt_surface_initial_change</td> </tr> <tr> <td><tt>q_vertical_gradient</tt></td> … … 543 540 pt_vertical_gradient_level</td> </tr> <tr> <td><tt>surface_pressure</tt></td> 544 541 <td>R</td> <td>1013.25</td> <td>Luftdruck 545 an der Erdoberfl äche</td> </tr> <tr> <td><tt>surface_waterflux</tt></td>546 <td>R</td> <td>0.0</td> <td>oberfl ächennaher542 an der Erdoberfläche</td> </tr> <tr> <td><tt>surface_waterflux</tt></td> 543 <td>R</td> <td>0.0</td> <td>oberflächennaher 547 544 Wasser/Feuchtefluss</td> </tr> </tbody> </table> 548 Fortsetzungsl äufe mit Daten, die mit Modellversionen kleiner545 Fortsetzungsläufe mit Daten, die mit Modellversionen kleiner 549 546 2.0 550 547 erzeugt … … 555 552 deshalb nach Abschluss eines kompletten Modelllaufs erfolgen. 556 553 </p> <p><b>Achtung:</b> diese 557 Änderungen haben Auswirkungen auf554 Änderungen haben Auswirkungen auf 558 555 den Prognoseverlauf. 559 Die Testergebnis-Datei wurde entsprechend ge ändert.</p> </td>556 Die Testergebnis-Datei wurde entsprechend geändert.</p> </td> 560 557 <td nosave="" valign="top">modules, parin, 561 558 read_var_list, … … 572 569 Umstellung auf Modellversion 2.0 wurden 573 570 noch kleine 574 Änderungen an den Unterprogrammen diffusion_e, production_e571 Änderungen an den Unterprogrammen diffusion_e, production_e 575 572 und 576 573 diffusivities 577 574 vorgenommen. In allen genannten <br> 578 575 Unterprogrammen wurde bisher ein strengeres Kriterium zur Bestimmung 579 der Schichtungsverh ältnisse verwendet. Ein Modellvergleich576 der Schichtungsverhältnisse verwendet. Ein Modellvergleich 580 577 zeigte 581 578 jedoch, dass dadurch die Diffusion im Bereich der Inversion zu gross 582 wird. <p>Eine weitere Erg änzung betrifft die579 wird. <p>Eine weitere Ergänzung betrifft die 583 580 Advektionsverfahren: <br> 584 Das Bott-Chlond Advektionsverfahren ist nun auch f ür581 Das Bott-Chlond Advektionsverfahren ist nun auch für 585 582 Rechnungen 586 mit Feuchte/Wolkenphysik verf ügbar. </p> <p><b>Achtung:</b>587 diese Änderungen haben Auswirkungen auf583 mit Feuchte/Wolkenphysik verfügbar. </p> <p><b>Achtung:</b> 584 diese Änderungen haben Auswirkungen auf 588 585 den Prognoseverlauf. 589 Die Testergebnis-Datei wurde entsprechend ge ändert.</p> </td>586 Die Testergebnis-Datei wurde entsprechend geändert.</p> </td> 590 587 <td nosave="" valign="top">diffusion_e, 591 588 production_e, 592 589 diffusivities, check_parameters</td> </tr> <tr nosave="" valign="top"> <td nosave="">26/04/00</td> 593 590 <td>SR</td> <td>2.0a</td> <td>C</td> 594 <td>Vollst ändige Umstellung von vtk-Grafik-Software auf591 <td>Vollständige Umstellung von vtk-Grafik-Software auf 595 592 dvrp-Software. 596 593 Grafik-Ausgaben im dvr-Format sind mittlerweile begrenzt nutzbar 597 (Isooberfl ächen594 (Isooberflächen 598 595 sind darstellbar, siehe neue Laufparameter <tt>dt_dvrp, 599 596 threshold</tt>), 600 597 es wird aber in der Anwendung in naher Zukunft noch diverse 601 Änderungen602 geben. Partikelausgabe vorerst nicht mehr m öglich (Verlagerung598 Änderungen 599 geben. Partikelausgabe vorerst nicht mehr möglich (Verlagerung 603 600 schon).</td> <td>header, init_particles, init_3d_model, 604 601 leap_frog, modules, … … 610 607 init_dvrp, plot_dvrp</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">26/04/00</td> 611 608 <td>SR</td> <td>2.0a</td> <td>C</td> 612 <td>Durch zus ätzlichen zyklischen Rand im609 <td>Durch zusätzlichen zyklischen Rand im 613 610 Bott-Chlond-Schema 614 611 sehr 615 612 viele SENDRECV-Aufrufe eingespart. <br> 616 Änderung des Namens einer eingelesenen Environment-Variablen613 Änderung des Namens einer eingelesenen Environment-Variablen 617 614 in <tt>check_open</tt> 618 615 von <tt>remote_addres</tt> nach <tt>return_addres</tt>. … … 621 618 Initialisierungsparameter <tt>npex</tt> und <tt>npey</tt>. 622 619 Durch 623 Änderungen624 am Header-Format wurde Testergebnis-Datei ge ändert.</td>620 Änderungen 621 am Header-Format wurde Testergebnis-Datei geändert.</td> 625 622 <td>advec_s_bc, check_open, header, init_pegrid, modules, parin, 626 623 read_var_list, 627 624 write_var_list</td> </tr> <tr nosave="" valign="top"> <td nosave="">26/04/00</td> 628 625 <td>SR</td> <td>2.0a</td> <td>N</td> 629 <td>Rechnungen mit geneigter Oberfl äche erlauben jetzt626 <td>Rechnungen mit geneigter Oberfläche erlauben jetzt 630 627 die 631 628 Vorgabe 632 629 eines stabil geschichteten Temperaturprofils mit konstantem Gradienten. 633 Hinzuschalten von Feuchte funktioniert f ür Hangrechnungen noch630 Hinzuschalten von Feuchte funktioniert für Hangrechnungen noch 634 631 nicht.</td> <td>advec_s_bc, buoyancy, init_3d_model, 635 632 modules, … … 637 634 init_slope</p> </td> </tr> <tr nosave="" valign="top"> <td>13/06/00</td> <td>MS</td> 638 635 <td>2.0a</td> <td nosave="">B</td> <td>Es 639 werden nun auch die seitlichen R änder für ql636 werden nun auch die seitlichen Ränder für ql 640 637 gesetzt. 641 Ein Nichtsetzen der Randbedingungen f ührte zu Fehlern bei der638 Ein Nichtsetzen der Randbedingungen führte zu Fehlern bei der 642 639 Ausgabe 643 640 von pt</td> <td nosave="">calc_liquid_water_content</td> 644 641 </tr> <tr nosave="" valign="top"> <td nosave="">03/07/00</td> <td>SR</td> <td>2.0b</td> 645 <td>C</td> <td>F ür eine Reihe von 642 <td>C</td> <td>Für eine Reihe von 646 643 Unterprogrammen, die innerhalb 647 644 von <tt>prognostic_equation</tt> 648 aufgerufen werden und an die Pointer-Variablen übergeben645 aufgerufen werden und an die Pointer-Variablen übergeben 649 646 werden, 650 647 sind 651 648 jetzt explizite Interfaces vereinbart (neues Modul <tt>pointer_interfaces</tt>). 652 649 Auch im Vereinbarungsteil der entsprechenden Unterprogramme sind diese 653 Übergabeparameter jetzt aus Konsistenzgründen (und650 Übergabeparameter jetzt aus Konsistenzgründen (und 654 651 weil es 655 652 sonst 656 653 Laufzeitfehler gibt) als Pointer-Variablen deklariert. Mit dieser 657 Änderung658 ist der gro ße Speicher- und CPU-Zeit-Bedarf beim659 Übersetzen660 von prognostic_equations wieder auf ein ertr ägliches661 Ma ß654 Änderung 655 ist der große Speicher- und CPU-Zeit-Bedarf beim 656 Übersetzen 657 von prognostic_equations wieder auf ein erträgliches 658 Maß 662 659 reduziert. <br> 663 Im Rahmen dieser Änderungen wurde auch an diffusion_e ein660 Im Rahmen dieser Änderungen wurde auch an diffusion_e ein 664 661 Hilfsfeld 665 weniger übergeben. <br>666 Aus nicht vollst ändig gekärten Gründen haben662 weniger übergeben. <br> 663 Aus nicht vollständig gekärten Gründen haben 667 664 sich die 668 Zahlen in der Testergebnis-Datei minimal ge ändert (und zwar665 Zahlen in der Testergebnis-Datei minimal geändert (und zwar 669 666 die 670 667 Spalte 671 mit der Gesamtdivergenz nach Aufruf des Druckl ösers).</td>668 mit der Gesamtdivergenz nach Aufruf des Drucklösers).</td> 672 669 <td>advec_s_bc, buoyancy, diffusion_e, diffusion_s, diffusion_u, 673 670 diffusion_v, 674 671 diffusion_w, disturb_field, modules, production_e, prognostic_equations</td> 675 672 </tr> <tr nosave="" valign="top"> <td nosave="">04/07/00</td> <td>SR</td> <td>2.0b</td> 676 <td>B</td> <td>Diriclet-Randbedingungen f ür673 <td>B</td> <td>Diriclet-Randbedingungen für 677 674 Temperatur und 678 Fl üssigwassergehalt679 werden gesetzt. Dies w äre bereits nach Einführung der675 Flüssigwassergehalt 676 werden gesetzt. Dies wäre bereits nach Einführung der 680 677 Pointer 681 678 zwingend notwendig gewesen. Rechnungen mit vorgegebener 682 Oberfl ächentemperatur679 Oberflächentemperatur 683 680 schlugen deshalb bislang fehl.</td> <td>boundary_conds</td> 684 681 </tr> <tr nosave="" valign="top"> <td nosave="">07/09/00</td> <td>MS</td> <td>2.0b</td> 685 682 <td>B</td> <td>Die virtuelle potenielle Temperatur 686 wird nun auch f ür683 wird nun auch für 687 684 k=nzb und 688 k=nzt+1 berechnet, um die Randwerte auch f ür Ausgabezwecke685 k=nzt+1 berechnet, um die Randwerte auch für Ausgabezwecke 689 686 bereichtzustellen. </td> <td>compute_vpt</td> 690 687 </tr> <tr nosave="" valign="top"> <td nosave="">28/12/00</td> <td>SR</td> <td>2.1</td> … … 692 689 optionaler Software-Pakete, 693 690 die zwar 694 zum Modell geh ören, standardmäßig aber691 zum Modell gehören, standardmäßig aber 695 692 nicht mit 696 übersetzt693 übersetzt 697 694 werden, um so unter anderem Kompilationszeit einzusparen. Siehe neues 698 695 Kapitel 699 696 3.7 in der Modelldokumentation. <p>Einsatz der 700 dvrp-Software wurde vollst ändig701 überarbeitet.697 dvrp-Software wurde vollständig 698 überarbeitet. 702 699 Sie ist nun als Software-Paket optional im Modell einsetzbar. Siehe 703 überarbeitetes700 überarbeitetes 704 701 Kapitel 4.5.6 in der Modelldokumentation. </p> <p>Die 705 702 Prognose der Partikelverlagerung ist ebenfalls in ein … … 719 716 <td>parles</td> </tr> <tr nosave="" valign="top"> <td nosave="">02/01/01</td> 720 717 <td>SR</td> <td>2.1a</td> <td>C</td> 721 <td>Offene Dateien werden sobald m öglich geschlossen.722 Wiederer öffnung723 bestimmter Dateien mit POSITION='APPEND' m öglich. In diesem718 <td>Offene Dateien werden sobald möglich geschlossen. 719 Wiedereröffnung 720 bestimmter Dateien mit POSITION='APPEND' möglich. In diesem 724 721 Zusammenhang 725 722 wurde das Unterprogramm close_files in close_file umbenannt und hat nun … … 745 742 <br> 746 743 calc_spectra, plot_spectra</p> </td> </tr> <tr valign="top"> <td>25/01/01</td> <td>SR</td> 747 <td>2.1c</td> <td>C/N</td> <td>Druckl öser744 <td>2.1c</td> <td>C/N</td> <td>Drucklöser 748 745 kann jetzt die Singleton-FFT benutzen. Modul <tt>singleton</tt> 749 746 wurde deshalb von <tt>calc_spectra.f90</tt> nach <tt>modules.f90</tt> 750 747 geschoben 751 748 (neuer Initialisierungsparameter <tt>fft_method</tt>). 752 <br>Fl üsse können jetzt bei k=1753 unabh ängig von der749 <br>Flüsse können jetzt bei k=1 750 unabhängig von der 754 751 Verwendung 755 752 einer Prandtl-Schicht vorgegeben werden (neuer … … 760 757 umbenannt. <br> 761 758 Modell kann auf DEC-Workstations eingesetzt werden (-D <tt>dec</tt>). 762 <br>Zus ätzliche Zeitmessungen mit dvrp-Software.759 <br>Zusätzliche Zeitmessungen mit dvrp-Software. 763 760 Schreiben von 764 761 Partikel-Informationen … … 773 770 parles, module_test</p> </td> </tr> <tr valign="top"> <td>25/01/01</td> <td>SR</td> 774 771 <td>2.1c</td> <td>B</td> <td>Fehler 775 beim öffnen von Unit 23 beseitigt. Modul <tt>interface</tt>772 beim öffnen von Unit 23 beseitigt. Modul <tt>interface</tt> 776 773 fehlte in disturb_field. <br> 777 Wertebereich von theta* in Prandtl-Fluxes eingeschr änkt, weil774 Wertebereich von theta* in Prandtl-Fluxes eingeschränkt, weil 778 775 sonst auf einigen Rechnern overflow droht, wenn u*=0.</td> <td>check_open, 779 776 disturb_field, prandtl_fluxes</td> </tr> <tr valign="top"> <td>30/01/01</td> <td>SR</td> 780 777 <td>2.1d</td> <td>N/C</td> <td>Prognose 781 eines passiven Skalars ist m öglich (anstatt778 eines passiven Skalars ist möglich (anstatt 782 779 Feuchte), 783 780 dazu neuer Initialisierungsparameter <tt>passive_scalar</tt>. 784 <p>Falls kein Pfad f ür ftpcopy existiert, werden in <tt>check_open</tt>785 die Dateinamen f ür AVS-Koordinaten- und Datendatei auf781 <p>Falls kein Pfad für ftpcopy existiert, werden in <tt>check_open</tt> 782 die Dateinamen für AVS-Koordinaten- und Datendatei auf 786 783 "unknown" 787 784 gesetzt.</p> </td> <td>asselin_filter, … … 794 791 write_var_list</td> </tr> <tr valign="top"> <td>30/01/01</td> 795 792 <td>SR</td> <td>2.1d</td> <td>B</td> 796 <td>String-Ausgabe f ür use_surface_fluxes korrigiert.793 <td>String-Ausgabe für use_surface_fluxes korrigiert. 797 794 Version vom 798 25/01/01 brach deswegen bei Fortsetzungsl äufen ab. <br>795 25/01/01 brach deswegen bei Fortsetzungsläufen ab. <br> 799 796 xy-Plotausgabe von q und ql korrigiert.</td> <td>plot_2d, 800 797 write_var_list</td> </tr> <tr nosave="" valign="top"> <td nosave="">30/03/01</td> 801 798 <td>SR</td> <td>2.2</td> <td>N/C</td> 802 <td>S ämtliche 3D-Felder der Zeitebene t+dt wurden799 <td>Sämtliche 3D-Felder der Zeitebene t+dt wurden 803 800 entfernt. 804 801 Ebenso 805 802 die Arbeitsfelder work und work1. Hilfsfelder wurden soweit 806 m öglich807 eingespart und werden grunds ätzlich lokal allokiert. Dadurch803 möglich 804 eingespart und werden grundsätzlich lokal allokiert. Dadurch 808 805 konnten 809 eine Reihe von Übergabeparametern gestrichen werden (z.B. beim806 eine Reihe von Übergabeparametern gestrichen werden (z.B. beim 810 807 Aufruf 811 808 von <tt>pres, advec_s_bc, poisfft, production_e, sor </tt>sowie … … 815 812 der 816 813 Zeitfilter jetzt direkt in der prognostischen Gleichung angewendet wird 817 (bisher geschah dies erst nach Aufruf des Druckl ösers), haben814 (bisher geschah dies erst nach Aufruf des Drucklösers), haben 818 815 sich 819 die Zahlen in der Testergebnis-Datei ge ändert. Die mittleren816 die Zahlen in der Testergebnis-Datei geändert. Die mittleren 820 817 Profile 821 818 sind aber unbeeinflusst. <p>Eine weitere kleinere 822 Änderung der Testergebnis-Daten819 Änderung der Testergebnis-Daten 823 820 wird dadurch 824 821 verursacht, dass jetzt der Parameter <tt>adjust_mixing_length</tt> 825 defaultm äßig <tt>.FALSE.</tt>826 ist und der Mischungsweg grunds ätzlich zusätzlich822 defaultmäßig <tt>.FALSE.</tt> 823 ist und der Mischungsweg grundsätzlich zusätzlich 827 824 auf <tt>0.7 828 825 * zu</tt> begrenzt wird. Falls <tt>adjust_mixing_length=T</tt>, … … 830 827 analog zu <tt>diffusion_e</tt> 831 828 an allen Gitterpunkten modifiziert. <b>Achtung: </b>Diese 832 Änderungen833 k önnen bei anderen Simulationen t.w. erhebliche Auswirkungen829 Änderungen 830 können bei anderen Simulationen t.w. erhebliche Auswirkungen 834 831 auf 835 832 die 836 oberfl ächennahen Ergebnisse haben. </p> <p><b>Kleinere837 Änderungen:</b> In <tt>init_cloud_physics</tt>833 oberflächennahen Ergebnisse haben. </p> <p><b>Kleinere 834 Änderungen:</b> In <tt>init_cloud_physics</tt> 838 835 wird 839 836 der Wert der Variablen <tt>surface_pressure</tt> in hPa … … 842 839 Umwandlung in Pa). <br> 843 840 Die Namelist-Namen der Software-Pakete sind t.w. 844 ge ändert. <br>845 Unit 14 (bin äre Ausgabe für846 Fortsetzungsl äufe) wird841 geändert. <br> 842 Unit 14 (binäre Ausgabe für 843 Fortsetzungsläufe) wird 847 844 jetzt im Hauptprogramm geschlossen, damit auf diese Datei noch 848 845 benutzergesteuerte 849 Ausgaben erfolgen k önnen. <br>846 Ausgaben erfolgen können. <br> 850 847 Laufparameter werden in keinem Fall mehr auf Unit 14 ausgegeben bzw. 851 von Unit 13 gelesen, d.h. sie gelten jetzt tats ächlich nur852 f ür848 von Unit 13 gelesen, d.h. sie gelten jetzt tatsächlich nur 849 für 853 850 den jeweils aktuellen Lauf. <br> 854 851 Partikeladvektion funktioniert jetzt auch zusammen mit der … … 856 853 deutsche Variablen-, Unterprogramm- und Modulnamen 857 854 wurden 858 ins Englische übersetzt. Davon sind fast alle Programmteile855 ins Englische übersetzt. Davon sind fast alle Programmteile 859 856 betroffen 860 (nicht in rechter Spalte aufgef ührt).</p> </td> <td>advec_particles,857 (nicht in rechter Spalte aufgeführt).</p> </td> <td>advec_particles, 861 858 advec_s_bc, advec_s_ups, advec_u_ups, 862 859 advec_v_ups, … … 874 871 </tr> <tr nosave="" valign="top"> <td nosave="">12/07/01</td> <td>SR</td> <td>2.2a</td> 875 872 <td>N/C</td> <td>Defaultwert der unteren 876 Randbedingung f ür die TKE ist ab873 Randbedingung für die TKE ist ab 877 874 sofort <tt>bc_e_b 878 875 = 'neumann'</tt>. <br> 879 Partikelquelle ist über Paketparameter steuerbar.876 Partikelquelle ist über Paketparameter steuerbar. 880 877 Partikeleigenschaften 881 878 sind mittels benutzereigener Software steuerbar. Der dvrp-Ausgabeweg 882 879 ist 883 jetzt ebenfalls über Paketparameter steuerbar. <br>880 jetzt ebenfalls über Paketparameter steuerbar. <br> 884 881 Allen Modulen (bis auf singleton) wurde eine SAVE-Anweisung 885 hinzugef ügt,886 damit die durch sie vereinbarten Variablen w ährend der882 hinzugefügt, 883 damit die durch sie vereinbarten Variablen während der 887 884 Rechnung in 888 keinem Fall undefiniert werden k önnen (diese Gefahr besteht885 keinem Fall undefiniert werden können (diese Gefahr besteht 889 886 z.B. 890 887 auf … … 894 891 plot_dvrp, user_interface</td> </tr> <tr> <td>12/07/01</td> 895 892 <td>SR</td> <td>2.2a</td> <td>B</td> 896 <td>Verschiebung der OPEN-Anweisung f ür Unit 33893 <td>Verschiebung der OPEN-Anweisung für Unit 33 897 894 (AVS-FLD-Datei).</td> <td>check_open</td> </tr> 898 895 <tr nosave="" valign="top"> <td nosave="">20/07/01</td> 899 896 <td>SR</td> <td>2.3</td> <td>N</td> 900 <td>Einbau des Mehrgitterverfahrens zur L ösung der897 <td>Einbau des Mehrgitterverfahrens zur Lösung der 901 898 Poisson-Gleichung.</td> <td>check_parameters, 902 899 exchange_horiz, header, init_grid, … … 905 902 poismg</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">21/08/01</td> 906 903 <td>SR</td> <td>2.3a</td> <td>N/C</td> 907 <td>Erg änzung der dvrp-Software: Partikel908 k önnen mit909 Schw änzen904 <td>Ergänzung der dvrp-Software: Partikel 905 können mit 906 Schwänzen 910 907 versehen werden. Die Ausgabe der dvrp-Plotdaten durch einen separaten 911 908 PE 912 ist m öglich (mrun-Option "<tt>-p dvrp_graphics+1PE</tt>").909 ist möglich (mrun-Option "<tt>-p dvrp_graphics+1PE</tt>"). 913 910 Der 914 911 Kommunikator <tt>MPI_COMM_WORLD</tt> … … 917 914 Wandbereich kann wahlweise 918 915 abgeschaltet 919 werden. F ür km und kh wird jetzt eine Neumann-Randbedingung am916 werden. Für km und kh wird jetzt eine Neumann-Randbedingung am 920 917 unteren 921 918 Rand verwendet. Die TKE-Energieproduktionsberechnung durch Scherung des 922 919 Grundstroms am unteren Rand wurde verbessert. </p> <p>Generelle 923 920 Verwendung der bodennahen vertikalen 924 Impulsfl üsse in921 Impulsflüsse in 925 922 den Diffusionstermen von u und v, immer wenn diese als Randbedingung 926 923 vorgegeben 927 sind (gem äß<tt>use_surface_fluxes</tt>,924 sind (gemäß <tt>use_surface_fluxes</tt>, 928 925 - bisher wurde 929 926 dies 930 927 nur bei eingeschalteter Prandtl-Schicht gemacht). Schubspannungen 931 k önnen928 können 932 929 jetzt als Randbedingungen vorgegeben werden (allerdings unter 933 930 zuhilfenahme 934 benutzereigener Software). </p> <p>Zus ätzliche935 Schnittstellen f ür benutzereigene931 benutzereigener Software). </p> <p>Zusätzliche 932 Schnittstellen für benutzereigene 936 933 Software in 937 934 allen Tendenztermen, in <tt>flow_statistics,</tt> und 938 f ür935 für 939 936 Partikel- 940 937 bzw. dvrp-Programmteile (Bestimmung von Partikeleigenschaften, 941 938 Festlegung 942 939 von Farbtabellen). </p> <p><b>Achtung: </b>Diese 943 Änderungen können bei940 Änderungen können bei 944 941 Simulationen mit 945 mittlerem Wind durch die Änderungen im bereich des unteren942 mittlerem Wind durch die Änderungen im bereich des unteren 946 943 Randes 947 t.w. erhebliche Auswirkungen auf die oberfl ächennahen944 t.w. erhebliche Auswirkungen auf die oberflächennahen 948 945 Ergebnisse 949 946 haben.</p> </td> <td>advec_particles, diffusion_e, … … 961 958 <tr nosave="" valign="top"> <td nosave="">04/09/01</td> 962 959 <td>SR</td> <td>2.3b</td> <td>N</td> 963 <td>Zus ätzliche Profilausgaben für964 Energieproduktionsterme m öglich.965 St ördruck wird bei der FFT-Methode nicht mehr durch960 <td>Zusätzliche Profilausgaben für 961 Energieproduktionsterme möglich. 962 Stördruck wird bei der FFT-Methode nicht mehr durch 966 963 Aufsummieren 967 964 der … … 969 966 flow_statistics, modules, pres</td> </tr> <tr nosave="" valign="top"> <td nosave="">04/09/01</td> 970 967 <td>SR</td> <td>2.3b</td> <td>B</td> 971 <td>Austausch der Geisterr änder für die968 <td>Austausch der Geisterränder für die 972 969 zeitgefilterten 973 970 Felder 974 direkt nach Durchf ührung der Filterung (war fehlerhaft seit971 direkt nach Durchführung der Filterung (war fehlerhaft seit 975 972 Entfernen 976 der Zeitebene t+dt (Version 2.2) und f ührte zu sehr kleinen977 St örungen978 an den seitlichen R ändern der Teilgebiete). <p><b>Achtung:</b>979 Der Inhalt der Testergebnis-Datei ändert973 der Zeitebene t+dt (Version 2.2) und führte zu sehr kleinen 974 Störungen 975 an den seitlichen Rändern der Teilgebiete). <p><b>Achtung:</b> 976 Der Inhalt der Testergebnis-Datei ändert 980 977 sich durch 981 978 diese Fehlerkorrektur.</p> </td> <td>prognostic_equations</td> 982 979 </tr> <tr nosave="" valign="top"> <td nosave="">09/11/01</td> <td>SR</td> <td>2.3c</td> 983 980 <td>N</td> <td>Farbe entlang der 984 Partikelschw änze kann sich985 ändern. Einschalten986 der Partikeladvektion bei Fortsetzungsl äufen möglich.</td>981 Partikelschwänze kann sich 982 ändern. Einschalten 983 der Partikeladvektion bei Fortsetzungsläufen möglich.</td> 987 984 <td>advec_particles, init_particles, modules, package_parin, 988 985 plot_dvrp</td> </tr> <tr nosave="" valign="top"> 989 986 <td nosave="">16/04/02</td> <td>SR</td> 990 987 <td>2.3d</td> <td>N</td> <td>Vorgabe 991 von Randbedingungen f ür Partikeladvektion992 m öglich.993 Partikeldaten k önnen für spätere Analysen988 von Randbedingungen für Partikeladvektion 989 möglich. 990 Partikeldaten können für spätere Analysen 994 991 auf Datei 995 992 geschrieben 996 993 werden (Unit 85). <tt>PARTICLE</tt>-Datentyp 997 enth ält994 enthält 998 995 Informationen 999 über Partikelgeschwindigkeitskomponenten und996 über Partikelgeschwindigkeitskomponenten und 1000 997 Startposition. <p>Skalartransport sowie entsprechende 1001 Datenausgaben k önnen998 Datenausgaben können 1002 999 durch 1003 eigene Parameter gesteuert werden (nicht mehr über die1000 eigene Parameter gesteuert werden (nicht mehr über die 1004 1001 Feuchteparameter). </p> <p>Im Fall von 1005 1002 pdims(1)=1 (eindimensionales virtuelles … … 1017 1014 batch_scp-Script 1018 1015 ersetzt (dient der Ermittelung von Dateizyklusnummern als Information 1019 f ür1016 für 1020 1017 die AVS-fld-Datei). </p> <p>Bei fehlerhafter 1021 Er öffnung einer CPU_MEASURES-Datei1018 Eröffnung einer CPU_MEASURES-Datei 1022 1019 werden erneute 1023 Er öffnungsversuche durchgeführt.</p> </td>1020 Eröffnungsversuche durchgeführt.</p> </td> 1024 1021 <td>advec_particles, check_open, check_parameters, 1025 1022 exchange_horiz, flow_statistics, … … 1030 1027 transpose_zy, write_var_list</td> </tr> <tr nosave="" valign="top"> <td nosave="">02/05/02</td> 1031 1028 <td>SR</td> <td>2.3e</td> <td>B</td> 1032 <td>Wiedereinf ührung der 3D-Felder für die1029 <td>Wiedereinführung der 3D-Felder für die 1033 1030 Zeitebene 1034 1031 t+dt, die … … 1036 1033 mit der korrekten Zeitebene t-dt gearbeitet wird. Seit Version 2.2 1037 1034 wurde 1038 f älschlicherweise die Zeitebene t+dt verwendet. Asselin-Filter1035 fälschlicherweise die Zeitebene t+dt verwendet. Asselin-Filter 1039 1036 ist 1040 jetzt wieder eigenst ändiges Unterprogramm. <br>1041 Die Wiedereinf ührung dieser dritten Zeitebene ist auch1037 jetzt wieder eigenständiges Unterprogramm. <br> 1038 Die Wiedereinführung dieser dritten Zeitebene ist auch 1042 1039 Voraussetzung 1043 f ür die in der nächsten Version geplante skalare1040 für die in der nächsten Version geplante skalare 1044 1041 Optimierung. <p>Fehler bei Berechnung von 1045 1042 Ausgabezeitpunkten eliminiert (trat 1046 1043 auf, wenn 1047 bei Fortsetzungsl äufen von Ausgabeabständen 0.0 auf1044 bei Fortsetzungsläufen von Ausgabeabständen 0.0 auf 1048 1045 von Null 1049 1046 verschiedene Werte gewechselt werden sollte). </p> <p><b>Achtung:</b> 1050 Der Inhalt der Testergebnis-Datei ändert1047 Der Inhalt der Testergebnis-Datei ändert 1051 1048 sich durch 1052 1049 diese Fehlerkorrektur.</p> </td> <td>init_3d_model, … … 1057 1054 Anpassungen an IBM-Regatta-Systeme in 1058 1055 check_parameters, cpu_log. <p>Modul singleton in 1059 eigenst ändige Datei1060 überführt. </p> <p>Global_min_max1056 eigenständige Datei 1057 überführt. </p> <p>Global_min_max 1061 1058 arbeitet mit REAL*4, um Komplikationen auf 1062 1059 32-bit-Rechnern 1063 zu vermeiden (Datentyp MPI_2REAL m üsste dort sonst auf1060 zu vermeiden (Datentyp MPI_2REAL müsste dort sonst auf 1064 1061 MPI_2DOUBLE_PRECISION 1065 ge ändert werden). <br>1066 Horizontale Geschwindigkeitskomponenten f ür Partikeladvektion1062 geändert werden). <br> 1063 Horizontale Geschwindigkeitskomponenten für Partikeladvektion 1067 1064 werden exakt zwischen den jeweils benachbarten horizontalen 1068 Gitterfl ächen1065 Gitterflächen 1069 1066 interpoliert.</p> </td> <td nosave="">advec_particles, 1070 1067 check_parameters, cpu_log, … … 1073 1070 singleton</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">11/06/02</td> 1074 1071 <td>SR</td> <td>2.4</td> <td>N/C</td> 1075 <td>Optimierung und Anpassungen f ür einzelne Knoten der1072 <td>Optimierung und Anpassungen für einzelne Knoten der 1076 1073 IBM-Regatta-Systeme. <p><b>Skalare (Cache) Optimierung:</b> 1077 1074 Tendenzterme werden 1078 1075 innerhalb einer 1079 gro ßen (i,j)-Schleife berechnet, die nun jeweils die gesamte1076 großen (i,j)-Schleife berechnet, die nun jeweils die gesamte 1080 1077 prognostische 1081 1078 Gleichung umfasst. Schleifenindices i und j werden an die 1082 1079 Tendenzunterprogramme 1083 als Argumente übergeben. Die Tendenzunterprogramme sind als1080 als Argumente übergeben. Die Tendenzunterprogramme sind als 1084 1081 Module 1085 geschrieben, die mittels Technik der überladenen Funktionen1082 geschrieben, die mittels Technik der überladenen Funktionen 1086 1083 auch 1087 1084 ohne 1088 diese Indices i und j aufgerufen werden k önnen, und dann wie1085 diese Indices i und j aufgerufen werden können, und dann wie 1089 1086 in 1090 fr üheren1087 früheren 1091 1088 Versionen funktionieren (d.h. in ihnen laufen die Schleifen 1092 über1089 über 1093 1090 alle 1094 1091 3 Dimensionen). Solche Teile der prognostischen Gleichungen, die 1095 1092 globale 1096 1093 Kommunikation beinhalten, sind vor die (i,j)-Schleife gezogen. Falls 1097 ausschlie ßlich1094 ausschließlich 1098 1095 mit Piascek-Williams-Advektion gearbeitet wird, werden alle 1099 1096 prognostischen 1100 1097 Gleichungen in einer einzigen (i,j)-Schleife gerechnet. </p> <p><b>Kommunikationsoptimierung:</b> 1101 Zus ätzliche1098 Zusätzliche 1102 1099 Implementierung einer 1103 1100 1D-Gebietszerlegung in x-Richtung. Mit dem neuen direkten 1104 Druckl öser1101 Drucklöser 1105 1102 (poisfft_hybrid) werden dadurch 4 von 6 notwendigen Transpositionen 1106 1103 eingespart. 1107 1104 Die 1D-Zerlegung wird automatisch bei Verwendung dieses neuen 1108 L ösers1109 aktiviert. Austausch von Geisterr ändern in y-Richtung wird bei1105 Lösers 1106 aktiviert. Austausch von Geisterrändern in y-Richtung wird bei 1110 1107 dieser 1111 1108 1D-Zerlegung durch einfache zyklische Randbedingungen ersetzt. </p> … … 1114 1111 IBM-Compiler 1115 1112 in der Regel den Einsatz der TRIM-Funktion, um 1116 überschüssige1113 überschüssige 1117 1114 Blanks am Ende des Strings abzuschneiden (sonst liefert der Vergleich 1118 1115 .F.). <br> 1119 1116 Lokal allokierte Felder t.w. in sogenannte automatische Felder 1120 1117 umgewandelt 1121 (g ünstiger, um Gefahr von eventuellen Speicherlecks zu1118 (günstiger, um Gefahr von eventuellen Speicherlecks zu 1122 1119 vermeiden). 1123 1120 Einige Feldoperationen wurden wegen schlechter Performance auf der IBM 1124 1121 entfernt (s. pres). </p> <p><b>Achtung: </b>Die 1125 Übersetzung des Modells erfordert1122 Übersetzung des Modells erfordert 1126 1123 aufgrund 1127 1124 der Vielzahl nunmehr eingesetzter Module die Verwendung des 1128 1125 make-Mechanismus, 1129 um die korrekte Reihenfolge bei der Übersetzung der1126 um die korrekte Reihenfolge bei der Übersetzung der 1130 1127 Programmteile 1131 zu gew ährleisten. Die Beschreibung dieser1132 Abh ängigkeiten1128 zu gewährleisten. Die Beschreibung dieser 1129 Abhängigkeiten 1133 1130 erfolgt 1134 1131 in der Datei<tt> Makefile</tt>, die zusammen mit den … … 1139 1136 ebenfalls 1140 1137 als Modul geschrieben und bedient sich der Methode der 1141 überladenen1138 überladenen 1142 1139 Funktionen. Existierende benutzereigene Software muss entsprechend 1143 1140 angepasst … … 1159 1156 </tr> <tr nosave="" valign="top"> <td nosave="">11/06/02</td> <td>SR</td> <td>2.4</td> 1160 1157 <td>B</td> <td>Fehlende Variablenvereinbarung in <tt>advec_particles</tt> 1161 erg änzt.</td> <td>advec_particles</td> </tr>1158 ergänzt.</td> <td>advec_particles</td> </tr> 1162 1159 <tr> <td valign="top">12/09/02</td> <td valign="top">SR</td> <td valign="top">2.4a</td> 1163 1160 <td valign="top">N</td> <td valign="top">Partikel 1164 k önnen mit Trägheit versehen1161 können mit Trägheit versehen 1165 1162 werden und 1166 1163 einen Dichteunterschied zum umgebenden Fluid haben. Zur Steuerung dient … … 1173 1170 und q entfernt 1174 1171 (obere Feldgrenze von 10 konnte in <tt>pt_vertical_gradient_level_ind</tt> 1175 überschritten werden). Fehler führte unter1176 Umst änden zu1177 "segmentation fault" bei Verwendung sehr gro ßer1172 überschritten werden). Fehler führte unter 1173 Umständen zu 1174 "segmentation fault" bei Verwendung sehr großer 1178 1175 Gitterpunktzahlen 1179 1176 in z-Richtung. <p>Fehler in Berechnung der v-Komponente am 1180 1177 unteren Rand 1181 1178 entfernt, die 1182 f ür den Scherungsproduktionsterm der TKE verwendet wird1179 für den Scherungsproduktionsterm der TKE verwendet wird 1183 1180 (Variable 1184 v_0). F älschlicherweise wurde zur Berechnung u(k=1) statt1181 v_0). Fälschlicherweise wurde zur Berechnung u(k=1) statt 1185 1182 v(k=1) 1186 1183 verwendet.</p> </td> <td valign="top">check_parameters, … … 2081 2078 </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> 2082 2079 </td></tr><tr><td align="undefined" valign="undefined">dt_plsp</td><td align="undefined" valign="undefined">dt_dosp</td><td align="undefined" valign="undefined"><br> 2083 </td></tr><tr><td>dt_plts</td><td>dt_dots</td><td></td></tr><tr><td align="undefined" valign="undefined">dt_pr1d</td><td align="undefined" valign="undefined">dt_dopr_listing</td><td align="undefined" valign="undefined"><br> 2080 </td></tr><tr><td>dt_plts</td><td>dt_dots</td><td><br> 2081 </td></tr><tr><td align="undefined" valign="undefined">dt_pr1d</td><td align="undefined" valign="undefined">dt_dopr_listing</td><td align="undefined" valign="undefined"><br> 2084 2082 </td></tr><tr><td align="undefined" valign="undefined">dt_pl2d_xy</td><td align="undefined" valign="undefined">dt_do2d_xy</td><td align="undefined" valign="undefined"><br> 2085 2083 </td></tr><tr><td align="undefined" valign="undefined">dt_pl2d_xz</td><td align="undefined" valign="undefined">dt_do2d_xz</td><td align="undefined" valign="undefined"><br> … … 2097 2095 </td></tr><tr><td align="undefined" valign="undefined">pl3d_compress</td><td align="undefined" valign="undefined">do3d_compress</td><td align="undefined" valign="undefined"><br> 2098 2096 </td></tr><tr><td align="undefined" valign="undefined">pl3d_precision</td><td align="undefined" valign="undefined">do3d_comp_prec</td><td align="undefined" valign="undefined"><br> 2099 </td></tr><tr><td>pl_spectra</td><td>data_output_sp</td><td></td></tr><tr><td align="undefined" valign="undefined">read_fields_from_prior_run</td><td align="undefined" valign="undefined">read_restart_data</td><td align="undefined" valign="undefined"><br> 2097 </td></tr><tr><td>pl_spectra</td><td>data_output_sp</td><td><br> 2098 </td></tr><tr><td align="undefined" valign="undefined">read_fields_from_prior_run</td><td align="undefined" valign="undefined">read_restart_data</td><td align="undefined" valign="undefined"><br> 2100 2099 </td></tr><tr><td align="undefined" valign="undefined">write_plot2d_on_each_pe</td><td align="undefined" valign="undefined">data_output_2d_on_each_pe</td><td align="undefined" valign="undefined"><br> 2101 2100 </td></tr><tr><td align="undefined" valign="undefined">z_max_pl1d</td><td align="undefined" valign="undefined">z_max_do1d</td><td align="undefined" valign="undefined"><br> … … 2152 2151 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 2153 2152 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, 2154 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-2153 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- 2155 2154 and 3d-datasets now contain the exact coordinates of the variables 2156 2155 along x and y with respect to the positions where they are defined on … … 2465 2464 user_init_grid, user_init_particles, user_init_plant_canopy, 2466 2465 user_last_actions, user_parin, user_parin, user_particle_attributes, 2467 user_read_restart_data, user_spectra, user_statistics</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;"> Bugfix: array <span style="font-family: Courier New,Courier,monospace;">d</span> is reallocated in case that multigrid is used. (calc_spectra)<br><br>Bugfixes for nonparallel execution. (check_for_restart, cpu_statistics, inflow_turbulence, timestep)<br><br>Size of <span style="font-family: Courier New,Courier,monospace;">pf3d </span>adjusted to the required output size (1 gridpoint less, along all three dimensions), because output of a subset of the data (<span style="font-family: Courier New,Courier,monospace;">pf3d(nxa:nxe...)</span>) in the NF90_PUT_VAR statement caused segmentation fault with the INTEL compiler. (combine_plot_fields)<br><br>Bugfix: error in <span style="font-family: Courier New,Courier,monospace;">zu</span> index in case of <span style="font-family: Courier New,Courier,monospace;">section_xy=-1</span>. (header)<br><br>Inconsistency2466 user_read_restart_data, user_spectra, user_statistics</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;">Bugfix: array <span style="font-family: Courier New,Courier,monospace;">d</span> is reallocated in case that multigrid is used. (calc_spectra)<br><br>Bugfixes for nonparallel execution. (check_for_restart, cpu_statistics, inflow_turbulence, timestep)<br><br>Size of <span style="font-family: Courier New,Courier,monospace;">pf3d </span>adjusted to the required output size (1 gridpoint less, along all three dimensions), because output of a subset of the data (<span style="font-family: Courier New,Courier,monospace;">pf3d(nxa:nxe...)</span>) in the NF90_PUT_VAR statement caused segmentation fault with the INTEL compiler. (combine_plot_fields)<br><br>Bugfix: error in <span style="font-family: Courier New,Courier,monospace;">zu</span> index in case of <span style="font-family: Courier New,Courier,monospace;">section_xy=-1</span>. (header)<br><br>Inconsistency 2468 2467 removed: as the thermal stratification is not taken into account for 2469 2468 the evaluation of the wall fluxes at vertical walls, the eddy viscosity … … 2471 2470 calc_spectra, check_for_restart, cpu_statistics, data_output_2d, 2472 2471 header, inflow_turbulence, init_cloud_physics, production_e, 2473 read_3d_binary, timestep<br><br>combine_plot_fields</td></tr> 2472 read_3d_binary, timestep<br><br>combine_plot_fields</td></tr><tr> 2473 <td style="vertical-align: top;">23/09/09<br> 2474 </td> 2475 <td style="vertical-align: top;">SR<br> 2476 </td> 2477 <td style="vertical-align: top;">3.7 (r392)<br> 2478 </td> 2479 <td style="vertical-align: top;">N<br> 2480 </td> 2481 <td style="vertical-align: top;"><span style="font-weight: bold;">Most important changes:</span> 2482 atmosphere-ocean coupling has been improved and steering of 2483 dvr-software is more user-friendly now. Beside that, the release 2484 contains a lot of optimizations, changes, and bugfixes which have 2485 accumulated over the last months.<br> 2486 <br> 2487 The atmosphere-ocean coupling allows to use independent precursor runs 2488 in order to account for different spin-up times. The time when coupling 2489 has to be started is given by new inipar parameter <span style="font-family: Courier New,Courier,monospace;">coupling_start_time</span>. The precursor ocean run has to be started using new <span style="font-weight: bold;">mrun</span> option "<span style="font-family: Courier New,Courier,monospace;">-y</span>" in order to add the appendix "<span style="font-family: Courier New,Courier,monospace;">_O</span>" to all output files. (<span style="font-family: Courier New,Courier,monospace;">check_for_restart</span>, <span style="font-family: Courier New,Courier,monospace;">check_parameters</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_3d</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_profiles</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_ptseries</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_spectra</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_tseries</span>, <span style="font-family: Courier New,Courier,monospace;">header</span>, <span style="font-family: Courier New,Courier,monospace;">init_coupling</span>, <span style="font-family: Courier New,Courier,monospace;">modules</span>, <span style="font-family: Courier New,Courier,monospace;">mrun</span>, <span style="font-family: Courier New,Courier,monospace;">parin</span>, <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>, <span style="font-family: Courier New,Courier,monospace;">surface_coupler</span>, <span style="font-family: Courier New,Courier,monospace;">time_integration</span>, <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br> 2490 2491 <br> 2492 Clipping of dvrp output implemented. Default colourtable for particles 2493 implemented, particle attributes (color, dvrp_size) can be set with new 2494 parameters <span style="font-family: Courier New,Courier,monospace;">particle_color</span>, <span style="font-family: Courier New,Courier,monospace;">particle_dvrpsize</span>, <span style="font-family: Courier New,Courier,monospace;">color_interval</span>, <span style="font-family: Courier New,Courier,monospace;">dvrpsize_interval</span>. Slicer attributes (dvrp) are set with new routine <span style="font-family: Courier New,Courier,monospace;">set_slicer_attributes_dvrp</span> and are controlled with existing parameters <span style="font-family: Courier New,Courier,monospace;">slicer_range_limits</span>. (<span style="font-family: Courier New,Courier,monospace;">init_dvrp</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_dvrp</span>, <span style="font-family: Courier New,Courier,monospace;">modules</span>, <span style="font-family: Courier New,Courier,monospace;">set_slicer_attributes_dvrp</span>, <span style="font-family: Courier New,Courier,monospace;">user_data_output_dvrp</span>)<br> 2495 2496 <br> 2497 Polygon reduction for topography and ground plate isosurface. Reduction 2498 level for buildings can be chosen with parameter cluster_size. (<span style="font-family: Courier New,Courier,monospace;">init_dvrp</span>)<br> 2499 2500 <br> 2501 Variables <span style="font-family: Courier New,Courier,monospace;">bc_lr</span> / <span style="font-family: Courier New,Courier,monospace;">bc_ns</span> in most subroutines replaced by LOGICAL variables <span style="font-family: Courier New,Courier,monospace;">bc_lr_cyc</span>,<br> 2502 <span style="font-family: Courier New,Courier,monospace;"> 2503 bc_ns_cyc</span> for speed optimization. This gives a significant speedup of more than 10%. (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>, <span style="font-family: Courier New,Courier,monospace;">diffusion_u</span>, <span style="font-family: Courier New,Courier,monospace;">diffusion_v</span>, <span style="font-family: Courier New,Courier,monospace;">diffusion_w</span>, <span style="font-family: Courier New,Courier,monospace;">modules</span>)<br> 2504 2505 <br> 2506 Scripts have been adapted for machine lck (Yonsei Univ.). (<span style="font-family: Courier New,Courier,monospace;">mrun</span>, <span style="font-family: Courier New,Courier,monospace;">mbuild</span>, <span style="font-family: Courier New,Courier,monospace;">subjob</span>)<br> 2507 <br> 2508 Additional timestep criterion in case of simulations with plant canopy. (<span style="font-family: Courier New,Courier,monospace;">timestep</span>)<br> 2509 <br> 2510 Check for illegal entries in <span style="font-family: Courier New,Courier,monospace;">section_xy</span>|<span style="font-family: Courier New,Courier,monospace;">xz</span>|<span style="font-family: Courier New,Courier,monospace;">yz</span> that exceed <span style="font-family: Courier New,Courier,monospace;">nz+1</span>|<span style="font-family: Courier New,Courier,monospace;">ny+1</span>|<span style="font-family: Courier New,Courier,monospace;">nx+1</span>. (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>)<br> 2511 <br> 2512 External pressure gradient can be used as driving force. (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>, <span style="font-family: Courier New,Courier,monospace;">header</span>, <span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>, <span style="font-family: Courier New,Courier,monospace;">modules</span>, <span style="font-family: Courier New,Courier,monospace;">parin</span>, <span style="font-family: Courier New,Courier,monospace;">prognostic_equations</span>, <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>, <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br> 2513 <br> 2514 New topography case '<span style="font-style: italic;">single_street_canyon</span>'. (<span style="font-family: Courier New,Courier,monospace;">header</span>, <span style="font-family: Courier New,Courier,monospace;">init_grid</span>, <span style="font-family: Courier New,Courier,monospace;">modules</span>, <span style="font-family: Courier New,Courier,monospace;">parin</span>, <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>, <span style="font-family: Courier New,Courier,monospace;">user_check_parameters</span>, <span style="font-family: Courier New,Courier,monospace;">user_header</span>, <span style="font-family: Courier New,Courier,monospace;">user_init_grid</span>, <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br> 2515 <br> 2516 Option to predefine a target bulk velocity for <span style="font-family: Courier New,Courier,monospace;">conserve_volume_flow</span>. (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>, <span style="font-family: Courier New,Courier,monospace;">header</span>, <span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>, <span style="font-family: Courier New,Courier,monospace;">modules</span>, <span style="font-family: Courier New,Courier,monospace;">parin</span>, <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>, <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br> 2517 <br> 2518 Option for user defined 2D data output in xy cross sections at z=<span style="font-family: Courier New,Courier,monospace;">nzb</span>+1. (<span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>, <span style="font-family: Courier New,Courier,monospace;">user_data_output_2d</span>)<br> 2519 <br> 2520 xy cross section output of surface heatfluxes (sensible and latent). (<span style="font-family: Courier New,Courier,monospace;">average_3d_data</span>, <span style="font-family: Courier New,Courier,monospace;">check_parameters</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>, <span style="font-family: Courier New,Courier,monospace;">modules</span>, <span style="font-family: Courier New,Courier,monospace;">read_3d_binary</span>, <span style="font-family: Courier New,Courier,monospace;">sum_up_3d_data</span>, <span style="font-family: Courier New,Courier,monospace;">write_3d_binary</span>)<br> 2521 </td> 2522 <td style="vertical-align: top;">average_3d_data, 2523 check_for_restart, check_parameters, data_output_2d, data_output_3d, 2524 data_output_dvrp, data_output_profiles, data_output_ptseries, 2525 data_output_spectra, data_output_tseries, diffusion_u, diffusion_v, 2526 diffusion_w, init_coupling, init_dvrp, init_grid, init_3d_model, 2527 header, mbuild, modules, mrun, package_parin, parin, 2528 prognostic_equations, read_3d_binary, read_var_list, subjob, 2529 sum_up_3d_data, surface_coupler, timestep, time_integration, 2530 user_check_parameters, user_data_output_2d, user_data_output_dvrp, 2531 user_header, user_init_grid, write_3d_binary, write_var_list<br> 2532 <br> 2533 <span style="font-weight: bold;">new:</span> set_particle_attributes, set_slicer_attributes_dvrp<br> 2534 <br> 2535 </td> 2536 </tr> 2537 <tr> 2538 <td style="vertical-align: top;"><br> 2539 </td> 2540 <td style="vertical-align: top;"><br> 2541 </td> 2542 <td style="vertical-align: top;"><br> 2543 </td> 2544 <td style="vertical-align: top;">C<br> 2545 </td> 2546 <td style="vertical-align: top;">Output of messages replaced by message handling routine. (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>, <span style="font-family: Courier New,Courier,monospace;">advec_s_bc</span>, <span style="font-family: Courier New,Courier,monospace;">buoyancy</span>, <span style="font-family: Courier New,Courier,monospace;">calc_spectra</span>, <span style="font-family: Courier New,Courier,monospace;">check_for_restart</span>, <span style="font-family: Courier New,Courier,monospace;">check_open</span>, <span style="font-family: Courier New,Courier,monospace;">coriolis</span>, <span style="font-family: Courier New,Courier,monospace;">cpu_log</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_3d</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_dvrp</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_profiles</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_spectra</span>, <span style="font-family: Courier New,Courier,monospace;">fft_xy</span>, <span style="font-family: Courier New,Courier,monospace;">flow_statistics</span>, <span style="font-family: Courier New,Courier,monospace;">header</span>, <span style="font-family: Courier New,Courier,monospace;">init_1d_model</span>, <span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>, <span style="font-family: Courier New,Courier,monospace;">init_dvrp</span>, <span style="font-family: Courier New,Courier,monospace;">init_grid</span>, <span style="font-family: Courier New,Courier,monospace;">init_particles</span>, <span style="font-family: Courier New,Courier,monospace;">init_pegrid</span>, <span style="font-family: Courier New,Courier,monospace;">netcdf</span>, <span style="font-family: Courier New,Courier,monospace;">parin</span>, <span style="font-family: Courier New,Courier,monospace;">plant_canopy_model</span>, <span style="font-family: Courier New,Courier,monospace;">poisfft_hybrid</span>, <span style="font-family: Courier New,Courier,monospace;">poismg</span>, <span style="font-family: Courier New,Courier,monospace;">read_3d_binary</span>, <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>, <span style="font-family: Courier New,Courier,monospace;">surface_coupler</span>, <span style="font-family: Courier New,Courier,monospace;">temperton_fft</span>, <span style="font-family: Courier New,Courier,monospace;">timestep</span>, <span style="font-family: Courier New,Courier,monospace;">user_actions</span>, <span style="font-family: Courier New,Courier,monospace;">user_data_output_dvrp</span>, <span style="font-family: Courier New,Courier,monospace;">user_dvrp_coltab</span>, <span style="font-family: Courier New,Courier,monospace;">user_init_grid</span>, <span style="font-family: Courier New,Courier,monospace;">user_init_plant_canopy</span>, <span style="font-family: Courier New,Courier,monospace;">user_parin</span>, <span style="font-family: Courier New,Courier,monospace;">user_read_restart_data</span>, <span style="font-family: Courier New,Courier,monospace;">user_spectra</span>)<br> 2547 2548 <br> 2549 Output of NetCDF messages with aid of message handling routine. (<span style="font-family: Courier New,Courier,monospace;">check_open</span>, <span style="font-family: Courier New,Courier,monospace;">close_file</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_3d</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_profiles</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_ptseries</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_spectra</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_tseries</span>, <span style="font-family: Courier New,Courier,monospace;">netcdf</span><span style="font-family: Courier New,Courier,monospace;"></span>)<br> 2550 2551 <br> 2552 2553 __lcmuk changed to __lc to avoid problems with Intel compiler on sgi-ice. (<span style="font-family: Courier New,Courier,monospace;">poisfft</span>)<br> 2554 <br> 2555 For extended NetCDF files, the updated title attribute includes an update of time_average_text where appropriate. (<span style="font-family: Courier New,Courier,monospace;">netcdf</span>)<br> 2556 <br> 2557 In case of restart runs without extension, initial profiles are not 2558 written to NetCDF-file anymore. So far, time axis always started at 2559 t=0. (<span style="font-family: Courier New,Courier,monospace;">data_output_profiles</span>, <span style="font-family: Courier New,Courier,monospace;">modules</span>, <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>, <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br> 2560 <br> 2561 <span style="font-family: Courier New,Courier,monospace;">initializing_actions</span>='<span style="font-style: italic;">read_data_for_recycling</span>' renamed to '<span style="font-style: italic;">cyclic_fill</span>'. A cyclic fill can now <br> 2562 be used independent of <span style="font-family: Courier New,Courier,monospace;">turbulent_inflow</span>. (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>, <span style="font-family: Courier New,Courier,monospace;">header</span>, <span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>)<br> 2563 <br> 2564 2 NetCDF error numbers changed. (<span style="font-family: Courier New,Courier,monospace;">data_output_3d</span>)<br> 2565 <br> 2566 A link to the webspage appendix_a.html is printed for further 2567 information about possible reasons of the error that appeared. (message)<br> 2568 <br> 2569 Temperature gradient criterion for estimating the boundary layer height 2570 replaced by the gradient criterion of Sullivan et al. (1998). (<span style="font-family: Courier New,Courier,monospace;">flow_statistics</span>)<br> 2571 <br> 2572 NetCDF unit attribute in timeseries output in case of statistic regions added. (<span style="font-family: Courier New,Courier,monospace;">netcdf</span>)<br> 2573 <br> 2574 Maximum number of tails is calculated from maximum number of particles and <span style="font-family: Courier New,Courier,monospace;">skip_particles_for_tail</span>. (<span style="font-family: Courier New,Courier,monospace;">init_particles</span>)<br> 2575 <br> 2576 Value of <span style="font-family: Courier New,Courier,monospace;">vertical_particle_advection</span> is allowed to differ for each particle group. (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>, <span style="font-family: Courier New,Courier,monospace;">header</span>, <span style="font-family: Courier New,Courier,monospace;">modules</span>) <br> 2577 <br> 2578 First constant in array <span style="font-family: Courier New,Courier,monospace;">den</span> also defined as type double. (<span style="font-family: Courier New,Courier,monospace;">eqn_state_seawater</span>)<br> 2579 <br> 2580 Parameter <span style="font-family: Courier New,Courier,monospace;">dvrp_psize</span> moved from <span style="font-family: Courier New,Courier,monospace;">particles_par</span> to <span style="font-family: Courier New,Courier,monospace;">dvrp_graphics_par</span>. (<span style="font-family: Courier New,Courier,monospace;">package_parin</span>)<br> 2581 <br> 2582 <span style="font-family: Courier New,Courier,monospace;">topography_grid_convention</span> moved from <span style="font-family: Courier New,Courier,monospace;">userpar</span> to <span style="font-family: Courier New,Courier,monospace;">inipar</span> (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>, <span style="font-family: Courier New,Courier,monospace;">header</span>, <span style="font-family: Courier New,Courier,monospace;">parin</span>, <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>, <span style="font-family: Courier New,Courier,monospace;">user_check_parameters</span>, <span style="font-family: Courier New,Courier,monospace;">user_header</span>, <span style="font-family: Courier New,Courier,monospace;">user_init_grid</span>, <span style="font-family: Courier New,Courier,monospace;">user_parin</span>, <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br> 2583 <br> 2584 Default value of <span style="font-family: Courier New,Courier,monospace;">grid_matching</span> changed to '<span style="font-style: italic;">strict</span>'. (<span style="font-family: Courier New,Courier,monospace;">modules</span>)<br> 2585 <br> 2586 Adjustments for runs on lcxt4 (necessary due to a software update on CRAY) and for coupled runs on ibmy. (<span style="font-family: Courier New,Courier,monospace;">mrun</span>, <span style="font-family: Courier New,Courier,monospace;">subjob</span>)<br> 2587 </td> 2588 <td style="vertical-align: top;">advec_particles, advec_s_bc, 2589 buoyancy, calc_spectra, check_for_restart, check_open, 2590 check_parameters, close_file, coriolis, cpu_log, data_output_2d, 2591 data_output_3d, data_output_dvrp, data_output_profiles, 2592 data_output_ptseries, data_output_spectra, data_output_tseries, 2593 eqn_state_seawater, fft_xy, flow_statistics, header, init_1d_model, 2594 init_3d_model, init_dvrp, init_grid, init_particles, init_pegrid, 2595 message, mrun, netcdf, package_parin, parin, plant_canopy_model, 2596 poisfft, poisfft_hybrid, poismg, read_3d_binary, read_var_list, subjob, 2597 surface_coupler, temperton_fft, timestep, user_actions,user_check_parameters, user_data_output_dvrp, user_dvrp_coltab, user_header, user_init_grid, user_init_plant_canopy, user_parin, user_read_restart_data, user_spectra, write_var_list <br> 2598 </td> 2599 </tr> 2600 <tr> 2601 <td style="vertical-align: top;"><br> 2602 </td> 2603 <td style="vertical-align: top;"><br> 2604 </td> 2605 <td style="vertical-align: top;"><br> 2606 </td> 2607 <td style="vertical-align: top;">E<br> 2608 </td> 2609 <td style="vertical-align: top;">Bugfix: Initial hydrostatic pressure profile in case of ocean runs is now calculated in 5 iteration steps. (<span style="font-family: Courier New,Courier,monospace;">init_ocean</span>)<br> 2610 <br> 2611 Bugfix: wrong sign in buoyancy production of ocean part in case of not using the reference density (only in 3D routine <span style="font-family: Courier New,Courier,monospace;">production_e</span>). (<span style="font-family: Courier New,Courier,monospace;">production_e</span>)<br> 2612 <br> 2613 Bugfix: output of averaged 2d/3d quantities requires that an avaraging 2614 interval has been set, respective error message is included. (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>)<br> 2615 <br> 2616 Bugfix: Output on unit 14 only if requested by <span style="font-family: Courier New,Courier,monospace;">write_binary</span>. (<span style="font-family: Courier New,Courier,monospace;">user_last_actions</span>)<br> 2617 <br> 2618 Bugfix to avoid zero division by <span style="font-family: Courier New,Courier,monospace;">km_neutral</span>. (<span style="font-family: Courier New,Courier,monospace;">production_e</span>)<br> 2619 <br> 2620 Bugfix for extended NetCDF files: In order to avoid 'data mode' errors 2621 if updated attributes are larger than their original size, <span style="font-family: Courier New,Courier,monospace;">NF90_PUT_ATT</span> is called in 'define mode' enclosed by <span style="font-family: Courier New,Courier,monospace;">NF90_REDEF</span> and <span style="font-family: Courier New,Courier,monospace;">NF90_ENDDEF</span> 2622 calls. This implies a possible performance loss; an alternative 2623 strategy would be to ensure equal attribute size in a job chain. (<span style="font-family: Courier New,Courier,monospace;">netcdf</span>)<br> 2624 <br> 2625 Bugfix: correction of initial volume flow for non-flat topography. (<span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>)<br> 2626 <br> 2627 Bugfix: zero initialization of arrays within buildings for '<span style="font-style: italic;">cyclic_fill</span>'. (<span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>)<br> 2628 <br> 2629 Bugfix: <span style="font-family: Courier New,Courier,monospace;">to_be_resorted</span> => <span style="font-family: Courier New,Courier,monospace;">s_av</span> for time-averaged scalars. (<span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>, <span style="font-family: Courier New,Courier,monospace;">data_output_3d</span>)<br> 2630 <br> 2631 Bugfix: avoid that <span style="font-family: Courier New,Courier,monospace;">ngp_2dh_s_inner</span> becomes zero. (<span style="font-family: Courier New,Courier,monospace;">init_3_model</span>)<br> 2632 <br> 2633 Typographical error: unit of wpt in dots_unit. (<span style="font-family: Courier New,Courier,monospace;">modules</span>)<br> 2634 <br> 2635 Bugfix: error in check, if particles moved further than one subdomain 2636 length. This check must not be applied for newly released particles. (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>)<br> 2637 <br> 2638 Bugfix: several tail counters are initialized, <span style="font-family: Courier New,Courier,monospace;">particle_tail_coordinates</span> 2639 is only written to file if its third index is > 0, arrays for tails 2640 are allocated with a minimum size of 10 tails if there is no tail 2641 initially. (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>, <span style="font-family: Courier New,Courier,monospace;">init_particles</span><span style="font-family: Courier New,Courier,monospace;"></span>)<br> 2642 <br> 2643 Bugfix: pressure included for profile output. (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>)<br> 2644 <br> 2645 Bugfix: Type of <span style="font-family: Courier New,Courier,monospace;">count</span> and <span style="font-family: Courier New,Courier,monospace;">count_rate</span> changed to default <span style="font-family: Courier New,Courier,monospace;">INTEGER</span> on NEC machines. (<span style="font-family: Courier New,Courier,monospace;">cpu_log</span>)<br> 2646 <br> 2647 Bugfix: output of particle time series only if particle advection is switched on. (<span style="font-family: Courier New,Courier,monospace;">time_integration</span>)<br> 2648 <br> 2649 Bugfix: qsws was calculated in case of <span style="font-family: Courier New,Courier,monospace;">constant heatflux</span> = <span style="font-family: Courier New,Courier,monospace;">.FALSE.</span>. (<span style="font-family: Courier New,Courier,monospace;">prandtl_fluxes</span>)<br> 2650 <br> 2651 Bugfix: averaging along z is not allowed for 2d quantities (e.g. u* and z0). (<span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>)<br> 2652 <br> 2653 Typographical errors. (<span style="font-family: Courier New,Courier,monospace;">netcdf</span>)<br> 2654 <br> 2655 If the inversion height calculated by the prerun is zero, inflow_damping_height must be explicitly specified. (<span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>)<br> 2656 <br> 2657 Small bugfix concerning 3d 64bit netcdf output format. (<span style="font-family: Courier New,Courier,monospace;">header</span>)<br> 2658 <br> 2659 Bugfix: <span style="font-family: Courier New,Courier,monospace;">dt_fixed</span> 2660 removed from the restart file, because otherwise, no change from a 2661 fixed to a variable timestep would be possible in restart runs. (<span style="font-family: Courier New,Courier,monospace;">read_var_list</span>, <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br> 2662 <br> 2663 Bugfix: initial setting of <span style="font-family: Courier New,Courier,monospace;">time_coupling</span> in coupled restart runs. (<span style="font-family: Courier New,Courier,monospace;">time_integration</span>)<br> 2664 </td> 2665 <td style="vertical-align: top;">advec_particles, 2666 check_parameters, cpu_log, data_output_2d, data_output_3d, header, 2667 init_3d_model, init_particles, init_ocean, modules, netcdf, 2668 prandtl_fluxes, production_e, read_var_list, time_integration, 2669 user_last_actions, write_var_list<br> 2670 <br> 2671 </td> 2672 </tr> 2673 2474 2674 </tbody> 2475 2675 </table> <b><blink>Attention:</blink></b> … … 2484 2684 <h2><a name="Kapitel2.0"></a>2.0 How 2485 2685 to change the model source 2486 code</h2> 2487 Änderungen, Erweiterungen oder Fehlerkorrekturen am Modellcode 2488 (gemeint 2489 sind alle Modellteile, die nicht zur benutzereigenen Software 2490 gehören) 2491 erfordern vom jeweiligen Benutzer Gruppenzugehörigkeit zur 2492 Gruppe 2493 <tt>"palm"</tt>. 2494 Zu dieser Gruppe gehören alle Benutzer, die aktiv an der 2495 Modellentwicklung 2496 mitwirken (im weiteren "Modellarbeitsgruppe" genannt). 2497 <p>Die einzelnen Quelltextdateien des Modells werden mit RCS 2498 (Revision 2499 Control System) verwaltet. Grundlegende Kenntnisse über dieses 2500 System 2501 werden im weiteren vorausgesetzt. 2502 </p><p>Änderungen am Modellcode erfordern die 2503 Durchführung 2504 bestimmter 2505 Aktionen in einer festgelegten Reihenfolge, die nun beschrieben werden 2506 sollen. 2507 </p><ol> <li>Konzeptplanung.</li> <li>Diskussion 2508 der geplanten Änderungen in der 2509 Modellarbeitsgruppe.</li> <li>Auschecken der 2510 benötigten Quelltextdateien, 2511 Durchführung der 2512 Änderungen und Test des Modells. Die Modul-Datei <tt>modules.f90</tt> 2513 sollte nur dann ausgecheckt werden, wenn dort aufgeführte 2514 Variablen 2515 verändert werden. Neu geschaffene Variablen sollten 2516 während 2517 der 2518 Testphase in ein ebenfalls temporär neu zu schaffendes Modul 2519 geschrieben 2520 werden, das von den jeweiligen Programmteilen benutzt werden muss 2521 (Einfügen 2522 entsprechender USE-Anweisungen). Die in den einzelnen Dateien 2523 durchgeführten 2524 Arbeiten müssen in den Kopfzeilen (unter "Aktuelle 2525 Aenderungen") 2526 der 2527 Dateien vermerkt werden. Dies tut man am besten noch <b>bevor</b> 2528 man 2529 die 2530 eigentlichen Änderungen durchführt (sonst droht das 2531 Vergessen...).</li> <li>Zum Abschluß der Tests 2532 sollte mit dem Modell eine 2533 Kontrollrechnung 2534 durchgeführt und die Ergebnisse mit denen des sogenannten 2535 Standardlaufes 2536 verglichen werden. Diese Rechnung ist mit der <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/INSTALL/.mrun.config">Konfigurationsdatei</a> 2537 und der <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/INSTALL/example_p3d">Parameterdatei</a> 2538 durchzuführen, die auch zum Testen der Modellinstallation 2539 verwendet 2540 werden (s. <a href="Anwendung/Kapitel_5.0.html">Kapitel 2541 5.0</a> der 2542 Anwenderdokumentation). 2543 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> 2544 ist dann mit der <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/INSTALL/example_rc">Beispiel-Ergebnisdatei</a> 2545 zu vergleichen (Hinweis: Kopien der Konfigurationsdatei, Parameterdatei 2546 und Beispiel-Ergebnisdatei erhält man durch anklicken der 2547 entsprechenden 2548 Links mit der linken Maustaste bei gleichzeitig gedrückter 2549 <SHIFT>-Taste). 2550 Werden Abweichungen zwischen den beiden Ergebnisdateien festgestellt, 2551 so 2552 ist unbedingt eine Diskussion in der Modellarbeitsgruppe notwendig.</li> 2553 <li>Nach Durchführung der Tests müssen bei 2554 Bedarf die 2555 vorläufig 2556 im temporären Modul eingetragenen Variablen in die 2557 eigentlichen 2558 Module 2559 des Modells übertragen werden. Dazu ist die Modul-Datei <tt>modules.f90</tt> 2560 auszuchecken. Nach Eintrag der Variablen sollte das Modell noch einmal 2561 testweise übersetzt werden, um Sytaxfehler 2562 auszuschließen 2563 (evtl. 2564 auch noch einmal einen Probelauf starten). Die Verwendung des 2565 temporären 2566 Moduls ist aus allen entsprechenden Programmteilen zu entfernen.</li> 2567 <li>Vergabe einer neuen Versionsnummer (Variable <tt>version</tt> 2568 im 2569 Modul <tt>steuer</tt>), 2570 Bsp.: 2.3c. Kleinere Änderungen werden durch 2571 angehängten 2572 kleinen 2573 Buchstaben kenntlich gemacht, größere 2574 Änderungen durch 2575 die Zahl hinter dem Punkt. Eine Inkrementierung der Zahl vor dem Punkt 2576 geschieht nur bei ganz grundsätzlichen Revisionen des Modells, 2577 die 2578 weite Teile des Modellcodes betreffen.</li> <li>Nun 2579 können die ausgecheckten Dateien wieder eingecheckt 2580 werden. 2581 Vor 2582 dem Einchecken sind in den Dateien die Vermerke unter "Aktuelle 2583 Aenderungen" 2584 zu entfernen und beim check-in - Kommando als Log-Text anzugeben. Sie 2585 erscheinen 2586 dann in der jeweiligen Datei unter "Fruehere Aenderungen" zusammen mit 2587 der vom RCS-System vergebenen Versionsnummer (dafür sorgt die 2588 $Log: 2589 ... $ - Zeile, die in jeder Quelltextdatei steht).</li> <li>Eine 2590 Zusammenfassung der durchgeführten Arbeiten ist im <a href="#Kapitel1.0">Kapitel 2591 1.0</a> einzutragen.</li> <li>Die 2592 Modellarbeitsgruppe ist über den Abschluß der 2593 durchgeführten 2594 Arbeiten zu informieren.</li> 2595 </ol><h2> 2686 code</h2>Currently 2687 (Sep 23rd, 2009), only selected users are allowed to commit changes to 2688 the repository. Some rules about how to check in modified code will be 2689 added later.<br> 2690 <h2> 2596 2691 <a name="Kapitel3.0"></a>3.0 Description 2597 2692 of selected parts of the 2598 2693 model source code</h2> 2599 <b>This chapter is still under construction (Feb. 04).</b> 2600 <p>Dieser Abschnitt verweist auf vorhandene Beschreibungen der 2601 genauen 2602 Funktionsweise einzelner, spezieller Modellteile. Dazu gehört 2603 z.B. 2604 die Beschreibung bestimmter numerischer Verfahren (auch ihr 2605 theoretischer 2606 Hintergrund), die Erklärung besonderer Prinzipien bei der 2607 Parallelisierung 2608 oder Erläuterungen zur Vorgehensweise bei der Lösung 2609 bestimmter 2610 Aufgaben (z.B. der Berechnung horizontal gemittelter 2611 Größen). 2612 <br>Die Beschreibungen liegen in jeweils getrennten Dokumenten, 2613 in der 2614 Regel als TeX- und/oder Postscript-Datei vor und sind über die 2615 folgenden 2616 Links zugänglich. 2694 <b></b> 2695 <p>This section gives links to existing descriptions about special 2696 parts of the PALM code, e.g. special numerical methods that are used 2697 (including mathematical background). Documents are stored in single 2698 files with different formats. Some are in German (sorry for the 2699 inconvenience). 2617 2700 <br> 2618 2701 </p><ul> <li> <a name="UPS"></a>upstream-spline
Note: See TracChangeset
for help on using the changeset viewer.