Index: palm/trunk/DOC/tec/technical_documentation.html
===================================================================
--- palm/trunk/DOC/tec/technical_documentation.html	(revision 54)
+++ palm/trunk/DOC/tec/technical_documentation.html	(revision 62)
@@ -1,102 +1,39 @@
 <!DOCTYPE html PUBLIC "-//w3c//dtd html 4.0 transitional//en">
 <html><head>
-
-  
-  <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&nbsp;&nbsp;
+<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&nbsp;&nbsp;
 </font>a
 <b>pa</b>rallelized
 <b>L</b>ES <b>m</b>odel
-<br>
-
-<hr width="100%">
-<br>
-
-<font size="+2">technical/numerical documentation</font>
-<br>
-
-&nbsp;
-<br>
-
-&nbsp;
-<br>
-
-<table nosave="" cellpadding="0" cellspacing="0">
-
-  <caption>&nbsp; <br>
-
-  </caption><tbody>
-
-  
-  </tbody><tbody>
-
-  
-  </tbody> <tbody>
-
-    <tr nosave="">
-
-      <td nosave=""><b>Current model version:</b></td>
-
-      <td><span style="font-weight: bold;">3.1c</span></td>
-
-    </tr>
-
-    <tr nosave="">
-
-      <td nosave=""><b>Last change of this document</b>:&nbsp;</td>
-
-      <td nosave=""><b>05/03/07</b></td>
-
-    </tr>
-
-  
-  </tbody>
-</table>
-
-<hr width="100%">
+<br><hr width="100%">
+<br><font size="+2">technical/numerical documentation</font>
+<br>&nbsp;
+<br>&nbsp;
+<br><table nosave="" cellpadding="0" cellspacing="0"> <caption>&nbsp; <br> </caption><tbody>
+</tbody><tbody> </tbody> <tbody> <tr nosave=""> <td nosave=""><b>Current
+model version:</b></td> <td><span style="font-weight: bold;">3.1c</span></td> </tr>
+<tr nosave=""> <td nosave=""><b>Last
+change of this document</b>:&nbsp;</td> <td nosave=""><b>$Id$</b></td> </tr>
+</tbody>
+</table><hr width="100%">
 <p><font size="+2">Table of contents</font>
-<br>
-
-&nbsp;
-<br>
-
-&nbsp;
-</p>
-
-<ul>
-
-  <a href="#Kapitel1.0">1.0</a>&nbsp; History of changes
-  
-  <p><a href="#Kapitel2.0">2.0</a>&nbsp; How to change the model source
-code </p>
-
-  
-  <p><a href="#Kapitel3.0">3.0</a>&nbsp; Description of selected parts
+<br>&nbsp;
+<br>&nbsp;
+</p><ul> <a href="#Kapitel1.0">1.0</a>&nbsp;
+History of changes <p><a href="#Kapitel2.0">2.0</a>&nbsp;
+How to change the model source
+code </p> <p><a href="#Kapitel3.0">3.0</a>&nbsp;
+Description of selected parts
 of
 the model source code <br>
-
 &nbsp;</p>
-
-</ul>
-
-<b>Please note: Parts of this document are still in German. We
+</ul><b>Please note: Parts of this document are still in
+German. We
 apologize
 for any inconvenience.</b>
-<br>
-
-&nbsp;
-<h2><a name="Kapitel1.0"></a>1.0&nbsp; History of changes</h2>
-
+<br>&nbsp;
+<h2><a name="Kapitel1.0"></a>1.0&nbsp;
+History of changes</h2>
 The following table documents all changes to the model since July, 7th
 1998. Entries in the fourth column specify the type of the change:
@@ -107,293 +44,110 @@
 respective
 routines.
-<br>
-
-&nbsp;
-<br>
-
-<table nosave="" style="width: 100%;" border="1">
-
-  <caption>&nbsp; <br>
-
-  </caption><tbody>
-
-  
-  </tbody><tbody>
-
-  
-  </tbody> <tbody>
-
-    <tr>
-
-      <td><b>date</b></td>
-
-      <td><b>author</b></td>
-
-      <td><b>model</b> <br>
-
-      <b>version</b></td>
-
-      <td><b>type</b></td>
-
-      <td><b>description</b></td>
-
-      <td><b>affected routines</b></td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td>06/07/98</td>
-
-      <td>SR</td>
-
-      <td>1.0</td>
-
-      <td>C</td>
-
-      <td nosave="">&Auml;nderung von Default-Werten
+<br>&nbsp;
+<br><table nosave="" style="width: 100%;" border="1"> <caption>&nbsp; <br> </caption><tbody>
+</tbody><tbody> </tbody> <tbody> <tr> <td><b>date</b></td>
+<td><b>author</b></td> <td><b>model</b>
+<br> <b>version</b></td> <td><b>type</b></td>
+<td><b>description</b></td> <td><b>affected
+routines</b></td> </tr> <tr nosave="" valign="top"> <td>06/07/98</td> <td>SR</td>
+<td>1.0</td> <td>C</td> <td nosave="">&Auml;nderung
+von Default-Werten
 (initializing_actions, end_time,
 fcl_factor, use_prior_plot1d_parameters), zus&auml;tzliche
 &Uuml;berpr&uuml;fung
-von Parametern auf Zul&auml;ssigkeit (prandtl_layer, psolver, nx, ny,
+von Parametern auf Zul&auml;ssigkeit (prandtl_layer, psolver, nx,
+ny,
 nz,
-fcl_factor), H&ouml;he der Schnittebenen der Horizontalschnitte wird
+fcl_factor), H&ouml;he der Schnittebenen der Horizontalschnitte
+wird
 auf
-Plots exakt ausgegeben (uv- bzw. w-Gitter).</td>
-
-      <td>modules, check_parameters, plot_2d</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td>"</td>
-
-      <td>SR</td>
-
-      <td>1.0</td>
-
-      <td>N</td>
-
-      <td>Neues Modul test_variables wird in (fast) allen
+Plots exakt ausgegeben (uv- bzw. w-Gitter).</td> <td>modules,
+check_parameters, plot_2d</td> </tr> <tr nosave="" valign="top"> <td>"</td> <td>SR</td>
+<td>1.0</td> <td>N</td> <td>Neues Modul
+test_variables wird in (fast) allen
 Programmteilen benutzt
 (USE). <b>Wichtig</b>: dieses Modul steht in der Datei <tt>module_test.f90.
 D</tt>iese Datei mu&szlig; in der mrun- und add_library -
 Konfigurationsdatei
-bei der Variablen <tt>module_files</tt> aufgef&uuml;hrt werden.</td>
-
-      <td nosave="">(fast) alle&nbsp; <br>
-
-Neu: module_test</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td>16/07/98</td>
-
-      <td>SR</td>
-
-      <td>1.0</td>
-
-      <td>N</td>
-
-      <td nosave="">Anpassung an die T3E-Rechner in Juelich (cpp -
-Direktiven)</td>
-
-      <td>advec_s_bc, check_open, check_parameters, cpu_zeitmessung,
+bei der Variablen <tt>module_files</tt>
+aufgef&uuml;hrt werden.</td> <td nosave="">(fast)
+alle&nbsp; <br>
+Neu: module_test</td> </tr> <tr nosave="" valign="top"> <td>16/07/98</td> <td>SR</td>
+<td>1.0</td> <td>N</td> <td nosave="">Anpassung
+an die T3E-Rechner in Juelich (cpp -
+Direktiven)</td> <td>advec_s_bc, check_open,
+check_parameters, cpu_zeitmessung,
 init_pegrid,
 local_getenv, local_system, local_tremain, local_tremain_ini, poisfft</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td>24/07/98</td>
-
-      <td>MS</td>
-
-      <td>1.0</td>
-
-      <td>C</td>
-
-      <td nosave="">Anpassung des Druckloesealgorithmusses an
+</tr> <tr nosave="" valign="top"> <td>24/07/98</td>
+<td>MS</td> <td>1.0</td> <td>C</td>
+<td nosave="">Anpassung des Druckloesealgorithmusses an
 Vorgehensweise von
-Schmidt, Schumann und Volkert (DFVLR-Mitteilung 84-15, S. 20-26)&nbsp;</td>
-
-      <td>poisfft</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td>27/07/98</td>
-
-      <td>MS</td>
-
-      <td>1.0</td>
-
-      <td>C</td>
-
-      <td nosave="">Berechnung von l in maketri geschieht nun
+Schmidt, Schumann und Volkert (DFVLR-Mitteilung 84-15, S.
+20-26)&nbsp;</td> <td>poisfft</td> </tr> <tr nosave="" valign="top"> <td>27/07/98</td>
+<td>MS</td> <td>1.0</td> <td>C</td>
+<td nosave="">Berechnung von l in maketri geschieht nun
 ausserhalb der k-Schleife,
-da unabhaengig von k, dadurch Rechenzeitersparnis&nbsp;</td>
-
-      <td>poisfft</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td>05/08/98</td>
-
-      <td>SR</td>
-
-      <td>1.0</td>
-
-      <td>C</td>
-
-      <td>Entfernung der Sonderbehandlung der Advektionsterme
+da unabhaengig von k, dadurch Rechenzeitersparnis&nbsp;</td> <td>poisfft</td>
+</tr> <tr nosave="" valign="top"> <td>05/08/98</td>
+<td>SR</td> <td>1.0</td> <td>C</td>
+<td>Entfernung der Sonderbehandlung der Advektionsterme
 (Piascek-Williams)
-bei k=1. Terme werden nun &uuml;berall gleich berechnet. Evtl. leichte
+bei k=1. Terme werden nun &uuml;berall gleich berechnet. Evtl.
+leichte
 Erh&ouml;hung der Varianzen in Oberfl&auml;chenn&auml;he.</td>
-
-      <td nosave="">advec_u_pw, advec_v_pw, advec_w_pw</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">"</td>
-
-      <td>SR</td>
-
-      <td>1.0</td>
-
-      <td>C</td>
-
-      <td>Mischungsweganpassung (s. adjust_mixing_length) erfolgt nur
+<td nosave="">advec_u_pw, advec_v_pw, advec_w_pw</td>
+</tr> <tr nosave="" valign="top"> <td nosave="">"</td> <td>SR</td> <td>1.0</td>
+<td>C</td> <td>Mischungsweganpassung (s.
+adjust_mixing_length) erfolgt nur
 noch bei
-k=1. Unrealistische sekund&auml;re Maxima in den Km-Profilen oberhalb
+k=1. Unrealistische sekund&auml;re Maxima in den Km-Profilen
+oberhalb
 der
-Prandtl-Schicht werden damit vermieden. Leichte &Auml;nderungen in den
-oberfl&auml;chennahen Turbulenzgr&ouml;&szlig;en sind die Folge.</td>
-
-      <td>diffusivities</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td>"</td>
-
-      <td>SR</td>
-
-      <td>1.0</td>
-
-      <td>B</td>
-
-      <td>Falsche Berechnung der ben&ouml;tigten CPU-Zeit pro Sekunde
+Prandtl-Schicht werden damit vermieden. Leichte &Auml;nderungen in
+den
+oberfl&auml;chennahen Turbulenzgr&ouml;&szlig;en sind die
+Folge.</td> <td>diffusivities</td> </tr> <tr nosave="" valign="top"> <td>"</td> <td>SR</td>
+<td>1.0</td> <td>B</td> <td>Falsche
+Berechnung der ben&ouml;tigten CPU-Zeit pro Sekunde
 simulierter
-Zeit bei Fortsetzungsl&auml;ufen korrigiert.&nbsp;
-      
-      <p>Berechnung von cross_uxmin(1) und cross_uxmax(1) aus den
+Zeit bei Fortsetzungsl&auml;ufen korrigiert.&nbsp; <p>Berechnung
+von cross_uxmin(1) und cross_uxmax(1) aus den
 Werten von
-ug und vg entfernt. Lieferte falsche Wertebereiche f&uuml;r Plots von
+ug und vg entfernt. Lieferte falsche Wertebereiche f&uuml;r Plots
+von
 Vertikalprofilen,
 wenn im ersten Koordinatenkreuz nicht die Windgeschwindigkeiten
 gezeichnet
-wurden.</p>
-
-      </td>
-
-      <td nosave="">header, check_parameters</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">"</td>
-
-      <td>SR</td>
-
-      <td>1.0</td>
-
-      <td>N</td>
-
-      <td>Das erlaubte Maximum und Minimum der Rif-Zahl ist nun mittels
+wurden.</p> </td> <td nosave="">header,
+check_parameters</td> </tr> <tr nosave="" valign="top"> <td nosave="">"</td> <td>SR</td>
+<td>1.0</td> <td>N</td> <td>Das
+erlaubte Maximum und Minimum der Rif-Zahl ist nun mittels
 der Initialparameter
 rif_max und rif_min einstellbar (Werte waren bisher fest auf -2.0, +1.0
-eingestellt).</td>
-
-      <td>check_parameters, header, init_1d_model, modules, parin,
+eingestellt).</td> <td>check_parameters, header,
+init_1d_model, modules, parin,
 prandtl_fluxes,
-read_var_list, write_var_list</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">23/09/98</td>
-
-      <td>SR</td>
-
-      <td>1.0</td>
-
-      <td>B</td>
-
-      <td>Initialisierung von Feldern im Fall von initializing_action =
+read_var_list, write_var_list</td> </tr> <tr nosave="" valign="top"> <td nosave="">23/09/98</td>
+<td>SR</td> <td>1.0</td> <td>B</td>
+<td>Initialisierung von Feldern im Fall von initializing_action =
 "<i>set_constant_profiles</i>"
-und "<i>set_1d-model_profiles</i>" korrigiert.&nbsp;
-      
-      <p>2D-Plotausgabe von u* und theta* korrigiert. Hier wurden
+und "<i>set_1d-model_profiles</i>" korrigiert.&nbsp; <p>2D-Plotausgabe
+von u* und theta* korrigiert. Hier wurden
 bisher falsche
-Felder herausgeschrieben.&nbsp; </p>
-
-      
-      <p>Bei SOR-Aufrufen war bisher dp nicht initialisiert.
-Entsprechend ge&auml;ndert.</p>
-
-      </td>
-
-      <td>init_3d_model, plot_2d, pres</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td>"</td>
-
-      <td>SR</td>
-
-      <td>1.0</td>
-
-      <td>C</td>
-
-      <td nosave="">Ausgabeformate f&uuml;r Statistik-Ausgaben auf
+Felder herausgeschrieben.&nbsp; </p> <p>Bei
+SOR-Aufrufen war bisher dp nicht initialisiert.
+Entsprechend ge&auml;ndert.</p> </td> <td>init_3d_model,
+plot_2d, pres</td> </tr> <tr nosave="" valign="top"> <td>"</td> <td>SR</td>
+<td>1.0</td> <td>C</td> <td nosave="">Ausgabeformate
+f&uuml;r Statistik-Ausgaben auf
 Datei RUN_CONTROL
-um einige Kommastellen erweitert.&nbsp;
-      
-      <p>Lieferung von ausf&uuml;hrlicheren Informationen bei
+um einige Kommastellen erweitert.&nbsp; <p>Lieferung von
+ausf&uuml;hrlicheren Informationen bei
 Programmabbruch
-durch Unterschreitung des minimal erlaubten Zeitschritts.</p>
-
-      </td>
-
-      <td>run_control, timestep</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">"</td>
-
-      <td>SR</td>
-
-      <td>1.0</td>
-
-      <td>N</td>
-
-      <td>Rechnungen mit geneigter Ebene zur Simulation von
+durch Unterschreitung des minimal erlaubten Zeitschritts.</p> </td>
+<td>run_control, timestep</td> </tr> <tr nosave="" valign="top"> <td nosave="">"</td>
+<td>SR</td> <td>1.0</td> <td>N</td>
+<td>Rechnungen mit geneigter Ebene zur Simulation von
 Kaltluftabfl&uuml;ssen
 sind erlaubt. Gestartet und gesteuert werden sie mit dem
@@ -401,125 +155,65 @@
 Realisiert wird die geneigte Ebene durch entsprechende
 Schr&auml;gstellung
-des Schwerkraftvektors. Dadurch erh&auml;lt auch die u-Komponente der
+des Schwerkraftvektors. Dadurch erh&auml;lt auch die u-Komponente
+der
 Windgeschwindigkeit
-einen Anteil der Auftriebskraft.&nbsp;
-      
-      <p>Diese Implementierung ist noch nicht abgeschlossen. Die
+einen Anteil der Auftriebskraft.&nbsp; <p>Diese
+Implementierung ist noch nicht abgeschlossen. Die
 Simulation von
 Kaltluftabfl&uuml;ssen erfordert bisher noch Eingriffe "von Hand".</p>
-
-      </td>
-
-      <td>buoyancy, check_parameters, header, leap_frog, modules,
+</td> <td>buoyancy, check_parameters, header, leap_frog,
+modules,
 parin, prandtl_fluxes,
-read_var_list, write_var_list</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td>15/12/98</td>
-
-      <td>SR</td>
-
-      <td>1.0</td>
-
-      <td>N</td>
-
-      <td>Vereinbarung einer Cache-Layer zur schnelleren bin&auml;ren
+read_var_list, write_var_list</td> </tr> <tr nosave="" valign="top"> <td>15/12/98</td> <td>SR</td>
+<td>1.0</td> <td>N</td> <td>Vereinbarung
+einer Cache-Layer zur schnelleren bin&auml;ren
 I/O auf
-T3E-Rechnern.&nbsp;
-      
-      <p>Ermittelung wahrer I/O-Zeiten f&uuml;r bin&auml;re
+T3E-Rechnern.&nbsp; <p>Ermittelung wahrer I/O-Zeiten
+f&uuml;r bin&auml;re
 Ein-/Ausgabe auf
-T3E-Rechnern mit TIMEF.</p>
-
-      </td>
-
-      <td nosave="">check_open,&nbsp; <br>
-
-cpu_zeitmessung, init_3d_model, write_3d_model</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">"</td>
-
-      <td>SR</td>
-
-      <td>1.0</td>
-
-      <td>B</td>
-
-      <td>Bekanntgabe von run_description_header an alle PE's.&nbsp;
-      
-      <p>Diese Programm&auml;nderung erfordert auf T3E-Rechnern
+T3E-Rechnern mit TIMEF.</p> </td> <td nosave="">check_open,&nbsp;
+<br>cpu_zeitmessung, init_3d_model, write_3d_model</td> </tr>
+<tr nosave="" valign="top"> <td nosave="">"</td>
+<td>SR</td> <td>1.0</td> <td>B</td>
+<td>Bekanntgabe von run_description_header an alle
+PE's.&nbsp; <p>Diese Programm&auml;nderung erfordert auf
+T3E-Rechnern
 zwingend den
 Einsatz der MPI-Version mpt.1.2.0.1. In dieser neuen Version tritt auch
-der alte Fehler in MPI_ALLREDUCE nicht mehr auf, so da&szlig; ab jetzt
-mpt.1.1.0.1 nicht mehr verwendet werden mu&szlig;.</p>
-
-      </td>
-
-      <td>parles</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">05/02/99</td>
-
-      <td>SR</td>
-
-      <td>1.1</td>
-
-      <td>C/N</td>
-
-      <td>Implementierung des Upstream-Spline- Verfahrens zur Advektion
+der alte Fehler in MPI_ALLREDUCE nicht mehr auf, so da&szlig; ab
+jetzt
+mpt.1.1.0.1 nicht mehr verwendet werden mu&szlig;.</p> </td>
+<td>parles</td> </tr> <tr nosave="" valign="top"> <td nosave="">05/02/99</td>
+<td>SR</td> <td>1.1</td> <td>C/N</td>
+<td>Implementierung des Upstream-Spline- Verfahrens zur Advektion
 von Impuls
-und skalaren Gr&ouml;&szlig;en (neue UP's, s. rechts). Berechnung von
+und skalaren Gr&ouml;&szlig;en (neue UP's, s. rechts).
+Berechnung von
 Variablen
-f&uuml;r spezielle Advektionsverfahren in eigenst&auml;ndiges UP
+f&uuml;r spezielle Advektionsverfahren in eigenst&auml;ndiges
+UP
 ausgelagert
-(<tt>init_advec</tt> + eigenes Modul <tt>advection</tt>). Reines
+(<tt>init_advec</tt> + eigenes Modul <tt>advection</tt>).
+Reines
 Euler-/Upstream-Verfahren
-kann gerechnet werden.&nbsp;
-      
-      <p>Schnelle bin&auml;re I/O (s. 15/12/98) wieder abgeschaltet, da
+kann gerechnet werden.&nbsp; <p>Schnelle bin&auml;re I/O
+(s. 15/12/98) wieder abgeschaltet, da
 Dateien
-nicht wieder einlesbar sind. Druckst&ouml;rung (p) wird zu Beginn auf 0
+nicht wieder einlesbar sind. Druckst&ouml;rung (p) wird zu Beginn
+auf 0
 gesetzt, damit sie keine Tendenzen im ersten Zeitschritt hervorruft.</p>
-
-      </td>
-
-      <td>advec_s_bc, check_open, check_parameters, flow_statistics,
+</td> <td>advec_s_bc, check_open, check_parameters,
+flow_statistics,
 header,
 init_3d_model, init_grid, leap_frog, modules, parin, plot_2d, plot_ts,
-read_var_list, timestep, write_3d_binary, write_var_list&nbsp;
-      
-      <p><b>Neu:</b> <br>
-
+read_var_list, timestep, write_3d_binary, write_var_list&nbsp; <p><b>Neu:</b>
+<br>
 advec_s_ups, advec_u_ups, advec_v_ups, advec_w_ups, init_advec,
 long_filter,
-spline_x, spline_y, spline_z</p>
-
-      </td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">24/02/99</td>
-
-      <td>SR</td>
-
-      <td>1.1a</td>
-
-      <td>C/N</td>
-
-      <td>Kleinere Erg&auml;nzungen / &Auml;nderungen f&uuml;r Upstream
+spline_x, spline_y, spline_z</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">24/02/99</td>
+<td>SR</td> <td>1.1a</td> <td>C/N</td>
+<td>Kleinere Erg&auml;nzungen / &Auml;nderungen
+f&uuml;r Upstream
 - Spline,&nbsp; <br>
-
 Diverse Verbesserungen bei Berechnung der Diffusionskoeffizienten
 (exaktere
@@ -529,234 +223,94 @@
 Wertebegrenzung
 bei Spline - Berechnung zur Verhinderung von "&Uuml;berschwingern"</td>
-
-      <td>check_parameters, diffusion_e, diffusion_pt, diffusivities,
+<td>check_parameters, diffusion_e, diffusion_pt, diffusivities,
 header,
 init_advec, long_filter, modules, production_e, spline_x, spline_y,
-spline_z</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td>"</td>
-
-      <td nosave="">SR</td>
-
-      <td>1.1a</td>
-
-      <td>N/B</td>
-
-      <td>Erweiterung des des Bott - Chlond - Schemas auf Advektion der
+spline_z</td> </tr> <tr nosave="" valign="top">
+<td>"</td> <td nosave="">SR</td> <td>1.1a</td>
+<td>N/B</td> <td>Erweiterung des des Bott - Chlond -
+Schemas auf Advektion der
 TKE,
 Korrektur der Dichteberechnung (war falsch bei Galilei-Transformation)</td>
-
-      <td>advec_s_bc, leap_frog</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">01/03/99</td>
-
-      <td>MS</td>
-
-      <td>1.1b</td>
-
-      <td>N</td>
-
-      <td>Kleinere Erg&auml;nzungen f&uuml;r Upstream - Spline:&nbsp;
-      
-      <ul>
-
-        <li>Gradientenkontrolle f&uuml;r den nicht-parallelen Teil</li>
-
-        <li>statistische Auswertung &uuml;ber den prozentualen Anteil
+<td>advec_s_bc, leap_frog</td> </tr> <tr nosave="" valign="top"> <td nosave="">01/03/99</td>
+<td>MS</td> <td>1.1b</td> <td>N</td>
+<td>Kleinere Erg&auml;nzungen f&uuml;r Upstream -
+Spline:&nbsp; <ul> <li>Gradientenkontrolle
+f&uuml;r den nicht-parallelen Teil</li> <li>statistische
+Auswertung &uuml;ber den prozentualen Anteil
 des
 Upstream-Verfahrens
-an der Gesamtadvektion fuer nicht-parallelen Teil</li>
-
-      
-      </ul>
-
-      </td>
-
-      <td>modules, spline_x, spline_y, spline_z</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td>"</td>
-
-      <td>SR</td>
-
-      <td>1.1b</td>
-
-      <td>N</td>
-
-      <td>Datenkompression f&uuml;r 3D-Plotausgaben</td>
-
-      <td nosave="">check_open, check_parameters, close_files, header,
+an der Gesamtadvektion fuer nicht-parallelen Teil</li> </ul>
+</td> <td>modules, spline_x, spline_y, spline_z</td>
+</tr> <tr nosave="" valign="top"> <td>"</td>
+<td>SR</td> <td>1.1b</td> <td>N</td>
+<td>Datenkompression f&uuml;r 3D-Plotausgaben</td> <td nosave="">check_open, check_parameters, close_files, header,
 modules,
-parin, plot_3d&nbsp;
-      
-      <p><b>Neu:</b> <br>
-
-write_compressed</p>
-
-      </td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">03/03/99</td>
-
-      <td>SR</td>
-
-      <td>1.1c</td>
-
-      <td>B/C</td>
-
-      <td>Schnelle bin&auml;re I/O auf T3E-Rechnern funktioniert jetzt.
+parin, plot_3d&nbsp; <p><b>Neu:</b> <br>
+write_compressed</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">03/03/99</td>
+<td>SR</td> <td>1.1c</td> <td>B/C</td>
+<td>Schnelle bin&auml;re I/O auf T3E-Rechnern funktioniert
+jetzt.
 Beim
 bin&auml;ren Einlesen d&uuml;rfen Records nicht mit dem Trick "<tt>READ
-(..) idum</tt>" &uuml;berlesen werden. Zeitmessungen auf T3E erfolgen
+(..) idum</tt>" &uuml;berlesen werden. Zeitmessungen auf T3E
+erfolgen
 jetzt
-grunds&auml;tzlich mit <tt>TIMEF.</tt></td>
-
-      <td>check_open, cpu_zeitmessung, init_3d_model, modules,
-write_3d_binary</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">25/03/99</td>
-
-      <td>SR</td>
-
-      <td>1.1d</td>
-
-      <td>N</td>
-
-      <td>Einbau einer Rayleigh-D&auml;mpfungsschicht sowie flexiblere
+grunds&auml;tzlich mit <tt>TIMEF.</tt></td> <td>check_open,
+cpu_zeitmessung, init_3d_model, modules,
+write_3d_binary</td> </tr> <tr nosave="" valign="top"> <td nosave="">25/03/99</td>
+<td>SR</td> <td>1.1d</td> <td>N</td>
+<td>Einbau einer Rayleigh-D&auml;mpfungsschicht sowie
+flexiblere
 Wahl der
 numerischen Begrenzer f&uuml;r das Upstream - Spline - Verfahren.</td>
-
-      <td>check_parameters, header, init_3d_model, leap_frog, modules,
+<td>check_parameters, header, init_3d_model, leap_frog, modules,
 parin,
-read_var_list, spline_x, spline_y, spline_z, write_var_list</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">"</td>
-
-      <td>SR</td>
-
-      <td>1.1d</td>
-
-      <td>B</td>
-
-      <td>Korrektur des Zeitschrittkriteriums bei z=zp im 1D-Modell</td>
-
-      <td>init_1d_model</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td>25/11/99</td>
-
-      <td>SR</td>
-
-      <td>1.1e</td>
-
-      <td>N</td>
-
-      <td nosave="">Partikelplots mit vtk-Graphiksoftware sind auf t3eh
+read_var_list, spline_x, spline_y, spline_z, write_var_list</td> </tr>
+<tr nosave="" valign="top"> <td nosave="">"</td>
+<td>SR</td> <td>1.1d</td> <td>B</td>
+<td>Korrektur des Zeitschrittkriteriums bei z=zp im 1D-Modell</td>
+<td>init_1d_model</td> </tr> <tr nosave="" valign="top"> <td>25/11/99</td> <td>SR</td>
+<td>1.1e</td> <td>N</td> <td nosave="">Partikelplots
+mit vtk-Graphiksoftware sind auf t3eh
 m&ouml;glich.
 Ausgabe erfolgt im dvr-Format. Bisher noch nicht nutzbar, da noch
 Einstellm&ouml;glichkeiten
-f&uuml;r die Partikelquelle(n) fehlen. <b>Wichtig:</b> Die Benutzung
+f&uuml;r die Partikelquelle(n) fehlen. <b>Wichtig:</b>
+Die Benutzung
 der
 vtk-Software erfordert die zus&auml;tzliche mrun-Option <tt>-g</tt>.
 Das
-Hauptprogramm <tt>parles</tt> wird in diesem Fall zum Unterprogramm
+Hauptprogramm <tt>parles</tt> wird in diesem Fall zum
+Unterprogramm
 und
-von einem c++-Programm aus aufgerufen.</td>
-
-      <td>header, init_3d_model, init_pegrid, leap_frog, modules,
+von einem c++-Programm aus aufgerufen.</td> <td>header,
+init_3d_model, init_pegrid, leap_frog, modules,
 parin, parles,
-read_var_list, write_var_list&nbsp;
-      
-      <p><b>Neu:</b> <br>
-
-advec_particles, init_particles, plot_particles</p>
-
-      </td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td>"</td>
-
-      <td>SR</td>
-
-      <td>1.1e</td>
-
-      <td nosave="">B/C</td>
-
-      <td>TIMEF-Zeitmessungen von Subroutine- auf Function-Aufruf
+read_var_list, write_var_list&nbsp; <p><b>Neu:</b>
+<br>
+advec_particles, init_particles, plot_particles</p> </td> </tr>
+<tr nosave="" valign="top"> <td>"</td>
+<td>SR</td> <td>1.1e</td> <td nosave="">B/C</td>
+<td>TIMEF-Zeitmessungen von Subroutine- auf Function-Aufruf
 umgestellt.
-Laufindexkorrektur in <tt>spline_y</tt>. Fehler trat bei nx /= ny auf.</td>
-
-      <td>cpu_zeitmessung, spline_y</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">28/12/99</td>
-
-      <td>SR</td>
-
-      <td>1.1f</td>
-
-      <td>N</td>
-
-      <td>Plot von Isooberfl&auml;chen mit vtk-Graphiksoftware sind auf
+Laufindexkorrektur in <tt>spline_y</tt>. Fehler trat bei
+nx /= ny auf.</td> <td>cpu_zeitmessung, spline_y</td>
+</tr> <tr nosave="" valign="top"> <td nosave="">28/12/99</td> <td>SR</td> <td>1.1f</td>
+<td>N</td> <td>Plot von Isooberfl&auml;chen mit
+vtk-Graphiksoftware sind auf
 t3eh
 m&ouml;glich. Ausgabe erfolgt im dvr-Format. Bisher noch nicht
 allgemein
 nutzbar, da noch diverse manuelle Eingriffe n&ouml;tig sind.</td>
-
-      <td>header, leap_frog, modules, parin, plot_particles,
-read_var_list, write_var_list&nbsp;
-      
-      <p><b>Neu:</b> <br>
-
-plot_isosurface</p>
-
-      </td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">10/01/00</td>
-
-      <td>SR</td>
-
-      <td>1.2</td>
-
-      <td>C</td>
-
-      <td>Randbedingungen f&uuml;r w (=0) werden explizit gesetzt, und
+<td>header, leap_frog, modules, parin, plot_particles,
+read_var_list, write_var_list&nbsp; <p><b>Neu:</b>
+<br>
+plot_isosurface</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">10/01/00</td>
+<td>SR</td> <td>1.2</td> <td>C</td>
+<td>Randbedingungen f&uuml;r w (=0) werden explizit gesetzt,
+und
 zwar nicht
-nur bei&nbsp; <tt>k = nzb, nzt+1</tt> sondern auch bei&nbsp; <tt>k =
+nur bei&nbsp; <tt>k = nzb, nzt+1</tt> sondern auch
+bei&nbsp; <tt>k =
 nzt</tt>,
 wo bislang der vom vorigen Zeitschritt durch den Druckl&ouml;ser
@@ -764,16 +318,16 @@
 Wert stand. Auch jetzt wird der Wert bei k = nzt anschlie&szlig;end
 durch
-den Druckl&ouml;ser modifiziert. Soweit bis jetzt abzusehen, hat diese
-&Auml;nderung nur minimale Auswirkungen auf den Prognoseverlauf. Evtl.
+den Druckl&ouml;ser modifiziert. Soweit bis jetzt abzusehen, hat
+diese
+&Auml;nderung nur minimale Auswirkungen auf den Prognoseverlauf.
+Evtl.
 gibt es st&auml;rkere Auswirkungen bei exzessivem Auftreten von
 Schwerewellen
 am oberen Rand (die dort aber eigentlich sowieso nicht
 hingeh&ouml;ren).&nbsp; <br>
-
-Variablen&uuml;bergabe an <tt>boundary_conds</tt> erfolgt&nbsp; <br>
-
-&uuml;ber Modul und nicht mehr per Parameterliste.&nbsp;
-      
-      <p>Bei Galilei-Transformation wird jetzt defaultm&auml;&szlig;ig
+Variablen&uuml;bergabe an <tt>boundary_conds</tt>
+erfolgt&nbsp; <br>
+&uuml;ber Modul und nicht mehr per Parameterliste.&nbsp; <p>Bei
+Galilei-Transformation wird jetzt defaultm&auml;&szlig;ig
 der geostrophische
 Wind als Translationsgeschwindigkeit verwendet (neuer
@@ -782,7 +336,5 @@
 h&ouml;heren
 geostrophischen Windgeschwindigkeiten unterdr&uuml;ckt.&nbsp; </p>
-
-      
-      <p>3D-Feldvariablen (u, v, w, pt, e, km, kh) sowie einige
+<p>3D-Feldvariablen (u, v, w, pt, e, km, kh) sowie einige
 2D-Variablen
 (usws, vsws, rif, shf) sind nun als Pointer deklariert, denen nach
@@ -792,94 +344,56 @@
 3 Zeitebenen existieren die Targets ..._1 (z.B. u_1), ..._2, ..._3).
 Durch
-die jeweilige Neuzuordnung der Targets entf&auml;llt das Umspeichern
+die jeweilige Neuzuordnung der Targets entf&auml;llt das
+Umspeichern
 der
-Zeitebenen komplett (Zeitersparung: ca. 5% der Gesamtrechenzeit).&nbsp;
-      <br>
-
+Zeitebenen komplett (Zeitersparung: ca. 5% der
+Gesamtrechenzeit).&nbsp; <br>
 Die Umordnung der Zeitebenen wurde ins neue UP <tt>swap_timelevel</tt>
-ausgelagert.&nbsp; <br>
-
-      <b>Wichtig:</b> Die &Uuml;bersetzungszeiten haben durch die
+ausgelagert.&nbsp; <br> <b>Wichtig:</b> Die
+&Uuml;bersetzungszeiten haben durch die
 Einf&uuml;hrung
 der Pointer t.w. sehr deutlich zugenommen (<tt>leap_frog</tt>
 ben&ouml;tigt
-jetzt z.B. ca. 60 Sekunden). Au&szlig;erdem ben&ouml;tigt der Compiler
+jetzt z.B. ca. 60 Sekunden). Au&szlig;erdem ben&ouml;tigt der
+Compiler
 jetzt wesentlich mehr Speicher. Der "memory"-Parameter in den
 Konfigurationsdateien
 (.mrun.config, .al.config) muss nun mindestens den Wert 60 (MB) (t3eh,
-t3ej2) bzw. 45 (t3eb) und 96 (vpp) haben.&nbsp; <br>
-
-      <b>Achtung:</b> bislang wurde diese &Auml;nderung nur f&uuml;r
+t3ej2) bzw. 45 (t3eb) und 96 (vpp) haben.&nbsp; <br> <b>Achtung:</b>
+bislang wurde diese &Auml;nderung nur f&uuml;r
 das
-reine Leapfrog-Verfahren getestet!&nbsp; </p>
-
-      
-      <p>Neue Laderoption "<tt>-D preset=nan</tt>" empfohlen (s.
+reine Leapfrog-Verfahren getestet!&nbsp; </p> <p>Neue
+Laderoption "<tt>-D preset=nan</tt>" empfohlen (s.
 Beispiel-Konfigurationsdatei).
-Ebenso empfehlenswert f&uuml;r Debug-L&auml;ufe: Compileroption "<tt>-ei</tt>".
+Ebenso empfehlenswert f&uuml;r Debug-L&auml;ufe: Compileroption
+"<tt>-ei</tt>".
 Diese bewirkt einen Programmabbruch, falls mit nicht initialisierten
 Variablen
 gearbeitet wird (ist aber zeitaufwendig in der Ausf&uuml;hrung!).</p>
-
-      </td>
-
-      <td>boundary_conds, check_parameters, header, init_3d_model,
+</td> <td>boundary_conds, check_parameters, header,
+init_3d_model,
 leap_frog,
-modules, parin, read_var_list, timestep, write_var_list&nbsp;
-      
-      <p><b>Neu:</b> <br>
-
-swap_timelevel</p>
-
-      </td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">10/01/00</td>
-
-      <td>SR</td>
-
-      <td>1.2</td>
-
-      <td>B</td>
-
-      <td>Der Asselin-Filter wird vor Umordnung der Zeitebenen
+modules, parin, read_var_list, timestep, write_var_list&nbsp; <p><b>Neu:</b>
+<br>
+swap_timelevel</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">10/01/00</td>
+<td>SR</td> <td>1.2</td> <td>B</td>
+<td>Der Asselin-Filter wird vor Umordnung der Zeitebenen
 aufgerufen. Dies
 geschah bisher f&auml;lschlicherweise erst danach, sodass der
 Zeitfilter
-asymmetrisch arbeitete.&nbsp; <br>
-
-      <b>Achtung:</b> diese &Auml;nderung hat Auswirkungen auf den
+asymmetrisch arbeitete.&nbsp; <br> <b>Achtung:</b>
+diese &Auml;nderung hat Auswirkungen auf den
 Prognoseverlauf.
 Die Testergebnis-Datei wurde entsprechend ge&auml;ndert.</td>
-
-      <td>leap_frog</td>
-
-    </tr>
-
-    <tr nosave="">
-
-      <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> <br>
-
+<td>leap_frog</td> </tr> <tr nosave="">
+<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>
+<br>
 Der Asselin-Filter im 1D-Modell wird nun auch vor der Umordnung der
-Zeitebenen aufgerufen (s.o.).&nbsp;
-      
-      <p><b>3D-Modell:</b> <br>
-
-TKE wurde bisher bei <tt>scalar_advec=bc-scheme</tt> einer zeitlichen
-Filterung unterzogen. Dies wird nun unterbunden.&nbsp; </p>
-
-      
-      <p>Die Divergenz der Impulsfl&uuml;sse an der Prandtl - Schicht -
+Zeitebenen aufgerufen (s.o.).&nbsp; <p><b>3D-Modell:</b>
+<br>
+TKE wurde bisher bei <tt>scalar_advec=bc-scheme</tt> einer
+zeitlichen
+Filterung unterzogen. Dies wird nun unterbunden.&nbsp; </p> <p>Die
+Divergenz der Impulsfl&uuml;sse an der Prandtl - Schicht -
 Obergrenze
 wird nun &uuml;ber eine ganze Gitterweite gebildet. Die bisherige
@@ -887,8 +401,6 @@
 lieferte im Vergleich mit Ergebnissen anderer LES-Modelle&nbsp; zu
 gro&szlig;e
-Impulsfl&uuml;sse.&nbsp; </p>
-
-      
-      <p>Bei der Berechnung der Diffusion der Temperatur werden die
+Impulsfl&uuml;sse.&nbsp; </p> <p>Bei der
+Berechnung der Diffusion der Temperatur werden die
 Diffusionskoeffizienten
 nun wieder durch Mittelung f&uuml;r die entsprechenden&nbsp;
@@ -898,348 +410,143 @@
 Upstream-Spline Advektionsverfahrens eingebaut wurde, bewirkte ein zu
 rasches
-Anwachsen der Grenzschicht mit der Zeit.&nbsp; </p>
-
-      
-      <p><b>Achtung:</b> diese &Auml;nderungen haben Auswirkungen auf
+Anwachsen der Grenzschicht mit der Zeit.&nbsp; </p> <p><b>Achtung:</b>
+diese &Auml;nderungen haben Auswirkungen auf
 den Prognoseverlauf.
-Die Testergebnis-Datei wurde entsprechend ge&auml;ndert.</p>
-
-      </td>
-
-      <td nosave="" valign="top">init_1d_model, asselin_filter,
+Die Testergebnis-Datei wurde entsprechend ge&auml;ndert.</p> </td>
+<td nosave="" valign="top">init_1d_model,
+asselin_filter,
 diffusion_u, diffusion_v,&nbsp; <br>
-
-diffusion_pt</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">17/02/00</td>
-
-      <td>MS</td>
-
-      <td>2.0</td>
-
-      <td>N</td>
-
-      <td><b>Allgemeine Beschreibung der &Auml;nderungen:</b> <br>
-
-vollst&auml;ndige Implementierung des wolkenphysikalischen Moduls&nbsp;
-      
-      <p>Zur Implementierung der Wolkenphysik waren umfangreiche
+diffusion_pt</td> </tr> <tr nosave="" valign="top"> <td nosave="">17/02/00</td>
+<td>MS</td> <td>2.0</td> <td>N</td>
+<td><b>Allgemeine Beschreibung der &Auml;nderungen:</b>
+<br>
+vollst&auml;ndige Implementierung des wolkenphysikalischen
+Moduls&nbsp; <p>Zur Implementierung der Wolkenphysik waren
+umfangreiche
 &Auml;nderungen
 in zahlreichen Programmteilen notwendig (Beschreibung <a href="#Cloud_physics">siehe
 unten</a> ). Im einzelnen wurden folgende Unprogramme
-ge&auml;ndert:&nbsp; <br>
-
-      <b>modules</b>:&nbsp; <br>
-
-1) Bereitstellung neuer 3D-, 2D- und 1D-Felder f&uuml;r die Prognose
+ge&auml;ndert:&nbsp; <br> <b>modules</b>:&nbsp;
+<br>1) Bereitstellung neuer 3D-, 2D- und 1D-Felder f&uuml;r
+die Prognose
 des Gesamtwassergehaltes&nbsp; <br>
