Index: palm/trunk/DOC/tec/technical_documentation.html
===================================================================
--- palm/trunk/DOC/tec/technical_documentation.html	(revision 197)
+++ palm/trunk/DOC/tec/technical_documentation.html	(revision 198)
@@ -2351,5 +2351,88 @@
 prognostic_equations</td></tr><tr><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">E</td><td style="vertical-align: top;">Bugfix: summation of<span style="font-family: Courier New,Courier,monospace;"> sums_l_l </span>(<span style="font-family: Courier New,Courier,monospace;">diffusivities</span>).<br><br>Several bugfixes in the ocean part: Initial density<span style="font-family: Courier New,Courier,monospace;"> rho </span>is calculated (<span style="font-family: Courier New,Courier,monospace;">init_ocean</span>). Error in initializing<span style="font-family: Courier New,Courier,monospace;"> u_init </span>and<span style="font-family: Courier New,Courier,monospace;"> v_init </span>removed (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>). Calculation of density flux now starts from<span style="font-family: Courier New,Courier,monospace;"> nzb+1 </span>(<span style="font-family: Courier New,Courier,monospace;">production_e</span>).<br><br>Bugfix: <span style="font-family: Courier New,Courier,monospace;">pleft</span>/<span style="font-family: Courier New,Courier,monospace;">pright</span> changed to <span style="font-family: Courier New,Courier,monospace;">pnorth</span>/<span style="font-family: Courier New,Courier,monospace;">psouth</span> in sendrecv of particle tail numbers along y, small bugfixes in the SGS part (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>).<br><br>Bugfix:<span style="font-family: Courier New,Courier,monospace;"> model_string </span>needed a default value (<span style="font-family: Courier New,Courier,monospace;">combine_plot_fields</span>).<br><br>Bugfix: wavenumber calculation for even<span style="font-family: Courier New,Courier,monospace;"> nx </span>in routines<span style="font-family: Courier New,Courier,monospace;"> maketri </span>(<span style="font-family: Courier New,Courier,monospace;">poisfft</span>).<br><br>Bugfix: assignment of fluxes at walls.<br><br>Bugfix: absolute value of<span style="font-family: Courier New,Courier,monospace;"> f </span>must be used when calculating the Blackadar mixing length (<span style="font-family: Courier New,Courier,monospace;">init_1d_model</span>).</td><td style="vertical-align: top;">advec_particles,
 check_parameters, combine_plot_fields, diffusion_s, diffusivities,
