- Timestamp:
- Dec 14, 2009 10:40:54 AM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/DOC/tec/technical_documentation.html
r392 r414 47 47 <td><b>description</b></td> <td><b>affected 48 48 routines</b></td> </tr> <tr nosave="" valign="top"> <td>06/07/98</td> <td>SR</td> 49 <td>1.0</td> <td>C</td> <td nosave=""> Änderung49 <td>1.0</td> <td>C</td> <td nosave="">Ãnderung 50 50 von Default-Werten 51 51 (initializing_actions, end_time, 52 fcl_factor, use_prior_plot1d_parameters), zus ätzliche53 Überprüfung54 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, 55 55 ny, 56 56 nz, 57 fcl_factor), H öhe der Schnittebenen der Horizontalschnitte57 fcl_factor), Höhe der Schnittebenen der Horizontalschnitte 58 58 wird 59 59 auf … … 64 64 Programmteilen benutzt 65 65 (USE). <b>Wichtig</b>: dieses Modul steht in der Datei <tt>module_test.f90. 66 D</tt>iese Datei mu ßin der mrun- und add_library -66 D</tt>iese Datei muà in der mrun- und add_library - 67 67 Konfigurationsdatei 68 68 bei der Variablen <tt>module_files</tt> 69 aufgef ührt werden.</td> <td nosave="">(fast)69 aufgefÃŒhrt werden.</td> <td nosave="">(fast) 70 70 alle <br> 71 71 Neu: module_test</td> </tr> <tr nosave="" valign="top"> <td>16/07/98</td> <td>SR</td> … … 90 90 <td>Entfernung der Sonderbehandlung der Advektionsterme 91 91 (Piascek-Williams) 92 bei k=1. Terme werden nun überall gleich berechnet. Evtl.92 bei k=1. Terme werden nun ÃŒberall gleich berechnet. Evtl. 93 93 leichte 94 Erh öhung der Varianzen in Oberflächennähe.</td>94 Erhöhung der Varianzen in OberflÀchennÀhe.</td> 95 95 <td nosave="">advec_u_pw, advec_v_pw, advec_w_pw</td> 96 96 </tr> <tr nosave="" valign="top"> <td nosave="">"</td> <td>SR</td> <td>1.0</td> … … 98 98 adjust_mixing_length) erfolgt nur 99 99 noch bei 100 k=1. Unrealistische sekund äre Maxima in den Km-Profilen100 k=1. Unrealistische sekundÀre Maxima in den Km-Profilen 101 101 oberhalb 102 102 der 103 Prandtl-Schicht werden damit vermieden. Leichte Änderungen in103 Prandtl-Schicht werden damit vermieden. Leichte Ãnderungen in 104 104 den 105 oberfl ächennahen Turbulenzgrößen sind die105 oberflÀchennahen TurbulenzgröÃen sind die 106 106 Folge.</td> <td>diffusivities</td> </tr> <tr nosave="" valign="top"> <td>"</td> <td>SR</td> 107 107 <td>1.0</td> <td>B</td> <td>Falsche 108 Berechnung der ben ötigten CPU-Zeit pro Sekunde108 Berechnung der benötigten CPU-Zeit pro Sekunde 109 109 simulierter 110 Zeit bei Fortsetzungsl äufen korrigiert. <p>Berechnung110 Zeit bei FortsetzungslÀufen korrigiert. <p>Berechnung 111 111 von cross_uxmin(1) und cross_uxmax(1) aus den 112 112 Werten von 113 ug und vg entfernt. Lieferte falsche Wertebereiche f ür Plots113 ug und vg entfernt. Lieferte falsche Wertebereiche fÃŒr Plots 114 114 von 115 115 Vertikalprofilen, … … 134 134 Felder herausgeschrieben. </p> <p>Bei 135 135 SOR-Aufrufen war bisher dp nicht initialisiert. 136 Entsprechend ge ändert.</p> </td> <td>init_3d_model,136 Entsprechend geÀndert.</p> </td> <td>init_3d_model, 137 137 plot_2d, pres</td> </tr> <tr nosave="" valign="top"> <td>"</td> <td>SR</td> 138 138 <td>1.0</td> <td>C</td> <td nosave="">Ausgabeformate 139 f ür Statistik-Ausgaben auf139 fÃŒr Statistik-Ausgaben auf 140 140 Datei RUN_CONTROL 141 141 um einige Kommastellen erweitert. <p>Lieferung von 142 ausf ührlicheren Informationen bei142 ausfÃŒhrlicheren Informationen bei 143 143 Programmabbruch 144 144 durch Unterschreitung des minimal erlaubten Zeitschritts.</p> </td> … … 146 146 <td>SR</td> <td>1.0</td> <td>N</td> 147 147 <td>Rechnungen mit geneigter Ebene zur Simulation von 148 Kaltluftabfl üssen148 KaltluftabflÃŒssen 149 149 sind erlaubt. Gestartet und gesteuert werden sie mit dem 150 150 Intialisierungsparameter <tt>alpha_surface</tt>. 151 151 Realisiert wird die geneigte Ebene durch entsprechende 152 Schr ägstellung153 des Schwerkraftvektors. Dadurch erh ält auch die u-Komponente152 SchrÀgstellung 153 des Schwerkraftvektors. Dadurch erhÀlt auch die u-Komponente 154 154 der 155 155 Windgeschwindigkeit … … 157 157 Implementierung ist noch nicht abgeschlossen. Die 158 158 Simulation von 159 Kaltluftabfl üssen erfordert bisher noch Eingriffe "von Hand".</p>159 KaltluftabflÃŒssen erfordert bisher noch Eingriffe "von Hand".</p> 160 160 </td> <td>buoyancy, check_parameters, header, leap_frog, 161 161 modules, … … 163 163 read_var_list, write_var_list</td> </tr> <tr nosave="" valign="top"> <td>15/12/98</td> <td>SR</td> 164 164 <td>1.0</td> <td>N</td> <td>Vereinbarung 165 einer Cache-Layer zur schnelleren bin ären165 einer Cache-Layer zur schnelleren binÀren 166 166 I/O auf 167 167 T3E-Rechnern. <p>Ermittelung wahrer I/O-Zeiten 168 f ür binäre168 fÃŒr binÀre 169 169 Ein-/Ausgabe auf 170 170 T3E-Rechnern mit TIMEF.</p> </td> <td nosave="">check_open, … … 173 173 <td>SR</td> <td>1.0</td> <td>B</td> 174 174 <td>Bekanntgabe von run_description_header an alle 175 PE's. <p>Diese Programm änderung erfordert auf175 PE's. <p>Diese ProgrammÀnderung erfordert auf 176 176 T3E-Rechnern 177 177 zwingend den 178 178 Einsatz der MPI-Version mpt.1.2.0.1. In dieser neuen Version tritt auch 179 der alte Fehler in MPI_ALLREDUCE nicht mehr auf, so da ßab179 der alte Fehler in MPI_ALLREDUCE nicht mehr auf, so daà ab 180 180 jetzt 181 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> 182 182 <td>parles</td> </tr> <tr nosave="" valign="top"> <td nosave="">05/02/99</td> 183 183 <td>SR</td> <td>1.1</td> <td>C/N</td> 184 184 <td>Implementierung des Upstream-Spline- Verfahrens zur Advektion 185 185 von Impuls 186 und skalaren Gr ößen (neue UP's, s. rechts).186 und skalaren GröÃen (neue UP's, s. rechts). 187 187 Berechnung von 188 188 Variablen 189 f ür spezielle Advektionsverfahren in eigenständiges189 fÃŒr spezielle Advektionsverfahren in eigenstÀndiges 190 190 UP 191 191 ausgelagert … … 193 193 Reines 194 194 Euler-/Upstream-Verfahren 195 kann gerechnet werden. <p>Schnelle bin äre I/O195 kann gerechnet werden. <p>Schnelle binÀre I/O 196 196 (s. 15/12/98) wieder abgeschaltet, da 197 197 Dateien 198 nicht wieder einlesbar sind. Druckst örung (p) wird zu Beginn198 nicht wieder einlesbar sind. Druckstörung (p) wird zu Beginn 199 199 auf 0 200 200 gesetzt, damit sie keine Tendenzen im ersten Zeitschritt hervorruft.</p> … … 209 209 spline_x, spline_y, spline_z</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">24/02/99</td> 210 210 <td>SR</td> <td>1.1a</td> <td>C/N</td> 211 <td>Kleinere Erg änzungen / Änderungen212 f ür Upstream211 <td>Kleinere ErgÀnzungen / Ãnderungen 212 fÃŒr Upstream 213 213 - Spline, <br> 214 214 Diverse Verbesserungen bei Berechnung der Diffusionskoeffizienten … … 218 218 auch unterhalb des entsprechenden Gitterpunktes stabil ist), 219 219 Wertebegrenzung 220 bei Spline - Berechnung zur Verhinderung von " Überschwingern"</td>220 bei Spline - Berechnung zur Verhinderung von "Ãberschwingern"</td> 221 221 <td>check_parameters, diffusion_e, diffusion_pt, diffusivities, 222 222 header, … … 230 230 <td>advec_s_bc, leap_frog</td> </tr> <tr nosave="" valign="top"> <td nosave="">01/03/99</td> 231 231 <td>MS</td> <td>1.1b</td> <td>N</td> 232 <td>Kleinere Erg änzungen für Upstream -232 <td>Kleinere ErgÀnzungen fÃŒr Upstream - 233 233 Spline: <ul> <li>Gradientenkontrolle 234 f ür den nicht-parallelen Teil</li> <li>statistische235 Auswertung über den prozentualen Anteil234 fÃŒr den nicht-parallelen Teil</li> <li>statistische 235 Auswertung ÃŒber den prozentualen Anteil 236 236 des 237 237 Upstream-Verfahrens … … 240 240 </tr> <tr nosave="" valign="top"> <td>"</td> 241 241 <td>SR</td> <td>1.1b</td> <td>N</td> 242 <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, 243 243 modules, 244 244 parin, plot_3d <p><b>Neu:</b> <br> 245 245 write_compressed</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">03/03/99</td> 246 246 <td>SR</td> <td>1.1c</td> <td>B/C</td> 247 <td>Schnelle bin äre I/O auf T3E-Rechnern funktioniert247 <td>Schnelle binÀre I/O auf T3E-Rechnern funktioniert 248 248 jetzt. 249 249 Beim 250 bin ären Einlesen dürfen Records nicht mit dem Trick "<tt>READ251 (..) 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 252 252 erfolgen 253 253 jetzt 254 grunds ätzlich mit <tt>TIMEF.</tt></td> <td>check_open,254 grundsÀtzlich mit <tt>TIMEF.</tt></td> <td>check_open, 255 255 cpu_zeitmessung, init_3d_model, modules, 256 256 write_3d_binary</td> </tr> <tr nosave="" valign="top"> <td nosave="">25/03/99</td> 257 257 <td>SR</td> <td>1.1d</td> <td>N</td> 258 <td>Einbau einer Rayleigh-D ämpfungsschicht sowie258 <td>Einbau einer Rayleigh-DÀmpfungsschicht sowie 259 259 flexiblere 260 260 Wahl der 261 numerischen Begrenzer f ür das Upstream - Spline - Verfahren.</td>261 numerischen Begrenzer fÃŒr das Upstream - Spline - Verfahren.</td> 262 262 <td>check_parameters, header, init_3d_model, leap_frog, modules, 263 263 parin, … … 269 269 <td>1.1e</td> <td>N</td> <td nosave="">Partikelplots 270 270 mit vtk-Graphiksoftware sind auf t3eh 271 m öglich.271 möglich. 272 272 Ausgabe erfolgt im dvr-Format. Bisher noch nicht nutzbar, da noch 273 Einstellm öglichkeiten274 f ür die Partikelquelle(n) fehlen. <b>Wichtig:</b>273 Einstellmöglichkeiten 274 fÃŒr die Partikelquelle(n) fehlen. <b>Wichtig:</b> 275 275 Die Benutzung 276 276 der 277 vtk-Software erfordert die zus ätzliche mrun-Option <tt>-g</tt>.277 vtk-Software erfordert die zusÀtzliche mrun-Option <tt>-g</tt>. 278 278 Das 279 279 Hauptprogramm <tt>parles</tt> wird in diesem Fall zum … … 293 293 nx /= ny auf.</td> <td>cpu_zeitmessung, spline_y</td> 294 294 </tr> <tr nosave="" valign="top"> <td nosave="">28/12/99</td> <td>SR</td> <td>1.1f</td> 295 <td>N</td> <td>Plot von Isooberfl ächen mit295 <td>N</td> <td>Plot von IsooberflÀchen mit 296 296 vtk-Graphiksoftware sind auf 297 297 t3eh 298 m öglich. Ausgabe erfolgt im dvr-Format. Bisher noch nicht298 möglich. Ausgabe erfolgt im dvr-Format. Bisher noch nicht 299 299 allgemein 300 nutzbar, da noch diverse manuelle Eingriffe n ötig sind.</td>300 nutzbar, da noch diverse manuelle Eingriffe nötig sind.</td> 301 301 <td>header, leap_frog, modules, parin, plot_particles, 302 302 read_var_list, write_var_list <p><b>Neu:</b> … … 304 304 plot_isosurface</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">10/01/00</td> 305 305 <td>SR</td> <td>1.2</td> <td>C</td> 306 <td>Randbedingungen f ür w (=0) werden explizit gesetzt,306 <td>Randbedingungen fÃŒr w (=0) werden explizit gesetzt, 307 307 und 308 308 zwar nicht … … 310 310 bei <tt>k = 311 311 nzt</tt>, 312 wo bislang der vom vorigen Zeitschritt durch den Druckl öser312 wo bislang der vom vorigen Zeitschritt durch den Drucklöser 313 313 ermittelte 314 Wert stand. Auch jetzt wird der Wert bei k = nzt anschlie ßend314 Wert stand. Auch jetzt wird der Wert bei k = nzt anschlieÃend 315 315 durch 316 den Druckl öser modifiziert. Soweit bis jetzt abzusehen, hat316 den Drucklöser modifiziert. Soweit bis jetzt abzusehen, hat 317 317 diese 318 Änderung nur minimale Auswirkungen auf den Prognoseverlauf.318 Ãnderung nur minimale Auswirkungen auf den Prognoseverlauf. 319 319 Evtl. 320 gibt es st ärkere Auswirkungen bei exzessivem Auftreten von320 gibt es stÀrkere Auswirkungen bei exzessivem Auftreten von 321 321 Schwerewellen 322 322 am oberen Rand (die dort aber eigentlich sowieso nicht 323 hingeh ören). <br>324 Variablen übergabe an <tt>boundary_conds</tt>323 hingehören). <br> 324 VariablenÃŒbergabe an <tt>boundary_conds</tt> 325 325 erfolgt <br> 326 über Modul und nicht mehr per Parameterliste. <p>Bei327 Galilei-Transformation wird jetzt defaultm äßig326 ÃŒber Modul und nicht mehr per Parameterliste. <p>Bei 327 Galilei-Transformation wird jetzt defaultmÀÃig 328 328 der geostrophische 329 329 Wind als Translationsgeschwindigkeit verwendet (neuer 330 330 Initialisierungsparameter <tt>use_ug_for_galilei_tr</tt>). 331 Dadurch wird das verst ärkte Auftreten von Schwerewellen bei332 h öheren333 geostrophischen Windgeschwindigkeiten unterdr ückt. </p>331 Dadurch wird das verstÀrkte Auftreten von Schwerewellen bei 332 höheren 333 geostrophischen Windgeschwindigkeiten unterdrÃŒckt. </p> 334 334 <p>3D-Feldvariablen (u, v, w, pt, e, km, kh) sowie einige 335 335 2D-Variablen 336 336 (usws, vsws, rif, shf) sind nun als Pointer deklariert, denen nach 337 337 jedem 338 Zeitschritt neue Targets zugeordnet werden (f ür die339 ben ötigten 338 Zeitschritt neue Targets zugeordnet werden (fÃŒr die 339 benötigten 340 340 3 Zeitebenen existieren die Targets ..._1 (z.B. u_1), ..._2, ..._3). 341 341 Durch 342 die jeweilige Neuzuordnung der Targets entf ällt das342 die jeweilige Neuzuordnung der Targets entfÀllt das 343 343 Umspeichern 344 344 der … … 347 347 Die Umordnung der Zeitebenen wurde ins neue UP <tt>swap_timelevel</tt> 348 348 ausgelagert. <br> <b>Wichtig:</b> Die 349 Übersetzungszeiten haben durch die350 Einf ührung349 Ãbersetzungszeiten haben durch die 350 EinfÃŒhrung 351 351 der Pointer t.w. sehr deutlich zugenommen (<tt>leap_frog</tt> 352 ben ötigt353 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 354 354 Compiler 355 355 jetzt wesentlich mehr Speicher. Der "memory"-Parameter in den … … 357 357 (.mrun.config, .al.config) muss nun mindestens den Wert 60 (MB) (t3eh, 358 358 t3ej2) bzw. 45 (t3eb) und 96 (vpp) haben. <br> <b>Achtung:</b> 359 bislang wurde diese Änderung nur für359 bislang wurde diese Ãnderung nur fÃŒr 360 360 das 361 361 reine Leapfrog-Verfahren getestet! </p> <p>Neue 362 362 Laderoption "<tt>-D preset=nan</tt>" empfohlen (s. 363 363 Beispiel-Konfigurationsdatei). 364 Ebenso empfehlenswert f ür Debug-Läufe: Compileroption364 Ebenso empfehlenswert fÃŒr Debug-LÀufe: Compileroption 365 365 "<tt>-ei</tt>". 366 366 Diese bewirkt einen Programmabbruch, falls mit nicht initialisierten 367 367 Variablen 368 gearbeitet wird (ist aber zeitaufwendig in der Ausf ührung!).</p>368 gearbeitet wird (ist aber zeitaufwendig in der AusfÃŒhrung!).</p> 369 369 </td> <td>boundary_conds, check_parameters, header, 370 370 init_3d_model, … … 376 376 <td>Der Asselin-Filter wird vor Umordnung der Zeitebenen 377 377 aufgerufen. Dies 378 geschah bisher f älschlicherweise erst danach, sodass der378 geschah bisher fÀlschlicherweise erst danach, sodass der 379 379 Zeitfilter 380 380 asymmetrisch arbeitete. <br> <b>Achtung:</b> 381 diese Änderung hat Auswirkungen auf den381 diese Ãnderung hat Auswirkungen auf den 382 382 Prognoseverlauf. 383 Die Testergebnis-Datei wurde entsprechend ge ändert.</td>383 Die Testergebnis-Datei wurde entsprechend geÀndert.</td> 384 384 <td>leap_frog</td> </tr> <tr nosave=""> 385 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> … … 391 391 zeitlichen 392 392 Filterung unterzogen. Dies wird nun unterbunden. </p> <p>Die 393 Divergenz der Impulsfl üsse an der Prandtl - Schicht -393 Divergenz der ImpulsflÃŒsse an der Prandtl - Schicht - 394 394 Obergrenze 395 wird nun über eine ganze Gitterweite gebildet. Die bisherige395 wird nun ÃŒber eine ganze Gitterweite gebildet. Die bisherige 396 396 Methode 397 397 lieferte im Vergleich mit Ergebnissen anderer LES-Modelle zu 398 gro ße399 Impulsfl üsse. </p> <p>Bei der398 groÃe 399 ImpulsflÃŒsse. </p> <p>Bei der 400 400 Berechnung der Diffusion der Temperatur werden die 401 401 Diffusionskoeffizienten 402 nun wieder durch Mittelung f ür die entsprechenden 402 nun wieder durch Mittelung fÃŒr die entsprechenden 403 403 Gitterpunkte 404 404 bereitgestllt. Das Maximumkriterium, das im Zuge der Implementierung … … 407 407 rasches 408 408 Anwachsen der Grenzschicht mit der Zeit. </p> <p><b>Achtung:</b> 409 diese Änderungen haben Auswirkungen auf409 diese Ãnderungen haben Auswirkungen auf 410 410 den Prognoseverlauf. 411 Die Testergebnis-Datei wurde entsprechend ge ändert.</p> </td>411 Die Testergebnis-Datei wurde entsprechend geÀndert.</p> </td> 412 412 <td nosave="" valign="top">init_1d_model, 413 413 asselin_filter, … … 415 415 diffusion_pt</td> </tr> <tr nosave="" valign="top"> <td nosave="">17/02/00</td> 416 416 <td>MS</td> <td>2.0</td> <td>N</td> 417 <td><b>Allgemeine Beschreibung der Änderungen:</b>417 <td><b>Allgemeine Beschreibung der Ãnderungen:</b> 418 418 <br> 419 vollst ändige Implementierung des wolkenphysikalischen419 vollstÀndige Implementierung des wolkenphysikalischen 420 420 Moduls <p>Zur Implementierung der Wolkenphysik waren 421 421 umfangreiche 422 Änderungen422 Ãnderungen 423 423 in zahlreichen Programmteilen notwendig (Beschreibung <a href="#Cloud_physics">siehe 424 424 unten</a> ). Im einzelnen wurden folgende Unprogramme 425 ge ändert: <br> <b>modules</b>: 426 <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 427 427 die Prognose 428 428 des Gesamtwassergehaltes <br> … … 433 433 write_var_list, 434 434 write_3d_binary: </b>In/Output 435 der neuen Initialisierungsparamter und Felder f ür436 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> 437 437 Initialisierung 438 438 des … … 442 442 Temperatur,</font> <br> <font color="#ff0000">lpt</font><font color="#000000">: 443 443 potentielle 444 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">: 445 445 Gesamtwassergehalt,</font> <br> <font color="#ff0000">qv</font><font color="#000000">: 446 446 spezifische 447 447 Feuchte,</font> <br> <font color="#ff0000">ql</font><font color="#000000">: 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>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> 450 450 <font color="#ff0000">w"q", w*q*, wq</font><font color="#000000">: 451 451 Wasserfluss (subs., aufg., ges.),</font> <br> <font color="#ff0000">w"qv", w*qv*, wqv</font><font color="#000000">: … … 454 454 Feuchte/Wasserfelder <br> <b>init_1d_model</b>: 455 455 Anpassung an Rechnungen mit Feuchte <br> <b>header</b>: 456 Headerausgaben f ür Feuchtebehandlung und456 Headerausgaben fÃŒr Feuchtebehandlung und 457 457 Wolkenphysik 458 eingef ügt, außerdem Formatnummern erweitert um Platz459 f ür458 eingefÃŒgt, auÃerdem Formatnummern erweitert um Platz 459 fÃŒr 460 460 weitere Ausgaben zu schaffen <br> <b>leap_frog:</b> 461 461 <br> … … 468 468 flow_statistics: </b>Anpassung 469 469 an Feuchtebehandlung <br> <b>plot2d, plot3d</b>: 470 Ausgabem öglichkeiten von q und ql <br> <b>buoyancy</b>:470 Ausgabemöglichkeiten von q und ql <br> <b>buoyancy</b>: 471 471 Anpassung des Autriebstermes an 472 472 Feuchtebehandlung 473 durch neue Übergabeparameter: Temperatur + Nummer des473 durch neue Ãbergabeparameter: Temperatur + Nummer des 474 474 entsprechenden 475 475 mittleren Temperaturprofils <br> <b>diffusion_e</b>: 476 476 Anpassung der Mischungswegberechnung an 477 477 Feuchtebehandlung 478 durch Variation der Übergabeparameter <br> <b>production_e</b>:479 TKE-Produktion durch Auftriebskr äfte an478 durch Variation der Ãbergabeparameter <br> <b>production_e</b>: 479 TKE-Produktion durch AuftriebskrÀfte an 480 480 Feuchtebehandlung 481 481 angepasst <br> <b>diffusivities</b>: … … 483 483 Prandtl-Schicht fuer den 484 484 Gesamtwassergehalt 485 eingef ügt </p> <p>Zusätzlich485 eingefÃŒgt </p> <p>ZusÀtzlich 486 486 wurden folgende Unterprogramme NEU 487 eingef ührt: <br> <b>inti_cloud_physics</b>:487 eingefÃŒhrt: <br> <b>inti_cloud_physics</b>: 488 488 Initialisierung wolkenphysikalischer 489 489 Parameter <br> <b>prognostic_equations</b>: 490 490 Auslagerung der prognostischen 491 491 Gleichungen 492 aus leap_frog, zus ätzliche porgnostische Gleichung493 f ür den492 aus leap_frog, zusÀtzliche porgnostische Gleichung 493 fÃŒr den 494 494 Gesamtwassergehalt <br> <b>diffusion_s</b>: 495 495 Berechnung der Diffusion skalarer 496 Gr ößen496 GröÃen 497 497 (Temperatur und Gesamtwassergehalt), ersetzt diffusion_pt <br> 498 498 <b>calc_liquid_water_content</b>: Berechnung des 499 Fl üssigwassergehalts <br> <b>calc_radiation</b>:499 FlÃŒssigwassergehalts <br> <b>calc_radiation</b>: 500 500 Parametrisierung langwelliger 501 501 Strahlungsprozesse 502 über des Schema der effektiven Emissivität <br>502 ÃŒber des Schema der effektiven EmissivitÀt <br> 503 503 <b>calc_precipitation</b>:Parametrisierung von 504 504 Niederschlagsprozessen 505 505 mit Hilfe eines vereinfachten Kesslerschemas <br> <b>impact_of_latent_heat</b>: 506 Ber ücksichtigung der durch506 BerÃŒcksichtigung der durch 507 507 Niederschlagsprozesse 508 zur ückbleibenden latenten Waerme innerhalb eines508 zurÃŒckbleibenden latenten Waerme innerhalb eines 509 509 Gittervolumens <br> <b>compute_vpt</b>: 510 510 Berechnung der virtuellen potentiellen … … 526 526 der Niderschlagsparametrisierung</td> </tr> <tr> <td><tt>bc_q_b</tt></td> 527 527 <td>C</td> <td>'dirichlet'</td> <td>untere 528 Randbedingung f ür q (siehe auch bc_pt_b)</td> </tr>528 Randbedingung fÃŒr q (siehe auch bc_pt_b)</td> </tr> 529 529 <tr> <td><tt>bc_q_t</tt></td> <td>C</td> 530 <td>'neumann'</td> <td>obere RB f ür q530 <td>'neumann'</td> <td>obere RB fÃŒr q 531 531 (siehe auch bc_pt_t)</td> </tr> <tr> <td><tt>q_surface</tt></td> 532 532 <td>R</td> <td>0.0</td> <td>Feuchtewert 533 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> 534 534 <td>R</td> <td>0.0</td> <td>vgl. 535 535 pt_surface_initial_change</td> </tr> <tr> <td><tt>q_vertical_gradient</tt></td> … … 540 540 pt_vertical_gradient_level</td> </tr> <tr> <td><tt>surface_pressure</tt></td> 541 541 <td>R</td> <td>1013.25</td> <td>Luftdruck 542 an der Erdoberfl äche</td> </tr> <tr> <td><tt>surface_waterflux</tt></td>543 <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 544 544 Wasser/Feuchtefluss</td> </tr> </tbody> </table> 545 Fortsetzungsl äufe mit Daten, die mit Modellversionen kleiner545 FortsetzungslÀufe mit Daten, die mit Modellversionen kleiner 546 546 2.0 547 547 erzeugt … … 552 552 deshalb nach Abschluss eines kompletten Modelllaufs erfolgen. 553 553 </p> <p><b>Achtung:</b> diese 554 Änderungen haben Auswirkungen auf554 Ãnderungen haben Auswirkungen auf 555 555 den Prognoseverlauf. 556 Die Testergebnis-Datei wurde entsprechend ge ändert.</p> </td>556 Die Testergebnis-Datei wurde entsprechend geÀndert.</p> </td> 557 557 <td nosave="" valign="top">modules, parin, 558 558 read_var_list, … … 569 569 Umstellung auf Modellversion 2.0 wurden 570 570 noch kleine 571 Änderungen an den Unterprogrammen diffusion_e, production_e571 Ãnderungen an den Unterprogrammen diffusion_e, production_e 572 572 und 573 573 diffusivities 574 574 vorgenommen. In allen genannten <br> 575 575 Unterprogrammen wurde bisher ein strengeres Kriterium zur Bestimmung 576 der Schichtungsverh ältnisse verwendet. Ein Modellvergleich576 der SchichtungsverhÀltnisse verwendet. Ein Modellvergleich 577 577 zeigte 578 578 jedoch, dass dadurch die Diffusion im Bereich der Inversion zu gross 579 wird. <p>Eine weitere Erg änzung betrifft die579 wird. <p>Eine weitere ErgÀnzung betrifft die 580 580 Advektionsverfahren: <br> 581 Das Bott-Chlond Advektionsverfahren ist nun auch f ür581 Das Bott-Chlond Advektionsverfahren ist nun auch fÃŒr 582 582 Rechnungen 583 mit Feuchte/Wolkenphysik verf ügbar. </p> <p><b>Achtung:</b>584 diese Änderungen haben Auswirkungen auf583 mit Feuchte/Wolkenphysik verfÃŒgbar. </p> <p><b>Achtung:</b> 584 diese Ãnderungen haben Auswirkungen auf 585 585 den Prognoseverlauf. 586 Die Testergebnis-Datei wurde entsprechend ge ändert.</p> </td>586 Die Testergebnis-Datei wurde entsprechend geÀndert.</p> </td> 587 587 <td nosave="" valign="top">diffusion_e, 588 588 production_e, 589 589 diffusivities, check_parameters</td> </tr> <tr nosave="" valign="top"> <td nosave="">26/04/00</td> 590 590 <td>SR</td> <td>2.0a</td> <td>C</td> 591 <td>Vollst ändige Umstellung von vtk-Grafik-Software auf591 <td>VollstÀndige Umstellung von vtk-Grafik-Software auf 592 592 dvrp-Software. 593 593 Grafik-Ausgaben im dvr-Format sind mittlerweile begrenzt nutzbar 594 (Isooberfl ächen594 (IsooberflÀchen 595 595 sind darstellbar, siehe neue Laufparameter <tt>dt_dvrp, 596 596 threshold</tt>), 597 597 es wird aber in der Anwendung in naher Zukunft noch diverse 598 Änderungen599 geben. Partikelausgabe vorerst nicht mehr m öglich (Verlagerung598 Ãnderungen 599 geben. Partikelausgabe vorerst nicht mehr möglich (Verlagerung 600 600 schon).</td> <td>header, init_particles, init_3d_model, 601 601 leap_frog, modules, … … 607 607 init_dvrp, plot_dvrp</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">26/04/00</td> 608 608 <td>SR</td> <td>2.0a</td> <td>C</td> 609 <td>Durch zus ätzlichen zyklischen Rand im609 <td>Durch zusÀtzlichen zyklischen Rand im 610 610 Bott-Chlond-Schema 611 611 sehr 612 612 viele SENDRECV-Aufrufe eingespart. <br> 613 Änderung des Namens einer eingelesenen Environment-Variablen613 Ãnderung des Namens einer eingelesenen Environment-Variablen 614 614 in <tt>check_open</tt> 615 615 von <tt>remote_addres</tt> nach <tt>return_addres</tt>. … … 618 618 Initialisierungsparameter <tt>npex</tt> und <tt>npey</tt>. 619 619 Durch 620 Änderungen621 am Header-Format wurde Testergebnis-Datei ge ändert.</td>620 Ãnderungen 621 am Header-Format wurde Testergebnis-Datei geÀndert.</td> 622 622 <td>advec_s_bc, check_open, header, init_pegrid, modules, parin, 623 623 read_var_list, 624 624 write_var_list</td> </tr> <tr nosave="" valign="top"> <td nosave="">26/04/00</td> 625 625 <td>SR</td> <td>2.0a</td> <td>N</td> 626 <td>Rechnungen mit geneigter Oberfl äche erlauben jetzt626 <td>Rechnungen mit geneigter OberflÀche erlauben jetzt 627 627 die 628 628 Vorgabe 629 629 eines stabil geschichteten Temperaturprofils mit konstantem Gradienten. 630 Hinzuschalten von Feuchte funktioniert f ür Hangrechnungen noch630 Hinzuschalten von Feuchte funktioniert fÃŒr Hangrechnungen noch 631 631 nicht.</td> <td>advec_s_bc, buoyancy, init_3d_model, 632 632 modules, … … 634 634 init_slope</p> </td> </tr> <tr nosave="" valign="top"> <td>13/06/00</td> <td>MS</td> 635 635 <td>2.0a</td> <td nosave="">B</td> <td>Es 636 werden nun auch die seitlichen R änder für ql636 werden nun auch die seitlichen RÀnder fÃŒr ql 637 637 gesetzt. 638 Ein Nichtsetzen der Randbedingungen f ührte zu Fehlern bei der638 Ein Nichtsetzen der Randbedingungen fÃŒhrte zu Fehlern bei der 639 639 Ausgabe 640 640 von pt</td> <td nosave="">calc_liquid_water_content</td> 641 641 </tr> <tr nosave="" valign="top"> <td nosave="">03/07/00</td> <td>SR</td> <td>2.0b</td> 642 <td>C</td> <td>F ür eine Reihe von 642 <td>C</td> <td>FÃŒr eine Reihe von 643 643 Unterprogrammen, die innerhalb 644 644 von <tt>prognostic_equation</tt> 645 aufgerufen werden und an die Pointer-Variablen übergeben645 aufgerufen werden und an die Pointer-Variablen ÃŒbergeben 646 646 werden, 647 647 sind 648 648 jetzt explizite Interfaces vereinbart (neues Modul <tt>pointer_interfaces</tt>). 649 649 Auch im Vereinbarungsteil der entsprechenden Unterprogramme sind diese 650 Übergabeparameter jetzt aus Konsistenzgründen (und650 Ãbergabeparameter jetzt aus KonsistenzgrÃŒnden (und 651 651 weil es 652 652 sonst 653 653 Laufzeitfehler gibt) als Pointer-Variablen deklariert. Mit dieser 654 Änderung655 ist der gro ße Speicher- und CPU-Zeit-Bedarf beim656 Übersetzen657 von prognostic_equations wieder auf ein ertr ägliches658 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à 659 659 reduziert. <br> 660 Im Rahmen dieser Änderungen wurde auch an diffusion_e ein660 Im Rahmen dieser Ãnderungen wurde auch an diffusion_e ein 661 661 Hilfsfeld 662 weniger übergeben. <br>663 Aus nicht vollst ändig gekärten Gründen haben662 weniger ÃŒbergeben. <br> 663 Aus nicht vollstÀndig gekÀrten GrÃŒnden haben 664 664 sich die 665 Zahlen in der Testergebnis-Datei minimal ge ändert (und zwar665 Zahlen in der Testergebnis-Datei minimal geÀndert (und zwar 666 666 die 667 667 Spalte 668 mit der Gesamtdivergenz nach Aufruf des Druckl ösers).</td>668 mit der Gesamtdivergenz nach Aufruf des Drucklösers).</td> 669 669 <td>advec_s_bc, buoyancy, diffusion_e, diffusion_s, diffusion_u, 670 670 diffusion_v, 671 671 diffusion_w, disturb_field, modules, production_e, prognostic_equations</td> 672 672 </tr> <tr nosave="" valign="top"> <td nosave="">04/07/00</td> <td>SR</td> <td>2.0b</td> 673 <td>B</td> <td>Diriclet-Randbedingungen f ür673 <td>B</td> <td>Diriclet-Randbedingungen fÃŒr 674 674 Temperatur und 675 Fl üssigwassergehalt676 werden gesetzt. Dies w äre bereits nach Einführung der675 FlÃŒssigwassergehalt 676 werden gesetzt. Dies wÀre bereits nach EinfÃŒhrung der 677 677 Pointer 678 678 zwingend notwendig gewesen. Rechnungen mit vorgegebener 679 Oberfl ächentemperatur679 OberflÀchentemperatur 680 680 schlugen deshalb bislang fehl.</td> <td>boundary_conds</td> 681 681 </tr> <tr nosave="" valign="top"> <td nosave="">07/09/00</td> <td>MS</td> <td>2.0b</td> 682 682 <td>B</td> <td>Die virtuelle potenielle Temperatur 683 wird nun auch f ür683 wird nun auch fÃŒr 684 684 k=nzb und 685 k=nzt+1 berechnet, um die Randwerte auch f ür Ausgabezwecke685 k=nzt+1 berechnet, um die Randwerte auch fÃŒr Ausgabezwecke 686 686 bereichtzustellen. </td> <td>compute_vpt</td> 687 687 </tr> <tr nosave="" valign="top"> <td nosave="">28/12/00</td> <td>SR</td> <td>2.1</td> … … 689 689 optionaler Software-Pakete, 690 690 die zwar 691 zum Modell geh ören, standardmäßig aber691 zum Modell gehören, standardmÀÃig aber 692 692 nicht mit 693 übersetzt693 ÃŒbersetzt 694 694 werden, um so unter anderem Kompilationszeit einzusparen. Siehe neues 695 695 Kapitel 696 696 3.7 in der Modelldokumentation. <p>Einsatz der 697 dvrp-Software wurde vollst ändig698 überarbeitet.697 dvrp-Software wurde vollstÀndig 698 ÃŒberarbeitet. 699 699 Sie ist nun als Software-Paket optional im Modell einsetzbar. Siehe 700 überarbeitetes700 ÃŒberarbeitetes 701 701 Kapitel 4.5.6 in der Modelldokumentation. </p> <p>Die 702 702 Prognose der Partikelverlagerung ist ebenfalls in ein … … 716 716 <td>parles</td> </tr> <tr nosave="" valign="top"> <td nosave="">02/01/01</td> 717 717 <td>SR</td> <td>2.1a</td> <td>C</td> 718 <td>Offene Dateien werden sobald m öglich geschlossen.719 Wiederer öffnung720 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 721 721 Zusammenhang 722 722 wurde das Unterprogramm close_files in close_file umbenannt und hat nun … … 742 742 <br> 743 743 calc_spectra, plot_spectra</p> </td> </tr> <tr valign="top"> <td>25/01/01</td> <td>SR</td> 744 <td>2.1c</td> <td>C/N</td> <td>Druckl öser744 <td>2.1c</td> <td>C/N</td> <td>Drucklöser 745 745 kann jetzt die Singleton-FFT benutzen. Modul <tt>singleton</tt> 746 746 wurde deshalb von <tt>calc_spectra.f90</tt> nach <tt>modules.f90</tt> 747 747 geschoben 748 748 (neuer Initialisierungsparameter <tt>fft_method</tt>). 749 <br>Fl üsse können jetzt bei k=1750 unabh ängig von der749 <br>FlÃŒsse können jetzt bei k=1 750 unabhÀngig von der 751 751 Verwendung 752 752 einer Prandtl-Schicht vorgegeben werden (neuer … … 757 757 umbenannt. <br> 758 758 Modell kann auf DEC-Workstations eingesetzt werden (-D <tt>dec</tt>). 759 <br>Zus ätzliche Zeitmessungen mit dvrp-Software.759 <br>ZusÀtzliche Zeitmessungen mit dvrp-Software. 760 760 Schreiben von 761 761 Partikel-Informationen … … 770 770 parles, module_test</p> </td> </tr> <tr valign="top"> <td>25/01/01</td> <td>SR</td> 771 771 <td>2.1c</td> <td>B</td> <td>Fehler 772 beim öffnen von Unit 23 beseitigt. Modul <tt>interface</tt>772 beim öffnen von Unit 23 beseitigt. Modul <tt>interface</tt> 773 773 fehlte in disturb_field. <br> 774 Wertebereich von theta* in Prandtl-Fluxes eingeschr änkt, weil774 Wertebereich von theta* in Prandtl-Fluxes eingeschrÀnkt, weil 775 775 sonst auf einigen Rechnern overflow droht, wenn u*=0.</td> <td>check_open, 776 776 disturb_field, prandtl_fluxes</td> </tr> <tr valign="top"> <td>30/01/01</td> <td>SR</td> 777 777 <td>2.1d</td> <td>N/C</td> <td>Prognose 778 eines passiven Skalars ist m öglich (anstatt778 eines passiven Skalars ist möglich (anstatt 779 779 Feuchte), 780 780 dazu neuer Initialisierungsparameter <tt>passive_scalar</tt>. 781 <p>Falls kein Pfad f ür ftpcopy existiert, werden in <tt>check_open</tt>782 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 783 783 "unknown" 784 784 gesetzt.</p> </td> <td>asselin_filter, … … 791 791 write_var_list</td> </tr> <tr valign="top"> <td>30/01/01</td> 792 792 <td>SR</td> <td>2.1d</td> <td>B</td> 793 <td>String-Ausgabe f ür use_surface_fluxes korrigiert.793 <td>String-Ausgabe fÃŒr use_surface_fluxes korrigiert. 794 794 Version vom 795 25/01/01 brach deswegen bei Fortsetzungsl äufen ab. <br>795 25/01/01 brach deswegen bei FortsetzungslÀufen ab. <br> 796 796 xy-Plotausgabe von q und ql korrigiert.</td> <td>plot_2d, 797 797 write_var_list</td> </tr> <tr nosave="" valign="top"> <td nosave="">30/03/01</td> 798 798 <td>SR</td> <td>2.2</td> <td>N/C</td> 799 <td>S ämtliche 3D-Felder der Zeitebene t+dt wurden799 <td>SÀmtliche 3D-Felder der Zeitebene t+dt wurden 800 800 entfernt. 801 801 Ebenso 802 802 die Arbeitsfelder work und work1. Hilfsfelder wurden soweit 803 m öglich804 eingespart und werden grunds ätzlich lokal allokiert. Dadurch803 möglich 804 eingespart und werden grundsÀtzlich lokal allokiert. Dadurch 805 805 konnten 806 eine Reihe von Übergabeparametern gestrichen werden (z.B. beim806 eine Reihe von Ãbergabeparametern gestrichen werden (z.B. beim 807 807 Aufruf 808 808 von <tt>pres, advec_s_bc, poisfft, production_e, sor </tt>sowie … … 812 812 der 813 813 Zeitfilter jetzt direkt in der prognostischen Gleichung angewendet wird 814 (bisher geschah dies erst nach Aufruf des Druckl ösers), haben814 (bisher geschah dies erst nach Aufruf des Drucklösers), haben 815 815 sich 816 die Zahlen in der Testergebnis-Datei ge ändert. Die mittleren816 die Zahlen in der Testergebnis-Datei geÀndert. Die mittleren 817 817 Profile 818 818 sind aber unbeeinflusst. <p>Eine weitere kleinere 819 Änderung der Testergebnis-Daten819 Ãnderung der Testergebnis-Daten 820 820 wird dadurch 821 821 verursacht, dass jetzt der Parameter <tt>adjust_mixing_length</tt> 822 defaultm äßig <tt>.FALSE.</tt>823 ist und der Mischungsweg grunds ätzlich zusätzlich822 defaultmÀÃig <tt>.FALSE.</tt> 823 ist und der Mischungsweg grundsÀtzlich zusÀtzlich 824 824 auf <tt>0.7 825 825 * zu</tt> begrenzt wird. Falls <tt>adjust_mixing_length=T</tt>, … … 827 827 analog zu <tt>diffusion_e</tt> 828 828 an allen Gitterpunkten modifiziert. <b>Achtung: </b>Diese 829 Änderungen830 k önnen bei anderen Simulationen t.w. erhebliche Auswirkungen829 Ãnderungen 830 können bei anderen Simulationen t.w. erhebliche Auswirkungen 831 831 auf 832 832 die 833 oberfl ächennahen Ergebnisse haben. </p> <p><b>Kleinere834 Ä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> 835 835 wird 836 836 der Wert der Variablen <tt>surface_pressure</tt> in hPa … … 839 839 Umwandlung in Pa). <br> 840 840 Die Namelist-Namen der Software-Pakete sind t.w. 841 ge ändert. <br>842 Unit 14 (bin äre Ausgabe für843 Fortsetzungsl äufe) wird841 geÀndert. <br> 842 Unit 14 (binÀre Ausgabe fÃŒr 843 FortsetzungslÀufe) wird 844 844 jetzt im Hauptprogramm geschlossen, damit auf diese Datei noch 845 845 benutzergesteuerte 846 Ausgaben erfolgen k önnen. <br>846 Ausgaben erfolgen können. <br> 847 847 Laufparameter werden in keinem Fall mehr auf Unit 14 ausgegeben bzw. 848 von Unit 13 gelesen, d.h. sie gelten jetzt tats ächlich nur849 f ür848 von Unit 13 gelesen, d.h. sie gelten jetzt tatsÀchlich nur 849 fÃŒr 850 850 den jeweils aktuellen Lauf. <br> 851 851 Partikeladvektion funktioniert jetzt auch zusammen mit der … … 853 853 deutsche Variablen-, Unterprogramm- und Modulnamen 854 854 wurden 855 ins Englische übersetzt. Davon sind fast alle Programmteile855 ins Englische ÃŒbersetzt. Davon sind fast alle Programmteile 856 856 betroffen 857 (nicht in rechter Spalte aufgef ührt).</p> </td> <td>advec_particles,857 (nicht in rechter Spalte aufgefÃŒhrt).</p> </td> <td>advec_particles, 858 858 advec_s_bc, advec_s_ups, advec_u_ups, 859 859 advec_v_ups, … … 871 871 </tr> <tr nosave="" valign="top"> <td nosave="">12/07/01</td> <td>SR</td> <td>2.2a</td> 872 872 <td>N/C</td> <td>Defaultwert der unteren 873 Randbedingung f ür die TKE ist ab873 Randbedingung fÃŒr die TKE ist ab 874 874 sofort <tt>bc_e_b 875 875 = 'neumann'</tt>. <br> 876 Partikelquelle ist über Paketparameter steuerbar.876 Partikelquelle ist ÃŒber Paketparameter steuerbar. 877 877 Partikeleigenschaften 878 878 sind mittels benutzereigener Software steuerbar. Der dvrp-Ausgabeweg 879 879 ist 880 jetzt ebenfalls über Paketparameter steuerbar. <br>880 jetzt ebenfalls ÃŒber Paketparameter steuerbar. <br> 881 881 Allen Modulen (bis auf singleton) wurde eine SAVE-Anweisung 882 hinzugef ügt,883 damit die durch sie vereinbarten Variablen w ährend der882 hinzugefÃŒgt, 883 damit die durch sie vereinbarten Variablen wÀhrend der 884 884 Rechnung in 885 keinem Fall undefiniert werden k önnen (diese Gefahr besteht885 keinem Fall undefiniert werden können (diese Gefahr besteht 886 886 z.B. 887 887 auf … … 891 891 plot_dvrp, user_interface</td> </tr> <tr> <td>12/07/01</td> 892 892 <td>SR</td> <td>2.2a</td> <td>B</td> 893 <td>Verschiebung der OPEN-Anweisung f ür Unit 33893 <td>Verschiebung der OPEN-Anweisung fÃŒr Unit 33 894 894 (AVS-FLD-Datei).</td> <td>check_open</td> </tr> 895 895 <tr nosave="" valign="top"> <td nosave="">20/07/01</td> 896 896 <td>SR</td> <td>2.3</td> <td>N</td> 897 <td>Einbau des Mehrgitterverfahrens zur L ösung der897 <td>Einbau des Mehrgitterverfahrens zur Lösung der 898 898 Poisson-Gleichung.</td> <td>check_parameters, 899 899 exchange_horiz, header, init_grid, … … 902 902 poismg</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">21/08/01</td> 903 903 <td>SR</td> <td>2.3a</td> <td>N/C</td> 904 <td>Erg änzung der dvrp-Software: Partikel905 k önnen mit906 Schw änzen904 <td>ErgÀnzung der dvrp-Software: Partikel 905 können mit 906 SchwÀnzen 907 907 versehen werden. Die Ausgabe der dvrp-Plotdaten durch einen separaten 908 908 PE 909 ist m öglich (mrun-Option "<tt>-p dvrp_graphics+1PE</tt>").909 ist möglich (mrun-Option "<tt>-p dvrp_graphics+1PE</tt>"). 910 910 Der 911 911 Kommunikator <tt>MPI_COMM_WORLD</tt> … … 914 914 Wandbereich kann wahlweise 915 915 abgeschaltet 916 werden. F ür km und kh wird jetzt eine Neumann-Randbedingung am916 werden. FÃŒr km und kh wird jetzt eine Neumann-Randbedingung am 917 917 unteren 918 918 Rand verwendet. Die TKE-Energieproduktionsberechnung durch Scherung des 919 919 Grundstroms am unteren Rand wurde verbessert. </p> <p>Generelle 920 920 Verwendung der bodennahen vertikalen 921 Impulsfl üsse in921 ImpulsflÃŒsse in 922 922 den Diffusionstermen von u und v, immer wenn diese als Randbedingung 923 923 vorgegeben 924 sind (gem äß<tt>use_surface_fluxes</tt>,924 sind (gemÀà <tt>use_surface_fluxes</tt>, 925 925 - bisher wurde 926 926 dies 927 927 nur bei eingeschalteter Prandtl-Schicht gemacht). Schubspannungen 928 k önnen928 können 929 929 jetzt als Randbedingungen vorgegeben werden (allerdings unter 930 930 zuhilfenahme 931 benutzereigener Software). </p> <p>Zus ätzliche932 Schnittstellen f ür benutzereigene931 benutzereigener Software). </p> <p>ZusÀtzliche 932 Schnittstellen fÃŒr benutzereigene 933 933 Software in 934 934 allen Tendenztermen, in <tt>flow_statistics,</tt> und 935 f ür935 fÃŒr 936 936 Partikel- 937 937 bzw. dvrp-Programmteile (Bestimmung von Partikeleigenschaften, 938 938 Festlegung 939 939 von Farbtabellen). </p> <p><b>Achtung: </b>Diese 940 Änderungen können bei940 Ãnderungen können bei 941 941 Simulationen mit 942 mittlerem Wind durch die Änderungen im bereich des unteren942 mittlerem Wind durch die Ãnderungen im bereich des unteren 943 943 Randes 944 t.w. erhebliche Auswirkungen auf die oberfl ächennahen944 t.w. erhebliche Auswirkungen auf die oberflÀchennahen 945 945 Ergebnisse 946 946 haben.</p> </td> <td>advec_particles, diffusion_e, … … 958 958 <tr nosave="" valign="top"> <td nosave="">04/09/01</td> 959 959 <td>SR</td> <td>2.3b</td> <td>N</td> 960 <td>Zus ätzliche Profilausgaben für961 Energieproduktionsterme m öglich.962 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 963 963 Aufsummieren 964 964 der … … 966 966 flow_statistics, modules, pres</td> </tr> <tr nosave="" valign="top"> <td nosave="">04/09/01</td> 967 967 <td>SR</td> <td>2.3b</td> <td>B</td> 968 <td>Austausch der Geisterr änder für die968 <td>Austausch der GeisterrÀnder fÃŒr die 969 969 zeitgefilterten 970 970 Felder 971 direkt nach Durchf ührung der Filterung (war fehlerhaft seit971 direkt nach DurchfÃŒhrung der Filterung (war fehlerhaft seit 972 972 Entfernen 973 der Zeitebene t+dt (Version 2.2) und f ührte zu sehr kleinen974 St örungen975 an den seitlichen R ändern der Teilgebiete). <p><b>Achtung:</b>976 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 977 977 sich durch 978 978 diese Fehlerkorrektur.</p> </td> <td>prognostic_equations</td> 979 979 </tr> <tr nosave="" valign="top"> <td nosave="">09/11/01</td> <td>SR</td> <td>2.3c</td> 980 980 <td>N</td> <td>Farbe entlang der 981 Partikelschw änze kann sich982 ändern. Einschalten983 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> 984 984 <td>advec_particles, init_particles, modules, package_parin, 985 985 plot_dvrp</td> </tr> <tr nosave="" valign="top"> 986 986 <td nosave="">16/04/02</td> <td>SR</td> 987 987 <td>2.3d</td> <td>N</td> <td>Vorgabe 988 von Randbedingungen f ür Partikeladvektion989 m öglich.990 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 991 991 auf Datei 992 992 geschrieben 993 993 werden (Unit 85). <tt>PARTICLE</tt>-Datentyp 994 enth ält994 enthÀlt 995 995 Informationen 996 über Partikelgeschwindigkeitskomponenten und996 ÃŒber Partikelgeschwindigkeitskomponenten und 997 997 Startposition. <p>Skalartransport sowie entsprechende 998 Datenausgaben k önnen998 Datenausgaben können 999 999 durch 1000 eigene Parameter gesteuert werden (nicht mehr über die1000 eigene Parameter gesteuert werden (nicht mehr ÃŒber die 1001 1001 Feuchteparameter). </p> <p>Im Fall von 1002 1002 pdims(1)=1 (eindimensionales virtuelles … … 1014 1014 batch_scp-Script 1015 1015 ersetzt (dient der Ermittelung von Dateizyklusnummern als Information 1016 f ür1016 fÃŒr 1017 1017 die AVS-fld-Datei). </p> <p>Bei fehlerhafter 1018 Er öffnung einer CPU_MEASURES-Datei1018 Eröffnung einer CPU_MEASURES-Datei 1019 1019 werden erneute 1020 Er öffnungsversuche durchgeführt.</p> </td>1020 Eröffnungsversuche durchgefÃŒhrt.</p> </td> 1021 1021 <td>advec_particles, check_open, check_parameters, 1022 1022 exchange_horiz, flow_statistics, … … 1027 1027 transpose_zy, write_var_list</td> </tr> <tr nosave="" valign="top"> <td nosave="">02/05/02</td> 1028 1028 <td>SR</td> <td>2.3e</td> <td>B</td> 1029 <td>Wiedereinf ührung der 3D-Felder für die1029 <td>WiedereinfÃŒhrung der 3D-Felder fÃŒr die 1030 1030 Zeitebene 1031 1031 t+dt, die … … 1033 1033 mit der korrekten Zeitebene t-dt gearbeitet wird. Seit Version 2.2 1034 1034 wurde 1035 f älschlicherweise die Zeitebene t+dt verwendet. Asselin-Filter1035 fÀlschlicherweise die Zeitebene t+dt verwendet. Asselin-Filter 1036 1036 ist 1037 jetzt wieder eigenst ändiges Unterprogramm. <br>1038 Die Wiedereinf ührung dieser dritten Zeitebene ist auch1037 jetzt wieder eigenstÀndiges Unterprogramm. <br> 1038 Die WiedereinfÃŒhrung dieser dritten Zeitebene ist auch 1039 1039 Voraussetzung 1040 f ür die in der nächsten Version geplante skalare1040 fÃŒr die in der nÀchsten Version geplante skalare 1041 1041 Optimierung. <p>Fehler bei Berechnung von 1042 1042 Ausgabezeitpunkten eliminiert (trat 1043 1043 auf, wenn 1044 bei Fortsetzungsl äufen von Ausgabeabständen 0.0 auf1044 bei FortsetzungslÀufen von AusgabeabstÀnden 0.0 auf 1045 1045 von Null 1046 1046 verschiedene Werte gewechselt werden sollte). </p> <p><b>Achtung:</b> 1047 Der Inhalt der Testergebnis-Datei ändert1047 Der Inhalt der Testergebnis-Datei Àndert 1048 1048 sich durch 1049 1049 diese Fehlerkorrektur.</p> </td> <td>init_3d_model, … … 1054 1054 Anpassungen an IBM-Regatta-Systeme in 1055 1055 check_parameters, cpu_log. <p>Modul singleton in 1056 eigenst ändige Datei1057 überführt. </p> <p>Global_min_max1056 eigenstÀndige Datei 1057 ÃŒberfÃŒhrt. </p> <p>Global_min_max 1058 1058 arbeitet mit REAL*4, um Komplikationen auf 1059 1059 32-bit-Rechnern 1060 zu vermeiden (Datentyp MPI_2REAL m üsste dort sonst auf1060 zu vermeiden (Datentyp MPI_2REAL mÃŒsste dort sonst auf 1061 1061 MPI_2DOUBLE_PRECISION 1062 ge ändert werden). <br>1063 Horizontale Geschwindigkeitskomponenten f ür Partikeladvektion1062 geÀndert werden). <br> 1063 Horizontale Geschwindigkeitskomponenten fÃŒr Partikeladvektion 1064 1064 werden exakt zwischen den jeweils benachbarten horizontalen 1065 Gitterfl ächen1065 GitterflÀchen 1066 1066 interpoliert.</p> </td> <td nosave="">advec_particles, 1067 1067 check_parameters, cpu_log, … … 1070 1070 singleton</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">11/06/02</td> 1071 1071 <td>SR</td> <td>2.4</td> <td>N/C</td> 1072 <td>Optimierung und Anpassungen f ür einzelne Knoten der1072 <td>Optimierung und Anpassungen fÃŒr einzelne Knoten der 1073 1073 IBM-Regatta-Systeme. <p><b>Skalare (Cache) Optimierung:</b> 1074 1074 Tendenzterme werden 1075 1075 innerhalb einer 1076 gro ßen (i,j)-Schleife berechnet, die nun jeweils die gesamte1076 groÃen (i,j)-Schleife berechnet, die nun jeweils die gesamte 1077 1077 prognostische 1078 1078 Gleichung umfasst. Schleifenindices i und j werden an die 1079 1079 Tendenzunterprogramme 1080 als Argumente übergeben. Die Tendenzunterprogramme sind als1080 als Argumente ÃŒbergeben. Die Tendenzunterprogramme sind als 1081 1081 Module 1082 geschrieben, die mittels Technik der überladenen Funktionen1082 geschrieben, die mittels Technik der ÃŒberladenen Funktionen 1083 1083 auch 1084 1084 ohne 1085 diese Indices i und j aufgerufen werden k önnen, und dann wie1085 diese Indices i und j aufgerufen werden können, und dann wie 1086 1086 in 1087 fr üheren1087 frÃŒheren 1088 1088 Versionen funktionieren (d.h. in ihnen laufen die Schleifen 1089 über1089 ÃŒber 1090 1090 alle 1091 1091 3 Dimensionen). Solche Teile der prognostischen Gleichungen, die 1092 1092 globale 1093 1093 Kommunikation beinhalten, sind vor die (i,j)-Schleife gezogen. Falls 1094 ausschlie ßlich1094 ausschlieÃlich 1095 1095 mit Piascek-Williams-Advektion gearbeitet wird, werden alle 1096 1096 prognostischen 1097 1097 Gleichungen in einer einzigen (i,j)-Schleife gerechnet. </p> <p><b>Kommunikationsoptimierung:</b> 1098 Zus ätzliche1098 ZusÀtzliche 1099 1099 Implementierung einer 1100 1100 1D-Gebietszerlegung in x-Richtung. Mit dem neuen direkten 1101 Druckl öser1101 Drucklöser 1102 1102 (poisfft_hybrid) werden dadurch 4 von 6 notwendigen Transpositionen 1103 1103 eingespart. 1104 1104 Die 1D-Zerlegung wird automatisch bei Verwendung dieses neuen 1105 L ösers1106 aktiviert. Austausch von Geisterr ändern in y-Richtung wird bei1105 Lösers 1106 aktiviert. Austausch von GeisterrÀndern in y-Richtung wird bei 1107 1107 dieser 1108 1108 1D-Zerlegung durch einfache zyklische Randbedingungen ersetzt. </p> … … 1111 1111 IBM-Compiler 1112 1112 in der Regel den Einsatz der TRIM-Funktion, um 1113 überschüssige1113 ÃŒberschÃŒssige 1114 1114 Blanks am Ende des Strings abzuschneiden (sonst liefert der Vergleich 1115 1115 .F.). <br> 1116 1116 Lokal allokierte Felder t.w. in sogenannte automatische Felder 1117 1117 umgewandelt 1118 (g ünstiger, um Gefahr von eventuellen Speicherlecks zu1118 (gÃŒnstiger, um Gefahr von eventuellen Speicherlecks zu 1119 1119 vermeiden). 1120 1120 Einige Feldoperationen wurden wegen schlechter Performance auf der IBM 1121 1121 entfernt (s. pres). </p> <p><b>Achtung: </b>Die 1122 Übersetzung des Modells erfordert1122 Ãbersetzung des Modells erfordert 1123 1123 aufgrund 1124 1124 der Vielzahl nunmehr eingesetzter Module die Verwendung des 1125 1125 make-Mechanismus, 1126 um die korrekte Reihenfolge bei der Übersetzung der1126 um die korrekte Reihenfolge bei der Ãbersetzung der 1127 1127 Programmteile 1128 zu gew ährleisten. Die Beschreibung dieser1129 Abh ängigkeiten1128 zu gewÀhrleisten. Die Beschreibung dieser 1129 AbhÀngigkeiten 1130 1130 erfolgt 1131 1131 in der Datei<tt> Makefile</tt>, die zusammen mit den … … 1136 1136 ebenfalls 1137 1137 als Modul geschrieben und bedient sich der Methode der 1138 überladenen1138 ÃŒberladenen 1139 1139 Funktionen. Existierende benutzereigene Software muss entsprechend 1140 1140 angepasst … … 1156 1156 </tr> <tr nosave="" valign="top"> <td nosave="">11/06/02</td> <td>SR</td> <td>2.4</td> 1157 1157 <td>B</td> <td>Fehlende Variablenvereinbarung in <tt>advec_particles</tt> 1158 erg änzt.</td> <td>advec_particles</td> </tr>1158 ergÀnzt.</td> <td>advec_particles</td> </tr> 1159 1159 <tr> <td valign="top">12/09/02</td> <td valign="top">SR</td> <td valign="top">2.4a</td> 1160 1160 <td valign="top">N</td> <td valign="top">Partikel 1161 k önnen mit Trägheit versehen1161 können mit TrÀgheit versehen 1162 1162 werden und 1163 1163 einen Dichteunterschied zum umgebenden Fluid haben. Zur Steuerung dient … … 1170 1170 und q entfernt 1171 1171 (obere Feldgrenze von 10 konnte in <tt>pt_vertical_gradient_level_ind</tt> 1172 überschritten werden). Fehler führte unter1173 Umst änden zu1174 "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 1175 1175 Gitterpunktzahlen 1176 1176 in z-Richtung. <p>Fehler in Berechnung der v-Komponente am 1177 1177 unteren Rand 1178 1178 entfernt, die 1179 f ür den Scherungsproduktionsterm der TKE verwendet wird1179 fÃŒr den Scherungsproduktionsterm der TKE verwendet wird 1180 1180 (Variable 1181 v_0). F älschlicherweise wurde zur Berechnung u(k=1) statt1181 v_0). FÀlschlicherweise wurde zur Berechnung u(k=1) statt 1182 1182 v(k=1) 1183 1183 verwendet.</p> </td> <td valign="top">check_parameters, … … 2151 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 2152 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, 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-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- 2154 2154 and 3d-datasets now contain the exact coordinates of the variables 2155 2155 along x and y with respect to the positions where they are defined on … … 2702 2702 advection scheme ( <a href="Upstream_Spline/ups.ps">.ps</a> 2703 2703 )</li> <li> <a name="Cloud_physics"></a>cloud 2704 physics module ( <a href="Cloud_physics/ wolken.pdf">.pdf</a>2704 physics module ( <a href="Cloud_physics/cloud_physics.pdf">.pdf</a> 2705 2705 )</li> <li> <a name="Bitkompression"></a>Datenkompression 2706 2706 mit dem Verfahren
Note: See TracChangeset
for help on using the changeset viewer.