-
 2) "cloud_parameters" als neues Modul mit wolkenphysikalischen
-Variablen&nbsp; <br>
-
-      <b>parin: </b>Initialisierungsparameterliste wurde erweitert
-(s.u.)&nbsp; <br>
-
-      <font color="#000000"><b>read_var_list, write_var_list,
+Variablen&nbsp; <br> <b>parin: </b>Initialisierungsparameterliste
+wurde erweitert
+(s.u.)&nbsp; <br> <font color="#000000"><b>read_var_list,
+write_var_list,
 write_3d_binary: </b>In/Output
 der neuen Initialisierungsparamter und Felder f&uuml;r
-Fortsetzungsl&auml;ufe&nbsp;</font> <br>
-
-      <font color="#000000"><b>check_parameters:</b> Initialisierung
+Fortsetzungsl&auml;ufe&nbsp;</font> <br> <font color="#000000"><b>check_parameters:</b>
+Initialisierung
 des
-Feuchteprofils und Ausgabe von neuen Vertikalprofilen:&nbsp;</font> <br>
-
-      <font color="#ff0000">vpt</font><font color="#000000">: virtuelle
+Feuchteprofils und Ausgabe von neuen Vertikalprofilen:&nbsp;</font>
+<br> <font color="#ff0000">vpt</font><font color="#000000">: virtuelle
 poteniellen
-Temperatur,</font> <br>
-
-      <font color="#ff0000">lpt</font><font color="#000000">:
+Temperatur,</font> <br> <font color="#ff0000">lpt</font><font color="#000000">:
 potentielle
-Fl&uuml;ssigwassertemperatur,</font> <br>
-
-      <font color="#ff0000">q</font><font color="#000000">:
-Gesamtwassergehalt,</font> <br>
-
-      <font color="#ff0000">qv</font><font color="#000000">:
+Fl&uuml;ssigwassertemperatur,</font> <br> <font color="#ff0000">q</font><font color="#000000">:
+Gesamtwassergehalt,</font> <br> <font color="#ff0000">qv</font><font color="#000000">:
 spezifische
-Feuchte,</font> <br>
-
-      <font color="#ff0000">ql</font><font color="#000000">:
-Fl&uuml;ssigwassergehalt,</font> <br>
-
-      <font color="#ff0000">w"vpt", w*vpt*, wvpt</font><font color="#000000">:
+Feuchte,</font> <br> <font color="#ff0000">ql</font><font color="#000000">:
+Fl&uuml;ssigwassergehalt,</font> <br> <font color="#ff0000">w"vpt", w*vpt*, wvpt</font><font color="#000000">:
 Auftriebsfluss (subskalig, aufgel&ouml;st, gesamt),</font> <br>
-
-      <font color="#ff0000">w"q", w*q*, wq</font><font color="#000000">:
-Wasserfluss (subs., aufg., ges.),</font> <br>
-
-      <font color="#ff0000">w"qv", w*qv*, wqv</font><font color="#000000">:
-Feuchtefluss (subs, aufg., ges.)</font> <br>
-
-      <b>init_3d_model</b>: Initialisierung aller
-Feuchte/Wasserfelder&nbsp; <br>
-
-      <b>init_1d_model</b>: Anpassung an Rechnungen mit Feuchte&nbsp; <br>
-
-      <b>header</b>: Headerausgaben f&uuml;r Feuchtebehandlung und
+<font color="#ff0000">w"q", w*q*, wq</font><font color="#000000">:
+Wasserfluss (subs., aufg., ges.),</font> <br> <font color="#ff0000">w"qv", w*qv*, wqv</font><font color="#000000">:
+Feuchtefluss (subs, aufg., ges.)</font> <br> <b>init_3d_model</b>:
+Initialisierung aller
+Feuchte/Wasserfelder&nbsp; <br> <b>init_1d_model</b>:
+Anpassung an Rechnungen mit Feuchte&nbsp; <br> <b>header</b>:
+Headerausgaben f&uuml;r Feuchtebehandlung und
 Wolkenphysik
 eingef&uuml;gt, au&szlig;erdem Formatnummern erweitert um Platz
 f&uuml;r
-weitere Ausgaben zu schaffen&nbsp; <br>
-
-      <b>leap_frog:</b> <br>
-
+weitere Ausgaben zu schaffen&nbsp; <br> <b>leap_frog:</b>
+<br>
 1) alle prognostischen Gleichungen wurden in prognostic_equations
 ausgelagert&nbsp; <br>
-
 2) Aufrufe der neuen Unterprogramme (UP) compute_vpt,
 calc_liquid_water_content&nbsp; <br>
-
-3) veraenderter Aufruf von diffusivities&nbsp; <br>
-
-      <b><tt>b</tt>oundary_conds, asselin_filter, swap_timelevel,
+3) veraenderter Aufruf von diffusivities&nbsp; <br> <b><tt>b</tt>oundary_conds,
+asselin_filter, swap_timelevel,
 flow_statistics: </b>Anpassung
-an Feuchtebehandlung&nbsp; <br>
-
-      <b>plot2d, plot3d</b>: Ausgabem&ouml;glichkeiten von q und ql <br>
-
-      <b>buoyancy</b>: Anpassung des Autriebstermes an
+an Feuchtebehandlung&nbsp; <br> <b>plot2d, plot3d</b>:
+Ausgabem&ouml;glichkeiten von q und ql <br> <b>buoyancy</b>:
+Anpassung des Autriebstermes an
 Feuchtebehandlung
 durch neue &Uuml;bergabeparameter: Temperatur + Nummer des
 entsprechenden
-mittleren Temperaturprofils&nbsp; <br>
-
-      <b>diffusion_e</b>: Anpassung der Mischungswegberechnung an
+mittleren Temperaturprofils&nbsp; <br> <b>diffusion_e</b>:
+Anpassung der Mischungswegberechnung an
 Feuchtebehandlung
-durch Variation der &Uuml;bergabeparameter&nbsp; <br>
-
-      <b>production_e</b>: TKE-Produktion durch Auftriebskr&auml;fte an
+durch Variation der &Uuml;bergabeparameter&nbsp; <br> <b>production_e</b>:
+TKE-Produktion durch Auftriebskr&auml;fte an
 Feuchtebehandlung
-angepasst&nbsp; <br>
-
-      <b>diffusivities</b>: siehe diffusion_e&nbsp; <br>
-
-      <b>prandtl-fluxes</b>: Prandtl-Schicht fuer den
+angepasst&nbsp; <br> <b>diffusivities</b>:
+siehe diffusion_e&nbsp; <br> <b>prandtl-fluxes</b>:
+Prandtl-Schicht fuer den
 Gesamtwassergehalt
-eingef&uuml;gt&nbsp; </p>
-
-      
-      <p>Zus&auml;tzlich wurden folgende Unterprogramme NEU
-eingef&uuml;hrt:&nbsp; <br>
-
-      <b>inti_cloud_physics</b>: Initialisierung wolkenphysikalischer
-Parameter&nbsp; <br>
-
-      <b>prognostic_equations</b>: Auslagerung der prognostischen
+eingef&uuml;gt&nbsp; </p> <p>Zus&auml;tzlich
+wurden folgende Unterprogramme NEU
+eingef&uuml;hrt:&nbsp; <br> <b>inti_cloud_physics</b>:
+Initialisierung wolkenphysikalischer
+Parameter&nbsp; <br> <b>prognostic_equations</b>:
+Auslagerung der prognostischen
 Gleichungen
-aus leap_frog, zus&auml;tzliche porgnostische Gleichung f&uuml;r den
-Gesamtwassergehalt&nbsp; <br>
-
-      <b>diffusion_s</b>: Berechnung der Diffusion skalarer
+aus leap_frog, zus&auml;tzliche porgnostische Gleichung
+f&uuml;r den
+Gesamtwassergehalt&nbsp; <br> <b>diffusion_s</b>:
+Berechnung der Diffusion skalarer
 Gr&ouml;&szlig;en
 (Temperatur und Gesamtwassergehalt), ersetzt diffusion_pt&nbsp; <br>
-
-      <b>calc_liquid_water_content</b>: Berechnung des
-Fl&uuml;ssigwassergehalts&nbsp; <br>
-
-      <b>calc_radiation</b>: Parametrisierung langwelliger
+<b>calc_liquid_water_content</b>: Berechnung des
+Fl&uuml;ssigwassergehalts&nbsp; <br> <b>calc_radiation</b>:
+Parametrisierung langwelliger
 Strahlungsprozesse
 &uuml;ber des Schema der effektiven Emissivit&auml;t&nbsp; <br>
-
-      <b>calc_precipitation</b>:Parametrisierung von
+<b>calc_precipitation</b>:Parametrisierung von
 Niederschlagsprozessen
-mit Hilfe eines vereinfachten Kesslerschemas&nbsp; <br>
-
-      <b>impact_of_latent_heat</b>: Ber&uuml;cksichtigung der durch
+mit Hilfe eines vereinfachten Kesslerschemas&nbsp; <br> <b>impact_of_latent_heat</b>:
+Ber&uuml;cksichtigung der durch
 Niederschlagsprozesse
 zur&uuml;ckbleibenden latenten Waerme innerhalb eines
-Gittervolumens&nbsp; <br>
-
-      <b>compute_vpt</b>: Berechnung der virtuellen potentiellen
-Temperatur&nbsp; </p>
-
-      
-      <p><b>Liste der neuen Initialisierungsparameter:</b> <br>
-
-      
-      <table nosave="" border="1">
-
-        <caption>&nbsp; <br>
-
-        </caption><tbody>
-
-        
-        </tbody><tbody>
-
-        
-        </tbody> <tbody>
-
-          <tr nosave="">
-
-            <td nosave=""><b>Parameter</b></td>
-
-            <td nosave=""><b>Typ</b></td>
-
-            <td nosave="" width="4%"><b>Default</b></td>
-
-            <td nosave=""><b>Beschreibung</b></td>
-
-          </tr>
-
-          <tr nosave="">
-
-            <td nosave=""><tt>moisture</tt></td>
-
-            <td>L</td>
-
-            <td>FALSE</td>
-
-            <td>Ab-/Zuschalten des Gesamtwassergerhaltes als neue
-progn. Variable</td>
-
-          </tr>
-
-          <tr>
-
-            <td><tt>cloud_physics</tt></td>
-
-            <td>L</td>
-
-            <td>FALSE</td>
-
-            <td>Ab-/Zuschalten des Kondensationsschemas</td>
-
-          </tr>
-
-          <tr nosave="">
-
-            <td><tt>radiation</tt></td>
-
-            <td>L</td>
-
-            <td nosave="">FALSE</td>
-
-            <td>Ab-/Zuschalten des Strahlungsschemas</td>
-
-          </tr>
-
-          <tr>
-
-            <td><tt>precipitation</tt></td>
-
-            <td>L</td>
-
-            <td>FALSE</td>
-
-            <td>Ab-/Zuschalten der Niderschlagsparametrisierung</td>
-
-          </tr>
-
-          <tr>
-
-            <td><tt>bc_q_b</tt></td>
-
-            <td>C</td>
-
-            <td>'dirichlet'</td>
-
-            <td>untere Randbedingung f&uuml;r q (siehe auch bc_pt_b)</td>
-
-          </tr>
-
-          <tr>
-
-            <td><tt>bc_q_t</tt></td>
-
-            <td>C</td>
-
-            <td>'neumann'</td>
-
-            <td>obere RB f&uuml;r q (siehe auch bc_pt_t)</td>
-
-          </tr>
-
-          <tr>
-
-            <td><tt>q_surface</tt></td>
-
-            <td>R</td>
-
-            <td>0.0</td>
-
-            <td>Feuchtewert an der Erdoberfl&auml;che</td>
-
-          </tr>
-
-          <tr>
-
-            <td><tt>q_surface_initial_change</tt></td>
-
-            <td>R</td>
-
-            <td>0.0</td>
-
-            <td>vgl. pt_surface_initial_change</td>
-
-          </tr>
-
-          <tr>
-
-            <td><tt>q_vertical_gradient</tt></td>
-
-            <td>R(10)</td>
-
-            <td>0.0</td>
-
-            <td>vgl. pt_vertical_gradient</td>
-
-          </tr>
-
-          <tr nosave="">
-
-            <td nosave=""><tt>q_vertical_gradient_level</tt></td>
-
-            <td>R(10)</td>
-
-            <td>10000.0</td>
-
-            <td>vgl. pt_vertical_gradient_level</td>
-
-          </tr>
-
-          <tr>
-
-            <td><tt>surface_pressure</tt></td>
-
-            <td>R</td>
-
-            <td>1013.25</td>
-
-            <td>Luftdruck an der Erdoberfl&auml;che</td>
-
-          </tr>
-
-          <tr>
-
-            <td><tt>surface_waterflux</tt></td>
-
-            <td>R</td>
-
-            <td>0.0</td>
-
-            <td>oberfl&auml;chennaher Wasser/Feuchtefluss</td>
-
-          </tr>
-
-        
-        </tbody>
-      
-      </table>
-
-Fortsetzungsl&auml;ufe mit Daten, die mit Modellversionen kleiner 2.0
+Gittervolumens&nbsp; <br> <b>compute_vpt</b>:
+Berechnung der virtuellen potentiellen
+Temperatur&nbsp; </p> <p><b>Liste der neuen
+Initialisierungsparameter:</b> <br> <table nosave="" border="1"> <caption>&nbsp; <br> </caption><tbody>
+</tbody><tbody> </tbody> <tbody> <tr nosave=""> <td nosave=""><b>Parameter</b></td>
+<td nosave=""><b>Typ</b></td> <td nosave="" width="4%"><b>Default</b></td>
+<td nosave=""><b>Beschreibung</b></td> </tr>
+<tr nosave=""> <td nosave=""><tt>moisture</tt></td>
+<td>L</td> <td>FALSE</td> <td>Ab-/Zuschalten
+des Gesamtwassergerhaltes als neue
+progn. Variable</td> </tr> <tr> <td><tt>cloud_physics</tt></td>
+<td>L</td> <td>FALSE</td> <td>Ab-/Zuschalten
+des Kondensationsschemas</td> </tr> <tr nosave="">
+<td><tt>radiation</tt></td> <td>L</td>
+<td nosave="">FALSE</td> <td>Ab-/Zuschalten
+des Strahlungsschemas</td> </tr> <tr> <td><tt>precipitation</tt></td>
+<td>L</td> <td>FALSE</td> <td>Ab-/Zuschalten
+der Niderschlagsparametrisierung</td> </tr> <tr> <td><tt>bc_q_b</tt></td>
+<td>C</td> <td>'dirichlet'</td> <td>untere
+Randbedingung f&uuml;r q (siehe auch bc_pt_b)</td> </tr>
+<tr> <td><tt>bc_q_t</tt></td> <td>C</td>
+<td>'neumann'</td> <td>obere RB f&uuml;r q
+(siehe auch bc_pt_t)</td> </tr> <tr> <td><tt>q_surface</tt></td>
+<td>R</td> <td>0.0</td> <td>Feuchtewert
+an der Erdoberfl&auml;che</td> </tr> <tr> <td><tt>q_surface_initial_change</tt></td>
+<td>R</td> <td>0.0</td> <td>vgl.
+pt_surface_initial_change</td> </tr> <tr> <td><tt>q_vertical_gradient</tt></td>
+<td>R(10)</td> <td>0.0</td> <td>vgl.
+pt_vertical_gradient</td> </tr> <tr nosave="">
+<td nosave=""><tt>q_vertical_gradient_level</tt></td>
+<td>R(10)</td> <td>10000.0</td> <td>vgl.
+pt_vertical_gradient_level</td> </tr> <tr> <td><tt>surface_pressure</tt></td>
+<td>R</td> <td>1013.25</td> <td>Luftdruck
+an der Erdoberfl&auml;che</td> </tr> <tr> <td><tt>surface_waterflux</tt></td>
+<td>R</td> <td>0.0</td> <td>oberfl&auml;chennaher
+Wasser/Feuchtefluss</td> </tr> </tbody> </table>
+Fortsetzungsl&auml;ufe mit Daten, die mit Modellversionen kleiner
+2.0
 erzeugt
 wurden, sind aufgrund der erweiterten Initialisierungsparameterliste
@@ -1247,14 +554,11 @@
 zu den Einleseroutinen der Version 2.0. Ein Update auf Version 2.0
 sollte
-deshalb nach Abschluss eines kompletten Modelllaufs erfolgen.&nbsp; </p>
-
-      
-      <p><b>Achtung:</b> diese &Auml;nderungen haben Auswirkungen auf
+deshalb nach Abschluss eines kompletten Modelllaufs erfolgen.&nbsp;
+</p> <p><b>Achtung:</b> diese
+&Auml;nderungen haben Auswirkungen auf
 den Prognoseverlauf.
-Die Testergebnis-Datei wurde entsprechend ge&auml;ndert.</p>
-
-      </td>
-
-      <td nosave="" valign="top">modules, parin, read_var_list,
+Die Testergebnis-Datei wurde entsprechend ge&auml;ndert.</p> </td>
+<td nosave="" valign="top">modules, parin,
+read_var_list,
 check_parameters,
 init_3d_model, init_1d_model, header, leap_frog. buoyancy, diffusion_e,
@@ -1262,320 +566,159 @@
 asselin_filter,
 swap_timelevels, flow_statistics, plot_2d, plot_3d, write_var_list,
-write_3d_binary&nbsp;
-      
-      <p><b>Neu:</b> <br>
-
+write_3d_binary&nbsp; <p><b>Neu:</b> <br>
 init_cloud_physics, prognostic_equations, diffusion_s,
 calc_liquid_water_content,
 calc_radiation, calc_precipitation, impact_of_latent_heat, comute_vpt</p>
-
-      </td>
-
-    </tr>
-
-    <tr nosave="">
-
-      <td nosave="" valign="top">18/02/00</td>
-
-      <td nosave="" valign="top">MS</td>
-
-      <td nosave="" valign="top">2.0</td>
-
-      <td nosave="" valign="top">B/N</td>
-
-      <td>Im Zuge der Umstellung auf Modellversion 2.0&nbsp; wurden
+</td> </tr> <tr nosave=""> <td nosave="" valign="top">18/02/00</td> <td nosave="" valign="top">MS</td> <td nosave="" valign="top">2.0</td> <td nosave="" valign="top">B/N</td> <td>Im Zuge der
+Umstellung auf Modellversion 2.0&nbsp; wurden
 noch kleine
-&Auml;nderungen an den Unterprogrammen diffusion_e, production_e und
+&Auml;nderungen an den Unterprogrammen diffusion_e, production_e
+und
 diffusivities
 vorgenommen. In allen genannten&nbsp; <br>
-
 Unterprogrammen wurde bisher ein strengeres Kriterium zur Bestimmung
-der Schichtungsverh&auml;ltnisse verwendet. Ein Modellvergleich zeigte
+der Schichtungsverh&auml;ltnisse verwendet. Ein Modellvergleich
+zeigte
 jedoch, dass dadurch die Diffusion im Bereich der Inversion zu gross
-wird.&nbsp;
-      
-      <p>Eine weitere Erg&auml;nzung betrifft die
+wird.&nbsp; <p>Eine weitere Erg&auml;nzung betrifft die
 Advektionsverfahren:&nbsp; <br>
-
-Das Bott-Chlond Advektionsverfahren ist nun auch f&uuml;r Rechnungen
-mit Feuchte/Wolkenphysik verf&uuml;gbar.&nbsp; </p>
-
-      
-      <p><b>Achtung:</b> diese &Auml;nderungen haben Auswirkungen auf
+Das Bott-Chlond Advektionsverfahren ist nun auch f&uuml;r
+Rechnungen
+mit Feuchte/Wolkenphysik verf&uuml;gbar.&nbsp; </p> <p><b>Achtung:</b>
+diese &Auml;nderungen haben Auswirkungen auf
 den Prognoseverlauf.