-init_ocean, init_1d_model, poisfft, production_e</td></tr><tr><td style="vertical-align: top;">16/09/08</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.5</td><td style="vertical-align: top;">N</td><td style="vertical-align: top;"></td><td style="vertical-align: top;"></td></tr><tr><td align="undefined" valign="undefined">&nbsp;</td><td align="undefined" valign="undefined">&nbsp;</td><td align="undefined" valign="undefined">&nbsp;</td><td style="vertical-align: top;">C</td><td style="vertical-align: top;"></td><td align="undefined" valign="undefined"></td></tr><tr><td align="undefined" valign="undefined">&nbsp;</td><td align="undefined" valign="undefined">&nbsp;</td><td align="undefined" valign="undefined">&nbsp;</td><td style="vertical-align: top;">E</td><td style="vertical-align: top;"></td><td style="vertical-align: top;"></td></tr>
+init_ocean, init_1d_model, poisfft, production_e</td></tr><tr><td style="vertical-align: top;">16/09/08</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.5</td><td style="vertical-align: top;">N</td><td style="vertical-align: top;">The
+turbulence recycling method of Kataoka and Mizuno (2002) has been
+implemented. A precursor run must be used to initialize the main run
+with a quasi-stationary turbulence state (see <span style="font-family: monospace;">initializing_actions</span> = <span style="font-style: italic;">'read_data_for_recycling'</span>).
+The precursor run may have a smaller domain than the main run, which
+helps to spare CPU time. New parameters for steering the recycling
+method are <span style="font-family: monospace;">turbulent_inflow</span>, <span style="font-family: monospace;">recycling_width</span>, <span style="font-family: monospace;">inflow_damping_height</span>, and <span style="font-family: monospace;">inflow_damping_width</span>.<br><br>One
+new feature related with the recycling method is that within a job
+chain, the number of processors or the virtual processor grid is now
+allowed to&nbsp;change. Such a change modifies the domain decomposition
+and requires that the PEs within a restart job have to know from which
+file(s) they have to read their restart data (up to now there has been
+a one to one relation between the restart files and the PEs, i.e.
+processor 4 only reads data from file "_0004"). Informations about the
+subdomain index limits of all restart files are now written by PE0 (to
+file "_0000"). Also, all variables containing global values which do
+not depend on the subdomain are written to file "_0000" only. In a
+restart run, this file is opened by all PEs to read these global values
+and the informations about the subdomain index limits. Depending on
+these index limit information it is decided which restart file(s)
+has/have to be opened in order to read the complete data necessary for
+the current subdomain.<br><br>The plant canopy model has been extended
+to the prognostic equation for&nbsp;the potential temperature
+and&nbsp;passive scalar, as proposed by Watanabe (2004) and Shaw and
+Schumann (1992).<br><br>Particle advection adjusted for ocean runs.<br><br>Grid
+coarsening (multigrid method) is limited to subdomains (i.e. no
+collection of data on PE0 for further coarsening beyond subdomains) if <span style="font-family: monospace;">mg_switch_to_pe0_level</span> = <span style="font-style: italic;">-1</span>.<br><br>User-defined spectra are allowed.<br><br>Output of <span style="font-family: monospace;">q*2</span> profile and <span style="font-family: monospace;">w"u"0</span>, <span style="font-family: monospace;">w"v"0</span>, <span style="font-family: monospace;">w"q"0</span> timeseries are added.<br><br>Argument<span style="font-family: monospace;"> -c </span>introduced to most of the subjob calls, which allows the user to choose his own job catalog by setting the environment variable<span style="font-family: monospace;"> job_catalog </span>in the configuration file (default is <span style="font-family: monospace;">~/job_queue</span>).<br><br>This release contains first adjustments for the SGI-ICE system at HLRN.<br><br><span style="font-weight: bold;">Attention: On the SGI-ICE,&nbsp;runs are still aborted or are hanging for unknown reasons from time to time.</span></td><td style="vertical-align: top;">advec_particles, calc_spectra, check_open, check_parameters, data_output_spectra,
+header, init_3d_model, init_particles, init_pegrid, modules, netcdf, parin, particle_boundary_conds, plant_canopy_model,
+
+prognostic_equations, read_3d_binary, read_var_list, time_integration,
+user_interface, write_3d_binary, write_var_list<br><br><span style="font-weight: bold;">new:</span><br>inflow_turbulence<br></td></tr><tr><td align="undefined" valign="undefined">&nbsp;</td><td align="undefined" valign="undefined">&nbsp;</td><td align="undefined" valign="undefined">&nbsp;</td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Allow 100 spectra levels instead of 10 for consistency with define_netcdf_header. (<span style="font-family: monospace;">calc_spectra</span>, <span style="font-family: monospace;">data_output_spectra</span>, <span style="font-family: monospace;">header</span>, <span style="font-family: monospace;">modules</span>, <span style="font-family: monospace;">netcdf</span>)<br><br>Modification
+of the integrated version of the profile function for momentum for
+unstable stratification; more consistent flux definitions. (<span style="font-family: monospace;">wall_fluxes</span>, <span style="font-family: monospace;">production_e</span>)<br><br>Strict grid matching along z is not needed for mg-solver. (<span style="font-family: monospace;">check_parameters</span>)<br><br>Sharp
+lateral edges of the plant canopy introduced. Before, smaller values of
+the leaf area density than those applied in the prognostic equation for
+the vertical velocity were used in the prognostic equations for u and v
+at the boundaries of the plant_canopy. (<span style="font-family: monospace;">init_3d_model</span>)<br><br>-<span style="font-family: monospace;">cross_ts_*</span> (<span style="font-family: monospace;">read_var_list</span>, <span style="font-family: monospace;">write_var_list</span>), -<span style="font-family: monospace;">dopr_time_count</span>, old profil parameters (<span style="font-family: monospace;">cross_...</span>, <span style="font-family: monospace;">dopr_crossindex</span>,&nbsp;<span style="font-family: monospace;">profile_***</span>) removed. (<span style="font-family: monospace;">read_3d_binary</span>, <span style="font-family: monospace;">write_3d_binary</span>)<br><br>Reading/writing of <span style="font-family: monospace;">hom_sum</span>, <span style="font-family: monospace;">volume_flow_area</span>, <span style="font-family: monospace;">volume_flow_initial</span> is moved from <span style="font-family: monospace;">..._3d_binary</span> to <span style="font-family: monospace;">..._var_list</span>, <span style="font-family: monospace;">hom_sum</span> is now allocated in <span style="font-family: monospace;">parin</span> instead of <span style="font-family: monospace;">init_3d_model</span>.<br><br><span style="font-family: monospace;">npex</span>, <span style="font-family: monospace;">npey</span> moved from inipar- to d3par-namelist. (<span style="font-family: monospace;">parin</span>)<br><br><span style="font-family: monospace;">myid_char14</span>
+removed. This was used, since on IBM machines the process rank could be
+changed when the final communicator was defined, in order to save the
+preliminary processor-id for opening the binary output file for
+restarts (unit 14), because otherwise a mismatch occured when reading
+this file in a restart job. This mechanism is not needed any more since
+the subdomain informations are now contained in file _0000. (<span style="font-family: monospace;">parin</span>, <span style="font-family: monospace;">check_open</span>)<br><br>Transpositions for the 2D domain decomposition have been optimized by using <span style="font-family: monospace;">f_inv</span>
+as an automatic array instead of providing the memory by a dummy
+argument. This spares one copy loop per transposition. Order of indices
+in the 3D loops in some of the transpose routines have been rearranged
+for better cache utilization. Both have been suggested by Roland
+Richter (SGI) as part of the HLRN-II benchmark process. (<span style="font-family: monospace;">transpose</span>)<br><br>Workaround for getting information about the coupling mode, because the<span style="font-family: monospace;"> -env </span>Option (mpiexec) is not available on many systems (SGI, lcxt4 at&nbsp; Bergen Center for Computational Science). (<span style="font-family: monospace;">palm</span>, <span style="font-family: monospace;">mrun</span>)<br><br>Format adjustments in order to allow CPU# &gt; 999. (<span style="font-family: monospace;">cpu_statistics</span>)<br><br>Provisional solution for run_control_1d output: add 'CALL check_open( 15 )'. (<span style="font-family: monospace;">init_1d_model</span>)<br><br>Subdomain data are read into temporary arrays <span style="font-family: monospace;">pf_tmp</span> / <span style="font-family: monospace;">pf3d_tmp</span> in order to avoid INTEL compiler warnings about (automatic) creation of temporary arrays. (<span style="font-family: monospace;">combine_plot_fields</span>)<br><br>Time coordinate t=0 stored on netcdf-file only if an output is required for this time for at least one of the profiles.<br><br>Output of spline timeseries <span style="font-family: monospace;">splptx</span> etc. removed. (<span style="font-family: monospace;">flow_statistics</span>)</td><td style="vertical-align: top;">advec_s_ups,
+advec_u_ups, advec_v_ups, advec_w_ups, calc_spectra, check_open,
+check_parameters, cpu_statistics, data_output_profiles,
+data_output_spectra, flow_statistics, header, init_1d_model,
+init_3d_model, modules, netcdf, palm, parin, poisfft, production_e,
+read_3d_binary, read_var_list, transpose, wall_fluxes, write_3d_binary, write_var_list<br></td></tr><tr><td align="undefined" valign="undefined">&nbsp;</td><td align="undefined" valign="undefined">&nbsp;</td><td align="undefined" valign="undefined">&nbsp;</td><td style="vertical-align: top;">E</td><td style="vertical-align: top;">Bugfix:
+modification of the calculation of the vertical turbulent momentum
+fluxes u'w' and v'w'; the first usws that is computed corresponds to
+-u'w'/u* and not as priorily assumed to (-u'w')**0.5, the first vsws
+that is computed corresponds to -v'w'/u* and not as priorily assumed to
+(-v'w')**0.5. Therefore, the intermediate result for usws has to be
+multiplied by -u* instead by itself in order to get u'w'. Accordingly,
+the intermediate result for vsws has to be multiplied by -u* instead by
+itself in order to get v'w'. As u* is calculated for the position of a
+scalar an additional interpolation of u* to the position of u and v,
+respectively, is necessary. As u* is not determined for the ghost
+points on each PE, an additional exchange of information from
+neighbouring PEs is necessary. Tests showed that this
+decreases u* by some 10% and increases TKE and momentum fluxes by some
+10% because friction was underestimated and momentum fluxes were
+wrongly calculated due to the bug. (<span style="font-family: monospace;">prandtl_fluxes</span>, <span style="font-family: monospace;">wall_fluxes</span>)<br><br><span style="font-weight: bold;">Attention:</span> <br><span style="font-weight: bold;">
+The content of the testresult-file (run control output) has changed
+due to the above bugfix.</span><br><br>Bugfix:
+change definition of us_wall from 1D to 2D.&nbsp;<br><br>Bugfix: exponent = 1.0 / ( ny + 1.0 ) in calc_spectra_x. (<span style="font-family: monospace;">calc_spectra</span>)<br><br>Bugfix: calculation of horizontal fluxes at vertical walls. (<span style="font-family: monospace;">diffusion_s</span>)<br><br>Bugfix: zero assignments to tendency arrays in case of restarts. (<span style="font-family: monospace;">init_3d_model</span>)<br><br>Bugfix: wrong grid level used in multigrid solver, routine <span style="font-family: monospace;">restrict</span>. (<span style="font-family: monospace;">poismg</span>)<br><br>Bugfix:
+Construction of vertical profiles when 10 gradients have been specified
+in the parameter list (ug, vg, pt, q, sa, lad). (<span style="font-family: monospace;">check_parameters</span>)<br><br>Bugfix: Determination of the number of vertical levels for which spectra shall be output - 100 replaced by 10. (<span style="font-family: monospace;">netcdf</span>)<br><br>Bugfix:
+Lengths of those characters that are needed for the output of the
+characteristic levels of potential temperature, the geostrophic wind,
+scalar concentration, humidity and leaf area density to the header file
+has been increased. (<span style="font-family: monospace;">header</span>)<br><br>Bugfix: add ')' in <span style="font-family: monospace;">user_read_restart_data</span>.<br><br>The leaf area density is explicitly set to its surface value at k=0. (<span style="font-family: monospace;">check_parameters</span>)<br><br>Introduce<span style="font-family: monospace;"> prefix_chr</span> to ensure unique dvrp_file path.<br><br>small bugfixes for <span style="font-family: monospace;">user_interface</span> sample code (comments):<br>- initialize ustvst with 0.0 as it is now computed only until <span style="font-family: monospace;">nxr</span> and <span style="font-family: monospace;">nyn</span><br>- two ALLOCATE statements moved from <span style="font-family: monospace;">user_read_restart_data</span> back to <span style="font-family: monospace;">user_init</span><br>- remove 'READ (13) u2_av' statement in <span style="font-family: monospace;">user_read_restart_data</span><br><br>Bugfix: remove IF statement in <span style="font-family: monospace;">plant_canopy_model_ij</span>.<br><br>Bugfix: divide sums(k,8) (e) and sums(k,34) (e*) by ngp_2dh_s_inner(k,sr) (like other scalars). (<span style="font-family: monospace;">flow_statistics</span>)<br><br>Bugfix: <span style="font-family: monospace;">dopr_time_count</span>
+was written on the binary file, which caused that NetCDF files newly
+created by restart files (no append of existing files!) contained
+uneccessary time levels. (<span style="font-family: monospace;">read_3d_binary</span>, <span style="font-family: monospace;">write_3d_binary</span>)<br><br>Bugfix: extra '*' removed in user_statistics sample code. (<span style="font-family: monospace;">user_interface</span>)<br><br>Bugfix: a stop command was missing in some cases of the parallel branch. (<span style="font-family: monospace;">local_stop</span>)<br><br>Bugfix in volume flow control for non-cyclic boundary conditions. (<span style="font-family: monospace;">pres</span>)<br><br>Bugfix: misplaced #endif directives. (<span style="font-family: monospace;">combine_plot_fields</span>)<br><br>Data are collected from PE0 in an ordered sequence which seems to avoid hanging of processes on SGI-ICE. (<span style="font-family: monospace;">cpu_statistics</span>)</td><td style="vertical-align: top;">calc_spectra,
+check_parameters, cpu_statistics, diffusion_s, flow_statistics,
+init_3d_model, init_dvrp, local_stop, plant_canopy_model, poismg,
+prandtl_fluxes, pres, read_3d_binary, user_interface, wall_fluxes,
+write_3d_binary<br></td></tr>
 </tbody>
 </table>&nbsp;<b><blink>Attention:</blink></b>
