Index: palm/trunk/DOC/tec/technical_documentation.html
===================================================================
--- palm/trunk/DOC/tec/technical_documentation.html	(revision 552)
+++ palm/trunk/DOC/tec/technical_documentation.html	(revision 553)
@@ -1,1461 +1,1542 @@
-<!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;
-</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 fixed release:</b></td> <td><span style="font-weight: bold;">3.7a<br>
-</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 source
-code and how to document modifications<br>
- </p> <p><a href="#Kapitel3.0">3.0</a>&nbsp;
-Description of selected model numerics<br>
-&nbsp;</p>
-</ul><b></b><br>
-
-<h2><a name="Kapitel1.0"></a>1.0&nbsp;
-History of changes</h2>
-The following table documents all changes to the PALM repository since 2009. Older changes can be found <a href="tec_doc_1998-2008.html">here</a>.<br>
-<br>
-Entries in the fifth column specify the type of the change:
-change&nbsp;
-of existing code (C), new code (N) or bugfix (B). More detailed
-information
-about the changes may sometimes be found in the header comment lines of the
-respective
-routines.
-<br>&nbsp;
-<br>
-<table style="text-align: left; width: 100%; vertical-align: top;" border="1" cellpadding="2" cellspacing="2">
-  <tbody>
-    <tr>
-      <td style="vertical-align: middle; width: 27px;"><big><span style="font-weight: bold;">date</span></big><br>
-      </td>
-      <td style="width: 38px;"><big><span style="font-weight: bold;">author</span></big><br>
-      </td>
-      <td style="width: 42px;"><big><span style="font-weight: bold;">svn revision</span></big><br>
-      </td>
-      <td style="width: 40px;"><big><span style="font-weight: bold;">last release</span></big><br>
-      </td>
-      <td style="width: 27px;"><big><span style="font-weight: bold;">type</span></big><br>
-      </td>
-      <td style="width: 497px;"><big><span style="font-weight: bold;">description</span></big><br>
-      </td>
-      <td style="width: 189px;"><big><span style="font-weight: bold;">affected routines</span></big><br>
-      </td>
-    </tr>
-    
-    
-    
-    
-    
-    
-    
-    
-    
-    <tr>
-      <td style="vertical-align: top;">17/08/10<br>
-      </td>
-      <td style="vertical-align: top;">BM<br>
-      </td>
-      <td style="vertical-align: top;">544<br>
-      </td>
-      <td style="vertical-align: top;">3.7a<br>
-      </td>
-      <td style="vertical-align: top;">B/C/N<br>
-      </td>
-      <td style="vertical-align: top;">adjustments for lcxt4; bugfix (mbuild): loading of user-defined modules; new .mrun.config.lcxt4 for Cray XT4 in Bergen/Norway<br>
-      </td>
-      <td style="vertical-align: top;">mrun, mbuild, subjob<br>
-      </td>
-    </tr>
-
-    <tr>
-      <td style="vertical-align: top;">15/06/10<br>
-      </td>
-      <td style="vertical-align: top;">RH<br>
-      </td>
-      <td style="vertical-align: top;">540<br>
-      </td>
-      <td style="vertical-align: top;">3.7a<br>
-      </td>
-      <td style="vertical-align: top;">C<br>
-      </td>
-      <td style="vertical-align: top;">ssh from compute nodes to login
-nodes on hosts lcsgih/lcsgib is done using the names
-hicegate0/bicegate0 instead of the IP addresses<br>
-      </td>
-      <td style="vertical-align: top;">mrun<br>
-      </td>
-    </tr>
-<tr>
-      <td style="vertical-align: top;">21/04/10<br>
-      </td>
-      <td style="vertical-align: top;">RH<br>
-      </td>
-      <td style="vertical-align: top;">532<br>
-      </td>
-      <td style="vertical-align: top;">3.7a<br>
-      </td>
-      <td style="vertical-align: top;">C<br>
-      </td>
-      <td style="vertical-align: top;">
-Large scale vertical motion (subsidence/ascent) can be added as
-additional source to the prognostic equation for humidity/passive scalar<br>
-&nbsp;<br>
-character length for dopr_unit enlarged<br>
-      <br>
-Adjustment of the NCL scripts<span style="font-family: Courier New,Courier,monospace;"> </span>to allow for output of png files <br>
- </td>
-      <td style="vertical-align: top;">modules, prognostic_equations<br>
-      <br>
-.ncl.config.default, cross_sections.ncl, profiles.ncl, timeseries.ncl, spectra.ncl<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;"><br>
-      </td>
-      <td style="vertical-align: top;">B<br>
-      </td>
-      <td style="vertical-align: top;">Bugfix: unit of hyp changed to dbar<br>
-      </td>
-      <td style="vertical-align: top;">check_parameters<br>
-      </td>
-    </tr>
-<tr>
-      <td style="vertical-align: top;">01/04/10<br>
-      </td>
-      <td style="vertical-align: top;">RH<br>
-      </td>
-      <td style="vertical-align: top;">530<br>
-      </td>
-      <td style="vertical-align: top;">3.7a<br>
-      </td>
-      <td style="vertical-align: top;">B<br>
-      </td>
-      <td style="vertical-align: top;">Bugfix: treatment of special character * in NCL variables containing lists corrected<br>
-      </td>
-      <td style="vertical-align: top;">palmplot<br>
-      </td>
-    </tr>
-<tr>
-      <td style="vertical-align: top;">31/03/10<br>
-      </td>
-      <td style="vertical-align: top;">RH<br>
-      </td>
-      <td style="vertical-align: top;">529<br>
-      </td>
-      <td style="vertical-align: top;">3.7a<br>
-      </td>
-      <td style="vertical-align: top;">C<br>
-      </td>
-      <td style="vertical-align: top;">List of variable names have to specified using blank separated lists and single quotes<br>
-      </td>
-      <td style="vertical-align: top;">palmplot<br>
-      </td>
-    </tr>
-<tr>
-      <td style="vertical-align: top;">30/03/10<br>
-      </td>
-      <td style="vertical-align: top;">RH<br>
-      </td>
-      <td style="vertical-align: top;">526<br>
-      </td>
-      <td style="vertical-align: top;">3.7a<br>
-      </td>
-      <td style="vertical-align: top;">C<br>
-      </td>
-      <td style="vertical-align: top;">Adjustment of the NCL scripts and <span style="font-family: Courier New,Courier,monospace;">palmplot </span>to allow for special characters in NetCDF variable names <br>
-      </td>
-      <td style="vertical-align: top;">palmplot, .ncl.config.default, cross_sections.ncl, profiles.ncl, timeseries.ncl, spectra.ncl<br>
-      </td>
-    </tr>
-<tr>
-      <td style="vertical-align: top;">30/03/10<br>
-      </td>
-      <td style="vertical-align: top;">SR<br>
-      </td>
-      <td style="vertical-align: top;">525<br>
-      </td>
-      <td style="vertical-align: top;">3.7a<br>
-      </td>
-      <td style="vertical-align: top;">B<br>
-      </td>
-      <td style="vertical-align: top;">Bugfix: "/" in netcdf profile variable names replaced by ":"<br>
-      </td>
-      <td style="vertical-align: top;">check_parameters<br>
-      </td>
-    </tr>
-<tr>
-      <td style="vertical-align: top;">19/03/10<br>
-      </td>
-      <td style="vertical-align: top;">SR<br>
-      </td>
-      <td style="vertical-align: top;">520<br>
-      </td>
-      <td style="vertical-align: top;">3.7a<br>
-      </td>
-      <td style="vertical-align: top;">N<br>
-      </td>
-      <td style="vertical-align: top;">NetCDF4 support for particle
-data. Total number of particles is allowed to increase during the run,
-when netcdf_data_format &gt; <span style="font-style: italic;">2</span> is used.<br>
-      </td>
-      <td style="vertical-align: top;">advec_particles, check_open, 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;"><br>
-      </td>
-      <td style="vertical-align: top;">C<br>
-      </td>
-      <td style="vertical-align: top;">Special characters like * and "
-are now allowed for NetCDF variable names. Replacement of these
-characters has been removed (routine clean_netcdf_varname removed)<br>
-      </td>
-      <td style="vertical-align: top;"><br>
-      </td>
-    </tr>
-<tr>
-      <td style="vertical-align: top;">18/03/10<br>
-      </td>
-      <td style="vertical-align: top;">SR<br>
-      </td>
-      <td style="vertical-align: top;">516<br>
-      </td>
-      <td style="vertical-align: top;">3.7a<br>
-      </td>
-      <td style="vertical-align: top;">N<br>
-      </td>
-      <td style="vertical-align: top;">Simple installation method has been documented (see chapter 5.0). Respective scripts and configuration files have been updated.<br>
-      </td>
-      <td style="vertical-align: top;">palm_simple_install, palm_simple_run<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;"><br>
-      </td>
-      <td style="vertical-align: top;">B<br>
-      </td>
-      <td style="vertical-align: top;">opt=1 is forced for PGI-compiler in routine<span style="font-family: Courier New,Courier,monospace;"> buoyancy </span>by compiler-directive. Otherwise, program aborts with SIGFPE.<br>
-      </td>
-      <td style="vertical-align: top;">buoyancy<br>
-      </td>
-    </tr>
-<tr>
-      <td style="vertical-align: top;">16/03/10<br>
-      </td>
-      <td style="vertical-align: top;">RH<br>
-      </td>
-      <td style="vertical-align: top;">513<br>
-      </td>
-      <td style="vertical-align: top;">3.7a<br>
-      </td>
-      <td style="vertical-align: top;">N<br>
-      </td>
-      <td style="vertical-align: top;">NCL scripts can be run with the shell script <span style="font-family: Courier New,Courier,monospace;">palmplot.</span><br>
-      </td>
-      <td style="vertical-align: top;">palmplot, .ncl.config.default, cross_sections.ncl, profiles.ncl<br>
-      </td>
-    </tr>
-<tr>
-      <td style="vertical-align: top;">12/03/10<br>
-      </td>
-      <td style="vertical-align: top;">SR<br>
-      </td>
-      <td style="vertical-align: top;">509<br>
-      </td>
-      <td style="vertical-align: top;">3.7a<br>
-      </td>
-      <td style="vertical-align: top;">B<br>
-      </td>
-      <td style="vertical-align: top;">bugfix for rev 504: remode_md was given a wrong value<br>
-      </td>
-      <td style="vertical-align: top;">mbuild<br>
-      </td>
-    </tr>
-<tr>
-      <td style="vertical-align: top;">08/03/10<br>
-      </td>
-      <td style="vertical-align: top;">SR<br>
-      </td>
-      <td style="vertical-align: top;">506<br>
-      </td>
-      <td style="vertical-align: top;">3.7a<br>
-      </td>
-      <td style="vertical-align: top;">C<br>
-      </td>
-      <td style="vertical-align: top;">hack for better usage of <span style="font-weight: bold;">mrun</span> on Kyoto computing center machine (<span style="font-family: Courier New,Courier,monospace;">lckyoto</span>)<br>
-      </td>
-      <td style="vertical-align: top;">mrun<br>
-      </td>
-    </tr>
-<tr>
-      <td style="vertical-align: top;">03/03/10<br>
-      </td>
-      <td style="vertical-align: top;">SR<br>
-      </td>
-      <td style="vertical-align: top;">504<br>
-      </td>
-      <td style="vertical-align: top;">3.7a<br>
-      </td>
-      <td style="vertical-align: top;">C<br>
-      </td>
-      <td style="vertical-align: top;"><span style="font-family: Courier New,Courier,monospace;">combine_plot_fields</span> is now directly called within <span style="font-weight: bold;">mrun</span> and must not be given by an output command in the configuration file any more. Please remove the respective output command (<span style="font-family: Courier New,Courier,monospace;">OC:combine_plot_fields.x</span>) from your configuration file!<br>
-      <br>
-Utility programs (<span style="font-family: Courier New,Courier,monospace;">combine_plot_fields</span>, <span style="font-family: Courier New,Courier,monospace;">interpret_config</span>) are created by <span style="font-weight: bold;">mbuild</span> separately for each block in the configuration file. Depending on option <span style="font-family: Courier New,Courier,monospace;">-K</span>, <span style="font-weight: bold;">mrun</span> is using the utility programs for the respective block.<br>
-      <br>
-Utility programs are also compiled in parallel, depending on variable <span style="font-family: Courier New,Courier,monospace;">%mopts</span> in the configuration file..<br>
-      <br>
-Default configuration files updated.<br>
-      </td>
-      <td style="vertical-align: top;">mbuild, mrun, UTIL/Makefile, SCRIPTS/.mrun.config.*<br>
-      </td>
-    </tr>
-<tr>
-      <td style="vertical-align: top;">02/03/10<br>
-      </td>
-      <td style="vertical-align: top;">SR<br>
-      </td>
-      <td style="vertical-align: top;">496<br>
-      </td>
-      <td style="vertical-align: top;">3.7a<br>
-      </td>
-      <td style="vertical-align: top;">C<br>
-      </td>
-      <td style="vertical-align: top;">Particle data for restart runs are only written if environment variable <span style="font-family: Courier New,Courier,monospace;">write_binary</span> = <span style="font-style: italic;">.T.</span>, i.e. particle data are now handled in the same way as the other restart data.<br>
-      </td>
-      <td style="vertical-align: top;">palm<br>
-      </td>
-    </tr>
-<tr>
-      <td style="vertical-align: top;">01/03/10<br>
-      </td>
-      <td style="vertical-align: top;">SR<br>
-      </td>
-      <td style="vertical-align: top;">494<br>
-      </td>
-      <td style="vertical-align: top;">3.7a<br>
-      </td>
-      <td style="vertical-align: top;">N<br>
-      </td>
-      <td style="vertical-align: top;">Output in NetCDF4-format.<br>
-      <br>
-NetCDF4 allows for parallel I/O. Also, variables are not limited in
-size any more (i.e. arrays can be larger than 2GB). The new
-d3par-parameter <span style="font-family: Courier New,Courier,monospace;">netcdf_data_format</span>.can be used to select among the different NetCDF file formats (<span style="font-style: italic;">1</span>: classic format, <span style="font-style: italic;">2</span>: 64-bit-offset format (this is now the default), <span style="font-style: italic;">3</span>: NetCDF4/HDF5 format, <span style="font-style: italic;">4</span>: NetCDF4 format with NF90_CLASSIC_MODEL bit set). Setting <span style="font-family: Courier New,Courier,monospace;">netcdf_data_format</span> &gt; <span style="font-style: italic;">2</span>
-only affects 2d- and 3d-output, i.e. there is no parallel output for
-profiles, time series, spectra, and particles. Parallel output of
-particle data may realized in the future.<br>
-      <br>
-In case of parallel output, all PEs write into one file, which means that postprocessing the data with <span style="font-family: Courier New,Courier,monospace;">combine_plot_fields</span>
-(merging data from single files - one per PE - into one file) is not
-neccessary. In this case, flag files are automatically set, telling <span style="font-family: Courier New,Courier,monospace;">combine_plot_fields</span> that nothing is to be done.<br>
-      <br>
-Using NetCDF4 format requires a NetCDF4/HDF library. Respective
-settings for -I, -L, and -l compiler/loader options for the SGI-ICE at
-HLRN can be found in the new example configuration file <span style="font-family: Courier New,Courier,monospace;">..../trunk/SCRIPTS/.mrun.config.imuk_ice2_netcdf4</span><br>
-      <br>
-      <span style="font-weight: bold;">Attention:</span><br>
-So far, NetCDF parallel I/O has only been tested a little. Further
-performance optimizations may be required, especially in case of large
-data sets (e.g. data chunking, or lustre file system settings with <span style="font-family: Courier New,Courier,monospace;">lfs setstripe</span>). <span style="font-weight: bold;">Please report bugs and performance problems immediately!</span><br>
-      <br>
-      <span style="font-family: Courier New,Courier,monospace;">combine_plot_fields</span> outputs cpu-time needed for file merging.<br>
-      <br>
-Modules to be loaded for compilation (mbuild) or job execution (mrun) can/should now be given in the configuration file (<span style="font-family: Courier New,Courier,monospace;">.mrun.config</span>) using variable <span style="font-family: Courier New,Courier,monospace;">modules</span>. Example:<br>
-      <br>
-      <span style="font-family: Courier New,Courier,monospace;">&nbsp;&nbsp; %modules&nbsp;&nbsp;&nbsp;&nbsp; ifort/11.0.069:netcdf&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; lcsgih parallel</span><br>
-      <br>
-This method replaces the (undocumented) <span style="font-family: Courier New,Courier,monospace;">mpilib</span>-variable.<br>
-      <br>
-      <span style="font-weight: bold;">WARNING:</span><br>
-      <span style="font-weight: bold;">All fixed settings of modules in the scripts mbuild, mrun, and subjob have been removed! Please set the <span style="font-family: Courier New,Courier,monospace;">modules</span> variable appropriately and remove the mpilib variable in your configuration file.<br>
-      <br>
-</span>Compilation time can be reduced be using parallel make. The
-number of allowed simultaneous compiler calls (threads) is set in the
-configuration file using variable <span style="font-family: Courier New,Courier,monospace;">mopt</span>, e.g.<br>
-      <br>
-      <span style="font-family: Courier New,Courier,monospace;">&nbsp;&nbsp; %mopts&nbsp; -j:4 &nbsp; lcsgih parallel</span><br>
-      <br>
-allows a maximum of 4 routines to be compiled simultaneously (provided
-that the makefile does not show any dependencies between these
-routines).<span style="font-weight: bold;"><br>
-      </span>
-      </td>
-      <td style="vertical-align: top;">check_open, check_parameters, close_file, data_output_2d, data_output_3d, header, modules, netcdf, parin<br>
-      <br>
-combine_plot_fields<br>
-      <br>
-mbuild, mrun, subjob<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;"><br>
-      </td>
-      <td style="vertical-align: top;">C<br>
-      </td>
-      <td style="vertical-align: top;">d3par-parameters <span style="font-family: Courier New,Courier,monospace;">netcdf_64bit</span> and <span style="font-family: Courier New,Courier,monospace;">netcdf_64bit_3d</span> have been removed. Use the new d3par-variable <span style="font-family: Courier New,Courier,monospace;">netcdf_data_format</span> = <span style="font-style: italic;">2</span>
-for choosing the NetCDF 64bit-offset format (this is the default from
-now on). The offset-format can not be set independently for the
-3d-output-data any more.<br>
-      <br>
-d3par-parameters <span style="font-family: Courier New,Courier,monospace;">netcdf_format_mask</span>, <span style="font-family: Courier New,Courier,monospace;">netcdf_format_mask_av</span>, and variables <span style="font-family: Courier New,Courier,monospace;">nc_format_mask</span>, <span style="font-family: Courier New,Courier,monospace;">format_parallel_io</span> removed. They are replaced by the new parameter <span style="font-family: Courier New,Courier,monospace;">netcdf_data_format</span>. This change concerns the masked data output, which is still not documented.<br>
-      </td>
-      <td style="vertical-align: top;">check_open, close_file,<br>
-data_output_mask, header, init_masks, modules, parin<br>
-      <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;"><br>
-      </td>
-      <td style="vertical-align: top;">B<br>
-      </td>
-      <td style="vertical-align: top;">Bugfix in trunk/UTIL/Makefile: forgot to compile for interpret_config.<br>
-      <br>
-Bugfix: timeseries data have to be collected by PE0 (user_statistics).<br>
-      </td>
-      <td style="vertical-align: top;">user_statistics<br>
-      <br>
-trunk/UTIL/Makefile<br>
-      </td>
-    </tr>
-<tr>
-      <td style="vertical-align: top;">05/02/10<br>
-      </td>
-      <td style="vertical-align: top;">SR<br>
-      </td>
-      <td style="vertical-align: top;">485<br>
-      </td>
-      <td style="vertical-align: top;">3.7a<br>
-      </td>
-      <td style="vertical-align: top;">B<br>
-      </td>
-      <td style="vertical-align: top;"><span style="font-family: Courier New,Courier,monospace;">ngp_3d</span>, <span style="font-family: Courier New,Courier,monospace;">ngp_3d_inner</span>
-changed to 64 bit. Otherwise, on large grids &gt; 2*10**9 grid points
-some global statistics (e.g. total kinetic energy) have wrong values.<br>
-      </td>
-      <td style="vertical-align: top;">init_3d_model, modules<br>
-      </td>
-    </tr>
-    <tr>
-      <td style="vertical-align: top;">05/02/10<br>
-      </td>
-      <td style="vertical-align: top;">SR<br>
-      </td>
-      <td style="vertical-align: top;">484<br>
-      </td>
-      <td style="vertical-align: top;">3.7a<br>
-      </td>
-      <td style="vertical-align: top;">B<br>
-      </td>
-      <td style="vertical-align: top;">Typo in file headers removed ("current" instead of "actual").<br>
-      </td>
-      <td style="vertical-align: top;">many<br>
-      </td>
-    </tr>
-<tr>
-      <td style="vertical-align: top; width: 27px;">05/02/10<br>
-      </td>
-      <td style="vertical-align: top; width: 38px;">SR<br>
-      </td>
-      <td style="vertical-align: top; width: 42px;">483<br>
-      </td>
-      <td style="vertical-align: top; width: 40px;">3.7a<br>
-      </td>
-      <td style="vertical-align: top; width: 27px;">N<br>
-      </td>
-      <td style="vertical-align: top; width: 497px;"><span style="font-family: Courier New,Courier,monospace;">compare_palm_logs </span>is additionally compiled with <span style="font-weight: bold;">mbuild</span> <span style="font-family: Courier New,Courier,monospace;">-u</span> . (<span style="font-family: Courier New,Courier,monospace;">Makefile</span> in <span style="font-family: Courier New,Courier,monospace;">trunk/UTIL</span>)<br>
-      <br>
-make options (<span style="font-family: Courier New,Courier,monospace;">mopts</span>) to be set by configuration file implemented. (<span style="font-weight: bold;">mrun</span>, <span style="font-weight: bold;">mbuild</span>)<br>
-      <br>
-      <span style="font-family: Courier New,Courier,monospace;">humidity</span> = <span style="font-style: italic;">.T.</span> can be used for runs with topography. <span style="font-family: Courier New,Courier,monospace;">wall_humidityflux</span> and<br><span style="font-family: Courier New,Courier,monospace;">
-wall_scalarflux</span> are the corresponding new inipar parameters.<br>
-(<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>, <span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>, <span style="font-family: Courier New,Courier,monospace;">parin</span>)<br>
-      <br>
-Large scale vertical motion (subsidence/ascent) can be added as
-additional source to the prognostic equation for potential temperature.
-See new parameters <span style="font-family: Courier New,Courier,monospace;">large_scale_subsidence</span>, <span style="font-family: Courier New,Courier,monospace;">ws_vertical_gradient</span>, and <span style="font-family: Courier New,Courier,monospace;">ws_vertical_gradient_level</span>. (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>, <span style="font-family: Courier New,Courier,monospace;">header</span>, <span style="font-family: Courier New,Courier,monospace;">Makefile</span>, <span style="font-family: Courier New,Courier,monospace;">modules</span>, <span style="font-family: Courier New,Courier,monospace;">parin</span>, <span style="font-family: Courier New,Courier,monospace;">prognostic_equations</span>, <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>, <span style="font-family: Courier New,Courier,monospace;">subsidence</span>, <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br>
-      <br>Scripts have been updated for Cray-XT5 machine (lcxt5m) and Hitachi-machine (lckyoto).<br>
-      <br>
-A simple method for installing and running palm (with limited features) has been added. (<span style="font-family: Courier New,Courier,monospace;">Makefile</span>, new scripts <span style="font-weight: bold;">palm_simple_install</span>, <span style="font-weight: bold;">palm_simple_run</span>)<br>
-      <br>
-Masked data output has been implemented but it is not documented in the
-user manual so far. This feature can already use NetCDF4 file format.
-(check_open, check_parameters, close_file, modules, netcdf, parin,
-time_integration, user_parin, new: data_output_mask, init_masks,
-user_data_output_mask)<br>
-      </td>
-      <td style="vertical-align: top; width: 189px;">check_open, check
-parameters, close_file, header, init_3d_model, Makefile, mbuild, modules, mrun,
-netcdf, parin, prognostic_equations, read_var_list, time_integration, user_parin, write_var_list<br>
-      <br><span style="font-weight: bold;">
-new:</span> data_output_mask,&nbsp; init_masks, palm_simple_install, palm_simple_run, subsidence, user_data_output_mask<br>
-<br>
-      </td>
-    </tr>
-    <tr>
-      <td style="vertical-align: top; width: 27px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 38px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 42px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 40px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 27px;">C<br>
-      </td>
-      <td style="vertical-align: top; width: 497px;">2d-decomposition is default for Cray-XT machines. (<span style="font-family: Courier New,Courier,monospace;">init_pegrid</span>)<br>
-      <br><span style="font-family: Courier New,Courier,monospace;">
-var_ts</span> is replaced by <span style="font-family: Courier New,Courier,monospace;">dots_max</span>. (<span style="font-family: Courier New,Courier,monospace;">modules</span>, <span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>)<br>
-      <br>
-Every cloud droplet has now an own weighting factor and can be deleted
-due to collisions. Condensation and collision of cloud droplets are
-adjusted accordingly. (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>)<br>
-      <br>
-Collision efficiency for large cloud droplets has changed according to table of Rogers and Yau. (<span style="font-family: Courier New,Courier,monospace;">collision_efficiency</span>)<br>
-      </td>
-      <td style="vertical-align: top; width: 189px;">advec_particles, collision_efficiency, init_3d_model, init_pegrid, modules<br>
-      </td>
-    </tr>
-    <tr>
-      <td style="vertical-align: top; width: 27px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 38px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 42px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 40px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 27px;">B<br>
-      </td>
-      <td style="vertical-align: top; width: 497px;">Bugfix for generating serial jobs. (<span style="font-weight: bold;">subjob</span>)<br>
-      <br>
-Bugfix: index problem concerning gradient_level indices removed. (<span style="font-family: Courier New,Courier,monospace;">header</span>)<br>
-      <br>
-Dimension of array <span style="font-family: Courier New,Courier,monospace;">stat</span> in cascade change to prevent type problems with mpi2 libraries. (<span style="font-family: Courier New,Courier,monospace;">poisfft_hybrid</span>)<br>
-      <br>
-Loop was split to make runs reproducible when using ifort compiler. (<span style="font-family: Courier New,Courier,monospace;">disturb_field</span>)<br>
-      <br>
-Bugfix: exchange of ghost points for <span style="font-family: Courier New,Courier,monospace;">prho</span> included. (<span style="font-family: Courier New,Courier,monospace;">time_integration</span>)<br>
-      <br>
-Bugfix in calculation of time-averaged surface heatfluxes. (<span style="font-family: Courier New,Courier,monospace;">sum_up_3d_data</span>)<br>
-      <br>
-Bugfix in calculation of precipitation_rate. (<span style="font-family: Courier New,Courier,monospace;">calc_precipitation</span>)<br>
-      <br>
-Bugfix: initial data assignments to some dvrp arrays changed due to error messages from gfortran compiler. (<span style="font-family: Courier New,Courier,monospace;">modules</span>)<br>
-      <br>
-Bugfix in calculation of cloud droplet velocity. (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>)<br>
-      <br>
-Bugfix in transfer of particles at south/left edge. (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>)<br>
-      <br>
-Bugfix in calculation of collision efficiency. (<span style="font-family: Courier New,Courier,monospace;">collision_efficiency</span>)<span style="font-family: Courier New,Courier,monospace;"></span><span style="font-family: Courier New,Courier,monospace;"></span><br>
-      </td>
-      <td style="vertical-align: top; width: 189px;">advec_particles,
-calc_precipitation, collision_efficiency, disturb_field, header,
-modules, poisfft_hybrid, subjob, sum_up_3d_data,
-time_integration<br>
-<br>
-      </td>
-    </tr>
-    <tr>
-      <td style="vertical-align: top; width: 27px;">23/09/09<br>
-      </td>
-      <td style="vertical-align: top; width: 38px;">SR<br>
-      </td>
-      <td style="vertical-align: top; width: 42px;">392<br>
-      </td>
-      <td style="vertical-align: top; width: 40px;">3.7<br>
-      </td>
-      <td style="vertical-align: top; width: 27px;">N<br>
-      </td>
-      <td style="vertical-align: top; width: 497px;"><span style="font-weight: bold;">Most
-important changes:</span>
-atmosphere-ocean coupling has been improved and steering of
-dvr-software is more user-friendly now. Beside that, the release
-contains a lot of optimizations, changes, and bugfixes which have
-accumulated over the last months.<br>
-
-      <br>
-
-The atmosphere-ocean coupling allows to use independent precursor runs
-in order to account for different spin-up times. The time when coupling
-has to be started is given by new inipar parameter <span style="font-family: Courier New,Courier,monospace;">coupling_start_time</span>.
-The precursor ocean run has to be started using new <span style="font-weight: bold;">mrun</span> option "<span style="font-family: Courier New,Courier,monospace;">-y</span>" in
-order to add the appendix "<span style="font-family: Courier New,Courier,monospace;">_O</span>" to all
-output files. (<span style="font-family: Courier New,Courier,monospace;">check_for_restart</span>,
-      <span style="font-family: Courier New,Courier,monospace;">check_parameters</span>,
-      <span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>,
-      <span style="font-family: Courier New,Courier,monospace;">data_output_3d</span>,
-      <span style="font-family: Courier New,Courier,monospace;">data_output_profiles</span>,
-      <span style="font-family: Courier New,Courier,monospace;">data_output_ptseries</span>,
-      <span style="font-family: Courier New,Courier,monospace;">data_output_spectra</span>,
-      <span style="font-family: Courier New,Courier,monospace;">data_output_tseries</span>,
-      <span style="font-family: Courier New,Courier,monospace;">header</span>,
-      <span style="font-family: Courier New,Courier,monospace;">init_coupling</span>,
-      <span style="font-family: Courier New,Courier,monospace;">modules</span>,
-      <span style="font-family: Courier New,Courier,monospace;">mrun</span>,
-      <span style="font-family: Courier New,Courier,monospace;">parin</span>,
-      <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>,
-      <span style="font-family: Courier New,Courier,monospace;">surface_coupler</span>,
-      <span style="font-family: Courier New,Courier,monospace;">time_integration</span>,
-      <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br>
-
-      <br>
-
-Clipping of dvrp output implemented. Default colourtable for particles
-implemented, particle attributes (color, dvrp_size) can be set with new
-parameters <span style="font-family: Courier New,Courier,monospace;">particle_color</span>,
-      <span style="font-family: Courier New,Courier,monospace;">particle_dvrpsize</span>,
-      <span style="font-family: Courier New,Courier,monospace;">color_interval</span>,
-      <span style="font-family: Courier New,Courier,monospace;">dvrpsize_interval</span>.
-Slicer attributes (dvrp) are set with new routine <span style="font-family: Courier New,Courier,monospace;">set_slicer_attributes_dvrp</span>
-and are controlled with existing parameters <span style="font-family: Courier New,Courier,monospace;">slicer_range_limits</span>.
-(<span style="font-family: Courier New,Courier,monospace;">init_dvrp</span>,
-      <span style="font-family: Courier New,Courier,monospace;">data_output_dvrp</span>,
-      <span style="font-family: Courier New,Courier,monospace;">modules</span>,
-      <span style="font-family: Courier New,Courier,monospace;">set_slicer_attributes_dvrp</span>,&nbsp;
-      <span style="font-family: Courier New,Courier,monospace;">user_data_output_dvrp</span>)<br>
-
-      <br>
-
-Polygon reduction for topography and ground plate isosurface. Reduction
-level for buildings can be chosen with parameter cluster_size. (<span style="font-family: Courier New,Courier,monospace;">init_dvrp</span>)<br>
-
-      <br>
-
-Variables <span style="font-family: Courier New,Courier,monospace;">bc_lr</span>
-/ <span style="font-family: Courier New,Courier,monospace;">bc_ns</span>&nbsp;
-in most subroutines replaced by LOGICAL variables <span style="font-family: Courier New,Courier,monospace;">bc_lr_cyc</span>,<br>
-
-      <span style="font-family: Courier New,Courier,monospace;">bc_ns_cyc</span>
-for speed optimization. This gives a significant speedup of more than
-10%. (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>,
-      <span style="font-family: Courier New,Courier,monospace;">diffusion_u</span>,
-      <span style="font-family: Courier New,Courier,monospace;">diffusion_v</span>,
-      <span style="font-family: Courier New,Courier,monospace;">diffusion_w</span>,
-      <span style="font-family: Courier New,Courier,monospace;">modules</span>)<br>
-
-      <br>
-
-Scripts have been adapted for machine lck (Yonsei Univ.). (<span style="font-family: Courier New,Courier,monospace;">mrun</span>, <span style="font-family: Courier New,Courier,monospace;">mbuild</span>, <span style="font-family: Courier New,Courier,monospace;">subjob</span>)<br>
-
-      <br>
-
-Additional timestep criterion in case of simulations with plant canopy.
-(<span style="font-family: Courier New,Courier,monospace;">timestep</span>)<br>
-
-      <br>
-
-Check for illegal entries in <span style="font-family: Courier New,Courier,monospace;">section_xy</span>|<span style="font-family: Courier New,Courier,monospace;">xz</span>|<span style="font-family: Courier New,Courier,monospace;">yz</span> that
-exceed <span style="font-family: Courier New,Courier,monospace;">nz+1</span>|<span style="font-family: Courier New,Courier,monospace;">ny+1</span>|<span style="font-family: Courier New,Courier,monospace;">nx+1</span>. (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>)<br>
-
-      <br>
-
-External pressure gradient can be used as driving force. (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>,
-      <span style="font-family: Courier New,Courier,monospace;">header</span>,
-      <span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>,
-      <span style="font-family: Courier New,Courier,monospace;">modules</span>,
-      <span style="font-family: Courier New,Courier,monospace;">parin</span>,
-      <span style="font-family: Courier New,Courier,monospace;">prognostic_equations</span>,
-      <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>,
-      <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br>
-
-      <br>
-
-New topography case '<span style="font-style: italic;">single_street_canyon</span>'.
-(<span style="font-family: Courier New,Courier,monospace;">header</span>,
-      <span style="font-family: Courier New,Courier,monospace;">init_grid</span>,
-      <span style="font-family: Courier New,Courier,monospace;">modules</span>,
-      <span style="font-family: Courier New,Courier,monospace;">parin</span>,
-      <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>,
-      <span style="font-family: Courier New,Courier,monospace;">user_check_parameters</span>,
-      <span style="font-family: Courier New,Courier,monospace;">user_header</span>,
-      <span style="font-family: Courier New,Courier,monospace;">user_init_grid</span>,
-      <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br>
-
-      <br>
-
-Option to predefine a target bulk velocity for <span style="font-family: Courier New,Courier,monospace;">conserve_volume_flow</span>.
-(<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>,
-      <span style="font-family: Courier New,Courier,monospace;">header</span>,
-      <span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>,
-      <span style="font-family: Courier New,Courier,monospace;">modules</span>,
-      <span style="font-family: Courier New,Courier,monospace;">parin</span>,
-      <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>,
-      <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br>
-
-      <br>
-
-Option for user defined 2D data output in xy cross sections at z=<span style="font-family: Courier New,Courier,monospace;">nzb</span>+1. (<span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>,
-      <span style="font-family: Courier New,Courier,monospace;">user_data_output_2d</span>)<br>
-
-      <br>
-
-xy cross section output of surface heatfluxes (sensible and latent). (<span style="font-family: Courier New,Courier,monospace;">average_3d_data</span>,
-      <span style="font-family: Courier New,Courier,monospace;">check_parameters</span>,
-      <span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>,
-      <span style="font-family: Courier New,Courier,monospace;">modules</span>,
-      <span style="font-family: Courier New,Courier,monospace;">read_3d_binary</span>,
-      <span style="font-family: Courier New,Courier,monospace;">sum_up_3d_data</span>,
-      <span style="font-family: Courier New,Courier,monospace;">write_3d_binary</span>)</td>
-      <td style="vertical-align: top; width: 189px;">average_3d_data,
-check_for_restart, check_parameters, data_output_2d, data_output_3d,
-data_output_dvrp, data_output_profiles, data_output_ptseries,
-data_output_spectra, data_output_tseries, diffusion_u, diffusion_v,
-diffusion_w, init_coupling, init_dvrp, init_grid, init_3d_model,
-header, mbuild, modules, mrun, package_parin, parin,
-prognostic_equations, read_3d_binary, read_var_list, subjob,
-sum_up_3d_data, surface_coupler, timestep, time_integration,
-user_check_parameters, user_data_output_2d, user_data_output_dvrp,
-user_header, user_init_grid, write_3d_binary, write_var_list<br>
-
-      <br>
-
-      <span style="font-weight: bold;">new:</span>
-set_particle_attributes, set_slicer_attributes_dvrp</td>
-    </tr>
-    <tr>
-      <td style="vertical-align: top; width: 27px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 38px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 42px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 40px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 27px;">C<br>
-      </td>
-      <td style="vertical-align: top; width: 497px;">Output of messages replaced by
-message handling routine. (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>,
-      <span style="font-family: Courier New,Courier,monospace;">advec_s_bc</span>,
-      <span style="font-family: Courier New,Courier,monospace;">buoyancy</span>,
-      <span style="font-family: Courier New,Courier,monospace;">calc_spectra</span>,
-      <span style="font-family: Courier New,Courier,monospace;">check_for_restart</span>,
-      <span style="font-family: Courier New,Courier,monospace;">check_open</span>,
-      <span style="font-family: Courier New,Courier,monospace;">coriolis</span>,
-      <span style="font-family: Courier New,Courier,monospace;">cpu_log</span>,
-      <span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>,
-      <span style="font-family: Courier New,Courier,monospace;">data_output_3d</span>,
-      <span style="font-family: Courier New,Courier,monospace;">data_output_dvrp</span>,
-      <span style="font-family: Courier New,Courier,monospace;">data_output_profiles</span>,
-      <span style="font-family: Courier New,Courier,monospace;">data_output_spectra</span>,
-      <span style="font-family: Courier New,Courier,monospace;">fft_xy</span>,
-      <span style="font-family: Courier New,Courier,monospace;">flow_statistics</span>,
-      <span style="font-family: Courier New,Courier,monospace;">header</span>,&nbsp;
-      <span style="font-family: Courier New,Courier,monospace;">init_1d_model</span>,
-      <span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>,
-      <span style="font-family: Courier New,Courier,monospace;">init_dvrp</span>,
-      <span style="font-family: Courier New,Courier,monospace;">init_grid</span>,
-      <span style="font-family: Courier New,Courier,monospace;">init_particles</span>,
-      <span style="font-family: Courier New,Courier,monospace;">init_pegrid</span>,
-      <span style="font-family: Courier New,Courier,monospace;">netcdf</span>,
-      <span style="font-family: Courier New,Courier,monospace;">parin</span>,
-      <span style="font-family: Courier New,Courier,monospace;">plant_canopy_model</span>,
-      <span style="font-family: Courier New,Courier,monospace;">poisfft_hybrid</span>,
-      <span style="font-family: Courier New,Courier,monospace;">poismg</span>,
-      <span style="font-family: Courier New,Courier,monospace;">read_3d_binary</span>,
-      <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>,
-      <span style="font-family: Courier New,Courier,monospace;">surface_coupler</span>,
-      <span style="font-family: Courier New,Courier,monospace;">temperton_fft</span>,
-      <span style="font-family: Courier New,Courier,monospace;">timestep</span>,
-      <span style="font-family: Courier New,Courier,monospace;">user_actions</span>,
-      <span style="font-family: Courier New,Courier,monospace;">user_data_output_dvrp</span>,
-      <span style="font-family: Courier New,Courier,monospace;">user_dvrp_coltab</span>,
-      <span style="font-family: Courier New,Courier,monospace;">user_init_grid</span>,
-      <span style="font-family: Courier New,Courier,monospace;">user_init_plant_canopy</span>,
-      <span style="font-family: Courier New,Courier,monospace;">user_parin</span>,
-      <span style="font-family: Courier New,Courier,monospace;">user_read_restart_data</span>,
-      <span style="font-family: Courier New,Courier,monospace;">user_spectra</span>)<br>
-
-      <br>
-
-Output of NetCDF messages with aid of message handling routine. (<span style="font-family: Courier New,Courier,monospace;">check_open</span>,
-      <span style="font-family: Courier New,Courier,monospace;">close_file</span>,
-      <span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>,
-      <span style="font-family: Courier New,Courier,monospace;">data_output_3d</span>,
-      <span style="font-family: Courier New,Courier,monospace;">data_output_profiles</span>,
-      <span style="font-family: Courier New,Courier,monospace;">data_output_ptseries</span>,
-      <span style="font-family: Courier New,Courier,monospace;">data_output_spectra</span>,&nbsp;
-      <span style="font-family: Courier New,Courier,monospace;">data_output_tseries</span>,
-      <span style="font-family: Courier New,Courier,monospace;">netcdf</span><span style="font-family: Courier New,Courier,monospace;"></span>)<br>
-
-      <br>
-
-__lcmuk changed to __lc to avoid problems with Intel compiler on
-sgi-ice. (<span style="font-family: Courier New,Courier,monospace;">poisfft</span>)<br>
-
-      <br>
-
-For extended NetCDF files, the updated title attribute includes an
-update of time_average_text where appropriate. (<span style="font-family: Courier New,Courier,monospace;">netcdf</span>)<br>
-
-      <br>
-
-In case of restart runs without extension, initial profiles are not
-written to NetCDF-file anymore. So far, time axis always started at
-t=0. (<span style="font-family: Courier New,Courier,monospace;">data_output_profiles</span>,
-      <span style="font-family: Courier New,Courier,monospace;">modules</span>,
-      <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>,
-      <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br>
-
-      <br>
-
-      <span style="font-family: Courier New,Courier,monospace;">initializing_actions</span>='<span style="font-style: italic;">read_data_for_recycling</span>' renamed to
-'<span style="font-style: italic;">cyclic_fill</span>'. A cyclic fill
-can now <br>
-
-be used independent of <span style="font-family: Courier New,Courier,monospace;">turbulent_inflow</span>.
-(<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>,
-      <span style="font-family: Courier New,Courier,monospace;">header</span>,
-      <span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>)<br>
-
-      <br>
-
-2 NetCDF error numbers changed. (<span style="font-family: Courier New,Courier,monospace;">data_output_3d</span>)<br>
-
-      <br>
-
-A link to the webspage appendix_a.html is printed for further
-information about possible reasons of the error that appeared. (message)<br>
-
-      <br>
-
-Temperature gradient criterion for estimating the boundary layer height
-replaced by the gradient criterion of Sullivan et al. (1998). (<span style="font-family: Courier New,Courier,monospace;">flow_statistics</span>)<br>
-
-      <br>
-
-NetCDF unit attribute in timeseries output in case of statistic regions
-added. (<span style="font-family: Courier New,Courier,monospace;">netcdf</span>)<br>
-
-      <br>
-
-Maximum number of tails is calculated from maximum number of particles
-and <span style="font-family: Courier New,Courier,monospace;">skip_particles_for_tail</span>.
-(<span style="font-family: Courier New,Courier,monospace;">init_particles</span>)<br>
-
-      <br>
-
-Value of <span style="font-family: Courier New,Courier,monospace;">vertical_particle_advection</span>
-is allowed to differ for each particle group. (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>,
-      <span style="font-family: Courier New,Courier,monospace;">header</span>,
-      <span style="font-family: Courier New,Courier,monospace;">modules</span>)
-      <br>
-
-      <br>
-
-First constant in array <span style="font-family: Courier New,Courier,monospace;">den</span> also
-defined as type double. (<span style="font-family: Courier New,Courier,monospace;">eqn_state_seawater</span>)<br>
-
-      <br>
-
-Parameter <span style="font-family: Courier New,Courier,monospace;">dvrp_psize</span>
-moved from <span style="font-family: Courier New,Courier,monospace;">particles_par</span>
-to <span style="font-family: Courier New,Courier,monospace;">dvrp_graphics_par</span>.
-(<span style="font-family: Courier New,Courier,monospace;">package_parin</span>)<br>
-
-      <br>
-
-      <span style="font-family: Courier New,Courier,monospace;">topography_grid_convention</span>
-moved from <span style="font-family: Courier New,Courier,monospace;">userpar</span>
-to <span style="font-family: Courier New,Courier,monospace;">inipar</span>
-(<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>,
-      <span style="font-family: Courier New,Courier,monospace;">header</span>,
-      <span style="font-family: Courier New,Courier,monospace;">parin</span>,
-      <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>,
-      <span style="font-family: Courier New,Courier,monospace;">user_check_parameters</span>,
-      <span style="font-family: Courier New,Courier,monospace;">user_header</span>,
-      <span style="font-family: Courier New,Courier,monospace;">user_init_grid</span>,
-      <span style="font-family: Courier New,Courier,monospace;">user_parin</span>,
-      <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br>
-
-      <br>
-
-Default value of <span style="font-family: Courier New,Courier,monospace;">grid_matching</span>
-changed to '<span style="font-style: italic;">strict</span>'. (<span style="font-family: Courier New,Courier,monospace;">modules</span>)<br>
-
-      <br>
-
-Adjustments for runs on lcxt4 (necessary due to a software update on
-CRAY) and for coupled runs on ibmy. (<span style="font-family: Courier New,Courier,monospace;">mrun</span>, <span style="font-family: Courier New,Courier,monospace;">subjob</span>)</td>
-      <td style="vertical-align: top; width: 189px;">advec_particles,
-advec_s_bc,
-buoyancy, calc_spectra, check_for_restart, check_open,
-check_parameters, close_file, coriolis, cpu_log, data_output_2d,
-data_output_3d, data_output_dvrp, data_output_profiles,
-data_output_ptseries, data_output_spectra, data_output_tseries,
-eqn_state_seawater, fft_xy, flow_statistics, header, init_1d_model,
-init_3d_model, init_dvrp, init_grid, init_particles, init_pegrid,
-message, mrun, netcdf, package_parin, parin, plant_canopy_model,
-poisfft, poisfft_hybrid, poismg, read_3d_binary, read_var_list, subjob,
-surface_coupler, temperton_fft, timestep,
-user_actions,user_check_parameters, user_data_output_dvrp,
-user_dvrp_coltab, user_header, user_init_grid, user_init_plant_canopy,
-user_parin, user_read_restart_data, user_spectra, write_var_list </td>
-    </tr>
-    <tr>
-      <td style="vertical-align: top; width: 27px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 38px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 42px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 40px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 27px;">B<br>
-      </td>
-      <td style="vertical-align: top; width: 497px;">Bugfix: Initial hydrostatic
-pressure profile in case of ocean runs is now calculated in 5 iteration
-steps. (<span style="font-family: Courier New,Courier,monospace;">init_ocean</span>)<br>
-
-      <br>
-
-Bugfix: wrong sign in buoyancy production of ocean part in case of not
-using the reference density (only in 3D routine <span style="font-family: Courier New,Courier,monospace;">production_e</span>).
-(<span style="font-family: Courier New,Courier,monospace;">production_e</span>)<br>
-
-      <br>
-
-Bugfix: output of averaged 2d/3d quantities requires that an avaraging
-interval has been set, respective error message is included. (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>)<br>
-
-      <br>
-
-Bugfix: Output on unit 14 only if requested by <span style="font-family: Courier New,Courier,monospace;">write_binary</span>.
-(<span style="font-family: Courier New,Courier,monospace;">user_last_actions</span>)<br>
-
-      <br>
-
-Bugfix to avoid zero division by <span style="font-family: Courier New,Courier,monospace;">km_neutral</span>.
-(<span style="font-family: Courier New,Courier,monospace;">production_e</span>)<br>
-
-      <br>
-
-Bugfix for extended NetCDF files: In order to avoid 'data mode' errors
-if updated attributes are larger than their original size, <span style="font-family: Courier New,Courier,monospace;">NF90_PUT_ATT</span>
-is called in 'define mode' enclosed by <span style="font-family: Courier New,Courier,monospace;">NF90_REDEF</span>
-and <span style="font-family: Courier New,Courier,monospace;">NF90_ENDDEF</span>
-calls. This implies a possible performance loss; an alternative
-strategy would be to ensure equal attribute size in a job chain. (<span style="font-family: Courier New,Courier,monospace;">netcdf</span>)<br>
-
-      <br>
-
-Bugfix: correction of initial volume flow for non-flat topography. (<span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>)<br>
-
-      <br>
-
-Bugfix: zero initialization of arrays within buildings for '<span style="font-style: italic;">cyclic_fill</span>'. (<span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>)<br>
-
-      <br>
-
-Bugfix: <span style="font-family: Courier New,Courier,monospace;">to_be_resorted</span>
-=&gt; <span style="font-family: Courier New,Courier,monospace;">s_av</span>
-for time-averaged scalars. (<span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>,
-      <span style="font-family: Courier New,Courier,monospace;">data_output_3d</span>)<br>
-
-      <br>
-
-Bugfix: avoid that <span style="font-family: Courier New,Courier,monospace;">ngp_2dh_s_inner</span>
-becomes zero. (<span style="font-family: Courier New,Courier,monospace;">init_3_model</span>)<br>
-
-      <br>
-
-Typographical error: unit of wpt in dots_unit. (<span style="font-family: Courier New,Courier,monospace;">modules</span>)<br>
-
-      <br>
-
-Bugfix: error in check, if particles moved further than one subdomain
-length. This check must not be applied for newly released particles. (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>)<br>
-
-      <br>
-
-Bugfix: several tail counters are initialized, <span style="font-family: Courier New,Courier,monospace;">particle_tail_coordinates</span>
-is only written to file if its third index is &gt; 0, arrays for tails
-are allocated with a minimum size of 10 tails if there is no tail
-initially. (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>,
-      <span style="font-family: Courier New,Courier,monospace;">init_particles</span><span style="font-family: Courier New,Courier,monospace;"></span>)<br>
-
-      <br>
-
-Bugfix: pressure included for profile output. (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>)<br>
-
-      <br>
-
-Bugfix: Type of <span style="font-family: Courier New,Courier,monospace;">count</span> and <span style="font-family: Courier New,Courier,monospace;">count_rate</span>
-changed to default <span style="font-family: Courier New,Courier,monospace;">INTEGER</span> on
-NEC machines. (<span style="font-family: Courier New,Courier,monospace;">cpu_log</span>)<br>
-
-      <br>
-
-Bugfix: output of particle time series only if particle advection is
-switched on. (<span style="font-family: Courier New,Courier,monospace;">time_integration</span>)<br>
-
-      <br>
-
-Bugfix: qsws was calculated in case of <span style="font-family: Courier New,Courier,monospace;">constant heatflux</span>
-= <span style="font-family: Courier New,Courier,monospace;">.FALSE.</span>.
-(<span style="font-family: Courier New,Courier,monospace;">prandtl_fluxes</span>)<br>
-
-      <br>
-
-Bugfix: averaging along z is not allowed for 2d quantities (e.g. u* and
-z0). (<span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>)<br>
-
-      <br>
-
-Typographical errors. (<span style="font-family: Courier New,Courier,monospace;">netcdf</span>)<br>
-
-      <br>
-
-If the inversion height calculated by the prerun is zero,
-inflow_damping_height must be explicitly specified. (<span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>)<br>
-
-      <br>
-
-Small bugfix concerning 3d 64bit netcdf output format. (<span style="font-family: Courier New,Courier,monospace;">header</span>)<br>
-
-      <br>
-
-Bugfix: <span style="font-family: Courier New,Courier,monospace;">dt_fixed</span>
-removed from the restart file, because otherwise, no change from a
-fixed to a variable timestep would be possible in restart runs. (<span style="font-family: Courier New,Courier,monospace;">read_var_list</span>,
-      <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br>
-
-      <br>
-
-Bugfix: initial setting of <span style="font-family: Courier New,Courier,monospace;">time_coupling</span>
-in coupled restart runs. (<span style="font-family: Courier New,Courier,monospace;">time_integration</span>)</td>
-      <td style="vertical-align: top; width: 189px;">advec_particles,
-check_parameters, cpu_log, data_output_2d, data_output_3d, header,
-init_3d_model, init_particles, init_ocean, modules, netcdf,
-prandtl_fluxes, production_e, read_var_list, time_integration,
-user_last_actions, write_var_list<br>
-
-      </td>
-    </tr>
-    <tr>
-      <td style="vertical-align: top; width: 27px;">02/02/09<br>
-      </td>
-      <td style="vertical-align: top; width: 38px;">SR<br>
-      </td>
-      <td style="vertical-align: top; width: 42px;">228<br>
-      </td>
-      <td style="vertical-align: top; width: 40px;">3.6<br>
-      </td>
-      <td style="vertical-align: top; width: 27px;">N<br>
-      </td>
-      <td style="vertical-align: top; width: 497px;">This
-release contains further additions, changes, and bugfixes for the
-SGI-ICE system, as well as important changes like the splitting of the
-user-interface into single files, a complete revision of the
-pre-compiling mechanism allowing to use different make-depositories
-simultaneously, and a revision of the reading-mechanism from restart
-files, which has no more restrictions in case that the previous and
-current run have different domain/subdomain sizes. A new handling of
-all kind of output messages is introduced with this release and will be
-applied to all existing messages within the next releases.<br>
-
-      <br>
-
-      <span style="font-weight: bold;">New features are:</span><br style="font-weight: bold;">
-
-      <br>
-
-Restart runs on SGI-ICE are working. (mrun)<br>
-
-2d-decomposition is default on SGI-ICE systems. (init_pegrid)<br>
-
-      <br>
-
-Ocean-atmosphere coupling realized with MPI-1. mrun adjusted for this
-case (-Y option). Check that PALM is called with <span style="font-family: Courier New,Courier,monospace;">mrun</span>-option
-"<span style="font-style: italic;">-K parallel</span>" for coupling.
-Adjustments in <span style="font-family: Courier New,Courier,monospace;">mrun</span>, <span style="font-family: Courier New,Courier,monospace;">mbuild</span>, and
-      <span style="font-family: Courier New,Courier,monospace;">subjob</span>
-for lcxt4.<br>
-
-      <br>
-
-DVRP arguments changed to single precision, mode pathlines added.<br>
-
-      <br>
-
-User can add additional routines in files <span style="font-family: Courier New,Courier,monospace;">user_additional_routines</span>.<br>
-
-      <br>
-
-User can check user parameters and deduce further quantities in <span style="font-family: Courier New,Courier,monospace;">user_check_parameters</span></td>
-      <td style="vertical-align: top; width: 189px;">check_for_restart,
-check_parameters, data_output_dvrp, init_dvrp, init_pegrid, local_stop,
-modules, package_parin, palm, surface_coupler, timestep<br>
-
-      <br>
-
-Makefile,&nbsp;mbuild, mrun, subjob<br>
-
-      <br>
-
-      <span style="font-weight: bold;">new:</span><br>
-
-init_coupling, user_additional_routines, user_check_parameters</td>
-    </tr>
-    <tr>
-      <td style="vertical-align: top; width: 27px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 38px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 42px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 40px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 27px;">C<br>
-      </td>
-      <td style="vertical-align: top; width: 497px;">User interface has been
-split into single files.<br>
-
-      <br>
-
-Type of variables <span style="font-family: Courier New,Courier,monospace;">count</span> and <span style="font-family: Courier New,Courier,monospace;">count_rate</span>
-changed to INTEGER(8) in order to avoid "out of range" problems, which
-result in measured negative time intervals. (cpu_log, local_tremain,
-local_tremain_ini)<br>
-
-      <br>
-
-Reading mechanism from restart files
-completely revised. The subdomain/total domain size is now allowed to
-vary arbitrarily between the current and previous run. (read_var_list,
-read_3d_binary, user_read_restart_data)<br>
-
-      <br>
-
-Precompilation mechanism (<span style="font-family: Courier New,Courier,monospace;">mbuild</span>)
-completely revised: the source-tarfile remains in the source directory (<span style="font-family: Courier New,Courier,monospace;">trunk/SOURCE</span>).
-One depository per block (given in the config-file) is created. Always
-all files from the source directory are copied to the respective
-depository. No additional file checks are done any more (mrun version
-2.1). The depository name used contains the conditions given by <span style="font-family: Courier New,Courier,monospace;">mrun</span>-option
-"<span style="font-family: Courier New,Courier,monospace;">-K</span>".
-(mbuild, mrun)<br>
-
-      <br>
-
-Output messages (including required stop of execution) can now be
-handled with the new subroutine&nbsp;<span style="font-family: Courier New,Courier,monospace;">message</span>.
-All output messages will be replaced by this routine step by step
-within the next revisions.<br>
-
-      <br>
-
-Neumann boundary condition at <span style="font-family: Courier New,Courier,monospace;">k=nzb</span> is
-explicitly set for better reading, although this has been already done
-in boundary_conds. (advec_s_bc)<br>
-
-      <br>
-
-Origin of the xy-coordinate system shifted from the center of the first
-grid cell (indices <span style="font-family: Courier New,Courier,monospace;">i=0</span>, <span style="font-family: Courier New,Courier,monospace;">j=0</span>) to the
-south-left corner of this cell. (netcdf) Topography definition
-according to new user parameter <span style="font-family: Courier New,Courier,monospace;">topography_grid_convention</span>.
-(init_grid, modules, user_header, user_init_grid, user_parin)</td>
-      <td style="vertical-align: top; width: 189px;">advec_s_bc,
-check_parameters, cpu_log, init_grid, local_stop, local_tremain,
-local_tremain_ini, modules, netcdf, read_3d_binary, read_var_list,
-user_read_restart_data<br>
-
-      <br>
-
-Makefile, mbuild, mrun<br>
-
-      <br>
-
-      <span style="font-weight: bold;">removed:</span><br>
-
-user_interface<br>
-
-      <br>
-
-      <span style="font-weight: bold;">new:</span><br>
-
-message,
-user_3d_data_averaging, user_actions, user_advec_particles,
-user_check_data_output, user_check_data_output_pr, user_data_output_2d,
-user_data_output_3d, user_data_output_dvrp, user_define_netcdf_grid,
-user_dvrp_coltab, user_header, user_init_3d_model, user_init,
-user_init_grid, user_init_particles, user_init_plant_canopy,
-user_last_actions, user_parin, user_parin, user_particle_attributes,
-user_read_restart_data, user_spectra, user_statistics</td>
-    </tr>
-    <tr>
-      <td style="vertical-align: top; width: 27px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 38px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 42px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 40px;"><br>
-      </td>
-      <td style="vertical-align: top; width: 27px;">B<br>
-      </td>
-      <td style="vertical-align: top; width: 497px;">Bugfix: array <span style="font-family: Courier New,Courier,monospace;">d</span> is
-reallocated in case that multigrid is used. (calc_spectra)<br>
-
-      <br>
-
-Bugfixes for nonparallel execution. (check_for_restart, cpu_statistics,
-inflow_turbulence,&nbsp;timestep)<br>
-
-      <br>
-
-Size of <span style="font-family: Courier New,Courier,monospace;">pf3d
-      </span>adjusted to the required output size (1 gridpoint less,
-along all three dimensions), because output of a subset of the data (<span style="font-family: Courier New,Courier,monospace;">pf3d(nxa:nxe...)</span>)
-in the NF90_PUT_VAR statement caused segmentation fault with the INTEL
-compiler. (combine_plot_fields)<br>
-
-      <br>
-
-Bugfix: error in <span style="font-family: Courier New,Courier,monospace;">zu</span> index in
-case of <span style="font-family: Courier New,Courier,monospace;">section_xy=-1</span>.
-(header)<br>
-
-      <br>
-
-Inconsistency
-removed: as the thermal stratification is not taken into account for
-the evaluation of the wall fluxes at vertical walls, the eddy viscosity
-      <span style="font-family: Courier New,Courier,monospace;">km</span>
-must not be used for the evaluation of the velocity gradients <span style="font-family: Courier New,Courier,monospace;">dudy</span>, <span style="font-family: Courier New,Courier,monospace;">dwdy</span>, <span style="font-family: Courier New,Courier,monospace;">dvdx</span> and <span style="font-family: Courier New,Courier,monospace;">dwdx</span>.
-(production_e)<br>
-
-      <br>
-
-Bugfix in calculating <span style="font-family: Courier New,Courier,monospace;">k</span> index in
-case of oceans runs. (sort_particles)<br>
-
-      <br>
-
-Bugfix: no output of particle concentration and radius unless particles
-have been started. (data_output_2d)<br>
-
-      <br>
-
-Bugfix: reading of <span style="font-family: Courier New,Courier,monospace;">spectrum_x|y</span>
-from restart files ignored if total numbers of grid points do not
-match. (read_3d_binary)<br>
-
-      <br>
-
-Bugfix: abort in case that absolute temperature is below zero.
-(init_cloud_physics)</td>
-      <td style="vertical-align: top; width: 189px;">advec_particles,
-calc_spectra, check_for_restart, cpu_statistics, data_output_2d,
-header, inflow_turbulence, init_cloud_physics, production_e,
-read_3d_binary, timestep<br>
-
-      <br>
-
-combine_plot_fields</td>
-    </tr>
-  </tbody>
-</table>
-<b><blink></blink></b><br>&nbsp;
-<br>&nbsp;
-<h2><a name="Kapitel2.0"></a>2.0&nbsp; How
-to change the source
-code and how to document modifications<br>
-</h2>Currently
-(Sep 23rd, 2009), only selected users are allowed to commit changes to
-the repository. Rules given below are preliminary, until we have switched to the "trac"-System.<br>
-<br>
-<ol>
-  <li>Change source code in your current working copy. Document changes in the respective file headers unter "<span style="font-family: Courier New,Courier,monospace;">Current revisions:</span>".
-Additionally, document all changes in a local file (e.g.
-~/palm/current_version/CURRENT_MODIFICATIONS), including the files that
-have been changed. Classify changes in three cases, new (N), changed
-(C), and bugfix (B). <span style="font-weight: bold;">Try to be as short and as precise as possible in describing the changes!</span><br>
-    <br>
-  </li>
-  <li>Test the changed version, unless it works sufficiently.<br>
-    <br>
-  </li>
-  
-  <li>Commit the changes to the repository<br>
-    <br>
-    <span style="font-family: Courier New,Courier,monospace;">&nbsp;&nbsp; svn commit -m "message string" trunk<br>
-    </span><br>
-where "<span style="font-family: Courier New,Courier,monospace;">message string</span>" is a short, <span style="font-weight: bold;">meaningful</span> summary of the changes.<br>
-    <br>
-Alternatively, in case that you have made quite a lot of changes, you
-can also give the contents of your file CURRENT_MODIFICATIONS as the
-log message:<br>
-    <br>
-    <span style="font-family: Courier New,Courier,monospace;">&nbsp;&nbsp; svn commit -F CURRENT_MODIFICATIONS trunk</span><br>
-    <br>
-  </li>
-
-  <li>Enter the changes into the technical documentation (under <span style="font-family: Courier New,Courier,monospace;">.../trunk/DOC/tec/technical_documentation.html</span>).<br>
-    <br>
-  </li>
-  <li>Move the change comments in the file headers from "Current
-revisions:" to "Former revisions:". Set the current svn id in front of
-that.<br>
-    <br>
-    <span style="font-weight: bold;">Example:</span><br>
-    <br>
-    <span style="font-weight: bold;">before move:</span><br>
-    <br>
-    <span style="font-family: Courier New,Courier,monospace;">! Current revisions:</span><br style="font-family: Courier New,Courier,monospace;">
-    <span style="font-family: Courier New,Courier,monospace;">! -----------------</span><br style="font-family: Courier New,Courier,monospace;">
-    <span style="font-family: Courier New,Courier,monospace;">! <span style="font-weight: bold;">new dummy argument abcd</span></span><span style="font-family: Courier New,Courier,monospace;"></span><br style="font-family: Courier New,Courier,monospace;">
-    <span style="font-family: Courier New,Courier,monospace;">!</span><br style="font-family: Courier New,Courier,monospace;">
-    <span style="font-family: Courier New,Courier,monospace;">! Former revisions:</span><br style="font-family: Courier New,Courier,monospace;">
-    <span style="font-family: Courier New,Courier,monospace;">! -----------------</span><br style="font-family: Courier New,Courier,monospace;">
-    <span style="font-family: Courier New,Courier,monospace;">! ï¿œId: init_3d_model.f90 <span style="color: rgb(51, 102, 255); font-weight: bold;">425 2010-02-06 13:37:48Z raasch</span> ï¿œ</span><br style="font-family: Courier New,Courier,monospace;">
-    <span style="font-family: Courier New,Courier,monospace;">!</span><br style="font-family: Courier New,Courier,monospace;">
-    <span style="font-family: Courier New,Courier,monospace;"></span><span style="font-family: Courier New,Courier,monospace;"></span><span style="font-family: Courier New,Courier,monospace;"></span><span style="font-family: Courier New,Courier,monospace;">! 407 2009-12-01 15:01:15Z maronga</span><br style="font-family: Courier New,Courier,monospace;">
-    <span style="font-family: Courier New,Courier,monospace;">! var_ts is replaced by dots_max </span><br>
-    <br>
-    <span style="font-weight: bold;">after move:</span><br>
-    <br>
-    <span style="font-family: Courier New,Courier,monospace;">! Current revisions:</span><br style="font-family: Courier New,Courier,monospace;">
-
-    <span style="font-family: Courier New,Courier,monospace;">! -----------------</span><br style="font-family: Courier New,Courier,monospace;">
-
-    <span style="font-family: Courier New,Courier,monospace;">! </span><br style="font-family: Courier New,Courier,monospace;">
-
-    <span style="font-family: Courier New,Courier,monospace;">!</span><br style="font-family: Courier New,Courier,monospace;">
-
-    <span style="font-family: Courier New,Courier,monospace;">! Former revisions:</span><br style="font-family: Courier New,Courier,monospace;">
-
-    <span style="font-family: Courier New,Courier,monospace;">! -----------------</span><br style="font-family: Courier New,Courier,monospace;">
-
-    <span style="font-family: Courier New,Courier,monospace;">! ï¿œId: </span><span style="font-family: Courier New,Courier,monospace;">init_3d_model.f90 425 2010-02-06 13:37:48Z raasch</span><span style="font-family: Courier New,Courier,monospace;"> ï¿œ<br>
-!</span><br style="font-family: Courier New,Courier,monospace;">
-    <span style="font-family: Courier New,Courier,monospace;">! </span><span style="font-family: Courier New,Courier,monospace; color: rgb(51, 102, 255); font-weight: bold;">425 2010-02-06 13:37:48Z raasch</span><span style="font-family: Courier New,Courier,monospace;"><span style="font-weight: bold; color: rgb(51, 51, 255);"></span><br>
-! </span><span style="font-family: Courier New,Courier,monospace; font-weight: bold;">new dummy argument abcd</span><br style="font-family: Courier New,Courier,monospace;">
-
-    <span style="font-family: Courier New,Courier,monospace;">!</span><br style="font-family: Courier New,Courier,monospace;">
-
-    <span style="font-family: Courier New,Courier,monospace;">! 407 2009-12-01 15:01:15Z maronga</span><br style="font-family: Courier New,Courier,monospace;">
-
-    <span style="font-family: Courier New,Courier,monospace;">! var_ts is replaced by dots_max<br>
-    <br>
-    <span style="font-family: Times New Roman,Times,serif;">Note that the "<span style="font-family: Courier New,Courier,monospace;">ï¿œ</span>"
-in this example should read "$", but "$" cannot be used here, because
-otherwise svn would replace these entries after every change to this
-technical documentation.</span><br>
- </span><br>
-  </li>
-  <li>Commit the changes in the technical documentation and the file headers:<br>
-    <br>
-    <span style="font-family: Courier New,Courier,monospace;">&nbsp;&nbsp; svn commit -m "last commit documented" trunk</span><br>
-    <br>
- <br>
-  </li>
-</ol>
-
-<h2>
-<a name="Kapitel3.0"></a>3.0&nbsp; Description
-of selected parts of the
-model source code</h2>
-<b></b>
-<p>This section gives links to existing descriptions about special
+<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
+<HTML>
+<HEAD>
+	<META HTTP-EQUIV="CONTENT-TYPE" CONTENT="text/html; charset=utf-8">
+	<TITLE>PALM technical/numerical documentation</TITLE>
+	<META NAME="GENERATOR" CONTENT="OpenOffice.org 3.0  (Unix)">
+	<META NAME="AUTHOR" CONTENT="Siegfried Raasch">
+	<META NAME="CREATED" CONTENT="0;0">
+	<META NAME="CHANGED" CONTENT="20100901;15562900">
+	<META NAME="KEYWORDS" CONTENT="parallel LES model">
+</HEAD>
+<BODY LANG="en-US" DIR="LTR">
+<P STYLE="margin-bottom: 0in"><FONT SIZE=7>PALM&nbsp;&nbsp; </FONT>a
+<B>pa</B>rallelized <B>L</B>ES <B>m</B>odel 
+</P>
+<HR>
+<P STYLE="margin-bottom: 0in"><BR><FONT SIZE=5>technical/numerical
+documentation</FONT> <BR>&nbsp; <BR>&nbsp; 
+</P>
+<P ALIGN=CENTER STYLE="margin-bottom: 0in">&nbsp; 
+</P>
+<TABLE CELLPADDING=0 CELLSPACING=0>
+	<TR>
+		<TD>
+			<P><B>Current fixed release:</B></P>
+		</TD>
+		<TD>
+			<P><B>3.7a</B></P>
+		</TD>
+	</TR>
+	<TR>
+		<TD>
+			<P><B>Last change of this document</B>:&nbsp;</P>
+		</TD>
+		<TD>
+			<P><B>$Id: technical_documentation.html 544 2010-08-17 10:44:58Z
+			maronga $</B></P>
+		</TD>
+	</TR>
+</TABLE>
+<HR>
+<P><FONT SIZE=5>Table of contents</FONT> <BR>&nbsp; <BR>&nbsp; 
+</P>
+<UL>
+	<P><A HREF="#Kapitel1.0">1.0</A>&nbsp; History of changes 
+	</P>
+	<P><A HREF="#Kapitel2.0">2.0</A>&nbsp; How to change the source code
+	and how to document modifications</P>
+	<P><A HREF="#Kapitel3.0">3.0</A>&nbsp; Description of selected model
+	numerics<BR>&nbsp;</P>
+</UL>
+<P><BR><BR>
+</P>
+<H2><A NAME="Kapitel1.0"></A>1.0&nbsp; History of changes</H2>
+<P STYLE="margin-bottom: 0in">The following table documents all
+changes to the PALM repository since 2009. Older changes can be found
+<A HREF="tec_doc_1998-2008.html">here</A>.<BR><BR>Entries in the
+fifth column specify the type of the change: change&nbsp; of existing
+code (C), new code (N) or bugfix (B). More detailed information about
+the changes may sometimes be found in the header comment lines of the
+respective routines. <BR>&nbsp; 
+</P>
+<TABLE WIDTH=1223 BORDER=1 CELLPADDING=2 CELLSPACING=3>
+	<COL WIDTH=37>
+	<COL WIDTH=56>
+	<COL WIDTH=68>
+	<COL WIDTH=62>
+	<COL WIDTH=36>
+	<COL WIDTH=518>
+	<COL WIDTH=393>
+	<TR>
+		<TD WIDTH=37>
+			<P><FONT SIZE=4><B>date</B></FONT></P>
+		</TD>
+		<TD WIDTH=56>
+			<P><FONT SIZE=4><B>author</B></FONT></P>
+		</TD>
+		<TD WIDTH=68>
+			<P><FONT SIZE=4><B>svn revision</B></FONT></P>
+		</TD>
+		<TD WIDTH=62>
+			<P><FONT SIZE=4><B>last release</B></FONT></P>
+		</TD>
+		<TD WIDTH=36>
+			<P><FONT SIZE=4><B>type</B></FONT></P>
+		</TD>
+		<TD WIDTH=518>
+			<P><FONT SIZE=4><B>description</B></FONT></P>
+		</TD>
+		<TD WIDTH=393>
+			<P><FONT SIZE=4><B>affected routines</B></FONT></P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37 SDVAL="40187" SDNUM="1033;0;MM/DD/YY">
+			<P STYLE="font-weight: normal"><FONT SIZE=3>01/09/10</FONT></P>
+		</TD>
+		<TD WIDTH=56>
+			<P STYLE="font-weight: normal"><FONT SIZE=3>CH</FONT></P>
+		</TD>
+		<TD WIDTH=68 SDVAL="553" SDNUM="1033;">
+			<P STYLE="font-weight: normal"><FONT SIZE=3>553</FONT></P>
+		</TD>
+		<TD WIDTH=62>
+			<P STYLE="font-weight: normal"><FONT SIZE=3>3.7a</FONT></P>
+		</TD>
+		<TD WIDTH=36>
+			<P STYLE="font-weight: normal"><FONT SIZE=3>C</FONT></P>
+		</TD>
+		<TD WIDTH=518>
+			<P STYLE="font-weight: normal"><FONT SIZE=3>d3dpar-parameters for
+			masked output <FONT FACE="Courier New, monospace">data_output_mask_*,
+			data_output_mask_user_*, mask_*_x, mask_*_y, mask_*_z,
+			mask_*_x_loop, mask_*_y_loop </FONT><FONT FACE="Times New Roman, serif">and
+			</FONT><FONT FACE="Courier New, monospace">mask_*_z_loop </FONT>have
+			been replaced by corresponding arrays <FONT FACE="Courier New, monospace">data_output_masks,
+			data_output_masks_user, mask_x, mask_y, mask_z, mask_x_loop,
+			mask_y_loop </FONT><FONT FACE="Times New Roman, serif">and</FONT><FONT FACE="Courier New, monospace">
+			mask_z_loop</FONT></FONT></P>
+			<P STYLE="font-weight: normal"><FONT SIZE=3>inipar-parameter <FONT FACE="Courier New, monospace">masks</FONT>
+			has been removed. </FONT>
+			</P>
+			<P STYLE="font-weight: normal"><FONT SIZE=3>The changes concern
+			the masked data output, which is still not documented.</FONT></P>
+		</TD>
+		<TD WIDTH=393>
+			<P STYLE="font-weight: normal"><FONT SIZE=3>check_parameters,
+			init_masks, modules, parin, user_parin</FONT></P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>17/08/10</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>BM</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>544</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.7a</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>B/C/N</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>adjustments for lcxt4; bugfix (mbuild): loading of user-defined
+			modules; new .mrun.config.lcxt4 for Cray XT4 in Bergen/Norway</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>mrun, mbuild, subjob</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>15/06/10</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>RH</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>540</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.7a</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>C</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>ssh from compute nodes to login nodes on hosts lcsgih/lcsgib is
+			done using the names hicegate0/bicegate0 instead of the IP
+			addresses</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>mrun</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>21/04/10</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>RH</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>532</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.7a</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>C</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>Large scale vertical motion (subsidence/ascent) can be added as
+			additional source to the prognostic equation for humidity/passive
+			scalar<BR>&nbsp;<BR>character length for dopr_unit
+			enlarged<BR><BR>Adjustment of the NCL scripts to allow for output
+			of png files 
+			</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>modules, prognostic_equations<BR><BR>.ncl.config.default,
+			cross_sections.ncl, profiles.ncl, timeseries.ncl, spectra.ncl</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=56>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=68>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=62>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>B</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>Bugfix: unit of hyp changed to dbar</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>check_parameters</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>01/04/10</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>RH</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>530</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.7a</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>B</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>Bugfix: treatment of special character * in NCL variables
+			containing lists corrected</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>palmplot</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>31/03/10</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>RH</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>529</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.7a</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>C</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>List of variable names have to specified using blank separated
+			lists and single quotes</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>palmplot</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>30/03/10</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>RH</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>526</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.7a</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>C</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>Adjustment of the NCL scripts and <FONT FACE="Courier New, Courier, monospace">palmplot
+			</FONT>to allow for special characters in NetCDF variable names 
+			</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>palmplot, .ncl.config.default, cross_sections.ncl,
+			profiles.ncl, timeseries.ncl, spectra.ncl</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>30/03/10</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>SR</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>525</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.7a</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>B</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>Bugfix: &quot;/&quot; in netcdf profile variable names replaced
+			by &quot;:&quot;</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>check_parameters</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>19/03/10</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>SR</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>520</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.7a</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>N</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>NetCDF4 support for particle data. Total number of particles is
+			allowed to increase during the run, when netcdf_data_format &gt; <I>2</I>
+			is used.</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>advec_particles, check_open, netcdf</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=56>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=68>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=62>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>C</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>Special characters like * and &quot; are now allowed for NetCDF
+			variable names. Replacement of these characters has been removed
+			(routine clean_netcdf_varname removed)</P>
+		</TD>
+		<TD WIDTH=393>
+			<P><BR>
+			</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>18/03/10</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>SR</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>516</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.7a</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>N</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>Simple installation method has been documented (see chapter
+			5.0). Respective scripts and configuration files have been
+			updated.</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>palm_simple_install, palm_simple_run</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=56>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=68>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=62>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>B</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>opt=1 is forced for PGI-compiler in routine <FONT FACE="Courier New, Courier, monospace">buoyancy
+			</FONT>by compiler-directive. Otherwise, program aborts with
+			SIGFPE.</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>buoyancy</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>16/03/10</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>RH</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>513</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.7a</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>N</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>NCL scripts can be run with the shell script <FONT FACE="Courier New, Courier, monospace">palmplot.</FONT></P>
+		</TD>
+		<TD WIDTH=393>
+			<P>palmplot, .ncl.config.default, cross_sections.ncl, profiles.ncl</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>12/03/10</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>SR</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>509</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.7a</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>B</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>bugfix for rev 504: remode_md was given a wrong value</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>mbuild</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>08/03/10</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>SR</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>506</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.7a</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>C</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>hack for better usage of <B>mrun</B> on Kyoto computing center
+			machine (<FONT FACE="Courier New, Courier, monospace">lckyoto</FONT>)</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>mrun</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>03/03/10</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>SR</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>504</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.7a</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>C</P>
+		</TD>
+		<TD WIDTH=518>
+			<P><FONT FACE="Courier New, Courier, monospace">combine_plot_fields</FONT>
+			is now directly called within <B>mrun</B> and must not be given by
+			an output command in the configuration file any more. Please
+			remove the respective output command (<FONT FACE="Courier New, Courier, monospace">OC:combine_plot_fields.x</FONT>)
+			from your configuration file!<BR><BR>Utility programs
+			(<FONT FACE="Courier New, Courier, monospace">combine_plot_fields</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">interpret_config</FONT>)
+			are created by <B>mbuild</B> separately for each block in the
+			configuration file. Depending on option <FONT FACE="Courier New, Courier, monospace">-K</FONT>,
+			<B>mrun</B> is using the utility programs for the respective
+			block.<BR><BR>Utility programs are also compiled in parallel,
+			depending on variable <FONT FACE="Courier New, Courier, monospace">%mopts</FONT>
+			in the configuration file..<BR><BR>Default configuration files
+			updated.</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>mbuild, mrun, UTIL/Makefile, SCRIPTS/.mrun.config.*</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>02/03/10</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>SR</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>496</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.7a</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>C</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>Particle data for restart runs are only written if environment
+			variable <FONT FACE="Courier New, Courier, monospace">write_binary</FONT>
+			= <I>.T.</I>, i.e. particle data are now handled in the same way
+			as the other restart data.</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>palm</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>01/03/10</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>SR</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>494</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.7a</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>N</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>Output in NetCDF4-format.<BR><BR>NetCDF4 allows for parallel
+			I/O. Also, variables are not limited in size any more (i.e. arrays
+			can be larger than 2GB). The new d3par-parameter
+			<FONT FACE="Courier New, Courier, monospace">netcdf_data_format</FONT>.can
+			be used to select among the different NetCDF file formats (<I>1</I>:
+			classic format, <I>2</I>: 64-bit-offset format (this is now the
+			default), <I>3</I>: NetCDF4/HDF5 format, <I>4</I>: NetCDF4 format
+			with NF90_CLASSIC_MODEL bit set). Setting <FONT FACE="Courier New, Courier, monospace">netcdf_data_format</FONT>
+			&gt; <I>2</I> only affects 2d- and 3d-output, i.e. there is no
+			parallel output for profiles, time series, spectra, and particles.
+			Parallel output of particle data may realized in the future.<BR><BR>In
+			case of parallel output, all PEs write into one file, which means
+			that postprocessing the data with <FONT FACE="Courier New, Courier, monospace">combine_plot_fields</FONT>
+			(merging data from single files - one per PE - into one file) is
+			not neccessary. In this case, flag files are automatically set,
+			telling <FONT FACE="Courier New, Courier, monospace">combine_plot_fields</FONT>
+			that nothing is to be done.<BR><BR>Using NetCDF4 format requires a
+			NetCDF4/HDF library. Respective settings for -I, -L, and -l
+			compiler/loader options for the SGI-ICE at HLRN can be found in
+			the new example configuration file
+			<FONT FACE="Courier New, Courier, monospace">..../trunk/SCRIPTS/.mrun.config.imuk_ice2_netcdf4</FONT><BR><BR><B>Attention:</B><BR>So
+			far, NetCDF parallel I/O has only been tested a little. Further
+			performance optimizations may be required, especially in case of
+			large data sets (e.g. data chunking, or lustre file system
+			settings with <FONT FACE="Courier New, Courier, monospace">lfs
+			setstripe</FONT>). <B>Please report bugs and performance problems
+			immediately!</B><BR><BR><FONT FACE="Courier New, Courier, monospace">combine_plot_fields</FONT>
+			outputs cpu-time needed for file merging.<BR><BR>Modules to be
+			loaded for compilation (mbuild) or job execution (mrun) can/should
+			now be given in the configuration file (<FONT FACE="Courier New, Courier, monospace">.mrun.config</FONT>)
+			using variable <FONT FACE="Courier New, Courier, monospace">modules</FONT>.
+			Example:<BR><BR>&nbsp;&nbsp; <FONT FACE="Courier New, Courier, monospace">%modules&nbsp;&nbsp;&nbsp;&nbsp;
+			ifort/11.0.069:netcdf&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; lcsgih
+			parallel</FONT><BR><BR>This method replaces the (undocumented)
+			<FONT FACE="Courier New, Courier, monospace">mpilib</FONT>-variable.<BR><BR><B>WARNING:</B><BR><B>All
+			fixed settings of modules in the scripts mbuild, mrun, and subjob
+			have been removed! Please set the </B><FONT FACE="Courier New, Courier, monospace"><B>modules</B></FONT><B>
+			variable appropriately and remove the mpilib variable in your
+			configuration file.<BR><BR></B>Compilation time can be reduced be
+			using parallel make. The number of allowed simultaneous compiler
+			calls (threads) is set in the configuration file using variable
+			<FONT FACE="Courier New, Courier, monospace">mopt</FONT>, e.g.<BR><BR>&nbsp;&nbsp;
+			<FONT FACE="Courier New, Courier, monospace">%mopts&nbsp; -j:4 &nbsp;
+			lcsgih parallel</FONT><BR><BR>allows a maximum of 4 routines to be
+			compiled simultaneously (provided that the makefile does not show
+			any dependencies between these routines).</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>check_open, check_parameters, close_file, data_output_2d,
+			data_output_3d, header, modules, netcdf,
+			parin<BR><BR>combine_plot_fields<BR><BR>mbuild, mrun, subjob</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=56>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=68>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=62>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>C</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>d3par-parameters <FONT FACE="Courier New, Courier, monospace">netcdf_64bit</FONT>
+			and <FONT FACE="Courier New, Courier, monospace">netcdf_64bit_3d</FONT>
+			have been removed. Use the new d3par-variable <FONT FACE="Courier New, Courier, monospace">netcdf_data_format</FONT>
+			= <I>2</I> for choosing the NetCDF 64bit-offset format (this is
+			the default from now on). The offset-format can not be set
+			independently for the 3d-output-data any more.<BR><BR>d3par-parameters
+			<FONT FACE="Courier New, Courier, monospace">netcdf_format_mask</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">netcdf_format_mask_av</FONT>,
+			and variables <FONT FACE="Courier New, Courier, monospace">nc_format_mask</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">format_parallel_io</FONT>
+			removed. They are replaced by the new parameter
+			<FONT FACE="Courier New, Courier, monospace">netcdf_data_format</FONT>.
+			This change concerns the masked data output, which is still not
+			documented.</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>check_open, close_file,<BR>data_output_mask, header,
+			init_masks, modules, parin</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=56>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=68>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=62>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>B</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>Bugfix in trunk/UTIL/Makefile: forgot to compile for
+			interpret_config.<BR><BR>Bugfix: timeseries data have to be
+			collected by PE0 (user_statistics).</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>user_statistics<BR><BR>trunk/UTIL/Makefile</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>05/02/10</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>SR</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>485</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.7a</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>B</P>
+		</TD>
+		<TD WIDTH=518>
+			<P><FONT FACE="Courier New, Courier, monospace">ngp_3d</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">ngp_3d_inner</FONT>
+			changed to 64 bit. Otherwise, on large grids &gt; 2*10**9 grid
+			points some global statistics (e.g. total kinetic energy) have
+			wrong values.</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>init_3d_model, modules</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>05/02/10</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>SR</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>484</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.7a</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>B</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>Typo in file headers removed (&quot;current&quot; instead of
+			&quot;actual&quot;).</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>many</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>05/02/10</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>SR</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>483</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.7a</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>N</P>
+		</TD>
+		<TD WIDTH=518>
+			<P><FONT FACE="Courier New, Courier, monospace">compare_palm_logs
+			</FONT>is additionally compiled with <B>mbuild</B> <FONT FACE="Courier New, Courier, monospace">-u</FONT>
+			. (<FONT FACE="Courier New, Courier, monospace">Makefile</FONT> in
+			<FONT FACE="Courier New, Courier, monospace">trunk/UTIL</FONT>)<BR><BR>make
+			options (<FONT FACE="Courier New, Courier, monospace">mopts</FONT>)
+			to be set by configuration file implemented. (<B>mrun</B>,
+			<B>mbuild</B>)<BR><BR><FONT FACE="Courier New, Courier, monospace">humidity</FONT>
+			= <I>.T.</I> can be used for runs with topography.
+			<FONT FACE="Courier New, Courier, monospace">wall_humidityflux</FONT>
+			and<BR><FONT FACE="Courier New, Courier, monospace">wall_scalarflux</FONT>
+			are the corresponding new inipar parameters.<BR>(<FONT FACE="Courier New, Courier, monospace">check_parameters</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">init_3d_model</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">parin</FONT>)<BR><BR>Large
+			scale vertical motion (subsidence/ascent) can be added as
+			additional source to the prognostic equation for potential
+			temperature. See new parameters <FONT FACE="Courier New, Courier, monospace">large_scale_subsidence</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">ws_vertical_gradient</FONT>,
+			and <FONT FACE="Courier New, Courier, monospace">ws_vertical_gradient_level</FONT>.
+			(<FONT FACE="Courier New, Courier, monospace">check_parameters</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">header</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">Makefile</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">modules</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">parin</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">prognostic_equations</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">read_var_list</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">subsidence</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">write_var_list</FONT>)<BR><BR>Scripts
+			have been updated for Cray-XT5 machine (lcxt5m) and
+			Hitachi-machine (lckyoto).<BR><BR>A simple method for installing
+			and running palm (with limited features) has been added.
+			(<FONT FACE="Courier New, Courier, monospace">Makefile</FONT>, new
+			scripts <B>palm_simple_install</B>, <B>palm_simple_run</B>)<BR><BR>Masked
+			data output has been implemented but it is not documented in the
+			user manual so far. This feature can already use NetCDF4 file
+			format. (check_open, check_parameters, close_file, modules,
+			netcdf, parin, time_integration, user_parin, new:
+			data_output_mask, init_masks, user_data_output_mask)</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>check_open, check parameters, close_file, header,
+			init_3d_model, Makefile, mbuild, modules, mrun, netcdf, parin,
+			prognostic_equations, read_var_list, time_integration, user_parin,
+			write_var_list<BR><BR><B>new:</B> data_output_mask,&nbsp;
+			init_masks, palm_simple_install, palm_simple_run, subsidence,
+			user_data_output_mask</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=56>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=68>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=62>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>C</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>2d-decomposition is default for Cray-XT machines.
+			(<FONT FACE="Courier New, Courier, monospace">init_pegrid</FONT>)<BR><BR><FONT FACE="Courier New, Courier, monospace">var_ts</FONT>
+			is replaced by <FONT FACE="Courier New, Courier, monospace">dots_max</FONT>.
+			(<FONT FACE="Courier New, Courier, monospace">modules</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">init_3d_model</FONT>)<BR><BR>Every
+			cloud droplet has now an own weighting factor and can be deleted
+			due to collisions. Condensation and collision of cloud droplets
+			are adjusted accordingly. (<FONT FACE="Courier New, Courier, monospace">advec_particles</FONT>)<BR><BR>Collision
+			efficiency for large cloud droplets has changed according to table
+			of Rogers and Yau. (<FONT FACE="Courier New, Courier, monospace">collision_efficiency</FONT>)</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>advec_particles, collision_efficiency, init_3d_model,
+			init_pegrid, modules</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=56>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=68>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=62>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>B</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>Bugfix for generating serial jobs. (<B>subjob</B>)<BR><BR>Bugfix:
+			index problem concerning gradient_level indices removed.
+			(<FONT FACE="Courier New, Courier, monospace">header</FONT>)<BR><BR>Dimension
+			of array <FONT FACE="Courier New, Courier, monospace">stat</FONT>
+			in cascade change to prevent type problems with mpi2 libraries.
+			(<FONT FACE="Courier New, Courier, monospace">poisfft_hybrid</FONT>)<BR><BR>Loop
+			was split to make runs reproducible when using ifort compiler.
+			(<FONT FACE="Courier New, Courier, monospace">disturb_field</FONT>)<BR><BR>Bugfix:
+			exchange of ghost points for <FONT FACE="Courier New, Courier, monospace">prho</FONT>
+			included. (<FONT FACE="Courier New, Courier, monospace">time_integration</FONT>)<BR><BR>Bugfix
+			in calculation of time-averaged surface heatfluxes.
+			(<FONT FACE="Courier New, Courier, monospace">sum_up_3d_data</FONT>)<BR><BR>Bugfix
+			in calculation of precipitation_rate.
+			(<FONT FACE="Courier New, Courier, monospace">calc_precipitation</FONT>)<BR><BR>Bugfix:
+			initial data assignments to some dvrp arrays changed due to error
+			messages from gfortran compiler. (<FONT FACE="Courier New, Courier, monospace">modules</FONT>)<BR><BR>Bugfix
+			in calculation of cloud droplet velocity.
+			(<FONT FACE="Courier New, Courier, monospace">advec_particles</FONT>)<BR><BR>Bugfix
+			in transfer of particles at south/left edge.
+			(<FONT FACE="Courier New, Courier, monospace">advec_particles</FONT>)<BR><BR>Bugfix
+			in calculation of collision efficiency. (<FONT FACE="Courier New, Courier, monospace">collision_efficiency</FONT>)</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>advec_particles, calc_precipitation, collision_efficiency,
+			disturb_field, header, modules, poisfft_hybrid, subjob,
+			sum_up_3d_data, time_integration</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>23/09/09</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>SR</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>392</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.7</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>N</P>
+		</TD>
+		<TD WIDTH=518>
+			<P><B>Most important changes:</B> atmosphere-ocean coupling has
+			been improved and steering of dvr-software is more user-friendly
+			now. Beside that, the release contains a lot of optimizations,
+			changes, and bugfixes which have accumulated over the last
+			months.<BR><BR>The atmosphere-ocean coupling allows to use
+			independent precursor runs in order to account for different
+			spin-up times. The time when coupling has to be started is given
+			by new inipar parameter <FONT FACE="Courier New, Courier, monospace">coupling_start_time</FONT>.
+			The precursor ocean run has to be started using new <B>mrun</B>
+			option &quot;<FONT FACE="Courier New, Courier, monospace">-y</FONT>&quot;
+			in order to add the appendix &quot;<FONT FACE="Courier New, Courier, monospace">_O</FONT>&quot;
+			to all output files. (<FONT FACE="Courier New, Courier, monospace">check_for_restart</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">check_parameters</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">data_output_2d</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">data_output_3d</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">data_output_profiles</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">data_output_ptseries</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">data_output_spectra</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">data_output_tseries</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">header</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">init_coupling</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">modules</FONT>, <FONT FACE="Courier New, Courier, monospace">mrun</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">parin</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">read_var_list</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">surface_coupler</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">time_integration</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">write_var_list</FONT>)<BR><BR>Clipping
+			of dvrp output implemented. Default colourtable for particles
+			implemented, particle attributes (color, dvrp_size) can be set
+			with new parameters <FONT FACE="Courier New, Courier, monospace">particle_color</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">particle_dvrpsize</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">color_interval</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">dvrpsize_interval</FONT>.
+			Slicer attributes (dvrp) are set with new routine
+			<FONT FACE="Courier New, Courier, monospace">set_slicer_attributes_dvrp</FONT>
+			and are controlled with existing parameters <FONT FACE="Courier New, Courier, monospace">slicer_range_limits</FONT>.
+			(<FONT FACE="Courier New, Courier, monospace">init_dvrp</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">data_output_dvrp</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">modules</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">set_slicer_attributes_dvrp</FONT>,&nbsp;
+			<FONT FACE="Courier New, Courier, monospace">user_data_output_dvrp</FONT>)<BR><BR>Polygon
+			reduction for topography and ground plate isosurface. Reduction
+			level for buildings can be chosen with parameter cluster_size.
+			(<FONT FACE="Courier New, Courier, monospace">init_dvrp</FONT>)<BR><BR>Variables
+			<FONT FACE="Courier New, Courier, monospace">bc_lr</FONT> / <FONT FACE="Courier New, Courier, monospace">bc_ns</FONT>&nbsp;
+			in most subroutines replaced by LOGICAL variables
+			<FONT FACE="Courier New, Courier, monospace">bc_lr_cyc</FONT>,<BR><FONT FACE="Courier New, Courier, monospace">bc_ns_cyc</FONT>
+			for speed optimization. This gives a significant speedup of more
+			than 10%. (<FONT FACE="Courier New, Courier, monospace">check_parameters</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">diffusion_u</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">diffusion_v</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">diffusion_w</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">modules</FONT>)<BR><BR>Scripts
+			have been adapted for machine lck (Yonsei Univ.). (<FONT FACE="Courier New, Courier, monospace">mrun</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">mbuild</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">subjob</FONT>)<BR><BR>Additional
+			timestep criterion in case of simulations with plant canopy.
+			(<FONT FACE="Courier New, Courier, monospace">timestep</FONT>)<BR><BR>Check
+			for illegal entries in <FONT FACE="Courier New, Courier, monospace">section_xy</FONT>|<FONT FACE="Courier New, Courier, monospace">xz</FONT>|<FONT FACE="Courier New, Courier, monospace">yz</FONT>
+			that exceed <FONT FACE="Courier New, Courier, monospace">nz+1</FONT>|<FONT FACE="Courier New, Courier, monospace">ny+1</FONT>|<FONT FACE="Courier New, Courier, monospace">nx+1</FONT>.
+			(<FONT FACE="Courier New, Courier, monospace">check_parameters</FONT>)<BR><BR>External
+			pressure gradient can be used as driving force. (<FONT FACE="Courier New, Courier, monospace">check_parameters</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">header</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">init_3d_model</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">modules</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">parin</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">prognostic_equations</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">read_var_list</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">write_var_list</FONT>)<BR><BR>New
+			topography case '<I>single_street_canyon</I>'. (<FONT FACE="Courier New, Courier, monospace">header</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">init_grid</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">modules</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">parin</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">read_var_list</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">user_check_parameters</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">user_header</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">user_init_grid</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">write_var_list</FONT>)<BR><BR>Option
+			to predefine a target bulk velocity for <FONT FACE="Courier New, Courier, monospace">conserve_volume_flow</FONT>.
+			(<FONT FACE="Courier New, Courier, monospace">check_parameters</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">header</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">init_3d_model</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">modules</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">parin</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">read_var_list</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">write_var_list</FONT>)<BR><BR>Option
+			for user defined 2D data output in xy cross sections at z=<FONT FACE="Courier New, Courier, monospace">nzb</FONT>+1.
+			(<FONT FACE="Courier New, Courier, monospace">data_output_2d</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">user_data_output_2d</FONT>)<BR><BR>xy
+			cross section output of surface heatfluxes (sensible and latent).
+			(<FONT FACE="Courier New, Courier, monospace">average_3d_data</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">check_parameters</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">data_output_2d</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">modules</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">read_3d_binary</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">sum_up_3d_data</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">write_3d_binary</FONT>)</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>average_3d_data, check_for_restart, check_parameters,
+			data_output_2d, data_output_3d, data_output_dvrp,
+			data_output_profiles, data_output_ptseries, data_output_spectra,
+			data_output_tseries, diffusion_u, diffusion_v, diffusion_w,
+			init_coupling, init_dvrp, init_grid, init_3d_model, header,
+			mbuild, modules, mrun, package_parin, parin, prognostic_equations,
+			read_3d_binary, read_var_list, subjob, sum_up_3d_data,
+			surface_coupler, timestep, time_integration,
+			user_check_parameters, user_data_output_2d, user_data_output_dvrp,
+			user_header, user_init_grid, write_3d_binary, write_var_list<BR><BR><B>new:</B>
+			set_particle_attributes, set_slicer_attributes_dvrp</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=56>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=68>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=62>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>C</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>Output of messages replaced by message handling routine.
+			(<FONT FACE="Courier New, Courier, monospace">advec_particles</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">advec_s_bc</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">buoyancy</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">calc_spectra</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">check_for_restart</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">check_open</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">coriolis</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">cpu_log</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">data_output_2d</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">data_output_3d</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">data_output_dvrp</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">data_output_profiles</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">data_output_spectra</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">fft_xy</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">flow_statistics</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">header</FONT>,&nbsp;
+			<FONT FACE="Courier New, Courier, monospace">init_1d_model</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">init_3d_model</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">init_dvrp</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">init_grid</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">init_particles</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">init_pegrid</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">netcdf</FONT>, <FONT FACE="Courier New, Courier, monospace">parin</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">plant_canopy_model</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">poisfft_hybrid</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">poismg</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">read_3d_binary</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">read_var_list</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">surface_coupler</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">temperton_fft</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">timestep</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">user_actions</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">user_data_output_dvrp</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">user_dvrp_coltab</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">user_init_grid</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">user_init_plant_canopy</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">user_parin</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">user_read_restart_data</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">user_spectra</FONT>)<BR><BR>Output
+			of NetCDF messages with aid of message handling routine.
+			(<FONT FACE="Courier New, Courier, monospace">check_open</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">close_file</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">data_output_2d</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">data_output_3d</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">data_output_profiles</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">data_output_ptseries</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">data_output_spectra</FONT>,&nbsp;
+			<FONT FACE="Courier New, Courier, monospace">data_output_tseries</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">netcdf</FONT>)<BR><BR>__lcmuk
+			changed to __lc to avoid problems with Intel compiler on sgi-ice.
+			(<FONT FACE="Courier New, Courier, monospace">poisfft</FONT>)<BR><BR>For
+			extended NetCDF files, the updated title attribute includes an
+			update of time_average_text where appropriate. (<FONT FACE="Courier New, Courier, monospace">netcdf</FONT>)<BR><BR>In
+			case of restart runs without extension, initial profiles are not
+			written to NetCDF-file anymore. So far, time axis always started
+			at t=0. (<FONT FACE="Courier New, Courier, monospace">data_output_profiles</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">modules</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">read_var_list</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">write_var_list</FONT>)<BR><BR><FONT FACE="Courier New, Courier, monospace">initializing_actions</FONT>='<I>read_data_for_recycling</I>'
+			renamed to '<I>cyclic_fill</I>'. A cyclic fill can now <BR>be used
+			independent of <FONT FACE="Courier New, Courier, monospace">turbulent_inflow</FONT>.
+			(<FONT FACE="Courier New, Courier, monospace">check_parameters</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">header</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">init_3d_model</FONT>)<BR><BR>2
+			NetCDF error numbers changed. (<FONT FACE="Courier New, Courier, monospace">data_output_3d</FONT>)<BR><BR>A
+			link to the webspage appendix_a.html is printed for further
+			information about possible reasons of the error that appeared.
+			(message)<BR><BR>Temperature gradient criterion for estimating the
+			boundary layer height replaced by the gradient criterion of
+			Sullivan et al. (1998). (<FONT FACE="Courier New, Courier, monospace">flow_statistics</FONT>)<BR><BR>NetCDF
+			unit attribute in timeseries output in case of statistic regions
+			added. (<FONT FACE="Courier New, Courier, monospace">netcdf</FONT>)<BR><BR>Maximum
+			number of tails is calculated from maximum number of particles and
+			<FONT FACE="Courier New, Courier, monospace">skip_particles_for_tail</FONT>.
+			(<FONT FACE="Courier New, Courier, monospace">init_particles</FONT>)<BR><BR>Value
+			of <FONT FACE="Courier New, Courier, monospace">vertical_particle_advection</FONT>
+			is allowed to differ for each particle group. (<FONT FACE="Courier New, Courier, monospace">advec_particles</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">header</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">modules</FONT>)
+			<BR><BR>First constant in array <FONT FACE="Courier New, Courier, monospace">den</FONT>
+			also defined as type double. (<FONT FACE="Courier New, Courier, monospace">eqn_state_seawater</FONT>)<BR><BR>Parameter
+			<FONT FACE="Courier New, Courier, monospace">dvrp_psize</FONT>
+			moved from <FONT FACE="Courier New, Courier, monospace">particles_par</FONT>
+			to <FONT FACE="Courier New, Courier, monospace">dvrp_graphics_par</FONT>.
+			(<FONT FACE="Courier New, Courier, monospace">package_parin</FONT>)<BR><BR><FONT FACE="Courier New, Courier, monospace">topography_grid_convention</FONT>
+			moved from <FONT FACE="Courier New, Courier, monospace">userpar</FONT>
+			to <FONT FACE="Courier New, Courier, monospace">inipar</FONT>
+			(<FONT FACE="Courier New, Courier, monospace">check_parameters</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">header</FONT>, <FONT FACE="Courier New, Courier, monospace">parin</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">read_var_list</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">user_check_parameters</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">user_header</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">user_init_grid</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">user_parin</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">write_var_list</FONT>)<BR><BR>Default
+			value of <FONT FACE="Courier New, Courier, monospace">grid_matching</FONT>
+			changed to '<I>strict</I>'. (<FONT FACE="Courier New, Courier, monospace">modules</FONT>)<BR><BR>Adjustments
+			for runs on lcxt4 (necessary due to a software update on CRAY) and
+			for coupled runs on ibmy. (<FONT FACE="Courier New, Courier, monospace">mrun</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">subjob</FONT>)</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>advec_particles, advec_s_bc, buoyancy, calc_spectra,
+			check_for_restart, check_open, check_parameters, close_file,
+			coriolis, cpu_log, data_output_2d, data_output_3d,
+			data_output_dvrp, data_output_profiles, data_output_ptseries,
+			data_output_spectra, data_output_tseries, eqn_state_seawater,
+			fft_xy, flow_statistics, header, init_1d_model, init_3d_model,
+			init_dvrp, init_grid, init_particles, init_pegrid, message, mrun,
+			netcdf, package_parin, parin, plant_canopy_model, poisfft,
+			poisfft_hybrid, poismg, read_3d_binary, read_var_list, subjob,
+			surface_coupler, temperton_fft, timestep,
+			user_actions,user_check_parameters, user_data_output_dvrp,
+			user_dvrp_coltab, user_header, user_init_grid,
+			user_init_plant_canopy, user_parin, user_read_restart_data,
+			user_spectra, write_var_list 
+			</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=56>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=68>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=62>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>B</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>Bugfix: Initial hydrostatic pressure profile in case of ocean
+			runs is now calculated in 5 iteration steps. (<FONT FACE="Courier New, Courier, monospace">init_ocean</FONT>)<BR><BR>Bugfix:
+			wrong sign in buoyancy production of ocean part in case of not
+			using the reference density (only in 3D routine <FONT FACE="Courier New, Courier, monospace">production_e</FONT>).
+			(<FONT FACE="Courier New, Courier, monospace">production_e</FONT>)<BR><BR>Bugfix:
+			output of averaged 2d/3d quantities requires that an avaraging
+			interval has been set, respective error message is included.
+			(<FONT FACE="Courier New, Courier, monospace">check_parameters</FONT>)<BR><BR>Bugfix:
+			Output on unit 14 only if requested by <FONT FACE="Courier New, Courier, monospace">write_binary</FONT>.
+			(<FONT FACE="Courier New, Courier, monospace">user_last_actions</FONT>)<BR><BR>Bugfix
+			to avoid zero division by <FONT FACE="Courier New, Courier, monospace">km_neutral</FONT>.
+			(<FONT FACE="Courier New, Courier, monospace">production_e</FONT>)<BR><BR>Bugfix
+			for extended NetCDF files: In order to avoid 'data mode' errors if
+			updated attributes are larger than their original size,
+			<FONT FACE="Courier New, Courier, monospace">NF90_PUT_ATT</FONT>
+			is called in 'define mode' enclosed by <FONT FACE="Courier New, Courier, monospace">NF90_REDEF</FONT>
+			and <FONT FACE="Courier New, Courier, monospace">NF90_ENDDEF</FONT>
+			calls. This implies a possible performance loss; an alternative
+			strategy would be to ensure equal attribute size in a job chain.
+			(<FONT FACE="Courier New, Courier, monospace">netcdf</FONT>)<BR><BR>Bugfix:
+			correction of initial volume flow for non-flat topography.
+			(<FONT FACE="Courier New, Courier, monospace">init_3d_model</FONT>)<BR><BR>Bugfix:
+			zero initialization of arrays within buildings for '<I>cyclic_fill</I>'.
+			(<FONT FACE="Courier New, Courier, monospace">init_3d_model</FONT>)<BR><BR>Bugfix:
+			<FONT FACE="Courier New, Courier, monospace">to_be_resorted</FONT>
+			=&gt; <FONT FACE="Courier New, Courier, monospace">s_av</FONT> for
+			time-averaged scalars. (<FONT FACE="Courier New, Courier, monospace">data_output_2d</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">data_output_3d</FONT>)<BR><BR>Bugfix:
+			avoid that <FONT FACE="Courier New, Courier, monospace">ngp_2dh_s_inner</FONT>
+			becomes zero. (<FONT FACE="Courier New, Courier, monospace">init_3_model</FONT>)<BR><BR>Typographical
+			error: unit of wpt in dots_unit. (<FONT FACE="Courier New, Courier, monospace">modules</FONT>)<BR><BR>Bugfix:
+			error in check, if particles moved further than one subdomain
+			length. This check must not be applied for newly released
+			particles. (<FONT FACE="Courier New, Courier, monospace">advec_particles</FONT>)<BR><BR>Bugfix:
+			several tail counters are initialized, <FONT FACE="Courier New, Courier, monospace">particle_tail_coordinates</FONT>
+			is only written to file if its third index is &gt; 0, arrays for
+			tails are allocated with a minimum size of 10 tails if there is no
+			tail initially. (<FONT FACE="Courier New, Courier, monospace">advec_particles</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">init_particles</FONT>)<BR><BR>Bugfix:
+			pressure included for profile output. (<FONT FACE="Courier New, Courier, monospace">check_parameters</FONT>)<BR><BR>Bugfix:
+			Type of <FONT FACE="Courier New, Courier, monospace">count</FONT>
+			and <FONT FACE="Courier New, Courier, monospace">count_rate</FONT>
+			changed to default <FONT FACE="Courier New, Courier, monospace">INTEGER</FONT>
+			on NEC machines. (<FONT FACE="Courier New, Courier, monospace">cpu_log</FONT>)<BR><BR>Bugfix:
+			output of particle time series only if particle advection is
+			switched on. (<FONT FACE="Courier New, Courier, monospace">time_integration</FONT>)<BR><BR>Bugfix:
+			qsws was calculated in case of <FONT FACE="Courier New, Courier, monospace">constant
+			heatflux</FONT> = <FONT FACE="Courier New, Courier, monospace">.FALSE.</FONT>.
+			(<FONT FACE="Courier New, Courier, monospace">prandtl_fluxes</FONT>)<BR><BR>Bugfix:
+			averaging along z is not allowed for 2d quantities (e.g. u* and
+			z0). (<FONT FACE="Courier New, Courier, monospace">data_output_2d</FONT>)<BR><BR>Typographical
+			errors. (<FONT FACE="Courier New, Courier, monospace">netcdf</FONT>)<BR><BR>If
+			the inversion height calculated by the prerun is zero,
+			inflow_damping_height must be explicitly specified.
+			(<FONT FACE="Courier New, Courier, monospace">init_3d_model</FONT>)<BR><BR>Small
+			bugfix concerning 3d 64bit netcdf output format. (<FONT FACE="Courier New, Courier, monospace">header</FONT>)<BR><BR>Bugfix:
+			<FONT FACE="Courier New, Courier, monospace">dt_fixed</FONT>
+			removed from the restart file, because otherwise, no change from a
+			fixed to a variable timestep would be possible in restart runs.
+			(<FONT FACE="Courier New, Courier, monospace">read_var_list</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">write_var_list</FONT>)<BR><BR>Bugfix:
+			initial setting of <FONT FACE="Courier New, Courier, monospace">time_coupling</FONT>
+			in coupled restart runs. (<FONT FACE="Courier New, Courier, monospace">time_integration</FONT>)</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>advec_particles, check_parameters, cpu_log, data_output_2d,
+			data_output_3d, header, init_3d_model, init_particles, init_ocean,
+			modules, netcdf, prandtl_fluxes, production_e, read_var_list,
+			time_integration, user_last_actions, write_var_list</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P>02/02/09</P>
+		</TD>
+		<TD WIDTH=56>
+			<P>SR</P>
+		</TD>
+		<TD WIDTH=68>
+			<P>228</P>
+		</TD>
+		<TD WIDTH=62>
+			<P>3.6</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>N</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>This release contains further additions, changes, and bugfixes
+			for the SGI-ICE system, as well as important changes like the
+			splitting of the user-interface into single files, a complete
+			revision of the pre-compiling mechanism allowing to use different
+			make-depositories simultaneously, and a revision of the
+			reading-mechanism from restart files, which has no more
+			restrictions in case that the previous and current run have
+			different domain/subdomain sizes. A new handling of all kind of
+			output messages is introduced with this release and will be
+			applied to all existing messages within the next releases.<BR><BR><B>New
+			features are:</B><BR><BR>Restart runs on SGI-ICE are working.
+			(mrun)<BR>2d-decomposition is default on SGI-ICE systems.
+			(init_pegrid)<BR><BR>Ocean-atmosphere coupling realized with
+			MPI-1. mrun adjusted for this case (-Y option). Check that PALM is
+			called with <FONT FACE="Courier New, Courier, monospace">mrun</FONT>-option
+			&quot;<I>-K parallel</I>&quot; for coupling. Adjustments in <FONT FACE="Courier New, Courier, monospace">mrun</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">mbuild</FONT>, and
+			<FONT FACE="Courier New, Courier, monospace">subjob</FONT> for
+			lcxt4.<BR><BR>DVRP arguments changed to single precision, mode
+			pathlines added.<BR><BR>User can add additional routines in files
+			<FONT FACE="Courier New, Courier, monospace">user_additional_routines</FONT>.<BR><BR>User
+			can check user parameters and deduce further quantities in
+			<FONT FACE="Courier New, Courier, monospace">user_check_parameters</FONT></P>
+		</TD>
+		<TD WIDTH=393>
+			<P>check_for_restart, check_parameters, data_output_dvrp,
+			init_dvrp, init_pegrid, local_stop, modules, package_parin, palm,
+			surface_coupler, timestep<BR><BR>Makefile,&nbsp;mbuild, mrun,
+			subjob<BR><BR><B>new:</B><BR>init_coupling,
+			user_additional_routines, user_check_parameters</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=56>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=68>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=62>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>C</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>User interface has been split into single files.<BR><BR>Type of
+			variables <FONT FACE="Courier New, Courier, monospace">count</FONT>
+			and <FONT FACE="Courier New, Courier, monospace">count_rate</FONT>
+			changed to INTEGER(8) in order to avoid &quot;out of range&quot;
+			problems, which result in measured negative time intervals.
+			(cpu_log, local_tremain, local_tremain_ini)<BR><BR>Reading
+			mechanism from restart files completely revised. The
+			subdomain/total domain size is now allowed to vary arbitrarily
+			between the current and previous run. (read_var_list,
+			read_3d_binary, user_read_restart_data)<BR><BR>Precompilation
+			mechanism (<FONT FACE="Courier New, Courier, monospace">mbuild</FONT>)
+			completely revised: the source-tarfile remains in the source
+			directory (<FONT FACE="Courier New, Courier, monospace">trunk/SOURCE</FONT>).
+			One depository per block (given in the config-file) is created.
+			Always all files from the source directory are copied to the
+			respective depository. No additional file checks are done any more
+			(mrun version 2.1). The depository name used contains the
+			conditions given by <FONT FACE="Courier New, Courier, monospace">mrun</FONT>-option
+			&quot;<FONT FACE="Courier New, Courier, monospace">-K</FONT>&quot;.
+			(mbuild, mrun)<BR><BR>Output messages (including required stop of
+			execution) can now be handled with the new subroutine&nbsp;<FONT FACE="Courier New, Courier, monospace">message</FONT>.
+			All output messages will be replaced by this routine step by step
+			within the next revisions.<BR><BR>Neumann boundary condition at
+			<FONT FACE="Courier New, Courier, monospace">k=nzb</FONT> is
+			explicitly set for better reading, although this has been already
+			done in boundary_conds. (advec_s_bc)<BR><BR>Origin of the
+			xy-coordinate system shifted from the center of the first grid
+			cell (indices <FONT FACE="Courier New, Courier, monospace">i=0</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">j=0</FONT>) to the
+			south-left corner of this cell. (netcdf) Topography definition
+			according to new user parameter <FONT FACE="Courier New, Courier, monospace">topography_grid_convention</FONT>.
+			(init_grid, modules, user_header, user_init_grid, user_parin)</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>advec_s_bc, check_parameters, cpu_log, init_grid, local_stop,
+			local_tremain, local_tremain_ini, modules, netcdf, read_3d_binary,
+			read_var_list, user_read_restart_data<BR><BR>Makefile, mbuild,
+			mrun<BR><BR><B>removed:</B><BR>user_interface<BR><BR><B>new:</B><BR>message,
+			user_3d_data_averaging, user_actions, user_advec_particles,
+			user_check_data_output, user_check_data_output_pr,
+			user_data_output_2d, user_data_output_3d, user_data_output_dvrp,
+			user_define_netcdf_grid, user_dvrp_coltab, user_header,
+			user_init_3d_model, user_init, user_init_grid,
+			user_init_particles, user_init_plant_canopy, user_last_actions,
+			user_parin, user_parin, user_particle_attributes,
+			user_read_restart_data, user_spectra, user_statistics</P>
+		</TD>
+	</TR>
+	<TR>
+		<TD WIDTH=37>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=56>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=68>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=62>
+			<P><BR>
+			</P>
+		</TD>
+		<TD WIDTH=36>
+			<P>B</P>
+		</TD>
+		<TD WIDTH=518>
+			<P>Bugfix: array <FONT FACE="Courier New, Courier, monospace">d</FONT>
+			is reallocated in case that multigrid is used.
+			(calc_spectra)<BR><BR>Bugfixes for nonparallel execution.
+			(check_for_restart, cpu_statistics,
+			inflow_turbulence,&nbsp;timestep)<BR><BR>Size of <FONT FACE="Courier New, Courier, monospace">pf3d
+			</FONT>adjusted to the required output size (1 gridpoint less,
+			along all three dimensions), because output of a subset of the
+			data (<FONT FACE="Courier New, Courier, monospace">pf3d(nxa:nxe...)</FONT>)
+			in the NF90_PUT_VAR statement caused segmentation fault with the
+			INTEL compiler. (combine_plot_fields)<BR><BR>Bugfix: error in <FONT FACE="Courier New, Courier, monospace">zu</FONT>
+			index in case of <FONT FACE="Courier New, Courier, monospace">section_xy=-1</FONT>.
+			(header)<BR><BR>Inconsistency removed: as the thermal
+			stratification is not taken into account for the evaluation of the
+			wall fluxes at vertical walls, the eddy viscosity <FONT FACE="Courier New, Courier, monospace">km</FONT>
+			must not be used for the evaluation of the velocity gradients
+			<FONT FACE="Courier New, Courier, monospace">dudy</FONT>, <FONT FACE="Courier New, Courier, monospace">dwdy</FONT>,
+			<FONT FACE="Courier New, Courier, monospace">dvdx</FONT> and <FONT FACE="Courier New, Courier, monospace">dwdx</FONT>.
+			(production_e)<BR><BR>Bugfix in calculating <FONT FACE="Courier New, Courier, monospace">k</FONT>
+			index in case of oceans runs. (sort_particles)<BR><BR>Bugfix: no
+			output of particle concentration and radius unless particles have
+			been started. (data_output_2d)<BR><BR>Bugfix: reading of
+			<FONT FACE="Courier New, Courier, monospace">spectrum_x|y</FONT>
+			from restart files ignored if total numbers of grid points do not
+			match. (read_3d_binary)<BR><BR>Bugfix: abort in case that absolute
+			temperature is below zero. (init_cloud_physics)</P>
+		</TD>
+		<TD WIDTH=393>
+			<P>advec_particles, calc_spectra, check_for_restart,
+			cpu_statistics, data_output_2d, header, inflow_turbulence,
+			init_cloud_physics, production_e, read_3d_binary,
+			timestep<BR><BR>combine_plot_fields</P>
+		</TD>
+	</TR>
+</TABLE>
+<P><BR>&nbsp; <BR>&nbsp; 
+</P>
+<H2><A NAME="Kapitel2.0"></A>2.0&nbsp; How to change the source code
+and how to document modifications</H2>
+<P>Currently (Sep 23rd, 2009), only selected users are allowed to
+commit changes to the repository. Rules given below are preliminary,
+until we have switched to the &quot;trac&quot;-System.</P>
+<OL>
+	<LI><P>Change source code in your current working copy. Document
+	changes in the respective file headers unter &quot;<FONT FACE="Courier New, Courier, monospace">Current
+	revisions:</FONT>&quot;. Additionally, document all changes in a
+	local file (e.g. ~/palm/current_version/CURRENT_MODIFICATIONS),
+	including the files that have been changed. Classify changes in
+	three cases, new (N), changed (C), and bugfix (B). <B>Try to be as
+	short and as precise as possible in describing the changes!</B></P>
+	<LI><P>Test the changed version, unless it works sufficiently.</P>
+	<LI><P>Commit the changes to the repository<BR><BR>&nbsp;&nbsp; <FONT FACE="Courier New, Courier, monospace">svn
+	commit -m &quot;message string&quot; trunk<BR></FONT><BR>where
+	&quot;<FONT FACE="Courier New, Courier, monospace">message string</FONT>&quot;
+	is a short, <B>meaningful</B> summary of the
+	changes.<BR><BR>Alternatively, in case that you have made quite a
+	lot of changes, you can also give the contents of your file
+	CURRENT_MODIFICATIONS as the log message:<BR><BR>&nbsp;&nbsp; <FONT FACE="Courier New, Courier, monospace">svn
+	commit -F CURRENT_MODIFICATIONS trunk</FONT></P>
+	<LI><P>Enter the changes into the technical documentation (under
+	<FONT FACE="Courier New, Courier, monospace">.../trunk/DOC/tec/technical_documentation.html</FONT>).</P>
+	<LI><P>Move the change comments in the file headers from &quot;Current
+	revisions:&quot; to &quot;Former revisions:&quot;. Set the current
+	svn id in front of that.<BR><BR><B>Example:</B><BR><BR><B>before
+	move:</B><BR><BR><FONT FACE="Courier New, Courier, monospace">!
+	Current revisions:</FONT><BR><FONT FACE="Courier New, Courier, monospace">!
+	-----------------</FONT><BR><FONT FACE="Courier New, Courier, monospace">!
+	</FONT><FONT FACE="Courier New, Courier, monospace"><B>new dummy
+	argument abcd</B></FONT><BR><FONT FACE="Courier New, Courier, monospace">!</FONT><BR><FONT FACE="Courier New, Courier, monospace">!
+	Former revisions:</FONT><BR><FONT FACE="Courier New, Courier, monospace">!
+	-----------------</FONT><BR><FONT FACE="Courier New, Courier, monospace">!
+	Ã¯Â¿ÂœId: init_3d_model.f90 </FONT><FONT COLOR="#3366ff"><FONT FACE="Courier New, Courier, monospace"><B>425
+	2010-02-06 13:37:48Z raasch</B></FONT></FONT><FONT FACE="Courier New, Courier, monospace">
+	Ã¯Â¿Âœ</FONT><BR><FONT FACE="Courier New, Courier, monospace">!</FONT><BR><FONT FACE="Courier New, Courier, monospace">!
+	407 2009-12-01 15:01:15Z maronga</FONT><BR><FONT FACE="Courier New, Courier, monospace">!
+	var_ts is replaced by dots_max </FONT><BR><BR><B>after move:</B><BR><BR><FONT FACE="Courier New, Courier, monospace">!
+	Current revisions:</FONT><BR><FONT FACE="Courier New, Courier, monospace">!
+	-----------------</FONT><BR><FONT FACE="Courier New, Courier, monospace">!
+	</FONT><BR><FONT FACE="Courier New, Courier, monospace">!</FONT><BR><FONT FACE="Courier New, Courier, monospace">!
+	Former revisions:</FONT><BR><FONT FACE="Courier New, Courier, monospace">!
+	-----------------</FONT><BR><FONT FACE="Courier New, Courier, monospace">!
+	Ã¯Â¿ÂœId: init_3d_model.f90 425 2010-02-06 13:37:48Z raasch</FONT>
+	<FONT FACE="Courier New, Courier, monospace">Ã¯Â¿Âœ<BR>!</FONT><BR><FONT FACE="Courier New, Courier, monospace">!
+	</FONT><FONT COLOR="#3366ff"><FONT FACE="Courier New, Courier, monospace"><B>425
+	2010-02-06 13:37:48Z raasch</B></FONT></FONT><FONT FACE="Courier New, Courier, monospace"><BR>!
+	</FONT><FONT FACE="Courier New, Courier, monospace"><B>new dummy
+	argument abcd</B></FONT><BR><FONT FACE="Courier New, Courier, monospace">!</FONT><BR><FONT FACE="Courier New, Courier, monospace">!
+	407 2009-12-01 15:01:15Z maronga</FONT><BR><FONT FACE="Courier New, Courier, monospace">!
+	var_ts is replaced by dots_max<BR><BR></FONT><FONT FACE="Times New Roman, Times, serif">Note
+	that the &quot;</FONT><FONT FACE="Courier New, Courier, monospace">Ã¯Â¿Âœ</FONT><FONT FACE="Times New Roman, Times, serif">&quot;
+	in this example should read &quot;$&quot;, but &quot;$&quot; cannot
+	be used here, because otherwise svn would replace these entries
+	after every change to this technical documentation.</FONT></P>
+	<LI><P>Commit the changes in the technical documentation and the
+	file headers:<BR><BR>&nbsp;&nbsp; <FONT FACE="Courier New, Courier, monospace">svn
+	commit -m &quot;last commit documented&quot; trunk</FONT><BR><BR><BR>
+	</P>
+</OL>
+<H2><A NAME="Kapitel3.0"></A>3.0&nbsp; Description of selected parts
+of the model source code</H2>
+<P>This section gives links to existing descriptions about special
 parts of the PALM code, e.g. special numerical methods that are used
 (including mathematical background). Documents are stored in single
 files with different formats. Some are in German (sorry for the
-inconvenience).
-<br>&nbsp;
-</p><ul> <li> <a name="UPS"></a>upstream-spline
-advection scheme ( <a href="methods/upstream_spline/ups.pdf">.pdf</a>
-, only in German)</li> <li> <a name="Cloud_physics"></a>cloud
-physics module ( <a href="methods/cloud_physics/cloud_physics.pdf">.pdf</a>
-)</li> 
-  <li> <a name="Bitkompression"></a>data compression with <a href="methods/bit_compression/bit_compression.html">bit shifting method</a> (only in German)<br>
-  </li>
-  <li>Runge-Kutta
-time integration scheme ( <a href="methods/runge_kutta/runge_kutta.pdf">.pdf</a>
-, only in German)</li>
-
-</ul><br>
-&nbsp;
-<br>&nbsp;
-<br>&nbsp;
-<br>&nbsp;
-<br>&nbsp;
-<br>&nbsp;
-<br>&nbsp;
-<br>&nbsp;
-</body></html>
+inconvenience). <BR>&nbsp; 
+</P>
+<UL>
+	<LI><P STYLE="margin-bottom: 0in"><A NAME="UPS"></A>upstream-spline
+	advection scheme ( <A HREF="methods/upstream_spline/ups.pdf">.pdf</A>
+	, only in German) 
+	</P>
+	<LI><P STYLE="margin-bottom: 0in"><A NAME="Cloud_physics"></A>cloud
+	physics module ( <A HREF="methods/cloud_physics/cloud_physics.pdf">.pdf</A>
+	) 
+	</P>
+	<LI><P STYLE="margin-bottom: 0in"><A NAME="Bitkompression"></A>data
+	compression with <A HREF="methods/bit_compression/bit_compression.html">bit
+	shifting method</A> (only in German)</P>
+	<LI><P>Runge-Kutta time integration scheme ( <A HREF="methods/runge_kutta/runge_kutta.pdf">.pdf</A>
+	, only in German) 
+	</P>
+</UL>
+<P><BR>&nbsp; <BR>&nbsp; <BR>&nbsp; <BR>&nbsp; <BR>&nbsp; <BR>&nbsp;
+<BR>&nbsp; <BR>&nbsp; 
+</P>
+</BODY>
+</HTML>
Index: palm/trunk/SOURCE/check_parameters.f90
===================================================================
--- palm/trunk/SOURCE/check_parameters.f90	(revision 552)
+++ palm/trunk/SOURCE/check_parameters.f90	(revision 553)
@@ -4,5 +4,5 @@
 ! Current revisions:
 ! -----------------
-! 
+! masks is calculated and removed from inipar
 !
 ! Former revisions:
@@ -2680,4 +2680,10 @@
 !
 !-- Check mask conditions
+    DO mid = 1, max_masks
+       IF ( data_output_masks(mid,1) /= ' ') THEN
+          masks = masks + 1
+       ENDIF
+    ENDDO
+    
     IF ( masks < 0 .OR. masks > max_masks )  THEN
        WRITE( message_string, * )  'illegal value: masks must be >= 0 and ', &
Index: palm/trunk/SOURCE/init_masks.f90
===================================================================
--- palm/trunk/SOURCE/init_masks.f90	(revision 552)
+++ palm/trunk/SOURCE/init_masks.f90	(revision 553)
@@ -4,5 +4,5 @@
 ! Current revisions:
 ! -----------------
-! 
+! parameters for masked output are replaced by arrays
 !
 ! Former revisions:
@@ -60,97 +60,20 @@
        CALL message( 'init_masks', 'PA9998', 1, 2, 0, 6, 0 )
     ENDIF
+
 !
 !-- Store data output parameters for masked data output in few shared arrays
-    do_mask( 1,:) = data_output_mask_01; do_mask( 2,:) = data_output_mask_02
-    do_mask( 3,:) = data_output_mask_03; do_mask( 4,:) = data_output_mask_04
-    do_mask( 5,:) = data_output_mask_05; do_mask( 6,:) = data_output_mask_06
-    do_mask( 7,:) = data_output_mask_07; do_mask( 8,:) = data_output_mask_08
-    do_mask( 9,:) = data_output_mask_09; do_mask(10,:) = data_output_mask_10
-    do_mask(11,:) = data_output_mask_11; do_mask(12,:) = data_output_mask_12
-    do_mask(13,:) = data_output_mask_13; do_mask(14,:) = data_output_mask_14
-    do_mask(15,:) = data_output_mask_15; do_mask(16,:) = data_output_mask_16
-    do_mask(17,:) = data_output_mask_17; do_mask(18,:) = data_output_mask_18
-    do_mask(19,:) = data_output_mask_19; do_mask(20,:) = data_output_mask_20
-    do_mask_user( 1,:) = data_output_mask_user_01
-    do_mask_user( 2,:) = data_output_mask_user_02
-    do_mask_user( 3,:) = data_output_mask_user_03
-    do_mask_user( 4,:) = data_output_mask_user_04
-    do_mask_user( 5,:) = data_output_mask_user_05
-    do_mask_user( 6,:) = data_output_mask_user_06
-    do_mask_user( 7,:) = data_output_mask_user_07
-    do_mask_user( 8,:) = data_output_mask_user_08
-    do_mask_user( 9,:) = data_output_mask_user_09
-    do_mask_user(10,:) = data_output_mask_user_10
-    do_mask_user(11,:) = data_output_mask_user_11
-    do_mask_user(12,:) = data_output_mask_user_12
-    do_mask_user(13,:) = data_output_mask_user_13
-    do_mask_user(14,:) = data_output_mask_user_14
-    do_mask_user(15,:) = data_output_mask_user_15
-    do_mask_user(16,:) = data_output_mask_user_16
-    do_mask_user(17,:) = data_output_mask_user_17
-    do_mask_user(18,:) = data_output_mask_user_18
-    do_mask_user(19,:) = data_output_mask_user_19
-    do_mask_user(20,:) = data_output_mask_user_20
-    mask( 1,1,:) = mask_01_x; mask( 2,1,:) = mask_02_x
-    mask( 3,1,:) = mask_03_x; mask( 4,1,:) = mask_04_x
-    mask( 5,1,:) = mask_05_x; mask( 6,1,:) = mask_06_x
-    mask( 7,1,:) = mask_07_x; mask( 8,1,:) = mask_08_x
-    mask( 9,1,:) = mask_09_x; mask(10,1,:) = mask_10_x
-    mask(11,1,:) = mask_11_x; mask(12,1,:) = mask_12_x
-    mask(13,1,:) = mask_13_x; mask(14,1,:) = mask_14_x
-    mask(15,1,:) = mask_15_x; mask(16,1,:) = mask_16_x
-    mask(17,1,:) = mask_17_x; mask(18,1,:) = mask_18_x
-    mask(19,1,:) = mask_19_x; mask(20,1,:) = mask_20_x
-    mask( 1,2,:) = mask_01_y; mask( 2,2,:) = mask_02_y
-    mask( 3,2,:) = mask_03_y; mask( 4,2,:) = mask_04_y
-    mask( 5,2,:) = mask_05_y; mask( 6,2,:) = mask_06_y
-    mask( 7,2,:) = mask_07_y; mask( 8,2,:) = mask_08_y
-    mask( 9,2,:) = mask_09_y; mask(10,2,:) = mask_10_y
-    mask(11,2,:) = mask_11_y; mask(12,2,:) = mask_12_y
-    mask(13,2,:) = mask_13_y; mask(14,2,:) = mask_14_y
-    mask(15,2,:) = mask_15_y; mask(16,2,:) = mask_16_y
-    mask(17,2,:) = mask_17_y; mask(18,2,:) = mask_18_y
-    mask(19,2,:) = mask_19_y; mask(20,2,:) = mask_20_y
-    mask( 1,3,:) = mask_01_z; mask( 2,3,:) = mask_02_z
-    mask( 3,3,:) = mask_03_z; mask( 4,3,:) = mask_04_z
-    mask( 5,3,:) = mask_05_z; mask( 6,3,:) = mask_06_z
-    mask( 7,3,:) = mask_07_z; mask( 8,3,:) = mask_08_z
-    mask( 9,3,:) = mask_09_z; mask(10,3,:) = mask_10_z
-    mask(11,3,:) = mask_11_z; mask(12,3,:) = mask_12_z
-    mask(13,3,:) = mask_13_z; mask(14,3,:) = mask_14_z
-    mask(15,3,:) = mask_15_z; mask(16,3,:) = mask_16_z
-    mask(17,3,:) = mask_17_z; mask(18,3,:) = mask_18_z
-    mask(19,3,:) = mask_19_z; mask(20,3,:) = mask_20_z
-    mask_loop( 1,1,:) = mask_01_x_loop; mask_loop( 2,1,:) = mask_02_x_loop
-    mask_loop( 3,1,:) = mask_03_x_loop; mask_loop( 4,1,:) = mask_04_x_loop
-    mask_loop( 5,1,:) = mask_05_x_loop; mask_loop( 6,1,:) = mask_06_x_loop
-    mask_loop( 7,1,:) = mask_07_x_loop; mask_loop( 8,1,:) = mask_08_x_loop
-    mask_loop( 9,1,:) = mask_09_x_loop; mask_loop(10,1,:) = mask_10_x_loop
-    mask_loop(11,1,:) = mask_11_x_loop; mask_loop(12,1,:) = mask_12_x_loop
-    mask_loop(13,1,:) = mask_13_x_loop; mask_loop(14,1,:) = mask_14_x_loop
-    mask_loop(15,1,:) = mask_15_x_loop; mask_loop(16,1,:) = mask_16_x_loop
-    mask_loop(17,1,:) = mask_17_x_loop; mask_loop(18,1,:) = mask_18_x_loop
-    mask_loop(19,1,:) = mask_19_x_loop; mask_loop(20,1,:) = mask_20_x_loop
-    mask_loop( 1,2,:) = mask_01_y_loop; mask_loop( 2,2,:) = mask_02_y_loop
-    mask_loop( 3,2,:) = mask_03_y_loop; mask_loop( 4,2,:) = mask_04_y_loop
-    mask_loop( 5,2,:) = mask_05_y_loop; mask_loop( 6,2,:) = mask_06_y_loop
-    mask_loop( 7,2,:) = mask_07_y_loop; mask_loop( 8,2,:) = mask_08_y_loop
-    mask_loop( 9,2,:) = mask_09_y_loop; mask_loop(10,2,:) = mask_10_y_loop
-    mask_loop(11,2,:) = mask_11_y_loop; mask_loop(12,2,:) = mask_12_y_loop
-    mask_loop(13,2,:) = mask_13_y_loop; mask_loop(14,2,:) = mask_14_y_loop
-    mask_loop(15,2,:) = mask_15_y_loop; mask_loop(16,2,:) = mask_16_y_loop
-    mask_loop(17,2,:) = mask_17_y_loop; mask_loop(18,2,:) = mask_18_y_loop
-    mask_loop(19,2,:) = mask_19_y_loop; mask_loop(20,2,:) = mask_20_y_loop
-    mask_loop( 1,3,:) = mask_01_z_loop; mask_loop( 2,3,:) = mask_02_z_loop
-    mask_loop( 3,3,:) = mask_03_z_loop; mask_loop( 4,3,:) = mask_04_z_loop
-    mask_loop( 5,3,:) = mask_05_z_loop; mask_loop( 6,3,:) = mask_06_z_loop
-    mask_loop( 7,3,:) = mask_07_z_loop; mask_loop( 8,3,:) = mask_08_z_loop
-    mask_loop( 9,3,:) = mask_09_z_loop; mask_loop(10,3,:) = mask_10_z_loop
-    mask_loop(11,3,:) = mask_11_z_loop; mask_loop(12,3,:) = mask_12_z_loop
-    mask_loop(13,3,:) = mask_13_z_loop; mask_loop(14,3,:) = mask_14_z_loop
-    mask_loop(15,3,:) = mask_15_z_loop; mask_loop(16,3,:) = mask_16_z_loop
-    mask_loop(17,3,:) = mask_17_z_loop; mask_loop(18,3,:) = mask_18_z_loop
-    mask_loop(19,3,:) = mask_19_z_loop; mask_loop(20,3,:) = mask_20_z_loop
+    DO mid = 1, masks
+       do_mask     (mid,:) = data_output_masks(mid,:)
+       do_mask_user(mid,:) = data_output_masks_user(mid,:)
+       mask      (mid,1,:) = mask_x(mid,:) 
+       mask      (mid,2,:) = mask_y(mid,:)
+       mask      (mid,3,:) = mask_z(mid,:) 
+       mask_loop (mid,1,:) = mask_x_loop(mid,:) 
+       mask_loop (mid,2,:) = mask_y_loop(mid,:) 
+       mask_loop (mid,3,:) = mask_z_loop(mid,:) 
+    ENDDO
+    
     mask_i = -1; mask_j = -1; mask_k = -1
+    
 !
 !-- Global arrays are required by define_netcdf_header.
Index: palm/trunk/SOURCE/modules.f90
===================================================================
--- palm/trunk/SOURCE/modules.f90	(revision 552)
+++ palm/trunk/SOURCE/modules.f90	(revision 553)
@@ -5,5 +5,5 @@
 ! Current revisions:
 ! -----------------
-! 
+! parameters for masked output are replaced by arrays
 !
 ! Former revisions:
@@ -346,25 +346,6 @@
     CHARACTER (LEN=11), DIMENSION(100) ::  data_output = ' ',    &
                                            data_output_user = ' ', doav = ' '
-    CHARACTER (LEN=10), DIMENSION(100) ::  &
-         data_output_mask_01      = ' ', data_output_mask_02      = ' ',  &
-         data_output_mask_03      = ' ', data_output_mask_04      = ' ',  &
-         data_output_mask_05      = ' ', data_output_mask_06      = ' ',  &
-         data_output_mask_07      = ' ', data_output_mask_08      = ' ',  &
-         data_output_mask_09      = ' ', data_output_mask_10      = ' ',  &
-         data_output_mask_11      = ' ', data_output_mask_12      = ' ',  &
-         data_output_mask_13      = ' ', data_output_mask_14      = ' ',  &
-         data_output_mask_15      = ' ', data_output_mask_16      = ' ',  &
-         data_output_mask_17      = ' ', data_output_mask_18      = ' ',  &
-         data_output_mask_19      = ' ', data_output_mask_20      = ' ',  &
-         data_output_mask_user_01 = ' ', data_output_mask_user_02 = ' ',  &
-         data_output_mask_user_03 = ' ', data_output_mask_user_04 = ' ',  &
-         data_output_mask_user_05 = ' ', data_output_mask_user_06 = ' ',  &
-         data_output_mask_user_07 = ' ', data_output_mask_user_08 = ' ',  &
-         data_output_mask_user_09 = ' ', data_output_mask_user_10 = ' ',  &
-         data_output_mask_user_11 = ' ', data_output_mask_user_12 = ' ',  &
-         data_output_mask_user_13 = ' ', data_output_mask_user_14 = ' ',  &
-         data_output_mask_user_15 = ' ', data_output_mask_user_16 = ' ',  &
-         data_output_mask_user_17 = ' ', data_output_mask_user_18 = ' ',  &
-         data_output_mask_user_19 = ' ', data_output_mask_user_20 = ' '
+    CHARACTER (LEN=10), DIMENSION(max_masks,100) ::  &
+         data_output_masks = ' ', data_output_masks_user = ' '
 
     CHARACTER (LEN=10), DIMENSION(300) ::  data_output_pr = ' '
@@ -565,81 +546,46 @@
     REAL, DIMENSION(:), ALLOCATABLE ::  dp_smooth_factor
 
-    REAL, DIMENSION(mask_xyz_dimension) ::  &
-       mask_01_x = -1.0, mask_02_x = -1.0, mask_03_x = -1.0, mask_04_x = -1.0, &
-       mask_05_x = -1.0, mask_06_x = -1.0, mask_07_x = -1.0, mask_08_x = -1.0, &
-       mask_09_x = -1.0, mask_10_x = -1.0, mask_11_x = -1.0, mask_12_x = -1.0, &
-       mask_13_x = -1.0, mask_14_x = -1.0, mask_15_x = -1.0, mask_16_x = -1.0, &
-       mask_17_x = -1.0, mask_18_x = -1.0, mask_19_x = -1.0, mask_20_x = -1.0, &
-       mask_01_y = -1.0, mask_02_y = -1.0, mask_03_y = -1.0, mask_04_y = -1.0, &
-       mask_05_y = -1.0, mask_06_y = -1.0, mask_07_y = -1.0, mask_08_y = -1.0, &
-       mask_09_y = -1.0, mask_10_y = -1.0, mask_11_y = -1.0, mask_12_y = -1.0, &
-       mask_13_y = -1.0, mask_14_y = -1.0, mask_15_y = -1.0, mask_16_y = -1.0, &
-       mask_17_y = -1.0, mask_18_y = -1.0, mask_19_y = -1.0, mask_20_y = -1.0, &
-       mask_01_z = -1.0, mask_02_z = -1.0, mask_03_z = -1.0, mask_04_z = -1.0, &
-       mask_05_z = -1.0, mask_06_z = -1.0, mask_07_z = -1.0, mask_08_z = -1.0, &
-       mask_09_z = -1.0, mask_10_z = -1.0, mask_11_z = -1.0, mask_12_z = -1.0, &
-       mask_13_z = -1.0, mask_14_z = -1.0, mask_15_z = -1.0, mask_16_z = -1.0, &
-       mask_17_z = -1.0, mask_18_z = -1.0, mask_19_z = -1.0, mask_20_z = -1.0
-    REAL, DIMENSION(3) ::  &
-       mask_01_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_02_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_03_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_04_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_05_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_06_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_07_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_08_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_09_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_10_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_11_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_12_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_13_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_14_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_15_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_16_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_17_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_18_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_19_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_20_x_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_01_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_02_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_03_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_04_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_05_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_06_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_07_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_08_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_09_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_10_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_11_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_12_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_13_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_14_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_15_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_16_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_17_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_18_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_19_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_20_y_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_01_z_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_02_z_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_03_z_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_04_z_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_05_z_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_06_z_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_07_z_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_08_z_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_09_z_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_10_z_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_11_z_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_12_z_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_13_z_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_14_z_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_15_z_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_16_z_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_17_z_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_18_z_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_19_z_loop = (/ -1.0, -1.0, 0.0 /), &
-       mask_20_z_loop = (/ -1.0, -1.0, 0.0 /)
+    REAL, DIMENSION(max_masks,mask_xyz_dimension) ::       &
+        mask_x = -1.0, mask_y = -1.0, mask_z = -1.0
+    REAL, DIMENSION(max_masks,3) ::                        &
+        mask_x_loop = (/ -1.0, -1.0, -1.0, -1.0, -1.0,     &  
+	                 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			  0.0,  0.0,  0.0,  0.0,  0.0,     &
+	                  0.0,  0.0,  0.0,  0.0,  0.0,     &
+			  0.0,  0.0,  0.0,  0.0,  0.0,     &
+			  0.0,  0.0,  0.0,  0.0,  0.0  /), &
+	
+	mask_y_loop = (/ -1.0, -1.0, -1.0, -1.0, -1.0,     &  
+	                 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			  0.0,  0.0,  0.0,  0.0,  0.0,     &
+	                  0.0,  0.0,  0.0,  0.0,  0.0,     &
+			  0.0,  0.0,  0.0,  0.0,  0.0,     &
+			  0.0,  0.0,  0.0,  0.0,  0.0  /), &
+			   
+	mask_z_loop = (/ -1.0, -1.0, -1.0, -1.0, -1.0,     &  
+	                 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			 -1.0, -1.0, -1.0, -1.0, -1.0,     &
+			  0.0,  0.0,  0.0,  0.0,  0.0,     &
+	                  0.0,  0.0,  0.0,  0.0,  0.0,     &
+			  0.0,  0.0,  0.0,  0.0,  0.0,     &
+			  0.0,  0.0,  0.0,  0.0,  0.0  /)  
+    
 !
 !--    internal mask arrays ("mask,dimension,selection")
Index: palm/trunk/SOURCE/parin.f90
===================================================================
--- palm/trunk/SOURCE/parin.f90	(revision 552)
+++ palm/trunk/SOURCE/parin.f90	(revision 553)
@@ -4,5 +4,5 @@
 ! Current revisions:
 ! -----------------
-! 
+! parameters for masked output are replaced by arrays
 !
 ! Former revisions:
@@ -117,5 +117,5 @@
              lad_vertical_gradient, lad_vertical_gradient_level, &
              leaf_surface_concentration, long_filter_factor, &
-             loop_optimization, mixing_length_1d, masks, momentum_advec, &
+             loop_optimization, mixing_length_1d, momentum_advec, &
              netcdf_precision, nsor_ini, nx, ny, &
              nz, ocean, omega, outflow_damping_width, overshoot_limit_e, &
@@ -151,14 +151,5 @@
                        cross_ts_uymin, cross_xtext, cycle_mg, data_output, &
                        data_output_format, &
-                       data_output_mask_01, data_output_mask_02, &
-                       data_output_mask_03, data_output_mask_04, &
-                       data_output_mask_05, data_output_mask_06, &
-                       data_output_mask_07, data_output_mask_08, &
-                       data_output_mask_09, data_output_mask_10, &
-                       data_output_mask_11, data_output_mask_12, &
-                       data_output_mask_13, data_output_mask_14, &
-                       data_output_mask_15, data_output_mask_16, &
-                       data_output_mask_17, data_output_mask_18, &
-                       data_output_mask_19, data_output_mask_20, &
+                       data_output_masks, &
                        data_output_pr, &
                        data_output_2d_on_each_pe, disturbance_amplitude, &
@@ -172,40 +163,6 @@
                        end_time, force_print_header, &
                        mask_scale_x, mask_scale_y, mask_scale_z, &
-                       mask_01_x, mask_02_x, mask_03_x, mask_04_x, &
-                       mask_05_x, mask_06_x, mask_07_x, mask_08_x, &
-                       mask_09_x, mask_10_x, mask_11_x, mask_12_x, &
-                       mask_13_x, mask_14_x, mask_15_x, mask_16_x, &
-                       mask_17_x, mask_18_x, mask_19_x, mask_20_x, &
-                       mask_01_y, mask_02_y, mask_03_y, mask_04_y, &
-                       mask_05_y, mask_06_y, mask_07_y, mask_08_y, &
-                       mask_09_y, mask_10_y, mask_11_y, mask_12_y, &
-                       mask_13_y, mask_14_y, mask_15_y, mask_16_y, &
-                       mask_17_y, mask_18_y, mask_19_y, mask_20_y, &
-                       mask_01_z, mask_02_z, mask_03_z, mask_04_z, &
-                       mask_05_z, mask_06_z, mask_07_z, mask_08_z, &
-                       mask_09_z, mask_10_z, mask_11_z, mask_12_z, &
-                       mask_13_z, mask_14_z, mask_15_z, mask_16_z, &
-                       mask_17_z, mask_18_z, mask_19_z, mask_20_z, &
-                       mask_01_x_loop, mask_02_x_loop, mask_03_x_loop, &
-                       mask_04_x_loop, mask_05_x_loop, mask_06_x_loop, &
-                       mask_07_x_loop, mask_08_x_loop, mask_09_x_loop, &
-                       mask_10_x_loop, mask_11_x_loop, mask_12_x_loop, &
-                       mask_13_x_loop, mask_14_x_loop, mask_15_x_loop, &
-                       mask_16_x_loop, mask_17_x_loop, mask_18_x_loop, &
-                       mask_19_x_loop, mask_20_x_loop, &
-                       mask_01_y_loop, mask_02_y_loop, mask_03_y_loop, &
-                       mask_04_y_loop, mask_05_y_loop, mask_06_y_loop, &
-                       mask_07_y_loop, mask_08_y_loop, mask_09_y_loop, &
-                       mask_10_y_loop, mask_11_y_loop, mask_12_y_loop, &
-                       mask_13_y_loop, mask_14_y_loop, mask_15_y_loop, &
-                       mask_16_y_loop, mask_17_y_loop, mask_18_y_loop, &
-                       mask_19_y_loop, mask_20_y_loop, &
-                       mask_01_z_loop, mask_02_z_loop, mask_03_z_loop, &
-                       mask_04_z_loop, mask_05_z_loop, mask_06_z_loop, &
-                       mask_07_z_loop, mask_08_z_loop, mask_09_z_loop, &
-                       mask_10_z_loop, mask_11_z_loop, mask_12_z_loop, &
-                       mask_13_z_loop, mask_14_z_loop, mask_15_z_loop, &
-                       mask_16_z_loop, mask_17_z_loop, mask_18_z_loop, &
-                       mask_19_z_loop, mask_20_z_loop, &
+                       mask_x, mask_y, mask_z,  &
+                       mask_x_loop, mask_y_loop, mask_z_loop, &
                        mg_cycles, mg_switch_to_pe0_level, &
                        netcdf_data_format, &
Index: palm/trunk/SOURCE/user_parin.f90
===================================================================
--- palm/trunk/SOURCE/user_parin.f90	(revision 552)
+++ palm/trunk/SOURCE/user_parin.f90	(revision 553)
@@ -4,5 +4,5 @@
 ! Current revisions:
 ! -----------------
-! 
+! data_output_mask_user_* replaced by array data_output_masks_user
 !
 ! Former revisions:
@@ -39,14 +39,5 @@
 
     NAMELIST /userpar/  data_output_pr_user, data_output_user, region,  &
-                        data_output_mask_user_01, data_output_mask_user_02, &
-                        data_output_mask_user_03, data_output_mask_user_04, &
-                        data_output_mask_user_05, data_output_mask_user_06, &
-                        data_output_mask_user_07, data_output_mask_user_08, &
-                        data_output_mask_user_09, data_output_mask_user_10, &
-                        data_output_mask_user_11, data_output_mask_user_12, &
-                        data_output_mask_user_13, data_output_mask_user_14, &
-                        data_output_mask_user_15, data_output_mask_user_16, &
-                        data_output_mask_user_17, data_output_mask_user_18, &
-                        data_output_mask_user_19, data_output_mask_user_20
+                        data_output_masks_user
 
 !