-Die Testergebnis-Datei wurde entsprechend ge&auml;ndert.</p>
-
-      </td>
-
-      <td nosave="" valign="top">diffusion_e, production_e,
-diffusivities, check_parameters</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">26/04/00</td>
-
-      <td>SR</td>
-
-      <td>2.0a</td>
-
-      <td>C</td>
-
-      <td>Vollst&auml;ndige Umstellung von vtk-Grafik-Software auf
+Die Testergebnis-Datei wurde entsprechend ge&auml;ndert.</p> </td>
+<td nosave="" valign="top">diffusion_e,
+production_e,
+diffusivities, check_parameters</td> </tr> <tr nosave="" valign="top"> <td nosave="">26/04/00</td>
+<td>SR</td> <td>2.0a</td> <td>C</td>
+<td>Vollst&auml;ndige Umstellung von vtk-Grafik-Software auf
 dvrp-Software.
 Grafik-Ausgaben im dvr-Format sind mittlerweile begrenzt nutzbar
 (Isooberfl&auml;chen
-sind darstellbar, siehe neue Laufparameter&nbsp; <tt>dt_dvrp, threshold</tt>),
+sind darstellbar, siehe neue Laufparameter&nbsp; <tt>dt_dvrp,
+threshold</tt>),
 es wird aber in der Anwendung in naher Zukunft noch diverse
 &Auml;nderungen
 geben. Partikelausgabe vorerst nicht mehr m&ouml;glich (Verlagerung
-schon).</td>
-
-      <td>header, init_particles, init_3d_model, leap_frog, modules,
+schon).</td> <td>header, init_particles, init_3d_model,
+leap_frog, modules,
 parin, parles,
-read_var_list, write_var_list&nbsp; <br>
-
-      <b>Gestrichen:</b> <br>
-
-plot_isosurface, plot_particles&nbsp;
-      
-      <p><b>Neu:</b> <br>
-
-init_dvrp, plot_dvrp</p>
-
-      </td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">26/04/00</td>
-
-      <td>SR</td>
-
-      <td>2.0a</td>
-
-      <td>C</td>
-
-      <td>Durch zus&auml;tzlichen zyklischen Rand im Bott-Chlond-Schema
+read_var_list, write_var_list&nbsp; <br> <b>Gestrichen:</b>
+<br>
+plot_isosurface, plot_particles&nbsp; <p><b>Neu:</b>
+<br>
+init_dvrp, plot_dvrp</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">26/04/00</td>
+<td>SR</td> <td>2.0a</td> <td>C</td>
+<td>Durch zus&auml;tzlichen zyklischen Rand im
+Bott-Chlond-Schema
 sehr
 viele SENDRECV-Aufrufe eingespart.&nbsp; <br>
-
-&Auml;nderung des Namens einer eingelesenen Environment-Variablen in <tt>check_open</tt>
-von <tt>remote_addres</tt> nach <tt>return_addres</tt>.&nbsp; <br>
-
-Prozessor-Topologie kann durch Benutzer vorgegeben werden. Dazu neue
-Initialisierungsparameter <tt>npex</tt> und <tt>npey</tt>. Durch
+&Auml;nderung des Namens einer eingelesenen Environment-Variablen
+in <tt>check_open</tt>
+von <tt>remote_addres</tt> nach <tt>return_addres</tt>.&nbsp;
+<br>Prozessor-Topologie kann durch Benutzer vorgegeben werden.
+Dazu neue
+Initialisierungsparameter <tt>npex</tt> und <tt>npey</tt>.
+Durch
 &Auml;nderungen
 am Header-Format wurde Testergebnis-Datei ge&auml;ndert.</td>
-
-      <td>advec_s_bc, check_open, header, init_pegrid, modules, parin,
+<td>advec_s_bc, check_open, header, init_pegrid, modules, parin,
 read_var_list,
-write_var_list</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">26/04/00</td>
-
-      <td>SR</td>
-
-      <td>2.0a</td>
-
-      <td>N</td>
-
-      <td>Rechnungen mit geneigter Oberfl&auml;che erlauben jetzt die
+write_var_list</td> </tr> <tr nosave="" valign="top"> <td nosave="">26/04/00</td>
+<td>SR</td> <td>2.0a</td> <td>N</td>
+<td>Rechnungen mit geneigter Oberfl&auml;che erlauben jetzt
+die
 Vorgabe
 eines stabil geschichteten Temperaturprofils mit konstantem Gradienten.
 Hinzuschalten von Feuchte funktioniert f&uuml;r Hangrechnungen noch
-nicht.</td>
-
-      <td>advec_s_bc, buoyancy, init_3d_model, modules,
-prognostic_equations&nbsp;
-      
-      <p><b>Neu:</b> <br>
-
-init_slope</p>
-
-      </td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td>13/06/00</td>
-
-      <td>MS</td>
-
-      <td>2.0a</td>
-
-      <td nosave="">B</td>
-
-      <td>Es werden nun auch die seitlichen R&auml;nder f&uuml;r ql
+nicht.</td> <td>advec_s_bc, buoyancy, init_3d_model,
+modules,
+prognostic_equations&nbsp; <p><b>Neu:</b> <br>
+init_slope</p> </td> </tr> <tr nosave="" valign="top"> <td>13/06/00</td> <td>MS</td>
+<td>2.0a</td> <td nosave="">B</td> <td>Es
+werden nun auch die seitlichen R&auml;nder f&uuml;r ql
 gesetzt.
 Ein Nichtsetzen der Randbedingungen f&uuml;hrte zu Fehlern bei der
 Ausgabe
-von pt</td>
-
-      <td nosave="">calc_liquid_water_content</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">03/07/00</td>
-
-      <td>SR</td>
-
-      <td>2.0b</td>
-
-      <td>C</td>
-
-      <td>F&uuml;r eine Reihe von&nbsp; Unterprogrammen, die innerhalb
+von pt</td> <td nosave="">calc_liquid_water_content</td>
+</tr> <tr nosave="" valign="top"> <td nosave="">03/07/00</td> <td>SR</td> <td>2.0b</td>
+<td>C</td> <td>F&uuml;r eine Reihe von&nbsp;
+Unterprogrammen, die innerhalb
 von <tt>prognostic_equation</tt>
-aufgerufen werden und an die Pointer-Variablen &uuml;bergeben werden,
+aufgerufen werden und an die Pointer-Variablen &uuml;bergeben
+werden,
 sind
 jetzt explizite Interfaces vereinbart (neues Modul <tt>pointer_interfaces</tt>).
 Auch im Vereinbarungsteil der entsprechenden Unterprogramme sind diese
-&Uuml;bergabeparameter jetzt aus Konsistenzgr&uuml;nden (und weil es
+&Uuml;bergabeparameter jetzt aus Konsistenzgr&uuml;nden (und
+weil es
 sonst
 Laufzeitfehler gibt) als Pointer-Variablen deklariert. Mit dieser
 &Auml;nderung
-ist der gro&szlig;e Speicher- und CPU-Zeit-Bedarf beim &Uuml;bersetzen
-von prognostic_equations wieder auf ein ertr&auml;gliches Ma&szlig;
+ist der gro&szlig;e Speicher- und CPU-Zeit-Bedarf beim
+&Uuml;bersetzen
+von prognostic_equations wieder auf ein ertr&auml;gliches
+Ma&szlig;
 reduziert.&nbsp; <br>
-
 Im Rahmen dieser &Auml;nderungen wurde auch an diffusion_e ein
 Hilfsfeld
 weniger &uuml;bergeben.&nbsp; <br>
-
-Aus nicht vollst&auml;ndig gek&auml;rten Gr&uuml;nden haben sich die
-Zahlen in der Testergebnis-Datei minimal ge&auml;ndert (und zwar die
+Aus nicht vollst&auml;ndig gek&auml;rten Gr&uuml;nden haben
+sich die
+Zahlen in der Testergebnis-Datei minimal ge&auml;ndert (und zwar
+die
 Spalte
 mit der Gesamtdivergenz nach Aufruf des Druckl&ouml;sers).</td>
-
-      <td>advec_s_bc, buoyancy, diffusion_e, diffusion_s, diffusion_u,
+<td>advec_s_bc, buoyancy, diffusion_e, diffusion_s, diffusion_u,
 diffusion_v,
 diffusion_w, disturb_field, modules, production_e, prognostic_equations</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">04/07/00</td>
-
-      <td>SR</td>
-
-      <td>2.0b</td>
-
-      <td>B</td>
-
-      <td>Diriclet-Randbedingungen f&uuml;r Temperatur und
+</tr> <tr nosave="" valign="top"> <td nosave="">04/07/00</td> <td>SR</td> <td>2.0b</td>
+<td>B</td> <td>Diriclet-Randbedingungen f&uuml;r
+Temperatur und
 Fl&uuml;ssigwassergehalt
-werden gesetzt. Dies w&auml;re bereits nach Einf&uuml;hrung der Pointer
+werden gesetzt. Dies w&auml;re bereits nach Einf&uuml;hrung der
+Pointer
 zwingend notwendig gewesen. Rechnungen mit vorgegebener
 Oberfl&auml;chentemperatur
-schlugen deshalb bislang fehl.</td>
-
-      <td>boundary_conds</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">07/09/00</td>
-
-      <td>MS</td>
-
-      <td>2.0b</td>
-
-      <td>B</td>
-
-      <td>Die virtuelle potenielle Temperatur wird nun auch f&uuml;r
+schlugen deshalb bislang fehl.</td> <td>boundary_conds</td>
+</tr> <tr nosave="" valign="top"> <td nosave="">07/09/00</td> <td>MS</td> <td>2.0b</td>
+<td>B</td> <td>Die virtuelle potenielle Temperatur
+wird nun auch f&uuml;r
 k=nzb und
 k=nzt+1 berechnet, um die Randwerte auch f&uuml;r Ausgabezwecke
-bereichtzustellen.&nbsp;</td>
-
-      <td>compute_vpt</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">28/12/00</td>
-
-      <td>SR</td>
-
-      <td>2.1</td>
-
-      <td>C/N</td>
-
-      <td>PALM erlaubt jetzt den Einsatz optionaler Software-Pakete,
+bereichtzustellen.&nbsp;</td> <td>compute_vpt</td>
+</tr> <tr nosave="" valign="top"> <td nosave="">28/12/00</td> <td>SR</td> <td>2.1</td>
+<td>C/N</td> <td>PALM erlaubt jetzt den Einsatz
+optionaler Software-Pakete,
 die zwar
-zum Modell geh&ouml;ren, standardm&auml;&szlig;ig aber nicht mit
+zum Modell geh&ouml;ren, standardm&auml;&szlig;ig aber
+nicht mit
 &uuml;bersetzt
 werden, um so unter anderem Kompilationszeit einzusparen. Siehe neues
 Kapitel
-3.7 in der Modelldokumentation.&nbsp;
-      
-      <p>Einsatz der dvrp-Software wurde vollst&auml;ndig
+3.7 in der Modelldokumentation.&nbsp; <p>Einsatz der
+dvrp-Software wurde vollst&auml;ndig
 &uuml;berarbeitet.
 Sie ist nun als Software-Paket optional im Modell einsetzbar. Siehe
 &uuml;berarbeitetes
-Kapitel 4.5.6 in der Modelldokumentation.&nbsp; </p>
-
-      
-      <p>Die Prognose der Partikelverlagerung ist ebenfalls in ein
+Kapitel 4.5.6 in der Modelldokumentation.&nbsp; </p> <p>Die
+Prognose der Partikelverlagerung ist ebenfalls in ein
 optionales
-Software-Paket ausgelagert.</p>
-
-      </td>
-
-      <td>advec_particles, header, init_3d_model, init_dvrp,
+Software-Paket ausgelagert.</p> </td> <td>advec_particles,
+header, init_3d_model, init_dvrp,
 init_particles,
 init_pegrid, modules, parin, parles, plot_dvrp, prognostic_equations,
 read_var_list,
-write_var_list&nbsp;
-      
-      <p><b>Neu:</b> <br>
-
-package_parin</p>
-
-      </td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">28/12/00</td>
-
-      <td>SR</td>
-
-      <td>2.1</td>
-
-      <td>B</td>
-
-      <td>MPI_FINALIZE ans Ende des Hauptprogramms geschoben (stand
+write_var_list&nbsp; <p><b>Neu:</b> <br>
+package_parin</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">28/12/00</td>
+<td>SR</td> <td>2.1</td> <td>B</td>
+<td>MPI_FINALIZE ans Ende des Hauptprogramms geschoben (stand
 vorher vor
 cpu_auswertung, wo noch ein Barrier-Aufruf steht). comm2d war zu Beginn
 undefiniert und wird nun erst einmal gleich MPI_COMM_WORLD gesetzt.</td>
-
-      <td>parles</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">02/01/01</td>
-
-      <td>SR</td>
-
-      <td>2.1a</td>
-
-      <td>C</td>
-
-      <td>Offene Dateien werden sobald m&ouml;glich geschlossen.
+<td>parles</td> </tr> <tr nosave="" valign="top"> <td nosave="">02/01/01</td>
+<td>SR</td> <td>2.1a</td> <td>C</td>
+<td>Offene Dateien werden sobald m&ouml;glich geschlossen.
 Wiederer&ouml;ffnung
 bestimmter Dateien mit POSITION='APPEND' m&ouml;glich. In diesem
@@ -1583,157 +726,66 @@
 wurde das Unterprogramm close_files in close_file umbenannt und hat nun
 1 Argument.&nbsp; <br>
-
 Anpassung der Diffusionsparametrisierung im 1D-Modell an das 3D-Modell.</td>
-
-      <td>advec_particles, check_cpu_time, check_open, cpu_auswertung,
+<td>advec_particles, check_cpu_time, check_open, cpu_auswertung,
 init_1d_model,
 init_3d_model, init_particles, modules.f90, parin, parles, plot_2d,
-write_3d_binary&nbsp;
-      
-      <p><b>Neu:</b> <br>
-
-close_file&nbsp; </p>
-
-      
-      <p><b>Eliminiert:</b> <br>
-
-close_files</p>
-
-      </td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">02/01/01</td>
-
-      <td>SR</td>
-
-      <td>2.1a</td>
-
-      <td>B</td>
-
-      <td>Beseitigung kleiner Fehler.</td>
-
-      <td>check_cpu_time, diffusion_e, diffusivities, long_filter,
-production_e</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td>05/01/01</td>
-
-      <td>SR</td>
-
-      <td>2.1b</td>
-
-      <td>N</td>
-
-      <td nosave="">Neues Software-Paket zur Berechnung von Spektren im
-Ortsraum</td>
-
-      <td>check_open, header, init_3d_model, leap_frog, modules,
+write_3d_binary&nbsp; <p><b>Neu:</b> <br>
+close_file&nbsp; </p> <p><b>Eliminiert:</b>
+<br>
+close_files</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">02/01/01</td>
+<td>SR</td> <td>2.1a</td> <td>B</td>
+<td>Beseitigung kleiner Fehler.</td> <td>check_cpu_time,
+diffusion_e, diffusivities, long_filter,
+production_e</td> </tr> <tr nosave="" valign="top"> <td>05/01/01</td> <td>SR</td>
+<td>2.1b</td> <td>N</td> <td nosave="">Neues
+Software-Paket zur Berechnung von Spektren im
+Ortsraum</td> <td>check_open, header, init_3d_model,
+leap_frog, modules,
 package_parin,
-parles, read_var_list, write_3d_binary, write_var_list&nbsp;
-      
-      <p><b>Neu:</b> <br>
-
-calc_spectra, plot_spectra</p>
-
-      </td>
-
-    </tr>
-
-    <tr valign="top">
-
-      <td>25/01/01</td>
-
-      <td>SR</td>
-
-      <td>2.1c</td>
-
-      <td>C/N</td>
-
-      <td>Druckl&ouml;ser kann jetzt die Singleton-FFT benutzen. Modul <tt>singleton</tt>
+parles, read_var_list, write_3d_binary, write_var_list&nbsp; <p><b>Neu:</b>
+<br>
+calc_spectra, plot_spectra</p> </td> </tr> <tr valign="top"> <td>25/01/01</td> <td>SR</td>
+<td>2.1c</td> <td>C/N</td> <td>Druckl&ouml;ser
+kann jetzt die Singleton-FFT benutzen. Modul <tt>singleton</tt>
 wurde deshalb von <tt>calc_spectra.f90</tt> nach <tt>modules.f90</tt>
 geschoben
-(neuer Initialisierungsparameter <tt>fft_method</tt>).&nbsp; <br>
-
-Fl&uuml;sse k&ouml;nnen jetzt bei k=1 unabh&auml;ngig von der
+(neuer Initialisierungsparameter <tt>fft_method</tt>).&nbsp;
+<br>Fl&uuml;sse k&ouml;nnen jetzt bei k=1
+unabh&auml;ngig von der
 Verwendung
 einer Prandtl-Schicht vorgegeben werden (neuer
-Initialisierungsparameter <tt>use_surface_fluxes</tt>).&nbsp; <br>
-
-Modul <tt>test_variables</tt> wurde aus allen Programmteilen entfernt.
-Hauptprogramm wurde von parles nach <tt>palm</tt> umbenannt.&nbsp; <br>
-
+Initialisierungsparameter <tt>use_surface_fluxes</tt>).&nbsp;
+<br>Modul <tt>test_variables</tt> wurde aus allen
+Programmteilen entfernt.
+Hauptprogramm wurde von parles nach <tt>palm</tt>
+umbenannt.&nbsp; <br>
 Modell kann auf DEC-Workstations eingesetzt werden (-D <tt>dec</tt>).&nbsp;
-      <br>
-
-Zus&auml;tzliche Zeitmessungen mit dvrp-Software. Schreiben von
+<br>Zus&auml;tzliche Zeitmessungen mit dvrp-Software.
+Schreiben von
 Partikel-Informationen
 ist optional (neuer Parameter <tt>write_particle_informations</tt>).</td>
-
-      <td>advec_particles, calc_spectra, check_parameters,
+<td>advec_particles, calc_spectra, check_parameters,
 cpu_zeitmessung, diffusion_s,
 header, init_3d_model, init_dvrp, init_particles, leap_frog, modules,
 package_parin,
 parin, poisfft, read_var_list, swap_timelevel, write_var_list&nbsp;
-      
-      <p><b>Neu:</b> <br>
-
-palm&nbsp; </p>
-
-      
-      <p><b>Eliminiert:</b> <br>
-
-parles, module_test</p>
-
-      </td>
-
-    </tr>
-
-    <tr valign="top">
-
-      <td>25/01/01</td>
-
-      <td>SR</td>
-
-      <td>2.1c</td>
-
-      <td>B</td>
-
-      <td>Fehler beim &ouml;ffnen von Unit 23 beseitigt. Modul <tt>interface</tt>
+<p><b>Neu:</b> <br>
+palm&nbsp; </p> <p><b>Eliminiert:</b> <br>
+parles, module_test</p> </td> </tr> <tr valign="top"> <td>25/01/01</td> <td>SR</td>
+<td>2.1c</td> <td>B</td> <td>Fehler
+beim &ouml;ffnen von Unit 23 beseitigt. Modul <tt>interface</tt>
 fehlte in disturb_field.&nbsp; <br>
-
 Wertebereich von theta* in Prandtl-Fluxes eingeschr&auml;nkt, weil
-sonst auf einigen Rechnern overflow droht, wenn u*=0.</td>
-
-      <td>check_open, disturb_field, prandtl_fluxes</td>
-
-    </tr>
-
-    <tr valign="top">
-
-      <td>30/01/01</td>
-
-      <td>SR</td>
-
-      <td>2.1d</td>
-
-      <td>N/C</td>
-
-      <td>Prognose eines passiven Skalars ist m&ouml;glich (anstatt
+sonst auf einigen Rechnern overflow droht, wenn u*=0.</td> <td>check_open,
+disturb_field, prandtl_fluxes</td> </tr> <tr valign="top"> <td>30/01/01</td> <td>SR</td>
+<td>2.1d</td> <td>N/C</td> <td>Prognose
+eines passiven Skalars ist m&ouml;glich (anstatt
 Feuchte),
 dazu neuer Initialisierungsparameter <tt>passive_scalar</tt>.&nbsp;
-      
-      <p>Falls kein Pfad f&uuml;r ftpcopy existiert, werden in <tt>check_open</tt>
-die Dateinamen f&uuml;r AVS-Koordinaten- und Datendatei auf "unknown"
-gesetzt.</p>
-
-      </td>
-
-      <td>asselin_filter, boundary_conds, check_open, check_parameters,
+<p>Falls kein Pfad f&uuml;r ftpcopy existiert, werden in <tt>check_open</tt>
+die Dateinamen f&uuml;r AVS-Koordinaten- und Datendatei auf
+"unknown"
+gesetzt.</p> </td> <td>asselin_filter,
+boundary_conds, check_open, check_parameters,
 flow_statistics,
 header, init_1d_model, init_3d_model, modules, parin, plot_2d,
@@ -1741,39 +793,14 @@
 prandtl_fluxes, prognostic_equations, read_var_list, swap_timelevel,
 write_3d_binary,
-write_var_list</td>
-
-    </tr>
-
-    <tr valign="top">
-
-      <td>30/01/01</td>
-
-      <td>SR</td>
-
-      <td>2.1d</td>
-
-      <td>B</td>
-
-      <td>String-Ausgabe f&uuml;r use_surface_fluxes korrigiert.
+write_var_list</td> </tr> <tr valign="top"> <td>30/01/01</td>
+<td>SR</td> <td>2.1d</td> <td>B</td>
+<td>String-Ausgabe f&uuml;r use_surface_fluxes korrigiert.
 Version vom
 25/01/01 brach deswegen bei Fortsetzungsl&auml;ufen ab.&nbsp; <br>
-
-xy-Plotausgabe von q und ql korrigiert.</td>
-
-      <td>plot_2d, write_var_list</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">30/03/01</td>
-
-      <td>SR</td>
-
-      <td>2.2</td>
-
-      <td>N/C</td>
-
-      <td>S&auml;mtliche 3D-Felder der Zeitebene t+dt wurden entfernt.
+xy-Plotausgabe von q und ql korrigiert.</td> <td>plot_2d,
+write_var_list</td> </tr> <tr nosave="" valign="top"> <td nosave="">30/03/01</td>
+<td>SR</td> <td>2.2</td> <td>N/C</td>
+<td>S&auml;mtliche 3D-Felder der Zeitebene t+dt wurden
+entfernt.
 Ebenso
 die Arbeitsfelder work und work1. Hilfsfelder wurden soweit
@@ -1783,60 +810,55 @@
 eine Reihe von &Uuml;bergabeparametern gestrichen werden (z.B. beim
 Aufruf
-von <tt>pres, advec_s_bc, poisfft, production_e, sor </tt>sowie den
+von <tt>pres, advec_s_bc, poisfft, production_e, sor </tt>sowie
+den
 Spline-Unterprogrammen).
 Der Long-Filter wurde in das Upstream-Spline-Verfahren integriert. Da
 der
 Zeitfilter jetzt direkt in der prognostischen Gleichung angewendet wird
-(bisher geschah dies erst nach Aufruf des Druckl&ouml;sers), haben sich
+(bisher geschah dies erst nach Aufruf des Druckl&ouml;sers), haben
+sich
 die Zahlen in der Testergebnis-Datei ge&auml;ndert. Die mittleren
 Profile
-sind aber unbeeinflusst.&nbsp;
-      
-      <p>Eine weitere kleinere &Auml;nderung der Testergebnis-Daten
+sind aber unbeeinflusst.&nbsp; <p>Eine weitere kleinere
+&Auml;nderung der Testergebnis-Daten
 wird dadurch
 verursacht, dass jetzt der Parameter <tt>adjust_mixing_length</tt>
 defaultm&auml;&szlig;ig <tt>.FALSE.</tt>
-ist und der Mischungsweg grunds&auml;tzlich zus&auml;tzlich auf&nbsp; <tt>0.7
+ist und der Mischungsweg grunds&auml;tzlich zus&auml;tzlich
+auf&nbsp; <tt>0.7
 * zu</tt>&nbsp;&nbsp; begrenzt wird. Falls <tt>adjust_mixing_length=T</tt>,
-wird jetzt der Mischungsweg in <tt>diffusivities</tt> analog zu <tt>diffusion_e</tt>
+wird jetzt der Mischungsweg in <tt>diffusivities</tt>
+analog zu <tt>diffusion_e</tt>
 an allen Gitterpunkten modifiziert. <b>Achtung: </b>Diese
 &Auml;nderungen
-k&ouml;nnen bei anderen Simulationen t.w. erhebliche Auswirkungen auf
+k&ouml;nnen bei anderen Simulationen t.w. erhebliche Auswirkungen
+auf
 die
-oberfl&auml;chennahen Ergebnisse haben.&nbsp; </p>
-
-      
-      <p><b>Kleinere &Auml;nderungen:</b> In <tt>init_cloud_physics</tt>
+oberfl&auml;chennahen Ergebnisse haben.&nbsp; </p> <p><b>Kleinere
+&Auml;nderungen:</b> In <tt>init_cloud_physics</tt>
 wird
-der Wert der Variablen <tt>surface_pressure</tt> in hPa belassen
+der Wert der Variablen <tt>surface_pressure</tt> in hPa
+belassen
 (bisher
 Umwandlung in Pa).&nbsp; <br>
-
 Die Namelist-Namen der Software-Pakete sind&nbsp; t.w.
 ge&auml;ndert.&nbsp; <br>
-
-Unit 14 (bin&auml;re Ausgabe f&uuml;r Fortsetzungsl&auml;ufe) wird
+Unit 14 (bin&auml;re Ausgabe f&uuml;r
+Fortsetzungsl&auml;ufe) wird
 jetzt im Hauptprogramm geschlossen, damit auf diese Datei noch
 benutzergesteuerte
 Ausgaben erfolgen k&ouml;nnen.&nbsp; <br>
-
 Laufparameter werden in keinem Fall mehr auf Unit 14 ausgegeben bzw.
 von Unit 13 gelesen, d.h. sie gelten jetzt tats&auml;chlich nur
 f&uuml;r
 den jeweils aktuellen Lauf.&nbsp; <br>
-
 Partikeladvektion funktioniert jetzt auch zusammen mit der
-Galilei-Transformation.&nbsp; </p>
-
-      
-      <p>Restliche deutsche Variablen-, Unterprogramm- und Modulnamen
+Galilei-Transformation.&nbsp; </p> <p>Restliche
+deutsche Variablen-, Unterprogramm- und Modulnamen
 wurden
 ins Englische &uuml;bersetzt. Davon sind fast alle Programmteile
 betroffen
-(nicht in rechter Spalte aufgef&uuml;hrt).</p>
-
-      </td>
-
-      <td>advec_particles, advec_s_bc, advec_s_ups, advec_u_ups,
+(nicht in rechter Spalte aufgef&uuml;hrt).</p> </td> <td>advec_particles,
+advec_s_bc, advec_s_ups, advec_u_ups,
 advec_v_ups,
 advec_w_ups, boundary_conds, calc_spectra, check_parameters,
@@ -1847,33 +869,13 @@
 prognostic_equations,
 read_var_list, sor, swap_timelevel, transpose_*, write_3d_binary,
-write_var_list&nbsp;
-      
-      <p><b>Neu:</b> <br>
-
-cpu_log, cpu_statistics&nbsp; </p>
-
-      
-      <p><b>Eliminiert:</b> <br>
-
-asselin_filter, cpu_auswertung, cpu_zeitmessung, long_filter</p>
-
-      </td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">12/07/01</td>
-
-      <td>SR</td>
-
-      <td>2.2a</td>
-
-      <td>N/C</td>
-
-      <td>Defaultwert der unteren Randbedingung f&uuml;r die TKE ist ab
+write_var_list&nbsp; <p><b>Neu:</b> <br>
+cpu_log, cpu_statistics&nbsp; </p> <p><b>Eliminiert:</b>
+<br>
+asselin_filter, cpu_auswertung, cpu_zeitmessung, long_filter</p> </td>
+</tr> <tr nosave="" valign="top"> <td nosave="">12/07/01</td> <td>SR</td> <td>2.2a</td>
+<td>N/C</td> <td>Defaultwert der unteren
+Randbedingung f&uuml;r die TKE ist ab
 sofort <tt>bc_e_b
 = 'neumann'</tt>.&nbsp; <br>
-
 Partikelquelle ist &uuml;ber Paketparameter steuerbar.
 Partikeleigenschaften
@@ -1881,91 +883,49 @@
 ist
 jetzt ebenfalls &uuml;ber Paketparameter steuerbar.&nbsp; <br>
-
 Allen Modulen (bis auf singleton) wurde eine SAVE-Anweisung
 hinzugef&uuml;gt,
-damit die durch sie vereinbarten Variablen w&auml;hrend der Rechnung in
-keinem Fall undefiniert werden k&ouml;nnen (diese Gefahr besteht z.B.
+damit die durch sie vereinbarten Variablen w&auml;hrend der
+Rechnung in
+keinem Fall undefiniert werden k&ouml;nnen (diese Gefahr besteht
+z.B.
 auf
-SGI-Origin-Maschinen).</td>
-
-      <td>advec_particles, header, init_dvrp, init_particles,
+SGI-Origin-Maschinen).</td> <td>advec_particles, header,
+init_dvrp, init_particles,
 modules.f90, package_parin,
-plot_dvrp, user_interface</td>
-
-    </tr>
-
-    <tr>
-
-      <td>12/07/01</td>
-
-      <td>SR</td>
-
-      <td>2.2a</td>
-
-      <td>B</td>
-
-      <td>Verschiebung der OPEN-Anweisung f&uuml;r Unit 33
-(AVS-FLD-Datei).</td>
-
-      <td>check_open</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">20/07/01</td>
-
-      <td>SR</td>
-
-      <td>2.3</td>
-
-      <td>N</td>
-
-      <td>Einbau des Mehrgitterverfahrens zur L&ouml;sung der
-Poisson-Gleichung.</td>
-
-      <td>check_parameters, exchange_horiz, header, init_grid,
+plot_dvrp, user_interface</td> </tr> <tr> <td>12/07/01</td>
+<td>SR</td> <td>2.2a</td> <td>B</td>
+<td>Verschiebung der OPEN-Anweisung f&uuml;r Unit 33
+(AVS-FLD-Datei).</td> <td>check_open</td> </tr>
+<tr nosave="" valign="top"> <td nosave="">20/07/01</td>
+<td>SR</td> <td>2.3</td> <td>N</td>
+<td>Einbau des Mehrgitterverfahrens zur L&ouml;sung der
+Poisson-Gleichung.</td> <td>check_parameters,
+exchange_horiz, header, init_grid,
 init_pegrid, modules,
-parin, pres, run_control&nbsp;
-      
-      <p><b>Neu:</b> <br>
-
-poismg</p>
-
-      </td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">21/08/01</td>
-
-      <td>SR</td>
-
-      <td>2.3a</td>
-
-      <td>N/C</td>
-
-      <td>Erg&auml;nzung der dvrp-Software: Partikel k&ouml;nnen mit
+parin, pres, run_control&nbsp; <p><b>Neu:</b> <br>
+poismg</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">21/08/01</td>
+<td>SR</td> <td>2.3a</td> <td>N/C</td>
+<td>Erg&auml;nzung der dvrp-Software: Partikel
+k&ouml;nnen mit
 Schw&auml;nzen
 versehen werden. Die Ausgabe der dvrp-Plotdaten durch einen separaten
 PE
-ist m&ouml;glich (mrun-Option "<tt>-p dvrp_graphics+1PE</tt>"). Der
+ist m&ouml;glich (mrun-Option "<tt>-p dvrp_graphics+1PE</tt>").
+Der
 Kommunikator <tt>MPI_COMM_WORLD</tt>&nbsp;
-wurde deshalb durch&nbsp; <tt>comm_palm</tt>&nbsp; ersetzt.&nbsp;
-      
-      <p>Reduzierung des Mischungsweges im Wandbereich kann wahlweise
+wurde deshalb durch&nbsp; <tt>comm_palm</tt>&nbsp;
+ersetzt.&nbsp; <p>Reduzierung des Mischungsweges im
+Wandbereich kann wahlweise
 abgeschaltet
 werden. F&uuml;r km und kh wird jetzt eine Neumann-Randbedingung am
 unteren
 Rand verwendet. Die TKE-Energieproduktionsberechnung durch Scherung des
-Grundstroms am unteren Rand wurde verbessert.&nbsp; </p>
-
-      
-      <p>Generelle Verwendung der bodennahen vertikalen
+Grundstroms am unteren Rand wurde verbessert.&nbsp; </p> <p>Generelle
+Verwendung der bodennahen vertikalen
 Impulsfl&uuml;sse in
 den Diffusionstermen von u und v, immer wenn diese als Randbedingung
 vorgegeben
-sind (gem&auml;&szlig; <tt>use_surface_fluxes</tt>, - bisher wurde
+sind (gem&auml;&szlig; <tt>use_surface_fluxes</tt>,
+- bisher wurde
 dies
 nur bei eingeschalteter Prandtl-Schicht gemacht). Schubspannungen
@@ -1973,25 +933,21 @@
 jetzt als Randbedingungen vorgegeben werden (allerdings unter
 zuhilfenahme
-benutzereigener Software).&nbsp; </p>
-
-      
-      <p>Zus&auml;tzliche Schnittstellen f&uuml;r benutzereigene
+benutzereigener Software).&nbsp; </p> <p>Zus&auml;tzliche
+Schnittstellen f&uuml;r benutzereigene
 Software in
-allen Tendenztermen, in <tt>flow_statistics,</tt> und f&uuml;r
+allen Tendenztermen, in <tt>flow_statistics,</tt> und
+f&uuml;r
 Partikel-
 bzw. dvrp-Programmteile (Bestimmung von Partikeleigenschaften,
 Festlegung
-von Farbtabellen).&nbsp; </p>
-
-      
-      <p><b>Achtung: </b>Diese &Auml;nderungen k&ouml;nnen bei
+von Farbtabellen).&nbsp; </p> <p><b>Achtung: </b>Diese
+&Auml;nderungen k&ouml;nnen bei
 Simulationen mit
-mittlerem Wind durch die &Auml;nderungen im bereich des unteren Randes
-t.w. erhebliche Auswirkungen auf die oberfl&auml;chennahen Ergebnisse
-haben.</p>
-
-      </td>
-
-      <td>advec_particles, diffusion_e, diffusion_u, diffusion_v,
+mittlerem Wind durch die &Auml;nderungen im bereich des unteren
+Randes
+t.w. erhebliche Auswirkungen auf die oberfl&auml;chennahen
+Ergebnisse
+haben.</p> </td> <td>advec_particles, diffusion_e,
+diffusion_u, diffusion_v,
 diffusivities,
 flow_statistics, header, init_dvrp, init_particles, init_pegrid,
@@ -2000,55 +956,20 @@
 prognostic_equations,
 read_var_list, swap_timelevel, user_interface, write_var_list</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">21/08/01</td>
-
-      <td>SR</td>
-
-      <td>2.3a</td>
-
-      <td>B</td>
-
-      <td>Falsche Positionierung von MPI_ALLREDUCE in check_cpu_time
-korrigiert.</td>
-
-      <td>check_cpu_time</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">04/09/01</td>
-
-      <td>SR</td>
-
-      <td>2.3b</td>
-
-      <td>N</td>
-
-      <td>Zus&auml;tzliche Profilausgaben f&uuml;r
+</tr> <tr nosave="" valign="top"> <td nosave="">21/08/01</td> <td>SR</td> <td>2.3a</td>
+<td>B</td> <td>Falsche Positionierung von
+MPI_ALLREDUCE in check_cpu_time
+korrigiert.</td> <td>check_cpu_time</td> </tr>
+<tr nosave="" valign="top"> <td nosave="">04/09/01</td>
+<td>SR</td> <td>2.3b</td> <td>N</td>
+<td>Zus&auml;tzliche Profilausgaben f&uuml;r
 Energieproduktionsterme m&ouml;glich.
-St&ouml;rdruck wird bei der FFT-Methode nicht mehr durch Aufsummieren
+St&ouml;rdruck wird bei der FFT-Methode nicht mehr durch
+Aufsummieren
 der
-Werte zu den einzelnen Zeitschritten ermittelt.</td>
-
-      <td>check_parameters, flow_statistics, modules, pres</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">04/09/01</td>
-
-      <td>SR</td>
-
-      <td>2.3b</td>
-
-      <td>B</td>
-
-      <td>Austausch der Geisterr&auml;nder f&uuml;r die zeitgefilterten
+Werte zu den einzelnen Zeitschritten ermittelt.</td> <td>check_parameters,
+flow_statistics, modules, pres</td> </tr> <tr nosave="" valign="top"> <td nosave="">04/09/01</td>
+<td>SR</td> <td>2.3b</td> <td>B</td>
+<td>Austausch der Geisterr&auml;nder f&uuml;r die
+zeitgefilterten
 Felder
 direkt nach Durchf&uuml;hrung der Filterung (war fehlerhaft seit
@@ -2056,84 +977,51 @@
 der Zeitebene t+dt (Version 2.2) und f&uuml;hrte zu sehr kleinen
 St&ouml;rungen
-an den seitlichen R&auml;ndern der Teilgebiete).&nbsp;
-      
-      <p><b>Achtung:</b> Der Inhalt der Testergebnis-Datei &auml;ndert
+an den seitlichen R&auml;ndern der Teilgebiete).&nbsp; <p><b>Achtung:</b>
+Der Inhalt der Testergebnis-Datei &auml;ndert
 sich durch
-diese Fehlerkorrektur.</p>
-
-      </td>
-
-      <td>prognostic_equations</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">09/11/01</td>
-
-      <td>SR</td>
-
-      <td>2.3c</td>
-
-      <td>N</td>
-
-      <td>Farbe entlang der Partikelschw&auml;nze kann sich
+diese Fehlerkorrektur.</p> </td> <td>prognostic_equations</td>
+</tr> <tr nosave="" valign="top"> <td nosave="">09/11/01</td> <td>SR</td> <td>2.3c</td>
+<td>N</td> <td>Farbe entlang der
+Partikelschw&auml;nze kann sich
 &auml;ndern. Einschalten
 der Partikeladvektion bei Fortsetzungsl&auml;ufen m&ouml;glich.</td>
-
-      <td>advec_particles, init_particles, modules, package_parin,
-plot_dvrp</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">16/04/02</td>
-
-      <td>SR</td>
-
-      <td>2.3d</td>
-
-      <td>N</td>
-
-      <td>Vorgabe von Randbedingungen f&uuml;r Partikeladvektion
+<td>advec_particles, init_particles, modules, package_parin,
+plot_dvrp</td> </tr> <tr nosave="" valign="top">
+<td nosave="">16/04/02</td> <td>SR</td>
+<td>2.3d</td> <td>N</td> <td>Vorgabe
+von Randbedingungen f&uuml;r Partikeladvektion
 m&ouml;glich.
-Partikeldaten k&ouml;nnen f&uuml;r sp&auml;tere Analysen auf Datei
+Partikeldaten k&ouml;nnen f&uuml;r sp&auml;tere Analysen
+auf Datei
 geschrieben
-werden (Unit 85). <tt>PARTICLE</tt>-Datentyp enth&auml;lt
+werden (Unit 85). <tt>PARTICLE</tt>-Datentyp
+enth&auml;lt
 Informationen
-&uuml;ber Partikelgeschwindigkeitskomponenten und Startposition.&nbsp;
-      
-      <p>Skalartransport sowie entsprechende Datenausgaben k&ouml;nnen
+&uuml;ber Partikelgeschwindigkeitskomponenten und
+Startposition.&nbsp; <p>Skalartransport sowie entsprechende
+Datenausgaben k&ouml;nnen
 durch
 eigene Parameter gesteuert werden (nicht mehr &uuml;ber die
-Feuchteparameter).&nbsp; </p>
-
-      
-      <p>Im Fall von pdims(1)=1 (eindimensionales virtuelles
+Feuchteparameter).&nbsp; </p> <p>Im Fall von
+pdims(1)=1 (eindimensionales virtuelles
 Prozessornetz in
 y-Richtung) werden die zyklischen Randbedingungen in x-Richtung durch
 direktes
-Umspeichern anstatt durch&nbsp; <tt>sendrecv </tt>realisiert. Ebenso
+Umspeichern anstatt durch&nbsp; <tt>sendrecv </tt>realisiert.
+Ebenso
 werden
-die Transponierungen xz, yz, zx&nbsp; und zy eingespart (es wird aber
+die Transponierungen xz, yz, zx&nbsp; und zy eingespart (es wird
+aber
 weiter
-umsortiert).&nbsp; </p>
-
-      
-      <p>Im Fall von 3D-Plotausgaben Aufruf von ftpcopy-Script durch
+umsortiert).&nbsp; </p> <p>Im Fall von
+3D-Plotausgaben Aufruf von ftpcopy-Script durch
 batch_scp-Script
 ersetzt (dient der Ermittelung von Dateizyklusnummern als Information
 f&uuml;r
-die AVS-fld-Datei).&nbsp; </p>
-
-      
-      <p>Bei fehlerhafter Er&ouml;ffnung einer CPU_MEASURES-Datei
+die AVS-fld-Datei).&nbsp; </p> <p>Bei fehlerhafter
+Er&ouml;ffnung einer CPU_MEASURES-Datei
 werden erneute