Index: palm/trunk/SOURCE/CURRENT_MODIFICATIONS
===================================================================
--- palm/trunk/SOURCE/CURRENT_MODIFICATIONS	(revision 197)
+++ palm/trunk/SOURCE/CURRENT_MODIFICATIONS	(revision 198)
@@ -1,163 +1,13 @@
 New:
 ---
-The plant canopy model has been extended to the prognostic equation for   
-the potential temperature and the prognostic equation for the concentration 
-of a passive scalar, as proposed by Watanabe (2004) and Shaw and Schumann 
-(1992). (prognostic_equations, plant_canopy_model)  
 
-Within a job chain, the number of processors or the virtual processor grid
-is now allowed to be changed. This modifies the domain decomposition and
-requires that the PEs within a restart job have to know from which file(s) they
-have to read their restart data (up to now there is a one to one relation
-between the restart files and the PEs, i.e. processor 4 only reads data
-from file "_0004"). Informations about the subdomain index limits of all
-restart files are now written by PE0 (to file "_0000"). Also, all variables
-containing global values which do not depend on the subdomain are written
-to file "_0000" only. In a restart run, this file is opened by all PEs to
-read these global values and the informations about the subdomain index
-limits. Depending on these index limit information it is decided which
-restart file(s) has/have to be opened in order to read the complete data
-necessary for the current subdomain.
-
-The turbulence recycling method of Kataoka and Mizuno (2002) has been
-implemented. A pre-run with smaller model domain can be used to initialize
-the real run, in order to spare the CPU time needed to get the turbulence
-into a quasi-stationary state (see initializing_actions =
-'read_data_for_recycling').
-
-Particle advection adjusted for ocean runs.
-
-Multigrid levels are limited by subdomains if mg_switch_to_pe0_level = -1.
-
-User-defined spectra.
-
-Output of q*2 profile added.
-
-Argument -c introduced to most of the subjob calls, which allows the user to
-choose his own job catalog by setting job_catalog in the configuration file
-(default is ~/job_queue). Workaround for mpixec with -env option.
-Adjustments for lcxt4 (Bergen Center for Computational Science) (mrun, subjob)
-
-advec_particles, calc_spectra, check_open, check_parameters, data_output_spectra, header, init_particles, init_pegrid, init_3d_model, modules, netcdf, parin, particle_boundary_conds, plant_canopy_model, prognostic_equations, read_var_list, read_3d_binary, time_integration, user_interface, write_var_list, write_3d_binary
-
-New: inflow_turbulence
 
 
 Changed:
 -------
-Allow 100 spectra levels instead of 10 for consistency with 
-define_netcdf_header (calc_spectra, data_output_spectra, header, modules, 
-netcdf)
 
-Modification of the integrated version of the profile function for momentum 
-for unstable stratification; more consistent flux definitions. (wall_fluxes, 
-production_e)
-
-Strict grid matching along z is not needed for mg-solver. (check_parameters)
-
-Sharp lateral edges of the plant canopy introduced. Before smaller 
-values of the leaf area density than those applied in the prognostic 
-equation for the vertical velocity were used in the prognostic equations 
-for u and v at the boundaries of the plant_canopy. (init_3d_model)  
-
--cross_ts_* (read_var_list, write_var_list)
--dopr_time_count, old profil parameters (cross_..., dopr_crossindex,
- profile_***) removed (read_3d_binary, write_3d_binary)
-
-hom_sum moved, volume_flow_area, volume_flow_initial moved from ..._3d_binary
-to ..._var_list, hom_sum is now allocated in parin instead of init_3d_model
-
-npex, npey moved from inipar- to d3par-namelist (parin)
-
-myid_char14 removed. This was used since on IBM machines the process rank may
-be changed when the final communicator is defined, in order to save the
-preliminary processor-id for opening the binary output file for restarts
-(unit 14), because otherwise a mismatch occured when reading this file in a
-restart job. This mechanism is not needed any more since the subdomain
-informations are now contained in file _0000. (parin, check_open)
-
-Transpositions for the 2D domain decomposition have been optimized by using
-f_inv as an automatic array instead of providing the memory by a dummy argument.
-This spares one copy loop per transposition. Order of indices in the 3D loops
-in some of the transpose routines have been rearranged for better cache utilization.
-Both have been suggested by Roland Richter (SGI) as part of the
-HLRN-II benchmark process. (transpose)
-
-Workaround for getting information about the coupling mode. (palm)
-
-Format adjustments in order to allow CPU# > 999. (cpu_statistics)
-
-Provisional solution for run_control_1d output: add 'CALL check_open( 15 )'
-(init_1d_model)
-
-Subdomain data are read into temporary arrays pf_tmp/pf3d_tmp in order to
-avoid INTEL compiler warnings about (automatic) creation of temporary arrays
-(combine_plot_fields)
-
-Time coordinate t=0 stored on netcdf-file only if an output is required for
-this time for at least one of the profiles
-
-Spline timeseries splptx etc. removed, timeseries w'u', w'v', w'q' (k=0) added
-(flow_statistics)
-
-advec_s_ups, advec_u_ups, advec_v_ups, advec_w_ups, calc_spectra, check_open, check_parameters, cpu_statistics, data_output_profiles, data_output_spectra, flow_statistics, header, init_1d_model, init_3d_model, modules, netcdf, palm, parin, poisfft, production_e, read_var_list, read_3d_binary, transpose, wall_fluxes, write_var_list, write_3d_binary
 
 
 Errors:
 ------