-Er&ouml;ffnungsversuche durchgef&uuml;hrt.</p>
-
-      </td>
-
-      <td>advec_particles, check_open, check_parameters,
+Er&ouml;ffnungsversuche durchgef&uuml;hrt.</p> </td>
+<td>advec_particles, check_open, check_parameters,
 exchange_horiz, flow_statistics,
 header, init_particles, modules, package_parin, parin, plot_2d,
@@ -2141,110 +1029,52 @@
 prognostic_equations, read_var_list, transpose_xz, transpose_yz,
 transpose_zx,
-transpose_zy, write_var_list</td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">02/05/02</td>
-
-      <td>SR</td>
-
-      <td>2.3e</td>
-
-      <td>B</td>
-
-      <td>Wiedereinf&uuml;hrung der 3D-Felder f&uuml;r die Zeitebene
+transpose_zy, write_var_list</td> </tr> <tr nosave="" valign="top"> <td nosave="">02/05/02</td>
+<td>SR</td> <td>2.3e</td> <td>B</td>
+<td>Wiedereinf&uuml;hrung der 3D-Felder f&uuml;r die
+Zeitebene
 t+dt, die
 notwendig ist, damit in den Diffusionstermen bei Leapfrog-Zeitschritten
 mit der korrekten Zeitebene t-dt gearbeitet wird. Seit Version 2.2
 wurde
-f&auml;lschlicherweise die Zeitebene t+dt verwendet. Asselin-Filter ist
+f&auml;lschlicherweise die Zeitebene t+dt verwendet. Asselin-Filter
+ist
 jetzt wieder eigenst&auml;ndiges Unterprogramm. <br>
-
 Die Wiedereinf&uuml;hrung dieser dritten Zeitebene ist auch
 Voraussetzung
 f&uuml;r die in der n&auml;chsten Version geplante skalare
-Optimierung.&nbsp;
-      
-      <p>Fehler bei Berechnung von Ausgabezeitpunkten eliminiert (trat
+Optimierung.&nbsp; <p>Fehler bei Berechnung von
+Ausgabezeitpunkten eliminiert (trat
 auf, wenn
-bei Fortsetzungsl&auml;ufen von Ausgabeabst&auml;nden 0.0 auf von Null
-verschiedene Werte gewechselt werden sollte).&nbsp; </p>
-
-      
-      <p><b>Achtung:</b> Der Inhalt der Testergebnis-Datei &auml;ndert
+bei Fortsetzungsl&auml;ufen von Ausgabeabst&auml;nden 0.0 auf
+von Null
+verschiedene Werte gewechselt werden sollte).&nbsp; </p> <p><b>Achtung:</b>
+Der Inhalt der Testergebnis-Datei &auml;ndert
 sich durch
-diese Fehlerkorrektur.</p>
-
-      </td>
-
-      <td>init_3d_model, leap_frog, modules, prognostic_equations,
-swap_timelevel&nbsp;
-      
-      <p><b>Neu:</b> <br>
-
-asselin_filter</p>
-
-      </td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td>02/05/02</td>
-
-      <td>SR</td>
-
-      <td>2.3e</td>
-
-      <td>N/C</td>
-
-      <td>Kleinere Anpassungen an IBM-Regatta-Systeme in
-check_parameters, cpu_log.&nbsp;
-      
-      <p>Modul singleton in eigenst&auml;ndige Datei
-&uuml;berf&uuml;hrt.&nbsp; </p>
-
-      
-      <p>Global_min_max arbeitet mit REAL*4, um Komplikationen auf
+diese Fehlerkorrektur.</p> </td> <td>init_3d_model,
+leap_frog, modules, prognostic_equations,
+swap_timelevel&nbsp; <p><b>Neu:</b> <br>
+asselin_filter</p> </td> </tr> <tr nosave="" valign="top"> <td>02/05/02</td> <td>SR</td>
+<td>2.3e</td> <td>N/C</td> <td>Kleinere
+Anpassungen an IBM-Regatta-Systeme in
+check_parameters, cpu_log.&nbsp; <p>Modul singleton in
+eigenst&auml;ndige Datei
+&uuml;berf&uuml;hrt.&nbsp; </p> <p>Global_min_max
+arbeitet mit REAL*4, um Komplikationen auf
 32-bit-Rechnern
 zu vermeiden (Datentyp MPI_2REAL m&uuml;sste dort sonst auf
 MPI_2DOUBLE_PRECISION
 ge&auml;ndert werden). <br>
-
 Horizontale Geschwindigkeitskomponenten f&uuml;r Partikeladvektion
 werden exakt zwischen den jeweils benachbarten horizontalen
 Gitterfl&auml;chen
-interpoliert.</p>
-
-      </td>
-
-      <td nosave="">advec_particles, check_parameters, cpu_log,
+interpoliert.</p> </td> <td nosave="">advec_particles,
+check_parameters, cpu_log,
 global_min_max,
-modules&nbsp;
-      
-      <p><b>Neu:</b> <br>
-
-singleton</p>
-
-      </td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">11/06/02</td>
-
-      <td>SR</td>
-
-      <td>2.4</td>
-
-      <td>N/C</td>
-
-      <td>Optimierung und Anpassungen f&uuml;r einzelne Knoten der
-IBM-Regatta-Systeme.
-      
-      <p><b>Skalare (Cache) Optimierung:</b> Tendenzterme werden
+modules&nbsp; <p><b>Neu:</b> <br>
+singleton</p> </td> </tr> <tr nosave="" valign="top"> <td nosave="">11/06/02</td>
+<td>SR</td> <td>2.4</td> <td>N/C</td>
+<td>Optimierung und Anpassungen f&uuml;r einzelne Knoten der
+IBM-Regatta-Systeme. <p><b>Skalare (Cache) Optimierung:</b>
+Tendenzterme werden
 innerhalb einer
 gro&szlig;en (i,j)-Schleife berechnet, die nun jeweils die gesamte
@@ -2252,10 +1082,14 @@
 Gleichung umfasst. Schleifenindices i und j werden an die
 Tendenzunterprogramme
-als Argumente &uuml;bergeben. Die Tendenzunterprogramme sind als Module
-geschrieben, die mittels Technik der &uuml;berladenen Funktionen auch
+als Argumente &uuml;bergeben. Die Tendenzunterprogramme sind als
+Module
+geschrieben, die mittels Technik der &uuml;berladenen Funktionen
+auch
 ohne
-diese Indices i und j aufgerufen werden k&ouml;nnen, und dann wie in
+diese Indices i und j aufgerufen werden k&ouml;nnen, und dann wie
+in
 fr&uuml;heren
-Versionen funktionieren (d.h. in ihnen laufen die Schleifen &uuml;ber
+Versionen funktionieren (d.h. in ihnen laufen die Schleifen
+&uuml;ber
 alle
 3 Dimensionen). Solche Teile der prognostischen Gleichungen, die
@@ -2265,8 +1099,6 @@
 mit Piascek-Williams-Advektion gearbeitet wird, werden alle
 prognostischen
-Gleichungen in einer einzigen (i,j)-Schleife gerechnet. </p>
-
-      
-      <p><b>Kommunikationsoptimierung:</b> Zus&auml;tzliche
+Gleichungen in einer einzigen (i,j)-Schleife gerechnet. </p> <p><b>Kommunikationsoptimierung:</b>
+Zus&auml;tzliche
 Implementierung einer
 1D-Gebietszerlegung in x-Richtung. Mit dem neuen direkten
@@ -2279,41 +1111,38 @@
 dieser
 1D-Zerlegung durch einfache zyklische Randbedingungen ersetzt. </p>
-
-      
-      <p><b>Kleinere Anpassungen:</b> String-Vergleiche erfordern beim
+<p><b>Kleinere Anpassungen:</b> String-Vergleiche
+erfordern beim
 IBM-Compiler
-in der Regel den Einsatz der TRIM-Funktion, um &uuml;bersch&uuml;ssige
+in der Regel den Einsatz der TRIM-Funktion, um
+&uuml;bersch&uuml;ssige
 Blanks am Ende des Strings abzuschneiden (sonst liefert der Vergleich
 .F.). <br>
-
 Lokal allokierte Felder t.w. in sogenannte automatische Felder
 umgewandelt
-(g&uuml;nstiger, um Gefahr von eventuellen Speicherlecks zu vermeiden).
+(g&uuml;nstiger, um Gefahr von eventuellen Speicherlecks zu
+vermeiden).
 Einige Feldoperationen wurden wegen schlechter Performance auf der IBM
-entfernt (s. pres). </p>
-
-      
-      <p><b>Achtung: </b>Die &Uuml;bersetzung des Modells erfordert
+entfernt (s. pres). </p> <p><b>Achtung: </b>Die
+&Uuml;bersetzung des Modells erfordert
 aufgrund
 der Vielzahl nunmehr eingesetzter Module die Verwendung des
 make-Mechanismus,
-um die korrekte Reihenfolge bei der &Uuml;bersetzung der Programmteile
-zu gew&auml;hrleisten. Die Beschreibung dieser Abh&auml;ngigkeiten
+um die korrekte Reihenfolge bei der &Uuml;bersetzung der
+Programmteile
+zu gew&auml;hrleisten. Die Beschreibung dieser
+Abh&auml;ngigkeiten
 erfolgt
-in der Datei<tt> Makefile</tt>, die zusammen mit den Programmdateien
+in der Datei<tt> Makefile</tt>, die zusammen mit den
+Programmdateien
 abgelegt
-ist und ebenfalls unter RCS-Verwaltung steht. </p>
-
-      
-      <p>Im User-Interface ist das Unterprogramm<tt> user_actions </tt>nun
+ist und ebenfalls unter RCS-Verwaltung steht. </p> <p>Im
+User-Interface ist das Unterprogramm<tt> user_actions </tt>nun
 ebenfalls
-als Modul geschrieben und bedient sich der Methode der &uuml;berladenen
+als Modul geschrieben und bedient sich der Methode der
+&uuml;berladenen
 Funktionen. Existierende benutzereigene Software muss entsprechend
 angepasst
-werden.</p>
-
-      </td>
-
-      <td>advec_s_pw, advec_s_up, advec_u_pw, advec_u_up, advec_v_pw,
+werden.</p> </td> <td>advec_s_pw, advec_s_up,
+advec_u_pw, advec_u_up, advec_v_pw,
 advec_v_up,
 advec_w_pw, advec_w_up, buoyancy, calc_precipitation, calc_radiation,
@@ -2327,93 +1156,40 @@
 plot_3d, plot_spectra, poisfft, pres, production_e,
 prognostic_equations,
-user_interface,
-      
-      <p><b>Neu:</b> <br>
-
-fft_for_1d_decomp, Makefile, poisfft_hybrid,&nbsp;</p>
-
-      </td>
-
-    </tr>
-
-    <tr nosave="" valign="top">
-
-      <td nosave="">11/06/02</td>
-
-      <td>SR</td>
-
-      <td>2.4</td>
-
-      <td>B</td>
-
-      <td>Fehlende Variablenvereinbarung in <tt>advec_particles</tt>
-erg&auml;nzt.</td>
-
-      <td>advec_particles</td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top">12/09/02</td>
-
-      <td valign="top">SR</td>
-
-      <td valign="top">2.4a</td>
-
-      <td valign="top">N</td>
-
-      <td valign="top">Partikel k&ouml;nnen mit Tr&auml;gheit versehen
+user_interface, <p><b>Neu:</b> <br>
+fft_for_1d_decomp, Makefile, poisfft_hybrid,&nbsp;</p> </td>
+</tr> <tr nosave="" valign="top"> <td nosave="">11/06/02</td> <td>SR</td> <td>2.4</td>
+<td>B</td> <td>Fehlende Variablenvereinbarung in <tt>advec_particles</tt>
+erg&auml;nzt.</td> <td>advec_particles</td> </tr>
+<tr> <td valign="top">12/09/02</td> <td valign="top">SR</td> <td valign="top">2.4a</td>
+<td valign="top">N</td> <td valign="top">Partikel
+k&ouml;nnen mit Tr&auml;gheit versehen
 werden und
 einen Dichteunterschied zum umgebenden Fluid haben. Zur Steuerung dient
 der neue Paketparameter <tt>density_ratio</tt>.</td>
-
-      <td valign="top">advec_particles, header, init_particles,
-modules, package_parin</td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top">12/09/02</td>
-
-      <td valign="top">SR</td>
-
-      <td valign="top">2.4a</td>
-
-      <td valign="top">B</td>
-
-      <td valign="top">Fehler in Berechnung der Anfangsprofile von pt
+<td valign="top">advec_particles, header, init_particles,
+modules, package_parin</td> </tr> <tr> <td valign="top">12/09/02</td> <td valign="top">SR</td>
+<td valign="top">2.4a</td> <td valign="top">B</td>
+<td valign="top">Fehler in Berechnung der Anfangsprofile
+von pt
 und q entfernt
 (obere Feldgrenze von 10 konnte in <tt>pt_vertical_gradient_level_ind</tt>
-&uuml;berschritten werden). Fehler f&uuml;hrte unter Umst&auml;nden zu
-"segmentation fault" bei Verwendung sehr gro&szlig;er Gitterpunktzahlen
-in z-Richtung.
-      
-      <p>Fehler in Berechnung der v-Komponente am unteren Rand
+&uuml;berschritten werden). Fehler f&uuml;hrte unter
+Umst&auml;nden zu
+"segmentation fault" bei Verwendung sehr gro&szlig;er
+Gitterpunktzahlen
+in z-Richtung. <p>Fehler in Berechnung der v-Komponente am
+unteren Rand
 entfernt, die
-f&uuml;r den Scherungsproduktionsterm der TKE verwendet wird (Variable
-v_0). F&auml;lschlicherweise wurde zur Berechnung u(k=1) statt v(k=1)
-verwendet.</p>
-
-      </td>
-
-      <td valign="top">check_parameters, modules, production_e</td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top">19/12/02</td>
-
-      <td valign="top">SR</td>
-
-      <td valign="top">2.5</td>
-
-      <td valign="top">N</td>
-
-      <td valign="top">Restart times can be set by user with new
+f&uuml;r den Scherungsproduktionsterm der TKE verwendet wird
+(Variable
+v_0). F&auml;lschlicherweise wurde zur Berechnung u(k=1) statt
+v(k=1)
+verwendet.</p> </td> <td valign="top">check_parameters,
+modules, production_e</td> </tr> <tr> <td valign="top">19/12/02</td> <td valign="top">SR</td>
+<td valign="top">2.5</td> <td valign="top">N</td>
+<td valign="top">Restart times can be set by user with new
 runtime (d3par)
-parameters <tt>restart_time</tt> and <tt>dt_restart</tt>. Run
+parameters <tt>restart_time</tt> and <tt>dt_restart</tt>.
+Run
 description
 header is written on file CONTINUE_RUN. Output of cpu statistics
@@ -2423,40 +1199,11 @@
 Output of warnings in the job protocol in case of negative measured
 cpu-times.
-Remaining cpu-time is also evaluated on IBM-Regatta.</td>
-
-      <td valign="top">cpu_log, cpu_statistics,
+Remaining cpu-time is also evaluated on IBM-Regatta.</td> <td valign="top">cpu_log, cpu_statistics,
 check_parameters,header, leap_frog,
-local_tremain, local_tremain_ini, modules,&nbsp; palm,&nbsp;
-      
-      <p><b>new:</b> <br>
-
-check_for_restart </p>
-
-      
-      <p><b>deleted:</b> <br>
-
-check_cpu_time</p>
-
-      </td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top">&nbsp;<br>
-
-      </td>
-
-      <td valign="top">&nbsp;<br>
-
-      </td>
-
-      <td valign="top">&nbsp;<br>
-
-      </td>
-
-      <td valign="top">C</td>
-
-      <td valign="top">Unit 14 (BINOUT) must be opened using a special
+local_tremain, local_tremain_ini, modules,&nbsp; palm,&nbsp; <p><b>new:</b>
+<br>
+check_for_restart </p> <p><b>deleted:</b> <br>
+check_cpu_time</p> </td> </tr> <tr> <td valign="top">&nbsp;<br> </td> <td valign="top">&nbsp;<br> </td> <td valign="top">&nbsp;<br> </td> <td valign="top">C</td> <td valign="top">Unit
+14 (BINOUT) must be opened using a special
 process
 id string, because on IBM, the PE rank differs between communicators
@@ -2466,20 +1213,16 @@
 written on file. As the consequence, these files could not be read by
 restart-jobs.
-Reading of array <tt>hom</tt> moved from <tt>init_3d_model </tt>to <tt>read_var_list</tt>
-(binary version number incremented to 2.2).
-      
-      <p><b>Further optimization for IBM-Regatta-systems:</b> <br>
-
+Reading of array <tt>hom</tt> moved from <tt>init_3d_model
+</tt>to <tt>read_var_list</tt>
+(binary version number incremented to 2.2). <p><b>Further
+optimization for IBM-Regatta-systems:</b> <br>
 Additional optimization of the hybrid-solver for multinode usage
 (overlapping
-of communication and computation). </p>
-
-      
-      <p>Further cache optimization by using strides and joining loops (<tt>hybrid_solver,
-pres</tt> and <tt>timestep</tt>). Joining of MPI_ALLREDUCE calls in
-timestep. </p>
-
-      
-      <p>In case of 1d-decomposition along x only a part of the
+of communication and computation). </p> <p>Further cache
+optimization by using strides and joining loops (<tt>hybrid_solver,
+pres</tt> and <tt>timestep</tt>). Joining of
+MPI_ALLREDUCE calls in
+timestep. </p> <p>In case of 1d-decomposition along x
+only a part of the
 integral divisor
 conditions is checked, on IBM hosts a 1d-decomposition along x is the
@@ -2489,75 +1232,40 @@
 1d-decomposition along x, the hybrid-solver does not force a
 1d-decomposition
-any more. </p>
-
-      
-      <p>Array notation changed to do-loop constructs due to better
-performance. </p>
-
-      
-      <p>Cyclic boundary conditions along <tt>y</tt> used instead of
+any more. </p> <p>Array notation changed to do-loop
+constructs due to better
+performance. </p> <p>Cyclic boundary conditions along <tt>y</tt>
+used instead of
 sendrecv
-in case of a 1d-decomposition along <tt>x</tt>. SENDRECV replaced by
+in case of a 1d-decomposition along <tt>x</tt>. SENDRECV
+replaced by
 nonblocking
 routines ISEND and IRECV in <tt>exchange_horiz.f90</tt>. </p>
-
-      
-      <p>Speed optimization by removing MINVAL/MAXVAL calls and by
+<p>Speed optimization by removing MINVAL/MAXVAL calls and by
 handling <br>
-
 the "abs" case in a different way than the min/max cases (routine <tt>global_min_max</tt>).</p>
-
-      </td>
-
-      <td valign="top">asselin_filter, check_open, exchange_horiz,
+</td> <td valign="top">asselin_filter, check_open,
+exchange_horiz,
 global_min_max,
 init_3d_model, init_pegrid, parin, poisfft_hybrid, read_var_list,
 timestep,
-write_3d_binary, write_var_list</td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top">B</td>
-
-      <td valign="top">Correction of mixing length term (l(k)/ll(k)).
+write_3d_binary, write_var_list</td> </tr> <tr> <td valign="top"><br> </td> <td valign="top"><br>
+</td> <td valign="top"><br> </td> <td valign="top">B</td> <td valign="top">Correction
+of mixing length term (l(k)/ll(k)).
 The condition
 kh=3*km in the unstable case is now also exactly met in the wall
 adjustment
 region. Factor 0.7 in wall adjustment part replaced by variable <tt>wall_adjustment_factor</tt>,
-which is set to 1.8 in <tt>modules.f90</tt>. The factor 0.7 was the
+which is set to 1.8 in <tt>modules.f90</tt>. The factor
+0.7 was the
 possible
-reason for 2-delta-x-waves, which were observed since version 2.1d.
-      
-      <p>Calculation of deformation tensor re-designed (<tt>production_e</tt>).
-      </p>
-
-      
-      <p>STOP statements replaced by call of new subroutine <tt>local_stop,
-      </tt>where
+reason for 2-delta-x-waves, which were observed since version 2.1d. <p>Calculation
+of deformation tensor re-designed (<tt>production_e</tt>). </p>
+<p>STOP statements replaced by call of new subroutine <tt>local_stop,
+</tt>where
 MPI_FINALIZE is called before STOP in case of a parallel environment. </p>
-
-      
-      <p>ISO2D parameter dp set to true for ibm hosts. </p>
-
-      
-      <p>tend=p added in <tt>pres</tt> after calling sor method.</p>
-
-      </td>
-
-      <td valign="top">advec_s_bc, buoyancy, check_open,
+<p>ISO2D parameter dp set to true for ibm hosts. </p> <p>tend=p
+added in <tt>pres</tt> after calling sor method.</p>
+</td> <td valign="top">advec_s_bc, buoyancy,
+check_open,
 check_parameters, close_file,
 coriolis, diffusion_e, diffusivities, fft_for_1d_decomp,
@@ -2565,110 +1273,46 @@
 init_3d_model, init_dvrp, init_grid, init_particles, init_pegrid,
 parin,
-plot_2d, poisfft, pres, read_var_list, user_interface,&nbsp;
-      
-      <p><b>new:</b> <br>
-
-local_stop, production_e</p>
-
-      </td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top">03/03/03</td>
-
-      <td valign="top">SR</td>
-
-      <td valign="top">2.5a</td>
-
-      <td valign="top">N</td>
-
-      <td valign="top">Particle velocities are also stored in array
+plot_2d, poisfft, pres, read_var_list, user_interface,&nbsp; <p><b>new:</b>
+<br>
+local_stop, production_e</p> </td> </tr> <tr>
+<td valign="top">03/03/03</td> <td valign="top">SR</td>
+<td valign="top">2.5a</td> <td valign="top">N</td>
+<td valign="top">Particle velocities are also stored in
+array
 particles in
 case of zero density ratio. Steering of variables by dvrp browser
-included.</td>
-
-      <td valign="top">advec_particles, init_dvrp, modules,</td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top">C</td>
-
-      <td valign="top">AVS data format changed from float to xdr_float
+included.</td> <td valign="top">advec_particles,
+init_dvrp, modules,</td> </tr> <tr> <td valign="top"><br> </td> <td valign="top"><br>
+</td> <td valign="top"><br> </td> <td valign="top">C</td> <td valign="top">AVS
+data format changed from float to xdr_float
 (needed
 on linux machines <br>
-
 due to the little/big endian problem). Updates for new version of dvrp
-software (e.g. using module dvrp is now mandatory).</td>
-
-      <td valign="top">close_file, init_dvrp, leap_frog, plot_dvrp</td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top">B</td>
-
-      <td valign="top">Error in particle inertia part removed (exp_arg
+software (e.g. using module dvrp is now mandatory).</td> <td valign="top">close_file, init_dvrp, leap_frog, plot_dvrp</td>
+</tr> <tr> <td valign="top"><br> </td>
+<td valign="top"><br> </td> <td valign="top"><br> </td> <td valign="top">B</td>
+<td valign="top">Error in particle inertia part removed
+(exp_arg
 must not
 contain the timestep). <br>
-
 Error in calculation of the vertical flux of resolved scale energy
 (profile 57) removed. Displacement for integers in mpi_particle_type
 reduced
-from 16 to 8 on ibm.</td>
-
-      <td valign="top">advec_particles, flow_statistics, init_particles</td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top">12/03/03</td>
-
-      <td valign="top">SR</td>
-
-      <td valign="top">2.6</td>
-
-      <td valign="top">N</td>
-
-      <td valign="top"><b>Version optimized for NEC-SX6 parallel-vector
-machines.&nbsp;</b>
-      
-      <p>There are two main changes. A new vectorizable routine <tt>prognostic_equations_vec</tt>
+from 16 to 8 on ibm.</td> <td valign="top">advec_particles,
+flow_statistics, init_particles</td> </tr> <tr> <td valign="top">12/03/03</td> <td valign="top">SR</td>
+<td valign="top">2.6</td> <td valign="top">N</td>
+<td valign="top"><b>Version optimized for NEC-SX6
+parallel-vector
+machines.&nbsp;</b> <p>There are two main changes. A
+new vectorizable routine <tt>prognostic_equations_vec</tt>
 is added, where cache optimizations are undone. The tendency
 subroutines
-called by <tt>prognostic_equations_vec</tt> contain the full 3d-loops
+called by <tt>prognostic_equations_vec</tt> contain the
+full 3d-loops
 (compared
 with the cache-optimized versions where only loops over k are carried
 out).
-Additionally, a new pressure solver <tt>poisfft_hybrid_vec</tt> is
+Additionally, a new pressure solver <tt>poisfft_hybrid_vec</tt>
+is
 added,
 which requires a 1d-domain-decomposition (like poisfft_hybrid). In this
@@ -2676,29 +1320,18 @@
 optionally
 available by choosing <tt>fft_method </tt>= <i>'temperton-algorithm'</i>.&nbsp;
-      </p>
-
-      
-      <p>Additional changes in <tt>flow_statistics</tt> to allow
-better vectorization. </p>
-
-      
-      <p>Small changes in routines <tt>diffusivities</tt> and <tt>pres</tt>
+</p> <p>Additional changes in <tt>flow_statistics</tt>
+to allow
+better vectorization. </p> <p>Small changes in routines <tt>diffusivities</tt>
+and <tt>pres</tt>
 which
 caused run time errors on IBM and NEC due to compiler problems. </p>
-
-      
-      <p>Reading of environment variable <tt>tasks_per_node</tt> moved
+<p>Reading of environment variable <tt>tasks_per_node</tt>
+moved
 from routine <tt>poisfft_hybrid</tt>
-to routine <tt>parin</tt>. </p>
-
-      
-      <p><b>Changes in makefile:</b> modules and user_interface now
+to routine <tt>parin</tt>. </p> <p><b>Changes
+in makefile:</b> modules and user_interface now
 depend on
 the f90 files, dependency of singleton added, LDFLAGS moved to the end
-of PROG rule, temperton_fft added.</p>
-
-      </td>
-
-      <td valign="top">advec_s_pw, advec_s_up, advec_u_pw, advec_u_up,
+of PROG rule, temperton_fft added.</p> </td> <td valign="top">advec_s_pw, advec_s_up, advec_u_pw, advec_u_up,
 advec_v_pw,
 advec_v_up, advec_w_pw, advec_w_up, buoyancy, calc_precipitation,
@@ -2710,195 +1343,69 @@
 init_pegrid, leap_frog, local_tremain, modules, parin, poisfft,
 poisfft_hybrid,
-pres, production_e, prognostic_equations, user_interface
-      
-      <p><b>new:</b> <br>
-
-temperton_fft</p>
-
-      </td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top">B</td>
-
-      <td valign="top">tend=p added in routine <tt>pres</tt> after
+pres, production_e, prognostic_equations, user_interface <p><b>new:</b>
+<br>
+temperton_fft</p> </td> </tr> <tr> <td valign="top"><br> </td> <td valign="top"><br>
+</td> <td valign="top"><br> </td> <td valign="top">B</td> <td valign="top">tend=p
+added in routine <tt>pres</tt> after
 calling sor
-method.</td>
-
-      <td valign="top">pres</td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top">14/03/03</td>
-
-      <td valign="top">SR</td>
-
-      <td valign="top">2.6a</td>
-
-      <td valign="top">N</td>
-
-      <td valign="top">Additional vector optimization for NEC-SX6.
-      
-      <p>Optional system-specific random number generator available
+method.</td> <td valign="top">pres</td> </tr>
+<tr> <td valign="top">14/03/03</td> <td valign="top">SR</td> <td valign="top">2.6a</td>
+<td valign="top">N</td> <td valign="top">Additional
+vector optimization for NEC-SX6. <p>Optional system-specific
+random number generator available
 (for equal
-distributed numbers). </p>
-
-      
-      <p>Define strings in all relevant subroutines changed from&nbsp; </p>
-
-      
-      <p><tt>&nbsp;&nbsp;&nbsp;&nbsp; #if defined( define_string )</tt>
-      </p>
-
-      
-      <p>to </p>
-
-      
-      <p><tt>&nbsp;&nbsp;&nbsp;&nbsp; #if defined( __define_string )</tt>
-      </p>
-
-      
-      <p><b>Caution: </b>The relevant subroutines are NOT listed in
+distributed numbers). </p> <p>Define strings in all
+relevant subroutines changed from&nbsp; </p> <p><tt>&nbsp;&nbsp;&nbsp;&nbsp;
+#if defined( define_string )</tt> </p> <p>to </p>
+<p><tt>&nbsp;&nbsp;&nbsp;&nbsp; #if
+defined( __define_string )</tt> </p> <p><b>Caution:
+</b>The relevant subroutines are NOT listed in
 the right
-column!</p>
-
-      </td>
-
-      <td valign="top">asselin_filter, check_parameters, diffusion_e,
+column!</p> </td> <td valign="top">asselin_filter,
+check_parameters, diffusion_e,
 diffusivities,
 disturb_field, global_min_max, header, local_tremain, modules, parin,
 poisfft_hybrid,
-pres, read_var_list, write_var_list</td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top">B</td>
-
-      <td valign="top">Error in particle boundary condition removed
+pres, read_var_list, write_var_list</td> </tr> <tr>
+<td valign="top"><br> </td> <td valign="top"><br> </td> <td valign="top"><br>
+</td> <td valign="top">B</td> <td valign="top">Error in particle boundary condition removed
 (velocity must
-be inverted in case of reflection)</td>
-
-      <td valign="top">advec_particles</td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top">16/04/03</td>
-
-      <td valign="top">SR</td>
-
-      <td valign="top">2.6b</td>
-
-      <td valign="top">N</td>
-
-      <td valign="top">Time series output of Monin Obukhov length.
-      
-      <p>Temperton fft can now be used for all hosts and every domain
+be inverted in case of reflection)</td> <td valign="top">advec_particles</td>
+</tr> <tr> <td valign="top">16/04/03</td>
+<td valign="top">SR</td> <td valign="top">2.6b</td>
+<td valign="top">N</td> <td valign="top">Time
+series output of Monin Obukhov length. <p>Temperton fft can now
+be used for all hosts and every domain
 decomposition.
 Abort in case of Temperton fft, if number of gridpoints along x and/or
-y contain illegal factors.&nbsp; </p>
-
-      
-      <p>Index values for the extrema found in global_min_max are
+y contain illegal factors.&nbsp; </p> <p>Index values
+for the extrema found in global_min_max are
 limited to
 the range 0..nx, 0..ny (on IBM machines -1 and nx+1/ny+1 occured which
 produced different RUN_CONTROL output compared to other machines). </p>
-
-      
-      <p>Output format of iteration count in routine run_control
-enlarged.</p>
-
-      </td>
-
-      <td valign="top">check_parameters, fft_for_1d_decomp,
+<p>Output format of iteration count in routine run_control
+enlarged.</p> </td> <td valign="top">check_parameters,
+fft_for_1d_decomp,
 global_min_max, init_pegrid,
-modules, plot_ts, poisfft, run_control, temperton_fft</td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top">B</td>
-
-      <td valign="top">Additional checks of variables <tt>hybrid_solver</tt>
+modules, plot_ts, poisfft, run_control, temperton_fft</td> </tr>
+<tr> <td valign="top"><br> </td> <td valign="top"><br> </td> <td valign="top"><br>
+</td> <td valign="top">B</td> <td valign="top">Additional checks of variables <tt>hybrid_solver</tt>
 and <tt>host</tt>
-(<tt>hybrid_solver</tt> must not be used for a 2d-domain-decomposition,
-      <tt>host</tt>
+(<tt>hybrid_solver</tt> must not be used for a
+2d-domain-decomposition, <tt>host</tt>
 should be set by setting environment variable <tt>localhost</tt>).
-      
-      <p>Header output for mixing length limitations revised.</p>
-
-      </td>
-
-      <td valign="top">check_parameters, header</td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top">09/05/03</td>
-
-      <td valign="top">SR</td>
-
-      <td valign="top">2.7</td>
-
-      <td valign="top">N</td>
-
-      <td valign="top"><b>Version running on Linux Clusters using MPICH
+<p>Header output for mixing length limitations revised.</p>
+</td> <td valign="top">check_parameters, header</td>
+</tr> <tr> <td valign="top">09/05/03</td>
+<td valign="top">SR</td> <td valign="top">2.7</td>
+<td valign="top">N</td> <td valign="top"><b>Version
+running on Linux Clusters using MPICH
 and Intel
-FORTRAN compiler (ifc)</b>
-      
-      <p>So far, only absolutely neccessary changes have been done
+FORTRAN compiler (ifc)</b> <p>So far, only absolutely
+neccessary changes have been done
 (possible
-optimizations will follow in a later version): </p>
-
-      
-      <p>New time measurements for Linux (ifc) environment added. </p>
-
-      
-      <p>Some parameters in MPI calls had to be modified (arrays had to
+optimizations will follow in a later version): </p> <p>New
+time measurements for Linux (ifc) environment added. </p> <p>Some
+parameters in MPI calls had to be modified (arrays had to
 be replaced
 by the first element of the regarding array) in order to fulfill f90
@@ -2906,5 +1413,6 @@
 and rank requirements, since on the MUK-cluster a FORTRAN90-version of
 MPI is used. This also required to replace "mpif.h" by using a special
-module (named <tt>mpi</tt>).&nbsp; I did not find out how character
+module (named <tt>mpi</tt>).&nbsp; I did not find out
+how character
 strings
 can be send with MPI. Therefore, these strings are transformed to
@@ -2912,357 +1420,156 @@
 before they are send, and transformed back to characters on the
 receiving
-PE. </p>
-
-      
-      <p>On the MUK-cluster, only PE0 is able to read the values of
+PE. </p> <p>On the MUK-cluster, only PE0 is able to read
+the values of
 environment
 variables. Therefore, these values are communicated via broadcast to
 the
-other PEs. </p>
-
-      
-      <p>1d-decomposition is set as the default on Linux Clusters.
+other PEs. </p> <p>1d-decomposition is set as the default
+on Linux Clusters.
 Character
-strings have to be transformed to integer </p>
-
-      
-      <p>Smaller changes: </p>
-
-      
-      <p>Batch_scp paths for IBM, NEC and Linux added</p>
-
-      </td>
-
-      <td valign="top">check_open, cpu_log, flow_statistics,
+strings have to be transformed to integer </p> <p>Smaller
+changes: </p> <p>Batch_scp paths for IBM, NEC and Linux
+added</p> </td> <td valign="top">check_open,
+cpu_log, flow_statistics,
 init_pegrid, local_getenv,
-modules, palm, poisfft_hybrid</td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top">B</td>
-
-      <td valign="top">Measurements on IBM are now using function <tt>irtc</tt>,
+modules, palm, poisfft_hybrid</td> </tr> <tr> <td valign="top"><br> </td> <td valign="top"><br>
+</td> <td valign="top"><br> </td> <td valign="top">B</td> <td valign="top">Measurements
+on IBM are now using function <tt>irtc</tt>,
 which allows correct measurements for jobs running over the 24:00
-timeline.</td>
-
-      <td valign="top">cpu_log, local_tremain, local_tremain_ini</td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top">01/08/03</td>
-
-      <td valign="top">SR</td>
-
-      <td valign="top">2.7a</td>
-
-      <td valign="top">B</td>
-
-      <td valign="top">Check that the number processors is also an
+timeline.</td> <td valign="top">cpu_log,
+local_tremain, local_tremain_ini</td> </tr> <tr> <td valign="top">01/08/03</td> <td valign="top">SR</td>
+<td valign="top">2.7a</td> <td valign="top">B</td>
+<td valign="top">Check that the number processors is also
+an
 integral divisor
-of the number&nbsp; of gridpoints along y in case of a 1d-decomposition
-along x.
-      
-      <p>Error concerning the multinode-version of poisfft_hybrid
+of the number&nbsp; of gridpoints along y in case of a
+1d-decomposition
+along x. <p>Error concerning the multinode-version of
+poisfft_hybrid
 removed. It
 was caused by the Linux changes of version 2.7 and led to program abort
-due to MPI errors in <tt>MPI_ALLTOALL</tt>. </p>
-
-      
-      <p>No abort on t3e-systems if system specific routines are used. </p>
-
-      
-      <p>Array<tt> sums </tt>is initialized in<tt> init_3d_model </tt>before
-the first call of<tt> pres</tt>. Call of<tt> init_cloud_physics </tt>moved
-before<tt> init_particles</tt>.</p>
-
-      </td>
-
-      <td valign="top">fft_1dd_init, init_3d_model, init_pegrid,
-poisfft_hybrid</td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top">N</td>
-
-      <td valign="top">Interactive steering of dvrp-graphic features by
+due to MPI errors in <tt>MPI_ALLTOALL</tt>. </p> <p>No
+abort on t3e-systems if system specific routines are used. </p> <p>Array<tt>
+sums </tt>is initialized in<tt> init_3d_model </tt>before
+the first call of<tt> pres</tt>. Call of<tt>
+init_cloud_physics </tt>moved
+before<tt> init_particles</tt>.</p> </td> <td valign="top">fft_1dd_init, init_3d_model, init_pegrid,
+poisfft_hybrid</td> </tr> <tr> <td valign="top"><br>
+</td> <td valign="top"><br> </td> <td valign="top"><br> </td> <td valign="top">N</td>
+<td valign="top">Interactive steering of dvrp-graphic
+features by
 dvrp-browser
 plugin extended (position of slicers, threshold values of isosurfaces)
 . New dvrp_graphics package parameter<tt> slicer_range_limits_dvrp</tt>.</td>
-
-      <td valign="top">init_dvrp, modules, package_parin, plot_dvrp</td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top">29/10/03</td>
-
-      <td valign="top">SR</td>
-
-      <td valign="top">2.7b</td>
-
-      <td valign="top">N</td>
-
-      <td valign="top">In the multigrid method, on a defined level,
+<td valign="top">init_dvrp, modules, package_parin,
+plot_dvrp</td> </tr> <tr> <td valign="top">29/10/03</td>
+<td valign="top">SR</td> <td valign="top">2.7b</td>
+<td valign="top">N</td> <td valign="top">In
+the multigrid method, on a defined level,
 data are gathered
 on PE0 and further calculations are carried out only on this PE. New<tt>
 d3par </tt>parameter <tt>mg_switch_to_pe0_level</tt>.
-      
-      <p>Particle groups implemented. New data type<tt>
+<p>Particle groups implemented. New data type<tt>
 particle_groups_type </tt>and data type<tt> particle_type </tt>modified.
 Version
 numbers are
 output on the particle files. New<tt> particles_par </tt>parameter<tt>
-diameter</tt>. Parameter name<tt> uniform_psize</tt> changed to<tt>
+diameter</tt>. Parameter name<tt> uniform_psize</tt>
+changed to<tt>
 dvrp_psize</tt>.
 Parameters<tt> density_ratio </tt>is now an array.</p>
-
-      </td>
-
-      <td valign="top">advec_particles, check_open, exchange_horiz,
+</td> <td valign="top">advec_particles, check_open,
+exchange_horiz,
 header, init_particles,
-init_pegrid, modules, package_parin, parin, poismg</td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top">C</td>
-
-      <td valign="top">Random number generator from numerical recipes
+init_pegrid, modules, package_parin, parin, poismg</td> </tr>
+<tr> <td valign="top"><br> </td> <td valign="top"><br> </td> <td valign="top"><br>
+</td> <td valign="top">C</td> <td valign="top">Random number generator from numerical recipes
 is now rewritten
 as a module. In restart runs the seed of this generator is set to the
 last
 values of the previous run in order to keep the sequence of the random
-numbers. New binary version 2.2.
-      
-      <p>Buffer is flushed for file RUN_CONTROL immediately after every
+numbers. New binary version 2.2. <p>Buffer is flushed for file
+RUN_CONTROL immediately after every
 output
-on IBM and Linux cluster.</p>
-
-      </td>
-
-      <td valign="top">advec_particles, disturb_field, init_3d_model,
+on IBM and Linux cluster.</p> </td> <td valign="top">advec_particles,
+disturb_field, init_3d_model,
 init_particles,
-random_function, random_gauss, run_control, write_3d_binary</td>
-
-    </tr>
-
-    <tr>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top"><br>
-
-      </td>
-
-      <td valign="top">E</td>
-
-      <td valign="top">In routine buoyance, horizontal mean temperature
+random_function, random_gauss, run_control, write_3d_binary</td> </tr>
+<tr> <td valign="top"><br> </td> <td valign="top"><br> </td> <td valign="top"><br>
+</td> <td valign="top">E</td> <td valign="top">In routine buoyance, horizontal mean
+temperature
 is now
 taken from array hom instead of array sums (otherwise inconsistence in
-case of using more than one statistical region).
-      
-      <p>Information about the time of the last timestep change added
+case of using more than one statistical region). <p>Information
+about the time of the last timestep change added
 to the
-restart file. </p>
-
-      
-      <p>System call of<tt> batch_scp </tt>on IBM corrected.</p>
-
-      </td>
-
-      <td valign="top">buoyancy, check_open, read_var_list,
-write_var_list</td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;">28/01/04</td>
-
-      <td style="vertical-align: top;">SR</td>
-
-      <td style="vertical-align: top;">2.8</td>
-
-      <td style="vertical-align: top;">N</td>
-
-      <td style="vertical-align: top;">Runge-Kutta schemes (2nd and 3rd
+restart file. </p> <p>System call of<tt> batch_scp </tt>on
+IBM corrected.</p> </td> <td valign="top">buoyancy,
+check_open, read_var_list,
+write_var_list</td> </tr> <tr> <td style="vertical-align: top;">28/01/04</td> <td style="vertical-align: top;">SR</td> <td style="vertical-align: top;">2.8</td> <td style="vertical-align: top;">N</td> <td style="vertical-align: top;">Runge-Kutta schemes (2nd and
+3rd
 order)
 for time integration implemented. The third order scheme is the new
 default.
 It allows the timestep to be 0.9*CFL, which is much larger than for the
-former default leapfrog scheme.&nbsp; The intermediate steps, which are
+former default leapfrog scheme.&nbsp; The intermediate steps, which
+are
 part of these schemes, are realized by an additional loop within routine<span style="font-family: monospace;">
 time_integration&nbsp;</span>(currently<span style="font-family: monospace;">
-leap_frog</span>, but this routine will be renamed soon). Steering
+leap_frog</span>, but this routine will be renamed soon).
+Steering
 variables<span style="font-family: monospace;">
-at&nbsp;</span>and<span style="font-family: monospace;"> bt&nbsp;</span>have
-been replaced by array<span style="font-family: monospace;"> sct</span>.
+at&nbsp;</span>and<span style="font-family: monospace;">
+bt&nbsp;</span>have
+been replaced by array<span style="font-family: monospace;">
+sct</span>.
 , which is particularly used in the prognostic equations. Values of
 this
 steering array are calculated within the new routine<span style="font-family: monospace;">
-timestep_scheme_steering</span>.
-      
-      <p>When using Runge-Kutta schemes,&nbsp; the timestep increment
+timestep_scheme_steering</span>. <p>When using Runge-Kutta
+schemes,&nbsp; the timestep increment
 is freely
 allowed to adjust after each timestep (the older schemes are using some
-restrictions, see routine<span style="font-family: monospace;"> timestep</span>).
-Also, routine<span style="font-family: monospace;"> run_control&nbsp;</span>is
+restrictions, see routine<span style="font-family: monospace;">
+timestep</span>).
+Also, routine<span style="font-family: monospace;">
+run_control&nbsp;</span>is
 not automatically called in case of timestep changes, when Runge-Kutta
-schemes are switched on. </p>
-
-      
-      <p>The old leapfrog scheme is still implemented and should
+schemes are switched on. </p> <p>The old leapfrog scheme
+is still implemented and should
 produce the
-same results as in the previous version(s)! </p>
-
-      
-      <p>Depending on the user experience, further adjustments to the
+same results as in the previous version(s)! </p> <p>Depending
+on the user experience, further adjustments to the
 Runge-Kutta
 schemes may be necessary within the next minor versions.&nbsp; </p>
-
-      
-      <p><span style="font-weight: bold;">Attention:</span> <br>
-
+<p><span style="font-weight: bold;">Attention:</span>
+<br>
 The content of the testresult-file (run control output) has
 significantly
-changed due to the new default timestep scheme!</p>
-
-      </td>
-
-      <td style="vertical-align: top;">check_parameters, header,
+changed due to the new default timestep scheme!</p> </td> <td style="vertical-align: top;">check_parameters, header,
 init_3d_model,
 init_rankine, leap_frog, modules, prognostic_equations, read_var_list,
-swap_timelevel, timestep, write_var_list
-      
-      <p><span style="font-weight: bold;">new:</span> <br>
-
-timestep_scheme_steering</p>
-
-      </td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;">C</td>
-
-      <td style="vertical-align: top;">Output of particle infos in
+swap_timelevel, timestep, write_var_list <p><span style="font-weight: bold;">new:</span> <br>
+timestep_scheme_steering</p> </td> </tr> <tr>
+<td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">C</td> <td style="vertical-align: top;">Output of particle infos in
 subroutine<span style="font-family: monospace;">
-allocate_prt_memory&nbsp;</span>on demand only.
-      
-      <p>Type<span style="font-family: monospace;"> log&nbsp;</span>(used
+allocate_prt_memory&nbsp;</span>on demand only. <p>Type<span style="font-family: monospace;"> log&nbsp;</span>(used
 for
 cpu time measurements) changed to<span style="font-family: monospace;">
-logpoint&nbsp;</span>due to name conflict with FORTRAN intrinsic<span style="font-family: monospace;">
-log.</span></p>
-
-      </td>
-
-      <td style="vertical-align: top;">advec_particles, cpu_log,&nbsp;</td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;">E</td>
-
-      <td style="vertical-align: top;">In case of opening unit 80, a
+logpoint&nbsp;</span>due to name conflict with FORTRAN
+intrinsic<span style="font-family: monospace;">
+log.</span></p> </td> <td style="vertical-align: top;">advec_particles,
+cpu_log,&nbsp;</td> </tr> <tr> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">E</td> <td style="vertical-align: top;">In case of opening unit 80, a
 barrier
 is set only for the first call of<span style="font-family: monospace;">
 check_open&nbsp;</span>(from routine<span style="font-family: monospace;">
-init_particles</span>), in order to avoid the possibility of hanging
+init_particles</span>), in order to avoid the possibility of
+hanging
 jobs,
 which may occur if unit 80 is opened within routine<span style="font-family: monospace;">