-Bugfix: exponent = 1.0 / ( ny + 1.0 ) in calc_spectra_x (calc_spectra)
 
-Bugfix: change definition of us_wall from 1D to 2D. Tests showed that this decreases 
-u* by some 10% and increases TKE and momentum fluxes by some 10% because friction 
-was underestimated and momentum fluxes were wrongly calculated due to the bug. 
-(prandtl_fluxes, wall_fluxes)
-
-Bugfix: calculation of horizontal fluxes at vertical walls (diffusion_s)
-
-Bugfix: zero assignments to tendency arrays in case of restarts (init_3d_model)
-
-Bugfix: wrong grid level used in multigrid solver, routine restrict (poismg)
-
-Bugfix: Construction of vertical profiles when 10 gradients have been 
-specified in the parameter list (ug, vg, pt, q, sa, lad). (check_parameters)
-
-Bugfix: Determination of the number of vertical levels for which spectra 
-shall be output - 100 replaced by 10. (netcdf)
-
-Bugfix: Lengths of those characters that are needed for the output of the 
-characteristic levels of potential temperature, the geostrophic wind, 
-scalar concentration, humidity and leaf area density to the header file 
-has been increased. (header) 
-
-Bugfix: add ')' in user_read_restart_data.
-
-The leaf area density is explicitly set to its surface value at k=0. (check_parameters)
-
-Introduce prefix_chr to ensure unique dvrp_file path.
-
-small bugfixes for user_interface sample code (comments):
-- initialize ustvst with 0.0 as it is now computed only until nxr and nyn
-- two ALLOCATE statements moved from user_read_restart_data back to user_init
-- remove 'READ (13) u2_av' statement in user_read_restart_data
-
-Bugfix: remove IF statement in plant_canopy_model_ij (plant_canopy_model)
-
-Bugfix: divide sums(k,8) (e) and sums(k,34) (e*) by ngp_2dh_s_inner(k,sr) 
-(like other scalars) (flow_statistics)
-
-Bugfix: dopr_time_count was written on the binary file, which caused that
-NetCDF files newly created by restart files (no append of existing files!)
-contained uneccessary time levels. (read_3d_binary, write_3d_binary)
-
-Bugfix: extra '*' removed in user_statistics sample code (user_interface)
-
-Bugfix: a stop command was missing in some cases of the parallel branch (local_stop)
-
-Bugfix in volume flow control for non-cyclic boundary conditions (pres)
-
-Bugfix: misplaced #endif directives (combine_plot_fields)
-
-data are collected from PE0 in an ordered sequence which seems to avoid hanging of processes on SGI-ICE (cpu_statistics)
-
-calc_spectra, check_parameters, cpu_statistics, diffusion_s, flow_statistics, init_dvrp, init_3d_model, local_stop, plant_canopy_model, poismg, prandtl_fluxes, pres, read_3d_binary, user_interface, wall_fluxes, write_3d_binary
-
Index: palm/trunk/SOURCE/advec_particles.f90
===================================================================
--- palm/trunk/SOURCE/advec_particles.f90	(revision 197)
+++ palm/trunk/SOURCE/advec_particles.f90	(revision 198)
@@ -4,6 +4,4 @@
 ! Actual revisions:
 ! -----------------
-! Bottom boundary condition and vertical index calculations adjusted for
-! ocean runs.
 ! TEST: PRINT statements on unit 9 (commented out)
 !
@@ -11,4 +9,8 @@
 ! -----------------
 ! $Id$
+!
+! 150 2008-02-29 08:19:58Z raasch
+! Bottom boundary condition and vertical index calculations adjusted for
+! ocean runs.
 !
 ! 119 2007-10-17 10:27:13Z raasch
Index: palm/trunk/SOURCE/advec_s_ups.f90
===================================================================
--- palm/trunk/SOURCE/advec_s_ups.f90	(revision 197)
+++ palm/trunk/SOURCE/advec_s_ups.f90	(revision 198)
@@ -4,9 +4,14 @@
 ! Actual revisions:
 ! -----------------
-! Arguments removed from transpose routines
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 164 2008-05-15 08:46:15Z raasch
+! Arguments removed from transpose routines
+!
+! February 2007
 ! RCS Log replace by Id keyword, revision history cleaned up
 !
Index: palm/trunk/SOURCE/advec_u_ups.f90
===================================================================
--- palm/trunk/SOURCE/advec_u_ups.f90	(revision 197)
+++ palm/trunk/SOURCE/advec_u_ups.f90	(revision 198)
@@ -4,9 +4,14 @@
 ! Actual revisions:
 ! -----------------
-! Arguments removed from transpose routines
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 164 2008-05-15 08:46:15Z raasch
+! Arguments removed from transpose routines
+!
+! February 2007
 ! RCS Log replace by Id keyword, revision history cleaned up
 !
Index: palm/trunk/SOURCE/advec_v_ups.f90
===================================================================
--- palm/trunk/SOURCE/advec_v_ups.f90	(revision 197)
+++ palm/trunk/SOURCE/advec_v_ups.f90	(revision 198)
@@ -4,9 +4,14 @@
 ! Actual revisions:
 ! -----------------
-! Arguments removed from transpose routines
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 164 2008-05-15 08:46:15Z raasch
+! Arguments removed from transpose routines
+!
+! February 2007
 ! RCS Log replace by Id keyword, revision history cleaned up
 !
Index: palm/trunk/SOURCE/advec_w_ups.f90
===================================================================
--- palm/trunk/SOURCE/advec_w_ups.f90	(revision 197)
+++ palm/trunk/SOURCE/advec_w_ups.f90	(revision 198)
@@ -4,9 +4,14 @@
 ! Actual revisions:
 ! -----------------
-! Arguments removed from transpose routines
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 164 2008-05-15 08:46:15Z raasch
+! Arguments removed from transpose routines
+!
+! February 2007
 ! RCS Log replace by Id keyword, revision history cleaned up
 !