-allocate_prt_memory</span>.</td>
-
-      <td style="vertical-align: top;">check_open</td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;">30/01/04</td>
-
-      <td style="vertical-align: top;">SR</td>
-
-      <td style="vertical-align: top;">2.8a</td>
-
-      <td style="vertical-align: top;">N</td>
-
-      <td style="vertical-align: top;">In order to prepare the code for
+allocate_prt_memory</span>.</td> <td style="vertical-align: top;">check_open</td> </tr>
+<tr> <td style="vertical-align: top;">30/01/04</td>
+<td style="vertical-align: top;">SR</td> <td style="vertical-align: top;">2.8a</td> <td style="vertical-align: top;">N</td> <td style="vertical-align: top;">In order to prepare the code
+for
 the simulation
 of flow around buildings (to be realized in one of the next major
@@ -3270,15 +1577,12 @@
 the lower<span style="font-family: monospace;"> k&nbsp;</span>index<span style="font-family: monospace;">
 nzb&nbsp;</span>has been replaced by a two-dimensional array<span style="font-family: monospace;">
-nzb_2d&nbsp;</span>in many of the three-dimensional loops. So far, all
+nzb_2d&nbsp;</span>in many of the three-dimensional loops. So
+far, all
 elements of this array are set to<span style="font-family: monospace;">
-nzb</span>.
-      
-      <p>This change should not effect simulation results and cpu time
+nzb</span>. <p>This change should not effect simulation
+results and cpu time
 in any
-way.</p>
-
-      </td>
-
-      <td style="vertical-align: top;">advec_s_pw, advec_s_up,
+way.</p> </td> <td style="vertical-align: top;">advec_s_pw,
+advec_s_up,
 advec_u_pw, advec_u_up,
 advec_v_pw, advec_v_up, advec_w_pw, advec_w_up, asselin_filter,
@@ -3289,87 +1593,25 @@
 diffusivities,
 impact_of_latent_heat, init_3d_model, modules, production_e,
-prognostic_equations,&nbsp;</td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;">C</td>
-
-      <td style="vertical-align: top;">Euler step informations are
+prognostic_equations,&nbsp;</td> </tr> <tr> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">C</td> <td style="vertical-align: top;">Euler step informations are
 removed in
 the run control output in case of Runge-Kutta schemes. Also, the
 timestep
-increment limitation is removed for the Runge-Kutta schemes.
-      
-      <p><span style="font-weight: bold;">Attention:</span> <br>
-
+increment limitation is removed for the Runge-Kutta schemes. <p><span style="font-weight: bold;">Attention:</span> <br>
 The content of the testresult-file (run control output) has changed
-due to this modification.</p>
-
-      </td>
-
-      <td style="vertical-align: top;">run_control, timestep</td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;">E</td>
-
-      <td style="vertical-align: top;">Velocity gradients at the
+due to this modification.</p> </td> <td style="vertical-align: top;">run_control, timestep</td>
+</tr> <tr> <td style="vertical-align: top;"><br>
+</td> <td style="vertical-align: top;"><br> </td>
+<td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">E</td> <td style="vertical-align: top;">Velocity gradients at the
 surface calculated
 in routine production_e are now limited. In case of inhomogeneous
 surface
 heating small diffusivities (<span style="font-family: monospace;">km</span>)
-sometimes caused very small timesteps due to this problem.</td>
-
-      <td style="vertical-align: top;">production_e</td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;">30/04/04<br>
-
-      </td>
-
-      <td style="vertical-align: top;">SR<br>
-
-      </td>
-
-      <td style="vertical-align: top;">2.8b<br>
-
-      </td>
-
-      <td style="vertical-align: top;">N<br>
-
-      </td>
-
-      <td style="vertical-align: top;">The number of processors along
+sometimes caused very small timesteps due to this problem.</td> <td style="vertical-align: top;">production_e</td> </tr>
+<tr> <td style="vertical-align: top;">30/04/04<br>
+</td> <td style="vertical-align: top;">SR<br>
+</td> <td style="vertical-align: top;">2.8b<br>
+</td> <td style="vertical-align: top;">N<br> </td>
+<td style="vertical-align: top;">The number of processors
+along
 the respective directions of the virtual processor grid must not be
 integral divisors of the number of gridpoints along x, y or z any more.
@@ -3379,28 +1621,28 @@
 Although this causes a load imbalance, only runs with very small
 numbers of processors will feel a significant decrease of performance.
-In order to be able to carry out the transpositions, array<span style="font-family: monospace;"> d </span>(and also working array<span style="font-family: monospace;"> tend</span>) have the same size on
+In order to be able to carry out the transpositions, array<span style="font-family: monospace;"> d </span>(and also
+working array<span style="font-family: monospace;"> tend</span>)
+have the same size on
 all processors. This new feature is switched on by setting the
-inipar-parameter<span style="font-family: monospace;"> grid_matching </span>=
-      <span style="font-style: italic;">'match'</span>, which is the
-default now.<br>
-
-      <br>
-
-Routine<span style="font-family: monospace;"> poisfft </span>is now
+inipar-parameter<span style="font-family: monospace;">
+grid_matching </span>= <span style="font-style: italic;">'match'</span>,
+which is the
+default now.<br> <br>
+Routine<span style="font-family: monospace;"> poisfft </span>is
+now
 the
 default for solving the Poisson-equation for the case of 2d- as well as
 1d-domain-decompositions and has been changed from a subroutine to a
 module. This routine now includes most of the optimization from the
-(Ketelsen) routine<span style="font-family: monospace;"> poisfft_hybrid
-      </span>and additionally allows 1d-decompositions along y.<br>
-
-      <span style="font-weight: bold;">ATTENTION:</span> Routine<span style="font-family: monospace;"> poisfft_hybrid </span>might still be
+(Ketelsen) routine<span style="font-family: monospace;">
+poisfft_hybrid </span>and additionally allows 1d-decompositions
+along y.<br> <span style="font-weight: bold;">ATTENTION:</span>
+Routine<span style="font-family: monospace;">
+poisfft_hybrid </span>might still be
 faster when running PALM on more than one node of the IBM-Regatta,
-although the difference is not very big. The <span style="font-family: monospace;">poisfft</span>-code for a
+although the difference is not very big. The <span style="font-family: monospace;">poisfft</span>-code
+for a
 decomposition along y has still a poor performance and needs further
-optimization.<br>
-
-      <br>
-
+optimization.<br> <br>
 The spectra-package is now able to use all the available FFT-methods
 (selected by inipar-parameter<span style="font-family: monospace;">
@@ -3408,156 +1650,55 @@
 Performance of the transpositions needed for calculating the spectra
 has been improved by using new direct transpositions from z to y and
-from y to x instead of two subsequent transpositions.<br>
-
-      <br>
-
-Runge-Kutta scheme implemented in the 1d-version (<span style="font-family: monospace;">init_1d_model</span>). Diffusion
+from y to x instead of two subsequent transpositions.<br> <br>
+Runge-Kutta scheme implemented in the 1d-version (<span style="font-family: monospace;">init_1d_model</span>).
+Diffusion
 quantity arrays are allocated for previous timelevels only in case of
-leapfrog scheme.<br>
-
-      <br>
-
+leapfrog scheme.<br> <br>
 This version is adapted for use on COMPAQ-DECALPHA systems (main
 changes: calling of SHAPE function within subroutine arguments (e.g.<span style="font-family: monospace;"> singleton</span>)
 causes compiler errors, particle data type modified for better
 alignment, DATA statements from modules removed because they caused
-internal compiler errors).<br>
-
-      </td>
-
-      <td style="vertical-align: top;">advec_s_ups, advec_u_ups,
+internal compiler errors).<br> </td> <td style="vertical-align: top;">advec_s_ups, advec_u_ups,
 advec_v_ups, advec_w_ups, calc_spectra, check_parameters, fft_xy,
 header, init_1d_model, init_3d_model, init_particles, init_pegrid,
 modules, parin, plot_2d, poisfft, prandtl_fluxes, pres, read_var_list,
 singleton, spline_x, spline_y, spline_z, time_integration, transpose,
-write_3d_binary, write_var_list<br>
-
-      </td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;">C<br>
-
-      </td>
-
-      <td style="vertical-align: top;">FFT routines renamed. Module<span style="font-family: monospace;"> fft_for_1d_decomp </span>renamed<span style="font-family: monospace;"> fft_xy</span>. Subroutine<span style="font-family: monospace;"> leap_frog </span>renamed<span style="font-family: monospace;"> time_integration</span>.<br>
-
-      <br>
-
+write_3d_binary, write_var_list<br> </td> </tr> <tr>
+<td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">C<br> </td> <td style="vertical-align: top;">FFT routines renamed. Module<span style="font-family: monospace;"> fft_for_1d_decomp </span>renamed<span style="font-family: monospace;"> fft_xy</span>.
+Subroutine<span style="font-family: monospace;"> leap_frog
+</span>renamed<span style="font-family: monospace;">
+time_integration</span>.<br> <br>
 For all transpositions, MPI_ALLTOALL is used instead of MPI_ALLTOALLV.
-The names of many transposition indices have been changed.<br>
-
-      <br>
-
-Inipar-parameter<span style="font-family: monospace;"> impulse_advec </span>renamed<span style="font-family: monospace;"> momentum_advec</span>.<br>
-
-      <br>
-
-Reading of restart data is moved from<span style="font-family: monospace;"> init_3d_model </span>to new routine <span style="font-family: monospace;">read_3d_binary</span>.<br>
-
-      <br>
-
-All former transpose routines collected into one file.<br>
-
-      </td>
-
-      <td style="vertical-align: top;">check_parameters,
+The names of many transposition indices have been changed.<br> <br>
+Inipar-parameter<span style="font-family: monospace;">
+impulse_advec </span>renamed<span style="font-family: monospace;"> momentum_advec</span>.<br>
+<br>
+Reading of restart data is moved from<span style="font-family: monospace;"> init_3d_model </span>to
+new routine <span style="font-family: monospace;">read_3d_binary</span>.<br>
+<br>
+All former transpose routines collected into one file.<br> </td>
+<td style="vertical-align: top;">check_parameters,
 flow_statistics, header, init_1d_model, init_3d_model, init_advec,
 init_pegrid, palm, parin, poisfft, poisfft_hybrid, production_e,
-prognostic_equations, read_var_list, sline_x, spline_y, write_var_list<span style="font-weight: bold;"><br>
-
-      <br>
-
+prognostic_equations, read_var_list, sline_x, spline_y, write_var_list<span style="font-weight: bold;"><br> <br>
 new:</span><br>
-
-fft_xy, read_3d_binary, time_integration, transpose<br>
-
-      <br>
-
-      <span style="font-weight: bold;">deleted:</span><br>
-
+fft_xy, read_3d_binary, time_integration, transpose<br> <br>
+<span style="font-weight: bold;">deleted:</span><br>
 fft_for_1d_decomp, leap_frog, transpose_xy, transpose_xz, transpose_yx,
-transpose_yz, transpose_zx, transpose_zy<br>
-
-      </td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;">E<br>
-
-      </td>
-
-      <td style="vertical-align: top;">Small errors in the
+transpose_yz, transpose_zx, transpose_zy<br> </td> </tr>
+<tr> <td style="vertical-align: top;"><br> </td>
+<td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">E<br> </td> <td style="vertical-align: top;">Small errors in the
 particle-package removed, which sometimes caused aborts of restart runs
-on IBM.<br>
-
-      <br>
-
+on IBM.<br> <br>
 Small error in calculating the time series profiles removed (array
 bounds of some arrays needed for profil-output had to be increased).<br>
-
-      <br>
-
+<br>
 Function fft (singleton method) replaced by subroutine fftn due to
-problems with 64-bit mode on IBM<br>
-
-      </td>
-
-      <td style="vertical-align: top;">advec_particles,
-check_parameters, <br>
-
-      </td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;">26/03/05<br>
-
-      </td>
-
-      <td style="vertical-align: top;">SR<br>
-
-      </td>
-
-      <td style="vertical-align: top;">2.9<br>
-
-      </td>
-
-      <td style="vertical-align: top;">N<br>
-
-      </td>
-
-      <td style="vertical-align: top;">Non-cyclic lateral boundary
+problems with 64-bit mode on IBM<br> </td> <td style="vertical-align: top;">advec_particles,
+check_parameters, <br> </td> </tr> <tr> <td style="vertical-align: top;">26/03/05<br> </td>
+<td style="vertical-align: top;">SR<br> </td>
+<td style="vertical-align: top;">2.9<br> </td>
+<td style="vertical-align: top;">N<br> </td> <td style="vertical-align: top;">Non-cyclic lateral boundary
 conditions implemented.<br>
-
 Dirichlet/Neumann conditions at the inflow/outflow can be used either
 along x or along y. Cyclic conditions have to be used along the
@@ -3565,28 +1706,28 @@
 upper bound nxl (nyn) of array u (v) is enlarged by one gridpoint on
 the rightmost (northest) processor. These array index informations are
-passed to routine<span style="font-family: monospace;"> exchange_horiz </span>and<span style="font-family: monospace;"> disturb_field </span>by two
+passed to routine<span style="font-family: monospace;">
+exchange_horiz </span>and<span style="font-family: monospace;">
+disturb_field </span>by two
 additional arguments. Boundary conditions are set with new
-inipar-parameters<span style="font-family: monospace;"> bc_lr </span>and<span style="font-family: monospace;"> bc_ns</span>. A damping layer is
-needed at the outflow in order to damp reflections (see new parameters<span style="font-family: monospace;"> km_damp_max </span>and<span style="font-family: monospace;"> outflow_damping_width</span>) and an
-additional disturbance is imposed near the inflow (see new parameters<span style="font-family: monospace;">&nbsp; inflow_disturbance_begin </span>and<span style="font-family: monospace;"> inflow_disturbance_end</span>).
+inipar-parameters<span style="font-family: monospace;">
+bc_lr </span>and<span style="font-family: monospace;">
+bc_ns</span>. A damping layer is
+needed at the outflow in order to damp reflections (see new parameters<span style="font-family: monospace;"> km_damp_max </span>and<span style="font-family: monospace;"> outflow_damping_width</span>)
+and an
+additional disturbance is imposed near the inflow (see new parameters<span style="font-family: monospace;">&nbsp;
+inflow_disturbance_begin </span>and<span style="font-family: monospace;"> inflow_disturbance_end</span>).
 Non-cyclic conditions require the use of the multigrid-method for
-solving the Poisson equation. Argument<span style="font-family: monospace;"> range </span>added to routine<span style="font-family: monospace;"> boundary_conds </span>(where the
-non-cyclic conditions are set). <br>
-
-      <br>
-
+solving the Poisson equation. Argument<span style="font-family: monospace;"> range </span>added
+to routine<span style="font-family: monospace;">
+boundary_conds </span>(where the
+non-cyclic conditions are set). <br> <br>
 Minor changes:<br>
-
-Output of parameter<span style="font-family: monospace;"> datform </span>on
-&amp;GLOBAL-Namelist (iso2d). This is set to <span style="font-style: italic;">'big_endian'</span> on IBM- and
-NEC-machines.<br>
-
-      <br>
-
-Number of particle groups are additionally calculated in<span style="font-family: monospace;"> header </span>(formerly unknown).<br>
-
-      </td>
-
-      <td style="vertical-align: top;">advec_u_pw,
+Output of parameter<span style="font-family: monospace;">
+datform </span>on
+&amp;GLOBAL-Namelist (iso2d). This is set to <span style="font-style: italic;">'big_endian'</span> on
+IBM- and
+NEC-machines.<br> <br>
+Number of particle groups are additionally calculated in<span style="font-family: monospace;"> header </span>(formerly
+unknown).<br> </td> <td style="vertical-align: top;">advec_u_pw,
 advec_u_up, advec_v_pw, advec_v_up, boundary_conds,
 calc_liquid_water_content, check_parameters, close_file, coriolis,
@@ -3595,604 +1736,239 @@
 init_rankine, modules, parin, poismg, pres, prognostic_equations,
 read_var_list, sor, time_integration, timestep, write_var_list<br>
-
-      </td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;">C<br>
-
-      </td>
-
-      <td style="vertical-align: top;"><span style="font-family: monospace;">pres </span>is called by default only
+</td> </tr> <tr> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">C<br> </td> <td style="vertical-align: top;"><span style="font-family: monospace;">pres </span>is called
+by default only
 at the last Runge-Kutta-substep, which spares a big amount of CPU-time
 (see also new parameter<span style="font-family: monospace;">
-call_psolver_at_all_substeps</span>). No pressure term in the momentum
-equations in case of Runge-Kutta-schemes.<br>
-
-      <span style="font-weight: bold;">Attention:</span> <br>
-
+call_psolver_at_all_substeps</span>). No pressure term in the
+momentum
+equations in case of Runge-Kutta-schemes.<br> <span style="font-weight: bold;">Attention:</span> <br>
 The content of the testresult-file (run control output) has changed
-due to this modification.<br>
-
-      <br>
-
-Program speedup by removing the exchange of ghost points in routines <span style="font-family: monospace;">calc_liquid_water_content</span>, <span style="font-family: monospace;">diffusivities</span>, .<br>
-
-      <br>
-
-Default value for residual limit increased from 1E-6 to 1E-4. Routine<span style="font-family: monospace;"> prolong </span>in<span style="font-family: monospace;"> poismg </span>simplified (one call of<span style="font-family: monospace;"> exchange_horiz </span>spared).<br>
-
-      </td>
-
-      <td style="vertical-align: top;">calc_liquid_water_content,
+due to this modification.<br> <br>
+Program speedup by removing the exchange of ghost points in routines <span style="font-family: monospace;">calc_liquid_water_content</span>,
+<span style="font-family: monospace;">diffusivities</span>,
+.<br> <br>
+Default value for residual limit increased from 1E-6 to 1E-4. Routine<span style="font-family: monospace;"> prolong </span>in<span style="font-family: monospace;"> poismg </span>simplified
+(one call of<span style="font-family: monospace;">
+exchange_horiz </span>spared).<br> </td> <td style="vertical-align: top;">calc_liquid_water_content,
 diffusivities, modules, poismg, pres, time_integration,
-timestep_scheme_steering<br>
-
-      </td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;">E<br>
-
-      </td>
-
-      <td style="vertical-align: top;">Calculation of vertical particle
-velocity (with inertia) corrected, exp_arg had a wrong sign.<br>
-
-      <br>
-
+timestep_scheme_steering<br> </td> </tr> <tr>
+<td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">E<br> </td> <td style="vertical-align: top;">Calculation of vertical
+particle
+velocity (with inertia) corrected, exp_arg had a wrong sign.<br> <br>
 (e)**1.5 replaced by e*SQRT(e) in init_1d_model because of wrong
 results on NEC machines (as a side effect, the new calculation is much
-faster).<br>
-
-      <br>
-
+faster).<br> <br>
 Two errors concerning switch to PE0 (multigrid method) removed from<span style="font-family: monospace;"> init_pegrid</span>.<br>
-
-      <br>
-
-Two errors in<span style="font-family: monospace;"> plot_2d </span>removed:
+<br>
+Two errors in<span style="font-family: monospace;">
+plot_2d </span>removed:
 local_2d is allocated with upper bound<span style="font-family: monospace;"> nzt</span> (former<span style="font-family: monospace;"> nzt+1 </span>was
 wrong), additional barrier in case that PE0 gathers data from the other
 PEs, this barrier had to be set at the end of the gathering because
-otherwise PE0 sometimes received wrong data on tag 0.<br>
-
-      <br>
-
-Use of module indices removed from routine<span style="font-family: monospace;"> split_1dd </span>in<span style="font-family: monospace;"> poisfft </span>since it caused
-errors in case of<span style="font-family: monospace;"> nx /= ny </span>and
-a 1D-decomposition along y.<br>
-
-      <br>
-
+otherwise PE0 sometimes received wrong data on tag 0.<br> <br>
+Use of module indices removed from routine<span style="font-family: monospace;"> split_1dd </span>in<span style="font-family: monospace;"> poisfft </span>since
+it caused
+errors in case of<span style="font-family: monospace;"> nx
+/= ny </span>and
+a 1D-decomposition along y.<br> <br>
 Wrong re-definition of dx2 and dy2 in<span style="font-family: monospace;"> sor </span>removed.<br>
-
-      </td>
-
-      <td style="vertical-align: top;">advec_particles, init_1d_model,
-init_pegrid, plot_2d, poisfft, sor<br>
-
-      </td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;">23/04/05<br>
-
-      </td>
-
-      <td style="vertical-align: top;">SR<br>
-
-      </td>
-
-      <td style="vertical-align: top;">2.9a<br>
-
-      </td>
-
-      <td style="vertical-align: top;">C<br>
-
-      </td>
-
-      <td style="vertical-align: top;">Additional check for
+</td> <td style="vertical-align: top;">advec_particles,
+init_1d_model,
+init_pegrid, plot_2d, poisfft, sor<br> </td> </tr> <tr>
+<td style="vertical-align: top;">23/04/05<br> </td>
+<td style="vertical-align: top;">SR<br> </td>
+<td style="vertical-align: top;">2.9a<br> </td>
+<td style="vertical-align: top;">C<br> </td> <td style="vertical-align: top;">Additional check for
 incompatibilities between non-cyclic lateral boundary conditions and
-other schemes.<br>
-
-      <br>
-
-fcl_factor renamed cfl_factor.<br>
-
-      <br>
-
+other schemes.<br> <br>
+fcl_factor renamed cfl_factor.<br> <br>
 New local array sums_ll declared in flow_statistics instead of
-temporarily using sums_l.<br>
-
-      <br>
-
+temporarily using sums_l.<br> <br>
 crmax (maximum number of crosses allowed in a plot of vertical
-profiles) increased from 20 to 100.<br>
-
-      </td>
-
-      <td style="vertical-align: top;">calc_spectra, check_parameters,
-flow_statistics, header, modules, parin, timestep<br>
-
-      </td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;">E<br>
-
-      </td>
-
-      <td style="vertical-align: top;">Revised calculation of output
+profiles) increased from 20 to 100.<br> </td> <td style="vertical-align: top;">calc_spectra, check_parameters,
+flow_statistics, header, modules, parin, timestep<br> </td>
+</tr> <tr> <td style="vertical-align: top;"><br>
+</td> <td style="vertical-align: top;"><br> </td>
+<td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">E<br> </td> <td style="vertical-align: top;">Revised calculation of output
 time counters regarding a possible decrease of the output time interval
-in case of restart runs.<br>
-
-      <br>
-
+in case of restart runs.<br> <br>
 Error removed in Dirichlet bottom boundary conditions for pt and q in
-case of Runge-Kutta schemes.<br>
-
-      <br>
-
-      <span style="font-weight: bold;">batch_scp</span> calls (needed
+case of Runge-Kutta schemes.<br> <br> <span style="font-weight: bold;">batch_scp</span> calls
+(needed
 for determining correct filenames for AVS-fld-files) is given the
-remote username as an additional argument.<br>
-
-      <br>
-
-Default setting of outflow_damping_width corrected.<br>
-
-      <br>
-
+remote username as an additional argument.<br> <br>
+Default setting of outflow_damping_width corrected.<br> <br>
 Initial horizontal velocities at the lowest grid levels in the 1d-model
 are set to a very small value in order to avoid wrong results and the
-resulting too small timesteps.<br>
-
-      <br>
-
+resulting too small timesteps.<br> <br>
 Implicit counters i renamed i9 in modules due to declaration conflicts
-with i in other subroutines (reported by Portland compiler).<br>
-
-      </td>
-
-      <td style="vertical-align: top;">advec_particles, boundary_conds,
+with i in other subroutines (reported by Portland compiler).<br> </td>
+<td style="vertical-align: top;">advec_particles,
+boundary_conds,
 check_open, check_parameters, init_1d_model, modules, time_integration<br>
-
-      </td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;">18/05/05<br>
-
-      </td>
-
-      <td style="vertical-align: top;">SR<br>
-
-      </td>
-
-      <td style="vertical-align: top;">2.10<br>
-
-      </td>
-
-      <td style="vertical-align: top;">N<br>
-
-      </td>
-
-      <td style="vertical-align: top;">NetCDF support implemented.<br>
-
+</td> </tr> <tr> <td style="vertical-align: top;">18/05/05<br> </td>
+<td style="vertical-align: top;">SR<br> </td>
+<td style="vertical-align: top;">2.10<br> </td>
+<td style="vertical-align: top;">N<br> </td> <td style="vertical-align: top;">NetCDF support implemented.<br>
 Profiles, time series, spectra, cross-sections, 3d-volume-data and
 particle data can now be output in NetCDF format (file handles
 101-108). This output can be switched on with the new d3par-parameter <span style="font-family: monospace;">data_output_format</span>.
 So far, this support is available on IBM systems at HLRN and on the NEC
-system at DKRZ. Output in format suitable for graphic software <span style="font-weight: bold;">profil</span>, <span style="font-weight: bold;">iso2d</span> and <span style="font-weight: bold;">avs</span> is still possible.<br>
-
-      <br>
-
+system at DKRZ. Output in format suitable for graphic software <span style="font-weight: bold;">profil</span>, <span style="font-weight: bold;">iso2d</span> and <span style="font-weight: bold;">avs</span> is still
+possible.<br> <br>
 Output of cross sections for qv, vpt and lwp (liquid water path).<br>
-
-      </td>
-
-      <td style="vertical-align: top;">advec_particles, check_open,
+</td> <td style="vertical-align: top;">advec_particles,
+check_open,
 check_parameters, close_file, header, modules, parin, plot_1d, plot_2d,
 plot_3d, plot_spectra, plot_ts, read_var_list, write_var_list<br>
-
-      <br>
-
-      <span style="font-weight: bold;">new:</span><br>
-
-netcdf<br>
-
-      </td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;">C<br>
-
-      </td>
-
-      <td style="vertical-align: top;">In <span style="font-family: monospace;">calc_spectra</span>, n is assigned the
+<br> <span style="font-weight: bold;">new:</span><br>
+netcdf<br> </td> </tr> <tr> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">C<br> </td> <td style="vertical-align: top;">In <span style="font-family: monospace;">calc_spectra</span>, n
+is assigned the
 number of spectra (formerly, the number of spectra was given by n-1)<br>
-
-      <br>
-
+<br>
 Abort if <span style="font-family: monospace;">poisfft_hybrid</span>
-is called in a non-parallel environment.<br>
-
-      <br>
-
-Default values of spectra package parameters <span style="font-family: monospace;">pl_spectra</span> and <span style="font-family: monospace;">spectra_direction</span> changed.<br>
-
-      </td>
-
-      <td style="vertical-align: top;">calc_spectra, check_parameters,
-modules, pres, read_3d_binary, write_3d_binary<br>
-
-      </td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;">E<br>
-
-      </td>
-
-      <td style="vertical-align: top;">Missing argument in ffty
-(non-parallel case) added.<br>
-
-      <br>
-
-Error in output of particle inertia information removed.<br>
-
-      </td>
-
-      <td style="vertical-align: top;">poisfft, header<br>
-
-      </td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;">05/07/05<br>
-
-      </td>
-
-      <td style="vertical-align: top;">SR<br>
-
-      </td>
-
-      <td style="vertical-align: top;">2.11<br>
-
-      </td>
-
-      <td style="vertical-align: top;">N/C<br>
-
-      </td>
-
-      <td style="vertical-align: top;">New cloud physics code
+is called in a non-parallel environment.<br> <br>
+Default values of spectra package parameters <span style="font-family: monospace;">pl_spectra</span> and <span style="font-family: monospace;">spectra_direction</span>
+changed.<br> </td> <td style="vertical-align: top;">calc_spectra,
+check_parameters,
+modules, pres, read_3d_binary, write_3d_binary<br> </td> </tr>
+<tr> <td style="vertical-align: top;"><br> </td>
+<td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">E<br> </td> <td style="vertical-align: top;">Missing argument in ffty
+(non-parallel case) added.<br> <br>
+Error in output of particle inertia information removed.<br> </td>
+<td style="vertical-align: top;">poisfft, header<br>
+</td> </tr> <tr> <td style="vertical-align: top;">05/07/05<br> </td>
+<td style="vertical-align: top;">SR<br> </td>
+<td style="vertical-align: top;">2.11<br> </td>
+<td style="vertical-align: top;">N/C<br> </td>
+<td style="vertical-align: top;">New cloud physics code
 implemented, which explicitly simulates single cloud droplets,
 including droplet growth by condensation and collision. Using this code
-requires to set the new inipar-parameter <span style="font-family: monospace;">cloud_droplets</span> = .TRUE. and to use the particle package (by using the <span style="font-weight: bold;">mrun</span>-option <span style="font-family: monospace;">-p particles</span>). Coupling between droplets and the thermodynamic quantities (potential
-temperature and specific humidity) is regarded in the new subroutine <span style="font-family: monospace;">interaction_droplets_ptq</span>. The real number of droplets in a grid cell can steered by the new inipar-parameter <span style="font-family: monospace;">initial_weighting_factor</span>.<br>
-
-
-      <br>
-
-
-      <span style="font-weight: bold;">This code will be further extended and probably changed in the near future, so the current status is that of a test version!</span><br>
-
-
-      <br>
-
-
-In connection with this new code implementation, several changes has been made to the existing code:<br>
-
-
-      <br>
-
-
-The particle (droplet) size is now given as radius (previously diameter). E.g., the name of the respective package parameter <span style="font-family: monospace;">diameter</span> has changed to <span style="font-family: monospace;">radius</span>.<br>
-
-
-      <br>
-
-
-The current weighting_factor and radius are stored as additional particle attributes in the type structure <span style="font-family: monospace;">particle_type</span>. The MPI data type <span style="font-family: monospace;">mpi_particle_type</span>, needed for exchanging particles between subdomains, has been changed correspondingly.<br>
-
-
-      <br>
-
-
-      <span style="font-family: monospace;">gas_constant</span> is renamed <span style="font-family: monospace;">r_d</span>, <span style="font-family: monospace;">latent_heat</span> is renamed <span style="font-family: monospace;">r_v</span>, the allowed string length of <span style="font-family: monospace;">pl2d</span> and <span style="font-family: monospace;">pl3d</span> is increased from 6 to 10.<br>
-
-
-      <br>
-
-
-New 3d arrays <span style="font-family: monospace;">ql_1</span>, <span style="font-family: monospace;">ql_2</span>, <span style="font-family: monospace;">ql_v</span> and <span style="font-family: monospace;">ql_vp</span>, <span style="font-family: monospace;">ql</span> is now a pointer, new pointer <span style="font-family: monospace;">ql_c</span>.<br>
-
-
-      <br>
-
-
-Data format for unit 85 (particle data output) is changed.</td>
-
-      <td style="vertical-align: top;">advec_particles.f90
+requires to set the new inipar-parameter <span style="font-family: monospace;">cloud_droplets</span>
+= .TRUE. and to use the particle package (by using the <span style="font-weight: bold;">mrun</span>-option <span style="font-family: monospace;">-p particles</span>).
+Coupling between droplets and the thermodynamic quantities (potential
+temperature and specific humidity) is regarded in the new subroutine <span style="font-family: monospace;">interaction_droplets_ptq</span>.
+The real number of droplets in a grid cell can steered by the new
+inipar-parameter <span style="font-family: monospace;">initial_weighting_factor</span>.<br>
+<br> <span style="font-weight: bold;">This code
+will be further extended and probably changed in the near future, so
+the current status is that of a test version!</span><br> <br>
+In connection with this new code implementation, several changes has
+been made to the existing code:<br> <br>
+The particle (droplet) size is now given as radius (previously
+diameter). E.g., the name of the respective package parameter <span style="font-family: monospace;">diameter</span> has
+changed to <span style="font-family: monospace;">radius</span>.<br>
+<br>
+The current weighting_factor and radius are stored as additional
+particle attributes in the type structure <span style="font-family: monospace;">particle_type</span>.
+The MPI data type <span style="font-family: monospace;">mpi_particle_type</span>,
+needed for exchanging particles between subdomains, has been changed
+correspondingly.<br> <br> <span style="font-family: monospace;">gas_constant</span> is
+renamed <span style="font-family: monospace;">r_d</span>,
+<span style="font-family: monospace;">latent_heat</span>
+is renamed <span style="font-family: monospace;">r_v</span>,
+the allowed string length of <span style="font-family: monospace;">pl2d</span>
+and <span style="font-family: monospace;">pl3d</span>
+is increased from 6 to 10.<br> <br>
+New 3d arrays <span style="font-family: monospace;">ql_1</span>,
+<span style="font-family: monospace;">ql_2</span>, <span style="font-family: monospace;">ql_v</span> and <span style="font-family: monospace;">ql_vp</span>, <span style="font-family: monospace;">ql</span> is now a
+pointer, new pointer <span style="font-family: monospace;">ql_c</span>.<br>
+<br>
+Data format for unit 85 (particle data output) is changed.</td> <td style="vertical-align: top;">advec_particles.f90
 check_open.f90, header.f90, init_3d_model.f90, init_cloud_physics.f90,
 init_particles.f90, modules.f90, package_parin.f90, parin.f90,
 plot_2d.f90, plot_3d.f90, read_var_list.f90, time_integration.f90,
-write_var_list.f90<br>
-
-
-      <br>
-
-
-      <span style="font-weight: bold;">new:</span><br>
-
-
-interaction_droplets_ptq.f90</td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;">C<br>
-
-      </td>
-
-      <td style="vertical-align: top;">Default value of <span style="font-family: monospace;">call_psolver_at_all_substeps</span> is <span style="font-style: italic;">.TRUE.</span>, because small scale waves occured with <span style="font-style: italic;">.FALSE.</span></td>
-
-      <td style="vertical-align: top;">modules.f90</td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;">E<br>
-
-      </td>
-
-      <td style="vertical-align: top;">2d- and 3d-NetCDF-files are not opened if there is no output for them (<span style="font-family: monospace;">pl2d</span> = <span style="font-family: monospace;">pl3d</span> = <span style="font-style: italic;">' '</span>).<br>
-
-
-      <br>
-
-
+write_var_list.f90<br> <br> <span style="font-weight: bold;">new:</span><br>
+interaction_droplets_ptq.f90</td> </tr> <tr> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">C<br> </td> <td style="vertical-align: top;">Default value of <span style="font-family: monospace;">call_psolver_at_all_substeps</span>
+is <span style="font-style: italic;">.TRUE.</span>,
+because small scale waves occured with <span style="font-style: italic;">.FALSE.</span></td> <td style="vertical-align: top;">modules.f90</td> </tr>
+<tr> <td style="vertical-align: top;"><br> </td>
+<td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">E<br> </td> <td style="vertical-align: top;">2d- and 3d-NetCDF-files are not
+opened if there is no output for them (<span style="font-family: monospace;">pl2d</span> = <span style="font-family: monospace;">pl3d</span> = <span style="font-style: italic;">' '</span>).<br> <br>
 No opening of units 101 - 108 in case of missing NetCDF support.</td>
-
-      <td style="vertical-align: top;">check_open.f90, netcdf.f90<br>
-
-      </td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;">30/06/05<br>
-
-      </td>
-
-      <td style="vertical-align: top;">GS<br>
-
-      </td>
-
-      <td style="vertical-align: top;">2.11a<br>
-
-      </td>
-
-      <td style="vertical-align: top;">C<br>
-
-      </td>
-
-      <td style="vertical-align: top;">Scalars ug and vg have been
+<td style="vertical-align: top;">check_open.f90, netcdf.f90<br>
+</td> </tr> <tr> <td style="vertical-align: top;">30/06/05<br> </td>
+<td style="vertical-align: top;">GS<br> </td>
+<td style="vertical-align: top;">2.11a<br> </td>
+<td style="vertical-align: top;">C<br> </td> <td style="vertical-align: top;">Scalars ug and vg have been
 changed into
 arrays in order to allow the specification of a geostrophic wind that
 depends on height (baroclinicity). The initial profiles of the the
 u- and v-component of the geostrophic wind are initialized by
-specifying the new initialization parameters <span style="font-family: monospace;">ug_surface</span>,
-      <span style="font-family: monospace;">ug_vertical_gradient</span>, <span style="font-family: monospace;">ug_vertical_gradient_level</span> and <span style="font-family: monospace;">vg_surface</span>,
-      <span style="font-family: monospace;">vg_vertical_gradient</span>, <span style="font-family: monospace;">vg_vertical_gradient_level</span>, respectively. <br>
-
-      <br>
-
-The former initialization parameters <span style="font-family: monospace;">ug</span> and <span style="font-family: monospace;">vg</span> have been removed. <br>
-
-      </td>
-
-      <td style="vertical-align: top;">boundary_conds,
+specifying the new initialization parameters <span style="font-family: monospace;">ug_surface</span>, <span style="font-family: monospace;">ug_vertical_gradient</span>,
+<span style="font-family: monospace;">ug_vertical_gradient_level</span>
+and <span style="font-family: monospace;">vg_surface</span>,
+<span style="font-family: monospace;">vg_vertical_gradient</span>,
+<span style="font-family: monospace;">vg_vertical_gradient_level</span>,
+respectively. <br> <br>
+The former initialization parameters <span style="font-family: monospace;">ug</span> and <span style="font-family: monospace;">vg</span> have been
+removed. <br> </td> <td style="vertical-align: top;">boundary_conds,
 check_parameters, coriolis, header, init_1d_model, modules, parin,
 prognostic_equations, read_var_list, spline_z, write_var_list<br>
-
-      </td>
-
-    </tr>
-
-    <tr>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;"><br>
-
-      </td>
-
-      <td style="vertical-align: top;">E<br>
-
-      </td>
-
-      <td style="vertical-align: top;">Error in the calculation of u_0
+</td> </tr> <tr> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">E<br> </td> <td style="vertical-align: top;">Error in the calculation of u_0
 and
-v_0 in production_e.f90 removed.<br>
-
-      <br>
-
-      <span style="font-weight: bold;">Attention:</span> <br>
-
+v_0 in production_e.f90 removed.<br> <br> <span style="font-weight: bold;">Attention:</span> <br>
 The content of the testresult-file (run control output) has changed
-due to this modification.<br>
-
-      </td>
-
-      <td style="vertical-align: top;">production_e<br>
-
-      </td>
-
-    </tr>
-    <tr>
-      <td style="text-align: left; vertical-align: top;">20/10/05</td>
-      <td style="vertical-align: top;">SR</td>
-      <td style="vertical-align: top;">2.11b</td>
-      <td style="vertical-align: top;">N/C</td>
-      <td style="vertical-align: top;">Collision process for droplet growth completed.<br>
-      <br>
-      <span style="font-weight: bold;">The droplet code is still under development and requires further speed optimization!<br>
-      <br>
-      </span>Number of particles really used is additionally output on the netcdf particle data file.<span style="font-weight: bold;"></span></td>
-      <td style="vertical-align: top;">advec_particles, modules, netcdf</td>
-    </tr>
-    <tr>
-      <td style="vertical-align: top;">&nbsp;</td>
-      <td style="vertical-align: top;">&nbsp;</td>
-      <td style="vertical-align: top;">&nbsp;</td>
-      <td style="vertical-align: top;">E</td>
-      <td style="vertical-align: top;">The
+due to this modification.<br> </td> <td style="vertical-align: top;">production_e<br> </td>
+</tr> <tr> <td style="text-align: left; vertical-align: top;">20/10/05</td>
+<td style="vertical-align: top;">SR</td> <td style="vertical-align: top;">2.11b</td> <td style="vertical-align: top;">N/C</td> <td style="vertical-align: top;">Collision process for droplet
+growth completed.<br> <br> <span style="font-weight: bold;">The droplet code is still under
+development and requires further speed optimization!<br> <br>
+</span>Number of particles really used is additionally output on
+the netcdf particle data file.<span style="font-weight: bold;"></span></td>
+<td style="vertical-align: top;">advec_particles, modules,
+netcdf</td> </tr> <tr> <td style="vertical-align: top;">&nbsp;</td> <td style="vertical-align: top;">&nbsp;</td> <td style="vertical-align: top;">&nbsp;</td> <td style="vertical-align: top;">E</td> <td style="vertical-align: top;">The
 last PE in a row is not allowed to have more grid points than the other
 PES (only less). Jobs crashed in these cases that the last PE has more
 grid points. Therefore, the number of gridpoints along x (and y) must
-now meet a special condition which is checked in case that <span style="font-family: monospace;">grid_matching</span> = <span style="font-style: italic;">'match'</span> is used.<br>
-      <br>
-Error removed in calculating y axis data for yz cross sections. Error in output of netcdf yz-cross-sections removed.<br>
-      <br>
-2*r replaced by r in the exponential term of the particle momentum equation.<br>
-Error in output of yz-slice information removed.</td>
-      <td style="vertical-align: top;">advec_particles, header, init_pegrid, netcdf, plot_2d</td>
-    </tr>
-    <tr>
-      <td style="vertical-align: top;">06/12/05</td>
-      <td style="vertical-align: top;">SR</td>
-      <td style="vertical-align: top;">2.11c</td>
-      <td style="vertical-align: top;">N</td>
-      <td style="vertical-align: top;">Output of ql profile is allowed in case of using cloud droplets.</td>
-      <td style="vertical-align: top;">check_parameters, flow_statistics, modules</td>
-    </tr>
-    <tr>
-      <td style="vertical-align: top;">&nbsp;</td>
-      <td style="vertical-align: top;">&nbsp;</td>
-      <td style="vertical-align: top;">&nbsp;</td>
-      <td style="vertical-align: top;">E</td>
-      <td style="vertical-align: top;">The horizontally averaged pt profile is calculated within<span style="font-family: monospace;"> buoyancy </span>only
+now meet a special condition which is checked in case that <span style="font-family: monospace;">grid_matching</span> =
+<span style="font-style: italic;">'match'</span> is
+used.<br> <br>
+Error removed in calculating y axis data for yz cross sections. Error
+in output of netcdf yz-cross-sections removed.<br> <br>
+2*r replaced by r in the exponential term of the particle momentum
+equation.<br>
+Error in output of yz-slice information removed.</td> <td style="vertical-align: top;">advec_particles, header,
+init_pegrid, netcdf, plot_2d</td> </tr> <tr> <td style="vertical-align: top;">06/12/05</td> <td style="vertical-align: top;">SR</td> <td style="vertical-align: top;">2.11c</td> <td style="vertical-align: top;">N</td> <td style="vertical-align: top;">Output of ql profile is allowed
+in case of using cloud droplets.</td> <td style="vertical-align: top;">check_parameters,
+flow_statistics, modules</td> </tr> <tr> <td style="vertical-align: top;">&nbsp;</td> <td style="vertical-align: top;">&nbsp;</td> <td style="vertical-align: top;">&nbsp;</td> <td style="vertical-align: top;">E</td> <td style="vertical-align: top;">The horizontally averaged pt
+profile is calculated within<span style="font-family: monospace;">
+buoyancy </span>only
 in case of the first respective intermediate timestep. This is done in
 order to spare communication time and to produce identical model
-results with jobs which are calling<span style="font-family: monospace;"> flow_statistics </span>at different time intervals.<br>
-      <br>
-Error in netcdf variable declaration for<span style="font-family: monospace;"> statistic_regions </span>/= <span style="font-style: italic;">0</span> removed.</td>
-      <td style="vertical-align: top;">buoyancy, netcdf</td>
-    </tr><tr><td style="vertical-align: top;">23/02/06</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.0</td><td style="vertical-align: top;">N</td><td><span style="font-weight: bold;">Version for simulating flows over topography (buildings, mountains, etc.)</span><br style="font-weight: bold;"><br>Besides
+results with jobs which are calling<span style="font-family: monospace;"> flow_statistics </span>at
+different time intervals.<br> <br>
+Error in netcdf variable declaration for<span style="font-family: monospace;"> statistic_regions </span>/=
+<span style="font-style: italic;">0</span> removed.</td>
+<td style="vertical-align: top;">buoyancy, netcdf</td>
+</tr><tr><td style="vertical-align: top;">23/02/06</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.0</td><td style="vertical-align: top;">N</td><td><span style="font-weight: bold;">Version for simulating flows over
+topography (buildings, mountains, etc.)</span><br style="font-weight: bold;"><br>Besides
 the implementation of topography, this version contains several other
-improvements. The most important ones are listed here:<br><ul><li>The cloud droplet code has been optimized (it now runs about 100 times faster than before).</li><li>Output of time-averaged 2d/3d data (sections / volume data) is possible.</li><li>Output of 2d-section data averaged along the direction normal to this section can be done.</li><li>The user-interface has been extended in order to allow 2d-/3d-output of user-defined quantities.</li></ul>The kind of topography to be used is controlled by the new initialization parameter <span style="font-family: monospace;">topography</span>. Allowed values are <span style="font-style: italic;">'flat'</span>, <span style="font-style: italic;">'single_building'</span>, and <span style="font-style: italic;">'read_from_file'</span>. The user can define his own special topography setting within the new user-interface routine <span style="font-family: monospace;">user_init_grid</span>. For the single building case, the&nbsp;building (rectangular) size and position can be controlled with parameters <span style="font-family: monospace;">building_height</span>, <span style="font-family: monospace;">building_length_x</span>, <span style="font-family: monospace;">building_length_y</span>, <span style="font-family: monospace;">building_wall_left</span>, and <span style="font-family: monospace;">building_wall_south</span>.<br>
-      <br>
-The topography code is already optimized (vectorized) for use on the NEC-SX6.<br><br>Output of time-averaged 2d-/3d-data can be enabled by appending the string <span style="font-style: italic;">'_av'</span> to the respective output quantities given with the new runtime-parameter <span style="font-family: monospace;">data_output</span>
+improvements. The most important ones are listed here:<br><ul><li>The
+cloud droplet code has been optimized (it now runs about 100 times
+faster than before).</li><li>Output of time-averaged 2d/3d
+data (sections / volume data) is possible.</li><li>Output
+of 2d-section data averaged along the direction normal to this section
+can be done.</li><li>The user-interface has been extended
+in order to allow 2d-/3d-output of user-defined quantities.</li></ul>The
+kind of topography to be used is controlled by the new initialization
+parameter <span style="font-family: monospace;">topography</span>.
+Allowed values are <span style="font-style: italic;">'flat'</span>,
+<span style="font-style: italic;">'single_building'</span>,
+and <span style="font-style: italic;">'read_from_file'</span>.
+The user can define his own special topography setting within the new
+user-interface routine <span style="font-family: monospace;">user_init_grid</span>.
+For the single building case, the&nbsp;building (rectangular) size
+and position can be controlled with parameters <span style="font-family: monospace;">building_height</span>,
+<span style="font-family: monospace;">building_length_x</span>,
+<span style="font-family: monospace;">building_length_y</span>,
+<span style="font-family: monospace;">building_wall_left</span>,
+and <span style="font-family: monospace;">building_wall_south</span>.<br>
+<br>
+The topography code is already optimized (vectorized) for use on the
+NEC-SX6.<br><br>Output of time-averaged 2d-/3d-data can be
+enabled by appending the string <span style="font-style: italic;">'_av'</span>
+to the respective output quantities given with the new
+runtime-parameter <span style="font-family: monospace;">data_output</span>
 (see further below for the (name)changes in the parameters for data
 output). Example: if a time-average of the vertical cross-section of
@@ -4202,28 +1978,58 @@
 the averaging interval and the temporal distance of the single
 instantaneous sections (volumes) entering into this averaging is
-controlled by the new runtime-parameters <span style="font-family: monospace;">averaging_interval</span>and <span style="font-family: monospace;">dt_averaging_input</span>. <br><br>An average normal to the direction of a cross-section can be output by setting the respective cross section index to <span style="font-style: italic;">-1</span> (example: <span style="font-family: monospace;">section_xz</span> = <span style="font-style: italic;">-1</span>).<br><br>New user-interface routines allowing the output of user-defined quantities are <span style="font-family: monospace;">user_check_data_output</span>, <span style="font-family: monospace;">user_define_netcdf_grid</span>, <span style="font-family: monospace;">user_data_output_2d</span>, <span style="font-family: monospace;">user_data_output_3d</span>, and <span style="font-family: monospace;">user_3d_data_averaging</span>. These quantities can be calculated using the new <span style="font-family: monospace;">CASE</span>-Entry <span style="font-style: italic;">'after integration'</span> in the user-interface routine <span style="font-family: monospace;">user_actions</span>.<br><br>Optimization
+controlled by the new runtime-parameters <span style="font-family: monospace;">averaging_interval</span>and
+<span style="font-family: monospace;">dt_averaging_input</span>.
+<br><br>An average normal to the direction of a
+cross-section can be output by setting the respective cross section
+index to <span style="font-style: italic;">-1</span>
+(example: <span style="font-family: monospace;">section_xz</span>
+= <span style="font-style: italic;">-1</span>).<br><br>New
+user-interface routines allowing the output of user-defined quantities
+are <span style="font-family: monospace;">user_check_data_output</span>,
+<span style="font-family: monospace;">user_define_netcdf_grid</span>,
+<span style="font-family: monospace;">user_data_output_2d</span>,
+<span style="font-family: monospace;">user_data_output_3d</span>,
+and <span style="font-family: monospace;">user_3d_data_averaging</span>.
+These quantities can be calculated using the new <span style="font-family: monospace;">CASE</span>-Entry <span style="font-style: italic;">'after integration'</span>
+in the user-interface routine <span style="font-family: monospace;">user_actions</span>.<br><br>Optimization
 of the cloud droplet code is mainly done by re-sorting all particles
 after every timestep in a way that now all particles which belong to
 the same grid-box (k,j,i) are stored consecutively (wthout any strides)
-in memory.<br>
-      <br>
-Code adaptation for the IBM at the department of atmospheric sciences, Yonsei university (ibmy). Affected subroutines: <span style="font-family: monospace;">fft_xy</span>, <span style="font-family: monospace;">netcdf</span>, <br><br><span style="font-weight: bold;">Further new features:</span><br><ul><li>Use of particle tails now requires setting of the new package-parameter <span style="font-family: monospace;">use_particle_tails</span>. The new parameter skip_particles_for_tail allows only every <span style="font-family: monospace;">skip_particle_for_tail</span>'th parameter having a tail.</li><li>Initialization-parameter <span style="font-family: monospace;">e_min</span> can be used to set a lower limit for the subgrid-scale TKE (i.e. to guarantee a lower limit for the flow's Reynolds-number).</li><li>A
+in memory.<br> <br>
+Code adaptation for the IBM at the department of atmospheric sciences,
+Yonsei university (ibmy). Affected subroutines: <span style="font-family: monospace;">fft_xy</span>, <span style="font-family: monospace;">netcdf</span>, <br><br><span style="font-weight: bold;">Further new features:</span><br><ul><li>Use
+of particle tails now requires setting of the new package-parameter <span style="font-family: monospace;">use_particle_tails</span>.
+The new parameter skip_particles_for_tail allows only every <span style="font-family: monospace;">skip_particle_for_tail</span>'th
+parameter having a tail.</li><li>Initialization-parameter <span style="font-family: monospace;">e_min</span> can be
+used to set a lower limit for the subgrid-scale TKE (i.e. to guarantee
+a lower limit for the flow's Reynolds-number).</li><li>A
 conservation of the volume flow (through the complete xz- and
 yz-sections &nbsp;of the total domain) can be enforced by the new
-initialization-parameter <span style="font-family: monospace;">conserve_volume_flow</span>.</li><li>Data output can be skipped for a given time interval from simulation start (t=0) using new runtime-parameters <span style="font-family: monospace;">skip_time_dopr</span>, <span style="font-family: monospace;">skip_time_dosp</span>, <span style="font-family: monospace;">skip_time_do2d_xy</span>, <span style="font-family: monospace;">skip_time_do2d_xz</span>, <span style="font-family: monospace;">skip_time_do2d_yz</span>, and <span style="font-family: monospace;">skip_time_do3d</span>.<br>
-</li><li>By default, NetCDF output is now using 64-bit offset format
+initialization-parameter <span style="font-family: monospace;">conserve_volume_flow</span>.</li><li>Data
+output can be skipped for a given time interval from simulation start
+(t=0) using new runtime-parameters <span style="font-family: monospace;">skip_time_dopr</span>,
+<span style="font-family: monospace;">skip_time_dosp</span>,
+<span style="font-family: monospace;">skip_time_do2d_xy</span>,
+<span style="font-family: monospace;">skip_time_do2d_xz</span>,
+<span style="font-family: monospace;">skip_time_do2d_yz</span>,
+and <span style="font-family: monospace;">skip_time_do3d</span>.<br>
+</li><li>By default, NetCDF output is now using 64-bit
+offset format
 (large file support) on most machines. The user can switch to the
-32-bit offset format with the new runtime-parameter <span style="font-family: monospace;">netcdf_64bit</span>. Units 111:113&nbsp; and 116 are opened for NetCDF output of time-averaged 2d-3d data.<br>
-        </li>
-        <li>Data
+32-bit offset format with the new runtime-parameter <span style="font-family: monospace;">netcdf_64bit</span>.
+Units 111:113&nbsp; and 116 are opened for NetCDF output of
+time-averaged 2d-3d data.<br> </li> <li>Data
 logging routines are added (see file data_log.f90), which can be used
 for debugging purposes. Output is done on unit 20 (local file
-name/directory name DATA_LOG).</li>
-        <li>New routine <span style="font-family: monospace;">exchange_horiz_2d_int</span> for&nbsp; ghostpoint exchange of 2d-integer arrays. Routine <span style="font-family: monospace;">exchange_horiz_2d</span> is extended for non-cyclic boundary conditions in the non-parallel case.</li>
-        <li>Two different methods for calculating the mixing length and
+name/directory name DATA_LOG).</li> <li>New routine <span style="font-family: monospace;">exchange_horiz_2d_int</span>
+for&nbsp; ghostpoint exchange of 2d-integer arrays. Routine <span style="font-family: monospace;">exchange_horiz_2d</span>
+is extended for non-cyclic boundary conditions in the non-parallel case.</li>
+<li>Two different methods for calculating the mixing length and
 the dissipation can be used in the 1d-model. These are steered by the
-new initializing-parameters <span style="font-family: monospace;">dissipation_1d</span> and <span style="font-family: monospace;">mixing_length_1d</span>.<br>
-        </li>
-</ul></td><td style="vertical-align: top;">advec_particles, advec_s_pw,
+new initializing-parameters <span style="font-family: monospace;">dissipation_1d</span>
+and <span style="font-family: monospace;">mixing_length_1d</span>.<br>
+</li>
+</ul></td><td style="vertical-align: top;">advec_particles,
+advec_s_pw,
 advec_s_up, advec_u_pw, advec_u_up, advec_v_pw, advec_v_up, advec_w_pw,
 advec_w_up, boundary_conds, buoyancy, check_open, check_parameters,
@@ -4235,17 +2041,18 @@
 prandtl_fluxes, pres, production_e, prognostic_equations,
 read_3d_binary, read_var_list, run_control, time_integration,
-user_interface, write_3d_binary, write_var_list<br>
-      <br>
-      <span style="font-weight: bold;">new:</span><br>
-average_3d_data, data_log, sum_up_3d_data<br>
-      <br>
-      <span style="font-weight: bold;">renamed:<br>
-      </span>plot_dvrp -&gt; data_output_dvrp<br>
+user_interface, write_3d_binary, write_var_list<br> <br> <span style="font-weight: bold;">new:</span><br>
+average_3d_data, data_log, sum_up_3d_data<br> <br> <span style="font-weight: bold;">renamed:<br> </span>plot_dvrp
+-&gt; data_output_dvrp<br>
 plot_spectra -&gt; data_output_spectra<br>
 plot_ts -&gt; data_output_tseries<br>
-plot_1d -&gt; data_output_profiles<br>
-      <span style="font-weight: bold;"></span>plot_2d -&gt; data_output_2d<br>
+plot_1d -&gt; data_output_profiles<br> <span style="font-weight: bold;"></span>plot_2d -&gt;
+data_output_2d<br>
 plot_3d -&gt; data_output_3d<br>
-</td></tr><tr><td style="vertical-align: top;">&nbsp;</td><td>&nbsp;</td><td>&nbsp;</td><td style="vertical-align: top;">C</td><td>In parallel mode, all PEs are opening the same copy of the NAMELIST-parameter file <span style="font-family: monospace;">PARIN</span>, i.e. from now on, the file attribute "<span style="font-family: monospace;">npe</span>" must be removed from the corresponding file connection statements in the <span style="font-weight: bold;">mrun</span> configuration file (.mrun.config).<br><br>A
+</td></tr><tr><td style="vertical-align: top;">&nbsp;</td><td>&nbsp;</td><td>&nbsp;</td><td style="vertical-align: top;">C</td><td>In
+parallel mode, all PEs are opening the same copy of the
+NAMELIST-parameter file <span style="font-family: monospace;">PARIN</span>,
+i.e. from now on, the file attribute "<span style="font-family: monospace;">npe</span>" must be
+removed from the corresponding file connection statements in the <span style="font-weight: bold;">mrun</span> configuration
+file (.mrun.config).<br><br>A
 large number of parameter names and local filenames have been changed
 (see&nbsp;list below). Most of them are parameters concerning data
@@ -4253,13 +2060,23 @@
 to "data_output" and "do", respectively. In addition to these changes,
 the names of many internal PALM parameters and variables have also
-changed.<br><br>Restart runs now require setting of the initialization-parameter <span style="font-family: monospace;">initializing_actions</span> = <span style="font-style: italic;">'read_restart_data'</span>.<br>
-      <br>
+changed.<br><br>Restart runs now require setting of the
+initialization-parameter <span style="font-family: monospace;">initializing_actions</span>
+= <span style="font-style: italic;">'read_restart_data'</span>.<br>
+<br>
 Concerning particles, the parameters <span style="font-family: monospace;">psl</span>, <span style="font-family: monospace;">psr</span>, <span style="font-family: monospace;">pdx</span>,
 etc., which are controlling the position of the particle source and the
 number of particles within this source, are now 1d-arrays. This allows
 to define different particle sources for different particle groups.<br>
-      <br>
-The number of particle groups has to be set by the user (no automatic calculation from the value of parameter <span style="font-family: monospace;">density_ratio</span> any more)<span style="font-style: italic;"></span>. Variable <span style="font-family: monospace;">idum</span> in type <span style="font-family: monospace;">particle_type</span> has been renamed <span style="font-family: monospace;">tail_id</span>.<br><br>Output for units 15 and 17 flushed in <span style="font-family: monospace;">init_1d_model</span>.<br><br>List of changed parameter names:<br><table style="text-align: left; width: 559px; height: 831px;" border="1" cellpadding="2" cellspacing="2"><tbody><tr><td style="font-weight: bold;" align="undefined" valign="undefined">Old name</td><td align="undefined" valign="undefined"><span style="font-weight: bold;">New name</span></td><td align="undefined" valign="undefined"><span style="font-weight: bold;">Comments</span></td></tr><tr><td align="undefined" valign="undefined">average_period_pl1d</td><td align="undefined" valign="undefined">averaging_interval_pr</td><td align="undefined" valign="undefined"><br>
-</td></tr><tr><td style="vertical-align: top;">average_period_pr1d</td><td style="vertical-align: top;">---</td><td style="vertical-align: top;">averaging of printed profiles is not possible any more</td></tr><tr><td align="undefined" valign="undefined">average_period_sp</td><td align="undefined" valign="undefined">averaging_interval_sp</td><td align="undefined" valign="undefined"><br>
+<br>
+The number of particle groups has to be set by the user (no automatic
+calculation from the value of parameter <span style="font-family: monospace;">density_ratio</span>
+any more)<span style="font-style: italic;"></span>.
+Variable <span style="font-family: monospace;">idum</span>
+in type <span style="font-family: monospace;">particle_type</span>
+has been renamed <span style="font-family: monospace;">tail_id</span>.<br><br>Output
+for units 15 and 17 flushed in <span style="font-family: monospace;">init_1d_model</span>.<br><br>List
+of changed parameter names:<br><table style="text-align: left; width: 559px; height: 831px;" border="1" cellpadding="2" cellspacing="2"><tbody><tr><td style="font-weight: bold;" align="undefined" valign="undefined">Old name</td><td align="undefined" valign="undefined"><span style="font-weight: bold;">New name</span></td><td align="undefined" valign="undefined"><span style="font-weight: bold;">Comments</span></td></tr><tr><td align="undefined" valign="undefined">average_period_pl1d</td><td align="undefined" valign="undefined">averaging_interval_pr</td><td align="undefined" valign="undefined"><br>
+</td></tr><tr><td style="vertical-align: top;">average_period_pr1d</td><td style="vertical-align: top;">---</td><td style="vertical-align: top;">averaging of printed profiles
+is not possible any more</td></tr><tr><td align="undefined" valign="undefined">average_period_sp</td><td align="undefined" valign="undefined">averaging_interval_sp</td><td align="undefined" valign="undefined"><br>
 </td></tr><tr><td align="undefined" valign="undefined">dt_average</td><td align="undefined" valign="undefined">dt_averaging_input_pr</td><td align="undefined" valign="undefined"><br>
 </td></tr><tr><td align="undefined" valign="undefined">dt_pl1d</td><td align="undefined" valign="undefined">dt_dopr</td><td align="undefined" valign="undefined"><br>
@@ -4275,5 +2092,6 @@
 </td></tr><tr><td align="undefined" valign="undefined">plts</td><td align="undefined" valign="undefined">data_output_ts</td><td align="undefined" valign="undefined"><br>
 </td></tr><tr><td align="undefined" valign="undefined">pl1d</td><td align="undefined" valign="undefined">data_output_pr</td><td align="undefined" valign="undefined"><br>
-</td></tr><tr><td style="vertical-align: top;">pl2d</td><td style="vertical-align: top;">data_output</td><td style="vertical-align: top;">old parameters pl2d and pl3d are joined to new parameter data_output</td></tr><tr><td align="undefined" valign="undefined">pl3d</td><td align="undefined" valign="undefined">data_output</td><td align="undefined" valign="undefined"><br>
+</td></tr><tr><td style="vertical-align: top;">pl2d</td><td style="vertical-align: top;">data_output</td><td style="vertical-align: top;">old parameters pl2d and pl3d
+are joined to new parameter data_output</td></tr><tr><td align="undefined" valign="undefined">pl3d</td><td align="undefined" valign="undefined">data_output</td><td align="undefined" valign="undefined"><br>
 </td></tr><tr><td align="undefined" valign="undefined">pl2d_at_begin</td><td align="undefined" valign="undefined">do2d_at_begin</td><td align="undefined" valign="undefined"><br>
 </td></tr><tr><td align="undefined" valign="undefined">pl3d_at_begin</td><td align="undefined" valign="undefined">do3d_at_begin</td><td align="undefined" valign="undefined"><br>
@@ -4292,17 +2110,33 @@
 run_control, time_integration, timestep, write_compressed,
 write_var_list<br>
-</td></tr><tr><td>&nbsp;</td><td>&nbsp;</td><td>&nbsp;</td><td style="vertical-align: top;">E</td><td>Error in the particle release (defaults of <span style="font-family: monospace;">psl</span>, <span style="font-family: monospace;">psr</span>, <span style="font-family: monospace;">psn</span>, <span style="font-family: monospace;">pss</span>) at the PE boundary has been fixed.<br>
-      <br>
-Error removed getting the variable ids for <span style="font-family: monospace;">_ext</span> for profiles and time series.<br><br>NetCDF calls in <span style="font-family: monospace;">data_output_tseries</span> only if NetCDF output is switched on.<br><br>Output of the initial geostrophic wind profile corrected.<br><br>Maximum possibel timestep for Runge-Kutta-schemes reduced, error concerning initialization of&nbsp; <span style="font-family: monospace;">l_black</span> removed (both in <span style="font-family: monospace;">init_1d_model</span>).<br>
-</td><td style="vertical-align: top;">advec_particles, data_output_tseries, header, init_particles, init_1d_model, netcdf<br>
-</td></tr><tr><td align="left" valign="top">26/04/06</td><td align="left" valign="top">SR</td><td align="left" valign="top">3.0a</td><td align="left" valign="top">N</td><td align="left" valign="top">OpenMP parallelization for the cache-optimized PALM code.<br><br>OpenMP
+</td></tr><tr><td>&nbsp;</td><td>&nbsp;</td><td>&nbsp;</td><td style="vertical-align: top;">E</td><td>Error in
+the particle release (defaults of <span style="font-family: monospace;">psl</span>, <span style="font-family: monospace;">psr</span>, <span style="font-family: monospace;">psn</span>, <span style="font-family: monospace;">pss</span>) at the PE
+boundary has been fixed.<br> <br>
+Error removed getting the variable ids for <span style="font-family: monospace;">_ext</span> for
+profiles and time series.<br><br>NetCDF calls in <span style="font-family: monospace;">data_output_tseries</span>
+only if NetCDF output is switched on.<br><br>Output of the
+initial geostrophic wind profile corrected.<br><br>Maximum
+possibel timestep for Runge-Kutta-schemes reduced, error concerning
+initialization of&nbsp; <span style="font-family: monospace;">l_black</span>
+removed (both in <span style="font-family: monospace;">init_1d_model</span>).<br>
+</td><td style="vertical-align: top;">advec_particles,
+data_output_tseries, header, init_particles, init_1d_model, netcdf<br>
+</td></tr><tr><td align="left" valign="top">26/04/06</td><td align="left" valign="top">SR</td><td align="left" valign="top">3.0a</td><td align="left" valign="top">N</td><td align="left" valign="top">OpenMP parallelization for
+the cache-optimized PALM code.<br><br>OpenMP
 directives have been added where necessary. Still further tuning will
 be necessary to get optimum performance. Other parts of the code than
-the cache-optimized code still need to be (OpenMP-) parallelized. <br><br>So far, threads are activated by setting the <span style="font-weight: bold;">mrun</span>-options <span style="font-family: monospace;">-T</span> <span style="font-style: italic;">&lt;number_of_threads&gt;</span> and <span style="font-family: monospace;">-O</span>, i.e. when the Option <span style="font-family: monospace;">-O</span> is given, the <span style="font-family: monospace;">-T</span> argument is interpreted as the number of threads to be used. When <span style="font-family: monospace;">-O</span> is given, the number of MPI tasks per node is automatically set to <span style="font-style: italic;">1</span>.<br><br>Beside
+the cache-optimized code still need to be (OpenMP-) parallelized. <br><br>So
+far, threads are activated by setting the <span style="font-weight: bold;">mrun</span>-options <span style="font-family: monospace;">-T</span> <span style="font-style: italic;">&lt;number_of_threads&gt;</span>
+and <span style="font-family: monospace;">-O</span>,
+i.e. when the Option <span style="font-family: monospace;">-O</span>
+is given, the <span style="font-family: monospace;">-T</span>
+argument is interpreted as the number of threads to be used. When <span style="font-family: monospace;">-O</span> is given,
+the number of MPI tasks per node is automatically set to <span style="font-style: italic;">1</span>.<br><br>Beside
 the additional compiler directives, the main code change required by
 the OpenMP parallelization is that global sums are now calculated by
 first calculating local sums on each thread. As the next step, local
 sums are calculated from these thread sums for each MPI task before the
-global sum is computed via MPI_ALLREDUCE.&nbsp; For storing the thread
+global sum is computed via MPI_ALLREDUCE.&nbsp; For storing the
+thread
 sums, arrays sums_l and sums_l_l now have three dimensions where the
 number of elements of the third dimension is equal to the number of
@@ -4310,90 +2144,174 @@
 cpu_statistics, diffusivities, flow_statistics, header, init_pegrid,
 init_3d_model, modules, prandtl_fluxes, pres, production_e</td></tr><tr><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">ML</td><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">C</td><td align="left" valign="top">Provisional
-correction for Piacsek &amp; Williams advection scheme: keep u and v
-zero one layer below the topography (in case of set_1d-model_profiles).</td><td align="left" valign="top">init_3d_model</td></tr><tr><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">SR/ML</td><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">E</td><td align="left" valign="top">In <span style="font-family: monospace;">flow_statistics.f90</span> nzb_s_inner(j,i) replaced by nzb in determination of z_i.<br><br>Errors removed in the computation of the diabatic mixing length (<span style="font-family: monospace;">init_1d_model.f90</span>).<br><br>Error removed in extend mode when checking whether the selected cross sections match those in the already existing NetCDF file (<span style="font-family: monospace;">netcdf.f90</span>).</td><td align="left" valign="top">flow_statistics, init_1d_model, netcdf</td></tr><tr><td align="left" valign="top">02/06/06</td><td align="left" valign="top">SR</td><td align="left" valign="top">3.0b</td><td align="left" valign="top">N</td><td align="left" valign="top">This version is adapted to the Sun Fire X4600 system at TIT (to be used by setting the <span style="font-weight: bold;">mrun</span> option <span style="font-family: monospace;">-h lctit</span>).</td><td align="left" valign="top">cpu_log, init_1d_model, local_tremain, loca&ouml;_tremain_ini, run_control</td></tr><tr><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">C</td><td align="left" valign="top">NetCDF&nbsp;2d-
+correction for Piacsek &amp; Williams advection scheme: keep u and
+v
+zero one layer below the topography (in case of set_1d-model_profiles).</td><td align="left" valign="top">init_3d_model</td></tr><tr><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">SR/ML</td><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">E</td><td align="left" valign="top">In <span style="font-family: monospace;">flow_statistics.f90</span>
+nzb_s_inner(j,i) replaced by nzb in determination of z_i.<br><br>Errors
+removed in the computation of the diabatic mixing length (<span style="font-family: monospace;">init_1d_model.f90</span>).<br><br>Error
+removed in extend mode when checking whether the selected cross
+sections match those in the already existing NetCDF file (<span style="font-family: monospace;">netcdf.f90</span>).</td><td align="left" valign="top">flow_statistics,
+init_1d_model, netcdf</td></tr><tr><td align="left" valign="top">02/06/06</td><td align="left" valign="top">SR</td><td align="left" valign="top">3.0b</td><td align="left" valign="top">N</td><td align="left" valign="top">This version is adapted to the Sun Fire X4600
+system at TIT (to be used by setting the <span style="font-weight: bold;">mrun</span> option <span style="font-family: monospace;">-h lctit</span>).</td><td align="left" valign="top">cpu_log, init_1d_model,
+local_tremain, loca&ouml;_tremain_ini, run_control</td></tr><tr><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">C</td><td align="left" valign="top">NetCDF&nbsp;2d-
 and 3d-datasets now contain the exact coordinates of the variables
 along x and y with respect to the positions where they are defined on