Index: palm/trunk/SOURCE/calc_spectra.f90
===================================================================
--- palm/trunk/SOURCE/calc_spectra.f90	(revision 197)
+++ palm/trunk/SOURCE/calc_spectra.f90	(revision 198)
@@ -4,4 +4,11 @@
 ! Actual revisions:
 ! -----------------
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 192 2008-08-27 16:51:49Z letzel
 ! bugfix in calc_spectra_x: exponent = 1.0 / ( ny + 1.0 )
 ! allow 100 spectra levels instead of 10 for consistency with
@@ -9,7 +16,5 @@
 ! user-defined spectra, arguments removed from transpose routines
 !
-! Former revisions:
-! -----------------
-! $Id$
+! February 2007
 ! RCS Log replace by Id keyword, revision history cleaned up
 !
Index: palm/trunk/SOURCE/check_open.f90
===================================================================
--- palm/trunk/SOURCE/check_open.f90	(revision 197)
+++ palm/trunk/SOURCE/check_open.f90	(revision 198)
@@ -4,11 +4,14 @@
 ! Actual revisions:
 ! -----------------
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 146 2008-01-17 13:08:34Z raasch
 ! First opening of unit 13 openes file _0000 on all PEs (parallel version)
 ! because only this file contains the global variables,
 ! myid_char_14 removed
-!
-! Former revisions:
-! -----------------
-! $Id$
 !
 ! 120 2007-10-17 11:54:43Z raasch
Index: palm/trunk/SOURCE/check_parameters.f90
===================================================================
--- palm/trunk/SOURCE/check_parameters.f90	(revision 197)
+++ palm/trunk/SOURCE/check_parameters.f90	(revision 198)
@@ -4,4 +4,11 @@
 ! Actual revisions:
 ! -----------------
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 197 2008-09-16 15:29:03Z raasch
 ! Bug fix: Construction of vertical profiles when 10 gradients have been 
 ! specified in the parameter list (ug, vg, pt, q, sa, lad)
@@ -12,8 +19,4 @@
 ! check of turbulent_inflow and calculation of recycling_plane.
 ! q*2 profile added
-!
-! Former revisions:
-! -----------------
-! $Id$
 !
 ! 138 2007-11-28 10:03:58Z letzel
Index: palm/trunk/SOURCE/cpu_statistics.f90
===================================================================
--- palm/trunk/SOURCE/cpu_statistics.f90	(revision 197)
+++ palm/trunk/SOURCE/cpu_statistics.f90	(revision 198)
@@ -4,11 +4,14 @@
 ! Actual revisions:
 ! -----------------
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 197 2008-09-16 15:29:03Z raasch
 ! Format adjustments in order to allow CPU# > 999,
 ! data are collected from PE0 in an ordered sequence which seems to avoid
 ! hanging of processes on SGI-ICE
-!
-! Former revisions:
-! -----------------
-! $Id$
 !
 ! 82 2007-04-16 15:40:52Z raasch
Index: palm/trunk/SOURCE/data_output_profiles.f90
===================================================================
--- palm/trunk/SOURCE/data_output_profiles.f90	(revision 197)
+++ palm/trunk/SOURCE/data_output_profiles.f90	(revision 198)
@@ -4,10 +4,15 @@
 ! Actual revisions:
 ! -----------------
-! Time coordinate t=0 stored on netcdf-file only if an output is required for
-! this time for at least one of the profiles
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 197 2008-09-16 15:29:03Z raasch
+! Time coordinate t=0 stored on netcdf-file only if an output is required for
+! this time for at least one of the profiles
+!
+! February 2007
 ! RCS Log replace by Id keyword, revision history cleaned up
 !
Index: palm/trunk/SOURCE/data_output_spectra.f90
===================================================================
--- palm/trunk/SOURCE/data_output_spectra.f90	(revision 197)
+++ palm/trunk/SOURCE/data_output_spectra.f90	(revision 198)
@@ -4,11 +4,15 @@
 ! Actual revisions:
 ! -----------------
-! allow 100 spectra levels instead of 10 for consistency with
-! define_netcdf_header
-! user-defined spectra
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 189 2008-08-13 17:09:26Z letzel
+! allow 100 spectra levels instead of 10 for consistency with
+! define_netcdf_header, +user-defined spectra
+!
+! February 2007
 ! RCS Log replace by Id keyword, revision history cleaned up
 !
Index: palm/trunk/SOURCE/diffusion_s.f90
===================================================================
--- palm/trunk/SOURCE/diffusion_s.f90	(revision 197)
+++ palm/trunk/SOURCE/diffusion_s.f90	(revision 198)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! bugfix: calculation of fluxes at vertical surfaces
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 183 2008-08-04 15:39:12Z letzel
+! bugfix: calculation of fluxes at vertical surfaces
 !
 ! 129 2007-10-30 12:12:24Z letzel
Index: palm/trunk/SOURCE/flow_statistics.f90
===================================================================
--- palm/trunk/SOURCE/flow_statistics.f90	(revision 197)
+++ palm/trunk/SOURCE/flow_statistics.f90	(revision 198)
@@ -4,12 +4,15 @@
 ! Actual revisions:
 ! -----------------
-! Spline timeseries splptx etc. removed, timeseries w'u', w'v', w'q' (k=0)
-! added
-! Bugfix: divide sums(k,8) (e) and sums(k,34) (e*) by ngp_2dh_s_inner(k,sr) 
-! (like other scalars)
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 197 2008-09-16 15:29:03Z raasch
+! Spline timeseries splptx etc. removed, timeseries w'u', w'v', w'q' (k=0)
+! added,
+! bugfix: divide sums(k,8) (e) and sums(k,34) (e*) by ngp_2dh_s_inner(k,sr) 
+! (like other scalars)
 !
 ! 133 2007-11-20 10:10:53Z letzel
Index: palm/trunk/SOURCE/header.f90
===================================================================
--- palm/trunk/SOURCE/header.f90	(revision 197)
+++ palm/trunk/SOURCE/header.f90	(revision 198)
@@ -4,14 +4,16 @@
 ! Actual revisions:
 ! -----------------
-! TEST: mg_switch_to_pe0_level = -1!!!!!!!!!
-! allow 100 spectra levels instead of 10 for consistency with 
-! define_netcdf_header
-! Bugfix in the output of the characteristic levels of potential temperature, 
-! geostrophic wind, scalar concentration, humidity and leaf area density
-! Output of turbulence recycling informations
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 197 2008-09-16 15:29:03Z raasch
+! allow 100 spectra levels instead of 10 for consistency with 
+! define_netcdf_header,
+! bugfix in the output of the characteristic levels of potential temperature, 
+! geostrophic wind, scalar concentration, humidity and leaf area density,
+! output of turbulence recycling informations
 !
 ! 138 2007-11-28 10:03:58Z letzel
Index: palm/trunk/SOURCE/inflow_turbulence.f90
===================================================================
--- palm/trunk/SOURCE/inflow_turbulence.f90	(revision 197)
+++ palm/trunk/SOURCE/inflow_turbulence.f90	(revision 198)
@@ -10,4 +10,5 @@
 ! $Id$
 !
+! Initial version (2008/03/07)
 !
 ! Description:
Index: palm/trunk/SOURCE/init_1d_model.f90
===================================================================
--- palm/trunk/SOURCE/init_1d_model.f90	(revision 197)
+++ palm/trunk/SOURCE/init_1d_model.f90	(revision 198)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! provisional solution for run_control_1d output: add 'CALL check_open( 15 )'
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 184 2008-08-04 15:53:39Z letzel
+! provisional solution for run_control_1d output: add 'CALL check_open( 15 )'
 !
 ! 135 2007-11-22 12:24:23Z raasch
Index: palm/trunk/SOURCE/init_3d_model.f90
===================================================================
--- palm/trunk/SOURCE/init_3d_model.f90	(revision 197)
+++ palm/trunk/SOURCE/init_3d_model.f90	(revision 198)
@@ -7,15 +7,18 @@
 ! Actual revisions:
 ! -----------------
-! Bugfix: zero assignments to tendency arrays in case of restarts
-! Further amendments to and modifications in the initialisation of the plant 
-! canopy model 
-! Allocation of hom_sum moved to parin, initialization of spectrum_x|y directly
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 181 2008-07-30 07:07:47Z raasch
+! bugfix: zero assignments to tendency arrays in case of restarts,
+! further extensions and modifications in the initialisation of the plant 
+! canopy model,
+! allocation of hom_sum moved to parin, initialization of spectrum_x|y directly
 ! after allocating theses arrays,
 ! read data for recycling added as new initialization option,
 ! dummy allocation for diss
-!
-! Former revisions:
-! -----------------
-! $Id$
 !
 ! 138 2007-11-28 10:03:58Z letzel
Index: palm/trunk/SOURCE/init_dvrp.f90
===================================================================
--- palm/trunk/SOURCE/init_dvrp.f90	(revision 197)
+++ palm/trunk/SOURCE/init_dvrp.f90	(revision 198)
@@ -5,5 +5,4 @@
 ! -----------------
 ! 
-! introduce prefix_chr to ensure unique dvrp_file path
 ! TEST: print* statements
 ! ToDo: checking of mode_dvrp for legal values is not correct
@@ -12,4 +11,7 @@
 ! -----------------
 ! $Id$
+!
+! 155 2008-03-28 10:56:30Z letzel
+! introduce prefix_chr to ensure unique dvrp_file path
 !
 ! 130 2007-11-13 14:08:40Z letzel
Index: palm/trunk/SOURCE/init_particles.f90
===================================================================
--- palm/trunk/SOURCE/init_particles.f90	(revision 197)
+++ palm/trunk/SOURCE/init_particles.f90	(revision 198)
@@ -4,9 +4,13 @@
 ! Actual revisions:
 ! -----------------
-! Set offset_ocean_* needed for calculating vertical indices within ocean runs
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 150 2008-02-29 08:19:58Z raasch
+! Setting offset_ocean_* needed for calculating vertical indices within ocean
+! runs
 !
 ! 117 2007-10-11 03:27:59Z raasch
Index: palm/trunk/SOURCE/init_pegrid.f90
===================================================================
--- palm/trunk/SOURCE/init_pegrid.f90	(revision 197)
+++ palm/trunk/SOURCE/init_pegrid.f90	(revision 198)
@@ -5,13 +5,15 @@
 ! -----------------
 ! ATTENTION: nnz_x undefined problem still has to be solved!!!!!!!!
+! TEST OUTPUT (TO BE REMOVED) logging mpi2 ierr values
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 197 2008-09-16 15:29:03Z raasch
 ! multigrid levels are limited by subdomains if mg_switch_to_pe0_level = -1,
 ! nz is used instead nnz for calculating mg-levels
 ! Collect on PE0 horizontal index bounds from all other PEs,
 ! broadcast the id of the inflow PE (using the respective communicator)
-! TEST OUTPUT (TO BE REMOVED) logging mpi2 ierr values
-!
-! Former revisions:
-! -----------------
-! $Id$
 !
 ! 114 2007-10-10 00:03:15Z raasch
Index: palm/trunk/SOURCE/local_stop.f90
===================================================================
--- palm/trunk/SOURCE/local_stop.f90	(revision 197)
+++ palm/trunk/SOURCE/local_stop.f90	(revision 198)
@@ -4,5 +4,5 @@
 ! Actual revisions:
 ! -----------------
-! Bugfix: a stop command was missing in some cases of the parallel branch
+! 
 ! 
 !
@@ -10,4 +10,7 @@
 ! -----------------
 ! $Id$
+!
+! 147 2008-02-01 12:41:46Z raasch
+! Bugfix: a stop command was missing in some cases of the parallel branch
 !
 ! 108 2007-08-24 15:10:38Z letzel
Index: palm/trunk/SOURCE/modules.f90
===================================================================
--- palm/trunk/SOURCE/modules.f90	(revision 197)
+++ palm/trunk/SOURCE/modules.f90	(revision 198)
@@ -5,7 +5,13 @@
 ! Actual revisions:
 ! -----------------
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 197 2008-09-16 15:29:03Z raasch
 ! allow 100 spectra levels instead of 10 for consistency with 
-! define_netcdf_header
-! +canopy_heat_flux, cthf, lai,
+! define_netcdf_header, +canopy_heat_flux, cthf, lai,
 ! +leaf_surface_concentration, scalar_exchange_coefficient, sec, sls
 ! +hor_index_bounds, hor_index_bounds_previous_run, id_inflow, id_recycling,
@@ -13,8 +19,4 @@
 ! ny_on_file, offset_ocean_*, recycling_plane, recycling_width, turbulent_inflow
 ! -myid_char_14
-!
-! Former revisions:
-! -----------------
-! $Id$
 !
 ! 138 2007-11-28 10:03:58Z letzel
Index: palm/trunk/SOURCE/netcdf.f90
===================================================================
--- palm/trunk/SOURCE/netcdf.f90	(revision 197)
+++ palm/trunk/SOURCE/netcdf.f90	(revision 198)
@@ -7,12 +7,14 @@
 ! Current revisions:
 ! ------------------
-! consistently allow 100 spectra levels instead of 10
-! bug fix in the determination of the number of output heights for spectra
-!
-! user-defined spectra
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 189 2008-08-13 17:09:26Z letzel
+! consistently allow 100 spectra levels instead of 10
+! bug fix in the determination of the number of output heights for spectra,
+! +user-defined spectra
 !
 ! 97 2007-06-21 08:23:15Z raasch
Index: palm/trunk/SOURCE/palm.f90
===================================================================
--- palm/trunk/SOURCE/palm.f90	(revision 197)
+++ palm/trunk/SOURCE/palm.f90	(revision 198)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! Workaround for getting information about the coupling mode
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 197 2008-09-16 15:29:03Z raasch
+! Workaround for getting information about the coupling mode
 !
 ! 108 2007-08-24 15:10:38Z letzel
Index: palm/trunk/SOURCE/parin.f90
===================================================================
--- palm/trunk/SOURCE/parin.f90	(revision 197)
+++ palm/trunk/SOURCE/parin.f90	(revision 198)
@@ -4,16 +4,17 @@
 ! Actual revisions:
 ! -----------------