-the staggered grid.<br><br>NetCDF datasets of cross sections (e.g. <span style="font-family: monospace;">DATA_2D_XY_NETCDF</span>) are only opened if the respective cross sections are really requested by the settings of parameter <span style="font-family: monospace;">data_output</span>.</td><td align="left" valign="top">check_parameters, data_output_2d, data_output_3d, modules, netcdf, user_interface</td></tr><tr><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">E</td><td align="left" valign="top">p is assigned to <span style="font-family: monospace;">to_be_resorted</span> instead of w</td><td align="left" valign="top">data_output_2d</td></tr><tr><td align="left" valign="top">04/08/06</td><td align="left" valign="top">SR</td><td align="left" valign="top">3.1</td><td align="left" valign="top">N</td><td align="left" valign="top">Subgrid
+the staggered grid.<br><br>NetCDF datasets of cross
+sections (e.g. <span style="font-family: monospace;">DATA_2D_XY_NETCDF</span>)
+are only opened if the respective cross sections are really requested
+by the settings of parameter <span style="font-family: monospace;">data_output</span>.</td><td align="left" valign="top">check_parameters,
+data_output_2d, data_output_3d, modules, netcdf, user_interface</td></tr><tr><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">E</td><td align="left" valign="top">p is assigned to <span style="font-family: monospace;">to_be_resorted</span>
+instead of w</td><td align="left" valign="top">data_output_2d</td></tr><tr><td align="left" valign="top">04/08/06</td><td align="left" valign="top">SR</td><td align="left" valign="top">3.1</td><td align="left" valign="top">N</td><td align="left" valign="top">Subgrid
 scale velocities can (optionally) be included for calculating the
 particle advection, using the method of Weil et al. (2004, JAS, 61,
 2877-2887). This method is switched on by the new particle package
-parameter <span style="font-family: monospace;">use_sgs_for_particles</span>. This also forces the Euler/upstream method to be used for time advancement of the TKE (see new parameter <span style="font-family: monospace;">use_upstream_for_tke</span>). The minimum timestep during the sub-timesteps is controlled by parameter <span style="font-family: monospace;">dt_min_part</span>. The data type <span style="font-family: monospace;">particle_type</span> has additional new attributes <span style="font-family: monospace;">e_m</span>, <span style="font-family: monospace;">dt_sum</span>, and <span style="font-family: monospace;">speed_x/y/z_sgs</span>.<br><br>Other new particle features: output of particle quantities as timeseries in NetCDF format (on local file <span style="font-family: monospace;">DATA_1D_PTS_NETCDF</span>). The output time interval is controlled by parameter <span style="font-family: monospace;">dt_dopts</span>. Particle advection can be switched off after some time using parameter <span style="font-family: monospace;">end_time_prel</span>. More than one particle per point can be started with parameter <span style="font-family: monospace;">particles_per_point</span>.<br><br>Additional parameter in function <span style="font-family: monospace;">random_gauss</span> which limits the range of the created random numbers to five times the standard deviation (=1). Seeds <span style="font-family: monospace;">iran</span> and <span style="font-family: monospace;">iran_part</span> are stored for restart runs.</td><td align="left" valign="top">advec_particles,
+parameter <span style="font-family: monospace;">use_sgs_for_particles</span>.
+This also forces the Euler/upstream method to be used for time
+advancement of the TKE (see new parameter <span style="font-family: monospace;">use_upstream_for_tke</span>).
+The minimum timestep during the sub-timesteps is controlled by
+parameter <span style="font-family: monospace;">dt_min_part</span>.
+The data type <span style="font-family: monospace;">particle_type</span>
+has additional new attributes <span style="font-family: monospace;">e_m</span>,
+<span style="font-family: monospace;">dt_sum</span>,
+and <span style="font-family: monospace;">speed_x/y/z_sgs</span>.<br><br>Other
+new particle features: output of particle quantities as timeseries in
+NetCDF format (on local file <span style="font-family: monospace;">DATA_1D_PTS_NETCDF</span>).
+The output time interval is controlled by parameter <span style="font-family: monospace;">dt_dopts</span>.
+Particle advection can be switched off after some time using parameter <span style="font-family: monospace;">end_time_prel</span>.
+More than one particle per point can be started with parameter <span style="font-family: monospace;">particles_per_point</span>.<br><br>Additional
+parameter in function <span style="font-family: monospace;">random_gauss</span>
+which limits the range of the created random numbers to five times the
+standard deviation (=1). Seeds <span style="font-family: monospace;">iran</span>
+and <span style="font-family: monospace;">iran_part</span>
+are stored for restart runs.</td><td align="left" valign="top">advec_particles,
 check_open, check_parameters, close_file, diffusion_e,
 disturb_heatflux, header, init_3d_model, init_particles, modules,
 netcdf, package_parin, parin, prognostic_equations, random_gauss,
 read_3d_binary, read_var_list, time_integration, write_3d_binary,
-write_var_list<br><br><span style="font-weight: bold;">new:</span><br>data_output_ptseries</td></tr><tr><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">C</td><td align="left" valign="top">Call of subroutine <span style="font-family: monospace;">header</span> is moved after call of subroutine <span style="font-family: monospace;">init_3d_model</span>. Generation of <span style="font-family: monospace;">run_description_header</span> is moved from routines <span style="font-family: monospace;">palm</span> and <span style="font-family: monospace;">header</span> to <span style="font-family: monospace;">check_parameters</span>. Determination of the number of particle groups is moved from <span style="font-family: monospace;">header</span> to <span style="font-family: monospace;">init_particles</span>.<br><br><span style="font-family: monospace;">izuf</span> renamed <span style="font-family: monospace;">iran</span>.</td><td align="left" valign="top">check_parameters, disturb_field, disturb_heatflux, header, init_particles, modules, palm</td></tr><tr><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">E</td><td align="left" valign="top">Variables
+write_var_list<br><br><span style="font-weight: bold;">new:</span><br>data_output_ptseries</td></tr><tr><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">C</td><td align="left" valign="top">Call of subroutine <span style="font-family: monospace;">header</span> is moved
+after call of subroutine <span style="font-family: monospace;">init_3d_model</span>.
+Generation of <span style="font-family: monospace;">run_description_header</span>
+is moved from routines <span style="font-family: monospace;">palm</span>
+and <span style="font-family: monospace;">header</span>
+to <span style="font-family: monospace;">check_parameters</span>.
+Determination of the number of particle groups is moved from <span style="font-family: monospace;">header</span> to <span style="font-family: monospace;">init_particles</span>.<br><br><span style="font-family: monospace;">izuf</span> renamed <span style="font-family: monospace;">iran</span>.</td><td align="left" valign="top">check_parameters,
+disturb_field, disturb_heatflux, header, init_particles, modules, palm</td></tr><tr><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">&nbsp;</td><td align="left" valign="top">E</td><td align="left" valign="top">Variables
 do2d_unit and do3d_unit now defined as 2d-arrays. Before, in some
 cases, the units of variables in the 2d- and 3d-NetCDF-files have been
-wrong.<br><br>In routine <span style="font-family: monospace;">poisfft</span> default setting of the thread number tn in case of not using OpenMP.<br><br>Error removed in the non-parallel part of routine <span style="font-family: monospace;">flow_statistics</span> (number of arguments of array <span style="font-family: monospace;">sums_l</span>).<br><br>Error removed which appeared if the user had decreased the value of <span style="font-family: monospace;">dt_dopr</span> within a job chain and if simulataneous time averaging of profiles was switched on.</td><td align="left" valign="top">check_parameters, data_output_profiles, flow_statistics, modules, netcdf, poisfft</td></tr><tr><td style="vertical-align: top;">22/08/06</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.1a</td><td style="vertical-align: top;">N</td><td style="vertical-align: top;">In case of more than one particle group, seperate output of timeseries for each of the groups.<br><br>New initial parameter <span style="font-family: monospace;">dz_max</span>, which limits the vertical gridspacing in case of a vertically stretched grid.</td><td style="vertical-align: top;">data_output_ptseries, header, init_grid, modules, netcdf, parin, read_var_list, write_var_list</td></tr><tr><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Disturbances are imposed only for the last Runge-Kutta-substep. <br><span style="font-weight: bold;">Attention:</span> <br>
-
+wrong.<br><br>In routine <span style="font-family: monospace;">poisfft</span> default
+setting of the thread number tn in case of not using OpenMP.<br><br>Error
+removed in the non-parallel part of routine <span style="font-family: monospace;">flow_statistics</span>
+(number of arguments of array <span style="font-family: monospace;">sums_l</span>).<br><br>Error
+removed which appeared if the user had decreased the value of <span style="font-family: monospace;">dt_dopr</span> within
+a job chain and if simulataneous time averaging of profiles was
+switched on.</td><td align="left" valign="top">check_parameters,
+data_output_profiles, flow_statistics, modules, netcdf, poisfft</td></tr><tr><td style="vertical-align: top;">22/08/06</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.1a</td><td style="vertical-align: top;">N</td><td style="vertical-align: top;">In case of more than one
+particle group, seperate output of timeseries for each of the groups.<br><br>New
+initial parameter <span style="font-family: monospace;">dz_max</span>,
+which limits the vertical gridspacing in case of a vertically stretched
+grid.</td><td style="vertical-align: top;">data_output_ptseries,
+header, init_grid, modules, netcdf, parin, read_var_list, write_var_list</td></tr><tr><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Disturbances are imposed only
+for the last Runge-Kutta-substep. <br><span style="font-weight: bold;">Attention:</span> <br>
 The content of the testresult-file (run control output) has changed
-due to this modification.<br><br>Output of xz and yz cross sections now up to gridpoint nzt+1.<br><br>Default settings of particle start positions changed.</td><td style="vertical-align: top;">check_open, close_file, data_output_2d, init_particles, netcdf, time_integration</td></tr><tr><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">E</td><td style="vertical-align: top;">Initialisation of all tendency arrays (t.._m) needed for the Runge-Kutta schemes.<br><br>Bugfix in index array initialization for line- or point-like topography structures.<br><br>Bugfix: yv coordinates are now used for yz cross sections, where neccessary.</td><td style="vertical-align: top;">init_3d_model, init_grid, netcdf</td></tr><tr><td style="vertical-align: top;">12/02/07</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.1b</td><td style="vertical-align: top;">&nbsp;N</td><td style="vertical-align: top;"><span style="font-weight: bold;">Attention:</span><br>This is the last revision before switching to the subversion revison control system!<br><br>Informative output to the job protocoll in case of restarts.</td><td style="vertical-align: top;">check_for_restart</td></tr><tr><td>&nbsp;&nbsp;</td><td>&nbsp;&nbsp;</td><td>&nbsp;&nbsp;</td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Values of environment variables are now read from local file <span style="font-family: monospace;">ENVPAR</span> instead of getting them by a system call.<br><br>For interpolation in <span style="font-family: monospace;">advec_particles</span>, allways level k is used.</td><td style="vertical-align: top;">advec_particles, local_tremain_ini, parin</td></tr><tr><td>&nbsp;</td><td>&nbsp;</td><td>&nbsp;</td><td style="vertical-align: top;">E</td><td style="vertical-align: top;">Several bugfixes in <span style="font-family: monospace;">advec_particles</span>: new particles are released only if <span style="font-family: monospace;">end_time_prel</span> &gt; <span style="font-family: monospace;">simulated_time</span>,&nbsp;transfer
+due to this modification.<br><br>Output of xz and yz cross
+sections now up to gridpoint nzt+1.<br><br>Default settings
+of particle start positions changed.</td><td style="vertical-align: top;">check_open, close_file,
+data_output_2d, init_particles, netcdf, time_integration</td></tr><tr><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">E</td><td style="vertical-align: top;">Initialisation of all tendency
+arrays (t.._m) needed for the Runge-Kutta schemes.<br><br>Bugfix
+in index array initialization for line- or point-like topography
+structures.<br><br>Bugfix: yv coordinates are now used for
+yz cross sections, where neccessary.</td><td style="vertical-align: top;">init_3d_model, init_grid, netcdf</td></tr><tr><td style="vertical-align: top;">12/02/07</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.1b</td><td style="vertical-align: top;">&nbsp;N</td><td style="vertical-align: top;"><span style="font-weight: bold;">Attention:</span><br>This
+is the last revision before switching to the subversion revison control
+system!<br><br>Informative output to the job protocoll in
+case of restarts.</td><td style="vertical-align: top;">check_for_restart</td></tr><tr><td>&nbsp;&nbsp;</td><td>&nbsp;&nbsp;</td><td>&nbsp;&nbsp;</td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Values of environment variables
+are now read from local file <span style="font-family: monospace;">ENVPAR</span>
+instead of getting them by a system call.<br><br>For
+interpolation in <span style="font-family: monospace;">advec_particles</span>,
+allways level k is used.</td><td style="vertical-align: top;">advec_particles,
+local_tremain_ini, parin</td></tr><tr><td>&nbsp;</td><td>&nbsp;</td><td>&nbsp;</td><td style="vertical-align: top;">E</td><td style="vertical-align: top;">Several bugfixes in <span style="font-family: monospace;">advec_particles</span>:
+new particles are released only if <span style="font-family: monospace;">end_time_prel</span>
+&gt; <span style="font-family: monospace;">simulated_time</span>,&nbsp;transfer
 of particles when x &lt; -0.5*dx (0.0 before), etc., index i,j used
 instead of cartesian (x,y) coordinate to check for transfer because
 this failed under very rare conditions, calculation of number of
-particles with same radius as the current particle (cloud droplet code).<br><br>Allocation of <span style="font-family: monospace;">tail_mask</span> and <span style="font-family: monospace;">new_tail_id</span> in case of restart-runs. "__" added in a cpp-directive. (both <span style="font-family: monospace;">init_particles</span>)</td><td style="vertical-align: top;">advec_particles, init_particles</td></tr><tr><td style="vertical-align: top;">02/03/07</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.1c</td><td style="vertical-align: top;">N</td><td style="vertical-align: top;">A heatflux can be prescribed at the top with new inipar parameters <span style="font-family: monospace;">top_heatflux</span> and <span style="font-family: monospace;">use_top_fluxes</span>. New 2d-arrays <span style="font-family: monospace;">qswst</span>, <span style="font-family: monospace;">qswst_m</span>, <span style="font-family: monospace;">tswst</span> and <span style="font-family: monospace;">tswst_m</span> are used to store this flux. Use of fluxes are controlled with new index variable <span style="font-family: monospace;">nzt_diff</span>. A Neumann boundary condition for temperature can be applied under these conditions.<br><br>Additionally, a Dirichlet condition for temperature can be used at the top.<br><br><span style="font-weight: bold;">Attention:</span> <br>
-
+particles with same radius as the current particle (cloud droplet code).<br><br>Allocation
+of <span style="font-family: monospace;">tail_mask</span>
+and <span style="font-family: monospace;">new_tail_id</span>
+in case of restart-runs. "__" added in a cpp-directive. (both <span style="font-family: monospace;">init_particles</span>)</td><td style="vertical-align: top;">advec_particles, init_particles</td></tr><tr><td style="vertical-align: top;">02/03/07</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.1c</td><td style="vertical-align: top;">N</td><td style="vertical-align: top;">A heatflux can be prescribed at
+the top with new inipar parameters <span style="font-family: monospace;">top_heatflux</span>
+and <span style="font-family: monospace;">use_top_fluxes</span>.
+New 2d-arrays <span style="font-family: monospace;">qswst</span>,
+<span style="font-family: monospace;">qswst_m</span>,
+<span style="font-family: monospace;">tswst</span>
+and <span style="font-family: monospace;">tswst_m</span>
+are used to store this flux. Use of fluxes are controlled with new
+index variable <span style="font-family: monospace;">nzt_diff</span>.
+A Neumann boundary condition for temperature can be applied under these
+conditions.<br><br>Additionally, a Dirichlet condition for
+temperature can be used at the top.<br><br><span style="font-weight: bold;">Attention:</span> <br>
 The content of the testresult-file (run control output) has changed
 due to this modification.</td><td style="vertical-align: top;">check_parameters,
 diffusion_s, flow_statistics, header, init_grid, init_3d_model,
 modules, parin, production_e, prognostic_equations, read_var_list,
-read_3d_binary, swap_timelevel, write_var_list, write_3d_binary</td></tr><tr><td>&nbsp;</td><td>&nbsp;</td><td>&nbsp;</td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Prognostic equations for all scalars are now solved up to gridpoint <span style="font-family: monospace;">nzt</span> (formerly <span style="font-family: monospace;">nzt-1</span>). Boundary conditions for scalars at top adjusted respectively (now applied only at <span style="font-family: monospace;">nzt+1</span>).<br><br>The default top boundary condition for temperature has been renamed to <span style="font-style: italic;">'initial_gradient'</span>.<br><br>Calls of <span style="font-family: monospace;">dvrp_output_local</span>, which were commented out for a long time, are now activated for all streams.</td><td style="vertical-align: top;">advec_s_pw,
+read_3d_binary, swap_timelevel, write_var_list, write_3d_binary</td></tr><tr><td>&nbsp;</td><td>&nbsp;</td><td>&nbsp;</td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Prognostic equations for all
+scalars are now solved up to gridpoint <span style="font-family: monospace;">nzt</span> (formerly <span style="font-family: monospace;">nzt-1</span>).
+Boundary conditions for scalars at top adjusted respectively (now
+applied only at <span style="font-family: monospace;">nzt+1</span>).<br><br>The
+default top boundary condition for temperature has been renamed to <span style="font-style: italic;">'initial_gradient'</span>.<br><br>Calls
+of <span style="font-family: monospace;">dvrp_output_local</span>,
+which were commented out for a long time, are now activated for all
+streams.</td><td style="vertical-align: top;">advec_s_pw,
 boundary_conds, calc_precipitation, check_parameters, diffusion_e,
 diffusion_s, impact_of_latent_heat, init_dvrp, init_pt_anomaly,
-modules, production_e, prognostic_equations, spline_z</td></tr><tr><td>&nbsp;</td><td>&nbsp;</td><td>&nbsp;</td><td style="vertical-align: top;">E</td><td style="vertical-align: top;">Bugfix: 3d-array <span style="font-family: monospace;">p</span> is not a pointer any more.<br><br>Bugfix in <span style="font-family: monospace;">init_particles</span>: <span style="font-family: monospace;">MPI_REAL</span> argument in <span style="font-family: monospace;">MPI_ALLREDUCE</span> replaced by<br><span style="font-family: monospace;">MPI_INTEGER</span> (caused error on NEC only).<br><br>Bugfix:&nbsp; <span style="font-family: monospace;">ddzw</span> now dimensioned <span style="font-family: monospace;">1:nzt"+1".</span></td><td style="vertical-align: top;">diffusion_e, diffusion_s, diffusion_u, diffusion_v, diffusion_w, init_particles, modules</td></tr>
-
-  
-  </tbody>
-</table>
-
-&nbsp;<b><blink>Attention:</blink></b> If<tt> make </tt>is not used
+modules, production_e, prognostic_equations, spline_z</td></tr><tr><td>&nbsp;</td><td>&nbsp;</td><td>&nbsp;</td><td style="vertical-align: top;">E</td><td style="vertical-align: top;">Bugfix: 3d-array <span style="font-family: monospace;">p</span> is not a
+pointer any more.<br><br>Bugfix in <span style="font-family: monospace;">init_particles</span>:
+<span style="font-family: monospace;">MPI_REAL</span>
+argument in <span style="font-family: monospace;">MPI_ALLREDUCE</span>
+replaced by<br><span style="font-family: monospace;">MPI_INTEGER</span>
+(caused error on NEC only).<br><br>Bugfix:&nbsp; <span style="font-family: monospace;">ddzw</span> now
+dimensioned <span style="font-family: monospace;">1:nzt"+1".</span></td><td style="vertical-align: top;">diffusion_e, diffusion_s,
+diffusion_u, diffusion_v, diffusion_w, init_particles, modules</td></tr>
+</tbody>
+</table>&nbsp;<b><blink>Attention:</blink></b>
+If<tt> make </tt>is not used
 for
 compiling,
-<b>all routines</b> have to be re-compiled after every change
+<b>all routines</b> have to be re-compiled after every
+change
 to <b>modules.f90!</b>
-<br>
-
-&nbsp;
-<br>
-
-&nbsp;
-<h2><a name="Kapitel2.0"></a>2.0&nbsp; How to change the model source
+<br>&nbsp;
+<br>&nbsp;
+<h2><a name="Kapitel2.0"></a>2.0&nbsp; How
+to change the model source
 code</h2>
-
 &Auml;nderungen, Erweiterungen oder Fehlerkorrekturen am Modellcode
 (gemeint
 sind alle Modellteile, die nicht zur benutzereigenen Software
 geh&ouml;ren)
-erfordern vom jeweiligen Benutzer Gruppenzugeh&ouml;rigkeit zur Gruppe
+erfordern vom jeweiligen Benutzer Gruppenzugeh&ouml;rigkeit zur
+Gruppe
 <tt>"palm"</tt>.
 Zu dieser Gruppe geh&ouml;ren alle Benutzer, die aktiv an der
 Modellentwicklung
 mitwirken (im weiteren "Modellarbeitsgruppe" genannt).
-<p>Die einzelnen Quelltextdateien des Modells werden mit RCS (Revision
+<p>Die einzelnen Quelltextdateien des Modells werden mit RCS
+(Revision
 Control System) verwaltet. Grundlegende Kenntnisse &uuml;ber dieses
 System
 werden im weiteren vorausgesetzt.
-</p>
-
-<p>&Auml;nderungen am Modellcode erfordern die Durchf&uuml;hrung
+</p><p>&Auml;nderungen am Modellcode erfordern die
+Durchf&uuml;hrung
 bestimmter
 Aktionen in einer festgelegten Reihenfolge, die nun beschrieben werden
 sollen.
-</p>
-
-<ol>
-
-  <li>Konzeptplanung.</li>
-
-  <li>Diskussion der geplanten &Auml;nderungen in der
-Modellarbeitsgruppe.</li>
-
-  <li>Auschecken der ben&ouml;tigten Quelltextdateien,
+</p><ol> <li>Konzeptplanung.</li> <li>Diskussion
+der geplanten &Auml;nderungen in der
+Modellarbeitsgruppe.</li> <li>Auschecken der
+ben&ouml;tigten Quelltextdateien,
 Durchf&uuml;hrung der
 &Auml;nderungen und Test des Modells. Die Modul-Datei <tt>modules.f90</tt>
 sollte nur dann ausgecheckt werden, wenn dort aufgef&uuml;hrte
 Variablen
-ver&auml;ndert werden. Neu geschaffene Variablen sollten w&auml;hrend
+ver&auml;ndert werden. Neu geschaffene Variablen sollten
+w&auml;hrend
 der
 Testphase in ein ebenfalls tempor&auml;r neu zu schaffendes Modul
@@ -4403,12 +2321,13 @@
 entsprechender USE-Anweisungen). Die in den einzelnen Dateien
 durchgef&uuml;hrten
-Arbeiten m&uuml;ssen in den Kopfzeilen (unter "Aktuelle Aenderungen")
+Arbeiten m&uuml;ssen in den Kopfzeilen (unter "Aktuelle
+Aenderungen")
 der
-Dateien vermerkt werden. Dies tut man am besten noch <b>bevor</b> man
+Dateien vermerkt werden. Dies tut man am besten noch <b>bevor</b>
+man
 die
 eigentlichen &Auml;nderungen durchf&uuml;hrt (sonst droht das
-Vergessen...).</li>
-
-  <li>Zum Abschlu&szlig; der Tests sollte mit dem Modell eine
+Vergessen...).</li> <li>Zum Abschlu&szlig; der Tests
+sollte mit dem Modell eine
 Kontrollrechnung
 durchgef&uuml;hrt und die Ergebnisse mit denen des sogenannten
@@ -4418,5 +2337,6 @@
 durchzuf&uuml;hren, die auch zum Testen der Modellinstallation
 verwendet
-werden (s. <a href="Anwendung/Kapitel_5.0.html">Kapitel 5.0</a> der
+werden (s. <a href="Anwendung/Kapitel_5.0.html">Kapitel
+5.0</a> der
 Anwenderdokumentation).
 Die Ergebnisdatei mit dem lokalen Namen <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM-1/Dokumentationen/Anwendung/Kapitel_3.4.html#RUN_CONTROL">RUN_CONTROL</a>
@@ -4430,27 +2350,31 @@
 so
 ist unbedingt eine Diskussion in der Modellarbeitsgruppe notwendig.</li>
-
-  <li>Nach Durchf&uuml;hrung der Tests m&uuml;ssen bei Bedarf die
+<li>Nach Durchf&uuml;hrung der Tests m&uuml;ssen bei
+Bedarf die
 vorl&auml;ufig
-im tempor&auml;ren Modul eingetragenen Variablen in die eigentlichen
+im tempor&auml;ren Modul eingetragenen Variablen in die
+eigentlichen
 Module
 des Modells &uuml;bertragen werden. Dazu ist die Modul-Datei <tt>modules.f90</tt>
 auszuchecken. Nach Eintrag der Variablen sollte das Modell noch einmal
-testweise &uuml;bersetzt werden, um Sytaxfehler auszuschlie&szlig;en
+testweise &uuml;bersetzt werden, um Sytaxfehler
+auszuschlie&szlig;en
 (evtl.
 auch noch einmal einen Probelauf starten). Die Verwendung des
 tempor&auml;ren
 Moduls ist aus allen entsprechenden Programmteilen zu entfernen.</li>
-
-  <li>Vergabe einer neuen Versionsnummer (Variable <tt>version</tt> im
+<li>Vergabe einer neuen Versionsnummer (Variable <tt>version</tt>
+im
 Modul <tt>steuer</tt>),
-Bsp.: 2.3c. Kleinere &Auml;nderungen werden durch angeh&auml;ngten
+Bsp.: 2.3c. Kleinere &Auml;nderungen werden durch
+angeh&auml;ngten
 kleinen
-Buchstaben kenntlich gemacht, gr&ouml;&szlig;ere &Auml;nderungen durch
+Buchstaben kenntlich gemacht, gr&ouml;&szlig;ere
+&Auml;nderungen durch
 die Zahl hinter dem Punkt. Eine Inkrementierung der Zahl vor dem Punkt
-geschieht nur bei ganz grunds&auml;tzlichen Revisionen des Modells, die
-weite Teile des Modellcodes betreffen.</li>
-
-  <li>Nun k&ouml;nnen die ausgecheckten Dateien wieder eingecheckt
+geschieht nur bei ganz grunds&auml;tzlichen Revisionen des Modells,
+die
+weite Teile des Modellcodes betreffen.</li> <li>Nun
+k&ouml;nnen die ausgecheckten Dateien wieder eingecheckt
 werden.
 Vor
@@ -4462,22 +2386,19 @@
 der vom RCS-System vergebenen Versionsnummer (daf&uuml;r sorgt die
 $Log:
-... $ - Zeile, die in jeder Quelltextdatei steht).</li>
-
-  <li>Eine Zusammenfassung der durchgef&uuml;hrten Arbeiten ist im <a href="#Kapitel1.0">Kapitel
-1.0</a> einzutragen.</li>
-
-  <li>Die Modellarbeitsgruppe ist &uuml;ber den Abschlu&szlig; der
+... $ - Zeile, die in jeder Quelltextdatei steht).</li> <li>Eine
+Zusammenfassung der durchgef&uuml;hrten Arbeiten ist im <a href="#Kapitel1.0">Kapitel
+1.0</a> einzutragen.</li> <li>Die
+Modellarbeitsgruppe ist &uuml;ber den Abschlu&szlig; der
 durchgef&uuml;hrten
 Arbeiten zu informieren.</li>
-
-</ol>
-
-<h2>
-<a name="Kapitel3.0"></a>3.0&nbsp; Description of selected parts of the
+</ol><h2>
+<a name="Kapitel3.0"></a>3.0&nbsp; Description
+of selected parts of the
 model source code</h2>
-
 <b>This chapter is still under construction (Feb. 04).</b>
-<p>Dieser Abschnitt verweist auf vorhandene Beschreibungen der genauen
-Funktionsweise einzelner, spezieller Modellteile. Dazu geh&ouml;rt z.B.
+<p>Dieser Abschnitt verweist auf vorhandene Beschreibungen der
+genauen
+Funktionsweise einzelner, spezieller Modellteile. Dazu geh&ouml;rt
+z.B.
 die Beschreibung bestimmter numerischer Verfahren (auch ihr
 theoretischer
@@ -4488,55 +2409,28 @@
 Aufgaben (z.B. der Berechnung horizontal gemittelter
 Gr&ouml;&szlig;en).
-<br>
-
-Die Beschreibungen liegen in jeweils getrennten Dokumenten, in der
+<br>Die Beschreibungen liegen in jeweils getrennten Dokumenten,
+in der
 Regel als TeX- und/oder Postscript-Datei vor und sind &uuml;ber die
 folgenden
 Links zug&auml;nglich.
-<br>
-
+<br>&nbsp;
+</p><ul> <li> <a name="UPS"></a>upstream-spline
+advection scheme ( <a href="Upstream_Spline/ups.ps">.ps</a>
+)</li> <li> <a name="Cloud_physics"></a>cloud
+physics module ( <a href="Cloud_physics/wolken.pdf">.pdf</a>
+)</li> <li> <a name="Bitkompression"></a>Datenkompression
+mit dem Verfahren
+der <a href="Bitkompression/kompression.html">Bitverschiebung</a></li>
+<li> <a name="Runge-kutta"></a>Runge-Kutta
+time integration scheme ( <a href="numerik.heiko/zeitschrittverfahren.pdf">.pdf</a>
+)</li>
+</ul><br>
 &nbsp;
-</p>
-
-<ul>
-
-  <li> <a name="UPS"></a>upstream-spline advection scheme ( <a href="Upstream_Spline/ups.ps">.ps</a>
-)</li>
-
-  <li> <a name="Cloud_physics"></a>cloud physics module ( <a href="Cloud_physics/wolken.pdf">.pdf</a>
-)</li>
-
-  <li> <a name="Bitkompression"></a>Datenkompression mit dem Verfahren
-der <a href="Bitkompression/kompression.html">Bitverschiebung</a></li>
-
-  <li> <a name="Runge-kutta"></a>Runge-Kutta time integration scheme (
-    <a href="numerik.heiko/zeitschrittverfahren.pdf">.pdf</a>
-)</li>
-
-</ul>
-
-<br>
-
-&nbsp;
-<br>
-
-&nbsp;
-<br>
-
-&nbsp;
-<br>
-
-&nbsp;
-<br>
-
-&nbsp;
-<br>
-
-&nbsp;
-<br>
-
-&nbsp;
-<br>
-
-&nbsp;
+<br>&nbsp;
+<br>&nbsp;
+<br>&nbsp;
+<br>&nbsp;
+<br>&nbsp;
+<br>&nbsp;
+<br>&nbsp;
 </body></html>