-! +cthf
-! +leaf_surface_concentration, scalar_exchange_coefficient
-! +inflow_damping_height, inflow_damping_width, recycling_width,
-! turbulent_inflow in inipar
-! -skip_time_dosp in d3par,
-! Allocation of hom_sum moved from init_3d_model to here,
-! npex, npey moved from inipar to d3par, setting of myid_char_14 removed,
-! lad is allways allocated
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 197 2008-09-16 15:29:03Z raasch
+! +cthf,leaf_surface_concentration, scalar_exchange_coefficient
+! +inflow_damping_height, inflow_damping_width, recycling_width,
+! turbulent_inflow in inipar, -skip_time_dosp in d3par,
+! allocation of hom_sum moved from init_3d_model to here,
+! npex, npey moved from inipar to d3par, setting of myid_char_14 removed,
+! lad is allways allocated
 !
 ! 138 2007-11-28 10:03:58Z letzel
Index: palm/trunk/SOURCE/particle_boundary_conds.f90
===================================================================
--- palm/trunk/SOURCE/particle_boundary_conds.f90	(revision 197)
+++ palm/trunk/SOURCE/particle_boundary_conds.f90	(revision 198)
@@ -4,9 +4,13 @@
 ! Actual revisions:
 ! -----------------
-! Vertical index calculations adjusted for ocean runs.
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 150 2008-02-29 08:19:58Z raasch
+! Vertical index calculations adjusted for ocean runs.
+!
 ! Initial version (2007/03/09)
 !
Index: palm/trunk/SOURCE/plant_canopy_model.f90
===================================================================
--- palm/trunk/SOURCE/plant_canopy_model.f90	(revision 197)
+++ palm/trunk/SOURCE/plant_canopy_model.f90	(revision 198)
@@ -4,9 +4,5 @@
 ! Actual revisions:
 ! -----------------
-! Consideration of heat sources within the forest canopy which represent the
-! rate of heat input into the air from the forest leaves
-!
-! Evaluation of sinks and sources for scalar concentration due to canopy 
-! elements
+! 
 !
 ! Bugfix: remove IF statement in plant_canopy_model_ij
@@ -15,4 +11,9 @@
 ! -----------------
 ! $Id$
+!
+! 153 2008-03-19 09:41:30Z steinfeld
+! heat sources within the forest canopy are added, which represent the
+! rate of heat input into the air from the forest leaves, evaluation of sinks
+! and sources for scalar concentration due to canopy elements
 !
 ! 138 2007-11-28 10:03:58Z letzel
Index: palm/trunk/SOURCE/poisfft.f90
===================================================================
--- palm/trunk/SOURCE/poisfft.f90	(revision 197)
+++ palm/trunk/SOURCE/poisfft.f90	(revision 198)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! Arguments removed from transpose routines
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 164 2008-05-15 08:46:15Z raasch
+! Arguments removed from transpose routines
 !
 ! 128 2007-10-26 13:11:14Z raasch
Index: palm/trunk/SOURCE/poismg.f90
===================================================================
--- palm/trunk/SOURCE/poismg.f90	(revision 197)
+++ palm/trunk/SOURCE/poismg.f90	(revision 198)
@@ -8,9 +8,12 @@
 ! Actual revisions:
 ! -----------------
-! Bugfix: grid_level+1 has to be used in restrict for flags-array
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 181 2008-07-30 07:07:47Z raasch
+! Bugfix: grid_level+1 has to be used in restrict for flags-array
 !
 ! 114 2007-10-10 00:03:15Z raasch
Index: palm/trunk/SOURCE/prandtl_fluxes.f90
===================================================================
--- palm/trunk/SOURCE/prandtl_fluxes.f90	(revision 197)
+++ palm/trunk/SOURCE/prandtl_fluxes.f90	(revision 198)
@@ -4,24 +4,16 @@
 ! Actual revisions:
 ! -----------------
-! Bugfix: change definition of us_wall from 1D to 2D:
-!         Modification of the evaluation of the vertical turbulent momentum 
-!         fluxes u'w' and v'w'; the first usws that is computed corresponds 
-!         to -u'w'/u* and not as priorily assumed to (-u'w')**0.5, the first 
-!         vsws that is computed corresponds to -v'w'/u* and not as priorily
-!         assumed to (-v'w')**0.5. Therefore, the intermediate result for 
-!         usws has to be multiplied by -u* instead by itself in order to 
-!         get u'w'. Accordingly, the intermediate result for vsws has to be 
-!         multiplied by -u* instead by itself in order to get v'w'. As u* 
-!         is calculated for the position of a scalar an additional 
-!         interpolation of u* to the position of u and v, respectively, 
-!         is necessary. As u* is not determined for the ghost points on 
-!         each PE, an additional exchange of information from neighbouring 
-!         PEs is necessary.!  
-! Change: Modification of the integrated version of the profile function for 
-!         momentum for unstable stratification
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 187 2008-08-06 16:25:09Z letzel
+! Bugfix: modification of the calculation of the vertical turbulent momentum
+! fluxes u'w' and v'w'
+! Bugfix: change definition of us_wall from 1D to 2D
+! Change: modification of the integrated version of the profile function for 
+! momentum for unstable stratification (does not effect results)
 !
 ! 108 2007-08-24 15:10:38Z letzel
Index: palm/trunk/SOURCE/pres.f90
===================================================================
--- palm/trunk/SOURCE/pres.f90	(revision 197)
+++ palm/trunk/SOURCE/pres.f90	(revision 198)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! Bugfix in volume flow control for non-cyclic boundary conditions
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 151 2008-03-07 13:42:18Z raasch
+! Bugfix in volume flow control for non-cyclic boundary conditions
 !
 ! 106 2007-08-16 14:30:26Z raasch
Index: palm/trunk/SOURCE/production_e.f90
===================================================================
--- palm/trunk/SOURCE/production_e.f90	(revision 197)
+++ palm/trunk/SOURCE/production_e.f90	(revision 198)
@@ -4,10 +4,13 @@
 ! Actual revisions:
 ! -----------------
-! Change: add 'minus' sign to fluxes obtained from subroutine wall_fluxes_e for
-!         consistency with subroutine wall_fluxes
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 187 2008-08-06 16:25:09Z letzel
+! Change: add 'minus' sign to fluxes obtained from subroutine wall_fluxes_e for
+! consistency with subroutine wall_fluxes
 !
 ! 124 2007-10-19 15:47:46Z raasch
Index: palm/trunk/SOURCE/prognostic_equations.f90
===================================================================
--- palm/trunk/SOURCE/prognostic_equations.f90	(revision 197)
+++ palm/trunk/SOURCE/prognostic_equations.f90	(revision 198)
@@ -4,12 +4,13 @@
 ! Actual revisions:
 ! -----------------
-! add call of plant_canopy_model in the prognostic equation for
-! the potential temperature
-! add call of plant_canopy_model in the prognostic equation for
-! the concentration of a passive scalar
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 153 2008-03-19 09:41:30Z steinfeld
+! add call of plant_canopy_model in the prognostic equation for
+! the potential temperature and for the passive scalar
 !
 ! 138 2007-11-28 10:03:58Z letzel
Index: palm/trunk/SOURCE/read_3d_binary.f90
===================================================================
--- palm/trunk/SOURCE/read_3d_binary.f90	(revision 197)
+++ palm/trunk/SOURCE/read_3d_binary.f90	(revision 198)
@@ -4,4 +4,11 @@
 ! Actual revisions:
 ! -----------------
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 150 2008-02-29 08:19:58Z raasch
 ! Files from which restart data are to be read are determined and subsequently
 ! opened. The total domain on the restart file is allowed to be smaller than
@@ -12,8 +19,4 @@
 ! reading of old profil parameters (cross_..., dopr_crossindex, profile_***)
 ! removed, initialization of spectrum_x|y removed
-!
-! Former revisions:
-! -----------------
-! $Id$
 !
 ! 102 2007-07-27 09:09:17Z raasch
Index: palm/trunk/SOURCE/read_var_list.f90
===================================================================
--- palm/trunk/SOURCE/read_var_list.f90	(revision 197)
+++ palm/trunk/SOURCE/read_var_list.f90	(revision 198)
@@ -4,4 +4,11 @@
 ! Actual revisions:
 ! -----------------
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 173 2008-05-23 20:39:38Z raasch
 ! +cthf, leaf_surface_concentration, scalar_exchange_coefficient
 ! +numprocs_previous_run, hor_index_bounds_previous_run, inflow_damping_factor,
@@ -12,8 +19,4 @@
 ! read_3d_binary to here,
 ! routines read_parts_of_var_list and skip_var_list added at the end
-!
-! Former revisions:
-! -----------------
-! $Id$
 !
 ! 138 2007-11-28 10:03:58Z letzel
Index: palm/trunk/SOURCE/time_integration.f90
===================================================================
--- palm/trunk/SOURCE/time_integration.f90	(revision 197)
+++ palm/trunk/SOURCE/time_integration.f90	(revision 198)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! inflow turbulence is imposed by calling new routine inflow_turbulence
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 151 2008-03-07 13:42:18Z raasch
+! inflow turbulence is imposed by calling new routine inflow_turbulence
 !
 ! 108 2007-08-24 15:10:38Z letzel
Index: palm/trunk/SOURCE/transpose.f90
===================================================================
--- palm/trunk/SOURCE/transpose.f90	(revision 197)
+++ palm/trunk/SOURCE/transpose.f90	(revision 198)
@@ -4,11 +4,16 @@
 ! Actual revisions:
 ! -----------------
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 164 2008-05-15 08:46:15Z raasch
 ! f_inv changed from subroutine argument to automatic array in order to do
 ! re-ordering from f_in to f_inv in one step, one array work is needed instead
 ! of work1 and work2
 !
-! Former revisions:
-! -----------------
-! $Id$
+! February 2007
 ! RCS Log replace by Id keyword, revision history cleaned up
 !
Index: palm/trunk/SOURCE/user_interface.f90
===================================================================
--- palm/trunk/SOURCE/user_interface.f90	(revision 197)
+++ palm/trunk/SOURCE/user_interface.f90	(revision 198)
@@ -4,4 +4,11 @@
 ! Actual revisions:
 ! -----------------
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 159 2008-04-09 08:06:03Z raasch
 ! bugfix in user_read_restart_data
 ! small bugfixes for sample code (comments):
@@ -10,8 +17,4 @@
 ! - remove 'READ (13) u2_av' statement in user_read_restart_data
 ! +routines user_read_restart_data, user_spectra
-!
-! Former revisions:
-! -----------------
-! $Id$
 !
 ! 138 2007-11-28 10:03:58Z letzel
Index: palm/trunk/SOURCE/wall_fluxes.f90
===================================================================
--- palm/trunk/SOURCE/wall_fluxes.f90	(revision 197)
+++ palm/trunk/SOURCE/wall_fluxes.f90	(revision 198)
@@ -3,24 +3,21 @@
 ! Actual revisions:
 ! -----------------
-! Bugfix: change definition of us_wall from 1D to 2D:
-!         Modification of the evaluation of the vertical turbulent momentum 
-!         fluxes u'w' and v'w'; the first usws that is computed corresponds 
-!         to -u'w'/u* and not as priorily assumed to (-u'w')**0.5, the first 
-!         vsws that is computed corresponds to -v'w'/u* and not as priorily
-!         assumed to (-v'w')**0.5. Therefore, the intermediate result for 
-!         usws has to be multiplied by -u* instead by itself in order to 
-!         get u'w'. Accordingly, the intermediate result for vsws has to be 
-!         multiplied by -u* instead by itself in order to get v'w'.
-!         This requires the calculation of us_wall (and vel_total, u_i, v_i, ws)
-!         also in wall_fluxes_e.
-! Bugfix: storage of rifs to rifs_wall in wall_fluxes_e removed
-! Change: add 'minus' sign to fluxes produced by subroutine wall_fluxes_e for 
-!         consistency with subroutine wall_fluxes
-! Change: Modification of the integrated version of the profile function for 
-!         momentum for unstable stratification
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 187 2008-08-06 16:25:09Z letzel
+! Bugfix: Modification of the evaluation of the vertical turbulent momentum 
+! fluxes u'w' and v'w (see prandtl_fluxes), this requires the calculation of
+! us_wall (and vel_total, u_i, v_i, ws) also in wall_fluxes_e.
+! Bugfix: change definition of us_wall from 1D to 2D
+! Bugfix: storage of rifs to rifs_wall in wall_fluxes_e removed
+! Change: add 'minus' sign to fluxes produced by subroutine wall_fluxes_e for 
+! consistency with subroutine wall_fluxes
+! Change: Modification of the integrated version of the profile function for 
+! momentum for unstable stratification
+!
 ! Initial version (2007/03/07)
 !
Index: palm/trunk/SOURCE/write_3d_binary.f90
===================================================================
--- palm/trunk/SOURCE/write_3d_binary.f90	(revision 197)
+++ palm/trunk/SOURCE/write_3d_binary.f90	(revision 198)
@@ -4,13 +4,15 @@
 ! Actual revisions:
 ! -----------------
-! write_var_list is only called for pe0
-! -dopr_time_count,
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 145 2008-01-09 08:17:38Z raasch
+! write_var_list is only called for pe0, -dopr_time_count,
 ! hom_sum moved, volume_flow_area, volume_flow_initial moved to write_var_list,
 ! writing of old profil parameters (cross_..., dopr_crossindex, profile_***)
 ! removed
-!
-! Former revisions:
-! -----------------
-! $Id$
 !
 ! 102 2007-07-27 09:09:17Z raasch
Index: palm/trunk/SOURCE/write_var_list.f90
===================================================================
--- palm/trunk/SOURCE/write_var_list.f90	(revision 197)
+++ palm/trunk/SOURCE/write_var_list.f90	(revision 198)
@@ -4,4 +4,11 @@
 ! Actual revisions:
 ! -----------------
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 153 2008-03-19 09:41:30Z steinfeld
 ! +cthf, leaf_surface_concentration, scalar_exchange_coefficient
 ! +numprocs, hor_index_bounds, inflow_damping_height, inflow_damping_width,
@@ -10,8 +17,4 @@
 ! hom_sum, volume_flow_area, volume_flow_initial moved from write_3d_binary
 ! to here
-!
-! Former revisions:
-! -----------------
-! $Id$
 !
 ! 138 2007-11-28 10:03:58Z letzel
