Index: /palm/trunk/DOC/tec/technical_documentation.html
===================================================================
--- /palm/trunk/DOC/tec/technical_documentation.html	(revision 109)
+++ /palm/trunk/DOC/tec/technical_documentation.html	(revision 110)
@@ -2314,5 +2314,27 @@
 density is calculated from the equation of state for seawater after
 each timestep, using the algorithm proposed by Jackett et al. (2006, J.
-Atmos. Oceanic Technol., <span style="font-weight: bold;">23</span>, 1709-1728).<br>So far, only the initial hydrostatic pressure is entered into this equation.</li><li>z=0 (sea surface) is assumed at the model top (vertical grid index <span style="font-family: Courier New,Courier,monospace;">k=nzt</span> on the w-grid), with negative values of z indicating the depth.</li><li>Initial profiles are constructed (e.g. from <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#pt_vertical_gradient">pt_vertical_gradient</a> / <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#pt_vertical_gradient_level">pt_vertical_gradient_level</a>) starting from the sea surface, using surface values&nbsp;given by <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#pt_surface">pt_surface</a>, <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#sa_surface">sa_surface</a>, <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#ug_surface">ug_surface</a>, and <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#vg_surface">vg_surface</a>.</li><li>Zero salinity flux is used as default boundary condition at the bottom of the sea.</li><li>If switched on, random perturbations are by default imposed to the upper model domain from zu(nzt*2/3) to zu(nzt-3).</li></ul>Relevant new inipar-parameters to be exclusively used for steering ocean runs are <span style="font-family: Courier New,Courier,monospace;">bc_sa_t</span>, <span style="font-family: Courier New,Courier,monospace;">bottom_salinityflux</span>, <span style="font-family: Courier New,Courier,monospace;">sa_surface</span>, <span style="font-family: Courier New,Courier,monospace;">sa_vertical_gradient</span>, <span style="font-family: Courier New,Courier,monospace;">sa_vertical_gradient_level</span>, and <span style="font-family: Courier New,Courier,monospace;">top_salinityflux</span>.<br><br>Salinity (<span style="font-family: Courier New,Courier,monospace;">sa</span>) and potential density (<span style="font-family: Courier New,Courier,monospace;">rho</span>) are included as new 2d/3d output quantities. Vertical profiles of salinity (<span style="font-family: Courier New,Courier,monospace;">sa</span>), salinity fluxes (<span style="font-family: Courier New,Courier,monospace;">w"sa"</span>, <span style="font-family: Courier New,Courier,monospace;">w*sa*</span>, <span style="font-family: Courier New,Courier,monospace;">wsa</span>), and potential density (<span style="font-family: Courier New,Courier,monospace;">rho</span>) can also be output.<span style="font-family: Courier New,Courier,monospace;"></span></td><td style="vertical-align: top;">advec_s_bc, average_3d_data,<br>boundary_conds, buoyancy,<br>check_parameters,&nbsp;data_output_2d,<br>data_output_3d,&nbsp;diffusion_e,<br>diffusivities,&nbsp;flow_statistics,&nbsp;header,<br>init_grid,&nbsp;init_3d_model,&nbsp;modules,<br>netcdf,&nbsp;parin,&nbsp;production_e,<br>prognostic_equations,<br>read_3d_binary,&nbsp;&nbsp;read_var_list,<br>sum_up_3d_data,&nbsp;swap_timelevel,<br>time_integration,&nbsp;user_interface,<br>write_var_list,&nbsp;write_3d_binary<br><br><span style="font-weight: bold;">new:</span><br>eqn_state_seawater,&nbsp;init_ocean</td></tr><tr><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Inipar-parameter <span style="font-family: Courier New,Courier,monospace;">use_pt_reference</span> renamed <span style="font-family: Courier New,Courier,monospace;">use_reference.</span><br>Internal variable <span style="font-family: Courier New,Courier,monospace;">hydro_press</span> renamed <span style="font-family: Courier New,Courier,monospace;">hyp</span>, routine <span style="font-family: Courier New,Courier,monospace;">calc_mean_pt_profile</span> renamed <span style="font-family: Courier New,Courier,monospace;">calc_mean_profile</span>.<br><br>The format of the <span style="font-family: Courier New,Courier,monospace;">RUN_CONTROL</span> file has been adjusted for ocean runs.<br><br>Defaults of <span style="font-family: Courier New,Courier,monospace;">.._vertical_gradient_levels</span> changed from <span style="font-style: italic;">-1.0</span> to <span style="font-style: italic;">-9999999.9</span>.</td><td style="vertical-align: top;">advec_particles,&nbsp;buoyancy,&nbsp;<br>calc_liquid_water_content, <br>check_parameters,&nbsp;diffusion_e,<br>diffusivities,&nbsp;header,<br>init_cloud_physics,&nbsp;modules,<br>production_e,&nbsp;prognostic_equations,<br>run_control</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: height above topography instead of height above level k=0 is used for calculating the mixing length (<span style="font-family: Courier New,Courier,monospace;">diffusion_e</span> and <span style="font-family: Courier New,Courier,monospace;">diffusivities</span>).<br><br>Bugfix: error in boundary condition for TKE removed (<span style="font-family: Courier New,Courier,monospace;">advec_s_bc</span>).</td><td style="vertical-align: top;">advec_s_bc,&nbsp;diffusion_e,&nbsp;diffusivities</td></tr>
+Atmos. Oceanic Technol., <span style="font-weight: bold;">23</span>, 1709-1728).<br>So far, only the initial hydrostatic pressure is entered into this equation.</li><li>z=0 (sea surface) is assumed at the model top (vertical grid index <span style="font-family: Courier New,Courier,monospace;">k=nzt</span> on the w-grid), with negative values of z indicating the depth.</li><li>Initial profiles are constructed (e.g. from <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#pt_vertical_gradient">pt_vertical_gradient</a> / <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#pt_vertical_gradient_level">pt_vertical_gradient_level</a>) starting from the sea surface, using surface values&nbsp;given by <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#pt_surface">pt_surface</a>, <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#sa_surface">sa_surface</a>, <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#ug_surface">ug_surface</a>, and <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#vg_surface">vg_surface</a>.</li><li>Zero salinity flux is used as default boundary condition at the bottom of the sea.</li><li>If switched on, random perturbations are by default imposed to the upper model domain from zu(nzt*2/3) to zu(nzt-3).</li></ul>Relevant new inipar-parameters to be exclusively used for steering ocean runs are <span style="font-family: Courier New,Courier,monospace;">bc_sa_t</span>, <span style="font-family: Courier New,Courier,monospace;">bottom_salinityflux</span>, <span style="font-family: Courier New,Courier,monospace;">sa_surface</span>, <span style="font-family: Courier New,Courier,monospace;">sa_vertical_gradient</span>, <span style="font-family: Courier New,Courier,monospace;">sa_vertical_gradient_level</span>, and <span style="font-family: Courier New,Courier,monospace;">top_salinityflux</span>.<br><br>Salinity (<span style="font-family: Courier New,Courier,monospace;">sa</span>) and potential density (<span style="font-family: Courier New,Courier,monospace;">rho</span>) are included as new 2d/3d output quantities. Vertical profiles of salinity (<span style="font-family: Courier New,Courier,monospace;">sa</span>), salinity fluxes (<span style="font-family: Courier New,Courier,monospace;">w"sa"</span>, <span style="font-family: Courier New,Courier,monospace;">w*sa*</span>, <span style="font-family: Courier New,Courier,monospace;">wsa</span>), and potential density (<span style="font-family: Courier New,Courier,monospace;">rho</span>) can also be output.<span style="font-family: Courier New,Courier,monospace;"></span></td><td style="vertical-align: top;">advec_s_bc, average_3d_data,<br>boundary_conds, buoyancy,<br>check_parameters,&nbsp;data_output_2d,<br>data_output_3d,&nbsp;diffusion_e,<br>diffusivities,&nbsp;flow_statistics,&nbsp;header,<br>init_grid,&nbsp;init_3d_model,&nbsp;modules,<br>netcdf,&nbsp;parin,&nbsp;production_e,<br>prognostic_equations,<br>read_3d_binary,&nbsp;&nbsp;read_var_list,<br>sum_up_3d_data,&nbsp;swap_timelevel,<br>time_integration,&nbsp;user_interface,<br>write_var_list,&nbsp;write_3d_binary<br><br><span style="font-weight: bold;">new:</span><br>eqn_state_seawater,&nbsp;init_ocean</td></tr><tr><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Inipar-parameter <span style="font-family: Courier New,Courier,monospace;">use_pt_reference</span> renamed <span style="font-family: Courier New,Courier,monospace;">use_reference.</span><br>Internal variable <span style="font-family: Courier New,Courier,monospace;">hydro_press</span> renamed <span style="font-family: Courier New,Courier,monospace;">hyp</span>, routine <span style="font-family: Courier New,Courier,monospace;">calc_mean_pt_profile</span> renamed <span style="font-family: Courier New,Courier,monospace;">calc_mean_profile</span>.<br><br>The format of the <span style="font-family: Courier New,Courier,monospace;">RUN_CONTROL</span> file has been adjusted for ocean runs.<br><br>Defaults of <span style="font-family: Courier New,Courier,monospace;">.._vertical_gradient_levels</span> changed from <span style="font-style: italic;">-1.0</span> to <span style="font-style: italic;">-9999999.9</span>.</td><td style="vertical-align: top;">advec_particles,&nbsp;buoyancy,&nbsp;<br>calc_liquid_water_content, <br>check_parameters,&nbsp;diffusion_e,<br>diffusivities,&nbsp;header,<br>init_cloud_physics,&nbsp;modules,<br>production_e,&nbsp;prognostic_equations,<br>run_control</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: height above topography instead of height above level k=0 is used for calculating the mixing length (<span style="font-family: Courier New,Courier,monospace;">diffusion_e</span> and <span style="font-family: Courier New,Courier,monospace;">diffusivities</span>).<br><br>Bugfix: error in boundary condition for TKE removed (<span style="font-family: Courier New,Courier,monospace;">advec_s_bc</span>).</td><td style="vertical-align: top;">advec_s_bc,&nbsp;diffusion_e,&nbsp;diffusivities</td></tr><tr><td style="vertical-align: top;">05/10/07</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.4</td><td style="vertical-align: top;">N</td><td style="vertical-align: top;">This version allows runs for a coupled atmosphere-ocean LES.<br><br>Coupling frequency is controlled by new d3par-parameter <span style="font-family: Courier New,Courier,monospace;">dt_coupling</span>, the coupling mode ('<span style="font-style: italic;">atmosphere_to_ocean</span>' or '<span style="font-style: italic;">ocean_to_atmosphere</span>') for the respective processes is read from environment variable <span style="font-family: Courier New,Courier,monospace;">coupling_mode</span>, which is set by the mpiexec-command in <span style="font-weight: bold;">mrun</span>. Communication between the two models is handled by the intercommunicator <span style="font-family: Courier New,Courier,monospace;">comm_inter</span>. Local files opened by the ocean model get the additional suffic "<span style="font-family: Courier New,Courier,monospace;">_O</span>".
+Saturation at k=nzb is assumed for the atmosphere coupled to ocean.
+Usage of the coupled model is described in the new section 3.8.<br><br>A momentum flux can be set as top boundary condition using the new inipar parameters <span style="font-family: Courier New,Courier,monospace;">top_momentumflux_u|v</span>.<br><br>Non-cyclic boundary conditions can be used along all horizontal directions.<br><br>Quantities w*p* and w"e can be output as vertical profiles.<br><br>Initial profiles are reset to constant profiles in case that <span style="font-family: Courier New,Courier,monospace;">initializing_actions</span> /= '<span style="font-style: italic;">set_constant_profiles</span>'. (init_rankine)<br><br>Optionally calculate km and kh from initial TKE e_init.<br><br><span style="font-weight: bold;">Attention:</span> <br>
+The format of the testresult-file (run control output) has changed
+due to this modification.</td><td style="vertical-align: top;">boundary_conds,
+check_open, check_parameters, diffusion_u, diffusion_v,
+flow_statistics, header, init_pegrid, init_rankine, init_3d_model,
+local_stop, modules, palm, parin, prandtl_fluxes, pres,
+prognostic_equations, read_var_list, read_3d_binary, swap_timelevel,
+timestep, time_integration, write_var_list, write_3d_binary<br><br><span style="font-weight: bold;">new:</span><br>surface_coupler</td></tr><tr><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Remaining variables <span style="font-family: Courier New,Courier,monospace;">iran</span> changed to <span style="font-family: Courier New,Courier,monospace;">iran_part</span> (advec_particles, init_particles).<br><br>In case that the presure solver is not called for every Runge-Kutta substep (<span style="font-family: Courier New,Courier,monospace;">call_psolver_at_all_substeps</span> = <span style="font-style: italic;">.F.</span>),
+it is called after the first substep instead of the last. In that case,
+random perturbations are also added to the velocity field after the
+first substep.<br><br>Initialization of <span style="font-family: Courier New,Courier,monospace;">km</span>, <span style="font-family: Courier New,Courier,monospace;">kh</span> = 0.00001 for <span style="font-family: Courier New,Courier,monospace;">ocean</span> = <span style="font-style: italic;">.T.</span> (for <span style="font-family: Courier New,Courier,monospace;">ocean</span> = <span style="font-style: italic;">.F.</span> it remains 0.01).<br><br>Allow <span style="font-family: Courier New,Courier,monospace;">data_output_pr</span> = '<span style="font-style: italic;">q</span>', '<span style="font-style: italic;">wq</span>', '<span style="font-style: italic;">w"q"</span>', '<span style="font-style: italic;">w*q*</span>' for <span style="font-family: Courier New,Courier,monospace;">humidity</span> = <span style="font-style: italic;">.T.</span> (instead of <span style="font-family: Courier New,Courier,monospace;">cloud_physics</span> = <span style="font-style: italic;">.T.</span>).</td><td style="vertical-align: top;">advec_particles, check_for_restart, check_parameters, init_particles, init_3d_model, time_integration<br></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;">Bugs from code parts for non-cyclic boundary conditions are removed:<br>loops for u and v are starting from index <span style="font-family: Courier New,Courier,monospace;">nxlu</span>, <span style="font-family: Courier New,Courier,monospace;">nysv</span>,
+respectively. The radiation boundary condition is used for every
+Runge-Kutta substep. Velocity phase speeds for the radiation boundary
+conditions are calculated for the first Runge-Kutta substep only and
+reused for the further substeps. New arrays <span style="font-family: Courier New,Courier,monospace;">c_u</span>, <span style="font-family: Courier New,Courier,monospace;">c_v</span>, and <span style="font-family: Courier New,Courier,monospace;">c_w</span>
+are defined for this purpose. Several index errors are removed from the
+radiation boundary condition code parts. Upper bounds for calculating <span style="font-family: Courier New,Courier,monospace;">u_0</span> and <span style="font-family: Courier New,Courier,monospace;">v_0</span> (in production_e) are <span style="font-style: italic;">nxr+1</span> and <span style="font-style: italic;">nyn+1</span> because otherwise these values are not available in case of non-cyclic boundary conditions.<br><br>+<span style="font-family: Courier New,Courier,monospace;">dots_num_palm</span> in module user, +module netcdf_control in user_init (both in user_interface)<br><br>Bugfix: wrong sign removed from the buoyancy production term in the case <span style="font-family: Courier New,Courier,monospace;">use_reference</span> = <span style="font-style: italic;">.T.</span> (production_e)<br><br>Bugfix: Error message concerning output of particle concentration (pc) modified (check_parameters).<br><br>Bugfix: Rayleigh damping for ocean fixed.</td><td style="vertical-align: top;">advec_u_pw,
+advec_u_up, advec_v_pw, advec_v_up, boundary_conds, buoyancy,
+check_parameters, coriolis, diffusion_u, diffusion_v, header,
+init_pegrid, init_3d_model, modules, production_e,
+prognostic_equations, user_interface<br></td></tr>
 </tbody>
 </table>&nbsp;<b><blink>Attention:</blink></b>
Index: /palm/trunk/SOURCE/CURRENT_MODIFICATIONS
===================================================================
--- /palm/trunk/SOURCE/CURRENT_MODIFICATIONS	(revision 109)
+++ /palm/trunk/SOURCE/CURRENT_MODIFICATIONS	(revision 110)
@@ -1,68 +1,10 @@
 New:
 ---
-Allows runs for a coupled atmosphere-ocean LES,
-coupling frequency is controlled by new d3par-parameter dt_coupling,
-the coupling mode (atmosphere_to_ocean or ocean_to_atmosphere) for the
-respective processes is read from environment variable coupling_mode,
-which is set by the mpiexec-command,
-communication between the two models is done using the intercommunicator
-comm_inter,
-local files opened by the ocean model get the additional suffic "_O".
-Assume saturation at k=nzb_s_inner(j,i) for atmosphere coupled to ocean.
-
-A momentum flux can be set as top boundary condition using the new
-inipar parameter top_momentumflux_u|v.
-
-Non-cyclic boundary conditions can be used along all horizontal directions.
-
-Quantities w*p* and w"e can be output as vertical profiles.
-
-Initial profiles are reset to constant profiles in case that initializing_actions /= 'set_constant_profiles'. (init_rankine)
-
-Optionally calculate km and kh from initial TKE e_init.
-
-boundary_conds, check_open, check_parameters, diffusion_u, diffusion_v, flow_statistics, header, init_pegrid, init_rankine, init_3d_model, modules, palm, parin, pres, prandtl_fluxes, prognostic_equations, read_var_list, read_3d_binary, swap_timelevel, time_integration, write_var_list, write_3d_binary
-
-New:
-surface_coupler
 
 
 Changed:
 -------
-Remaining variables iran changed to iran_part (advec_particles, init_particles).
-
-In case that the presure solver is not called for every Runge-Kutta substep
-(call_psolver_at_all_substeps = .F.), it is called after the first substep
-instead of the last. In that case, random perturbations are also added to the
-velocity field after the first substep.
-
-Initialization of km,kh = 0.00001 for ocean = .T. (for ocean = .F. it remains 0.01).
-
-Allow data_output_pr= q, wq, w"q", w*q* for humidity = .T. (instead of cloud_physics = .T.).
-
-Modifications to terminate coupled runs.
-
-advec_particles, check_for_restart, check_parameters, init_particles, init_3d_model, local_stop, timestep, time_integration
 
 
 Errors:
 ------
-Bugs from code parts for non-cyclic boundary conditions are removed: loops for
-u and v are starting from index nxlu, nysv, respectively. The radiation boundary
-condition is used for every Runge-Kutta substep. Velocity phase speeds for
-the radiation boundary conditions are calculated for the first Runge-Kutta
-substep only and reused for the further substeps. New arrays c_u, c_v, and c_w
-are defined for this purpose. Several index errors are removed from the
-radiation boundary condition code parts. Upper bounds for calculating
-u_0 and v_0 (in production_e) are nxr+1 and nyn+1 because otherwise these
-values are not available in case of non-cyclic boundary conditions.
-
-+dots_num_palm in module user, +module netcdf_control in user_init (both in user_interface)
-
-Bugfix: wrong sign removed from the buoyancy production term in the case use_reference = .T. (production_e)
-
-Bugfix: Error message concerning output of particle concentration (pc) modified (check_parameters).
-
-Bugfix: Rayleigh damping for ocean fixed.
-
-advec_u_pw, advec_u_up, advec_v_pw, advec_v_up, boundary_conds, buoyancy, check_parameters, coriolis, diffusion_u, diffusion_v, header, init_pegrid, init_3d_model, modules, production_e, prognostic_equations, user_interface
Index: /palm/trunk/SOURCE/advec_particles.f90
===================================================================
--- /palm/trunk/SOURCE/advec_particles.f90	(revision 109)
+++ /palm/trunk/SOURCE/advec_particles.f90	(revision 110)
@@ -4,5 +4,4 @@
 ! Actual revisions:
 ! -----------------
-! remaining variables iran changed to iran_part
 ! TEST: PRINT statements on unit 9 (commented out)
 !
@@ -10,4 +9,7 @@
 ! -----------------
 ! $Id$
+!
+! 106 2007-08-16 14:30:26Z raasch
+! remaining variables iran changed to iran_part
 !
 ! 95 2007-06-02 16:48:38Z raasch
Index: /palm/trunk/SOURCE/advec_u_pw.f90
===================================================================
--- /palm/trunk/SOURCE/advec_u_pw.f90	(revision 109)
+++ /palm/trunk/SOURCE/advec_u_pw.f90	(revision 110)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! i loop is starting from nxlu (needed for non-cyclic boundary conditions)
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 106 2007-08-16 14:30:26Z raasch
+! i loop is starting from nxlu (needed for non-cyclic boundary conditions)
 !
 ! 75 2007-03-22 09:54:05Z raasch
Index: /palm/trunk/SOURCE/advec_u_up.f90
===================================================================
--- /palm/trunk/SOURCE/advec_u_up.f90	(revision 109)
+++ /palm/trunk/SOURCE/advec_u_up.f90	(revision 110)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! i loop is starting from nxlu (needed for non-cyclic boundary conditions)
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 106 2007-08-16 14:30:26Z raasch
+! i loop is starting from nxlu (needed for non-cyclic boundary conditions)
 !
 ! 75 2007-03-22 09:54:05Z raasch
Index: /palm/trunk/SOURCE/advec_v_pw.f90
===================================================================
--- /palm/trunk/SOURCE/advec_v_pw.f90	(revision 109)
+++ /palm/trunk/SOURCE/advec_v_pw.f90	(revision 110)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! j loop is starting from nysv (needed for non-cyclic boundary conditions)
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 106 2007-08-16 14:30:26Z raasch
+! j loop is starting from nysv (needed for non-cyclic boundary conditions)
 !
 ! 75 2007-03-22 09:54:05Z raasch
Index: /palm/trunk/SOURCE/advec_v_up.f90
===================================================================
--- /palm/trunk/SOURCE/advec_v_up.f90	(revision 109)
+++ /palm/trunk/SOURCE/advec_v_up.f90	(revision 110)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! j loop is starting from nysv (needed for non-cyclic boundary conditions)
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 106 2007-08-16 14:30:26Z raasch
+! j loop is starting from nysv (needed for non-cyclic boundary conditions)
 !
 ! 75 2007-03-22 09:54:05Z raasch
Index: /palm/trunk/SOURCE/boundary_conds.f90
===================================================================
--- /palm/trunk/SOURCE/boundary_conds.f90	(revision 109)
+++ /palm/trunk/SOURCE/boundary_conds.f90	(revision 110)
@@ -4,13 +4,15 @@
 ! Actual revisions:
 ! -----------------
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 107 2007-08-17 13:54:45Z raasch
 ! Boundary conditions for temperature adjusted for coupled runs,
 ! bugfixes for the radiation boundary conditions at the outflow: radiation
 ! conditions are used for every substep, phase speeds are calculated for the
 ! first Runge-Kutta substep only and then reused, several index values changed
-! 
-!
-! Former revisions:
-! -----------------
-! $Id$
 !
 ! 95 2007-06-02 16:48:38Z raasch
Index: /palm/trunk/SOURCE/buoyancy.f90
===================================================================
--- /palm/trunk/SOURCE/buoyancy.f90	(revision 109)
+++ /palm/trunk/SOURCE/buoyancy.f90	(revision 110)
@@ -4,6 +4,5 @@
 ! Actual revisions:
 ! -----------------
-! i loop for u-component (sloping surface) is starting from nxlu (needed for
-! non-cyclic boundary conditions)
+! 
 !
 ! Former revisions:
@@ -11,4 +10,8 @@
 ! $Id$
 !
+! 106 2007-08-16 14:30:26Z raasch
+! i loop for u-component (sloping surface) is starting from nxlu (needed for
+! non-cyclic boundary conditions)
+! 
 ! 97 2007-06-21 08:23:15Z raasch
 ! Routine reneralized to be used with temperature AND density:
Index: /palm/trunk/SOURCE/check_for_restart.f90
===================================================================
--- /palm/trunk/SOURCE/check_for_restart.f90	(revision 109)
+++ /palm/trunk/SOURCE/check_for_restart.f90	(revision 110)
@@ -4,5 +4,4 @@
 ! Actual revisions:
 ! -----------------
-! modifications to terminate coupled runs
 ! 
 !
@@ -10,4 +9,8 @@
 ! -----------------
 ! $Id$
+!
+! 108 2007-08-24 15:10:38Z letzel
+! modifications to terminate coupled runs
+!
 ! RCS Log replace by Id keyword, revision history cleaned up
 !
@@ -61,8 +64,9 @@
 !-- Output that job will be terminated
     IF ( terminate_run  .AND.  myid == 0 )  THEN
-       PRINT*, '*** WARNING: run will be terminated because it is running out of job cpu', &
-                    'limit'
+       PRINT*, '*** WARNING: run will be terminated because it is running out', &
+                    ' of job cpu limit'
        PRINT*, '             remaining time:         ', remaining_time, ' s'
-       PRINT*, '             termination time needed:', termination_time_needed, ' s'
+       PRINT*, '             termination time needed:', termination_time_needed,&
+                    ' s'
     ENDIF
 
@@ -72,12 +76,11 @@
 !-- informed of another termination reason (terminate_coupled > 0) before, 
 !-- or vice versa (terminate_coupled_remote > 0).
-    IF ( terminate_run .AND. TRIM( coupling_mode ) /= 'uncoupled'  &
-         .AND. terminate_coupled == 0  &
-         .AND. terminate_coupled_remote == 0 )  THEN
+    IF ( terminate_run .AND. TRIM( coupling_mode ) /= 'uncoupled'  .AND. &
+         terminate_coupled == 0  .AND.  terminate_coupled_remote == 0 )  THEN
+
        terminate_coupled = 3
-       CALL MPI_SENDRECV( &
-            terminate_coupled,        1, MPI_INTEGER, myid,  0, &
-            terminate_coupled_remote, 1, MPI_INTEGER, myid,  0, &
-            comm_inter, status, ierr )
+       CALL MPI_SENDRECV( terminate_coupled,        1, MPI_INTEGER, myid,  0, &
+                          terminate_coupled_remote, 1, MPI_INTEGER, myid,  0, &
+                          comm_inter, status, ierr )
     ENDIF
 
@@ -101,5 +104,6 @@
 
           IF ( myid == 0 )  THEN
-             PRINT*, '*** INFORMATIVE: run will be terminated due to user settings of'
+             PRINT*, '*** INFORMATIVE: run will be terminated due to user ', &
+                                       'settings of'
              PRINT*, '                 restart_time / dt_restart'
              PRINT*, '                 new restart time is: ', time_restart, ' s'
@@ -110,6 +114,7 @@
 !--       informed of another termination reason (terminate_coupled > 0) before,
 !--       or vice versa (terminate_coupled_remote > 0).
-          IF ( coupling_mode /= 'uncoupled' .AND. terminate_coupled == 0  &
-               .AND. terminate_coupled_remote == 0)  THEN
+          IF ( coupling_mode /= 'uncoupled' .AND. terminate_coupled == 0  .AND. &
+               terminate_coupled_remote == 0)  THEN
+
              IF ( dt_restart /= 9999999.9 )  THEN
                 terminate_coupled = 4
@@ -117,8 +122,8 @@
                 terminate_coupled = 5
              ENDIF
-             CALL MPI_SENDRECV( &
-                  terminate_coupled,        1, MPI_INTEGER, myid,  0, &
-                  terminate_coupled_remote, 1, MPI_INTEGER, myid,  0, &
-                  comm_inter, status, ierr )
+             CALL MPI_SENDRECV(                                                 &
+                            terminate_coupled,        1, MPI_INTEGER, myid,  0, &
+                            terminate_coupled_remote, 1, MPI_INTEGER, myid,  0, &
+                            comm_inter, status, ierr )
           ENDIF
        ELSE
Index: /palm/trunk/SOURCE/check_open.f90
===================================================================
--- /palm/trunk/SOURCE/check_open.f90	(revision 109)
+++ /palm/trunk/SOURCE/check_open.f90	(revision 110)
@@ -4,10 +4,13 @@
 ! Actual revisions:
 ! -----------------
-! Different filenames are used in case of a coupled simulation,
-! coupling_char added to all relevant filenames
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 105 2007-08-08 07:12:55Z raasch
+! Different filenames are used in case of a coupled simulation,
+! coupling_char added to all relevant filenames
 !
 ! 82 2007-04-16 15:40:52Z raasch
Index: /palm/trunk/SOURCE/check_parameters.f90
===================================================================
--- /palm/trunk/SOURCE/check_parameters.f90	(revision 109)
+++ /palm/trunk/SOURCE/check_parameters.f90	(revision 110)
@@ -4,4 +4,11 @@
 ! Actual revisions:
 ! -----------------
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 109 2007-08-28 15:26:47Z letzel
 ! Check coupling_mode and set default (obligatory) values (like boundary
 ! conditions for temperature and fluxes) in case of coupled runs.
@@ -14,8 +21,4 @@
 ! Rayleigh damping for ocean fixed.
 ! Check and, if necessary, set default value for dt_coupling
-!
-! Former revisions:
-! -----------------
-! $Id$
 !
 ! 97 2007-06-21 08:23:15Z raasch
Index: /palm/trunk/SOURCE/coriolis.f90
===================================================================
--- /palm/trunk/SOURCE/coriolis.f90	(revision 109)
+++ /palm/trunk/SOURCE/coriolis.f90	(revision 110)
@@ -4,10 +4,13 @@
 ! Actual revisions:
 ! -----------------
-! loops for u and v are starting from index nxlu, nysv, respectively (needed
-! for non-cyclic boundary conditions)
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 106 2007-08-16 14:30:26Z raasch
+! loops for u and v are starting from index nxlu, nysv, respectively (needed
+! for non-cyclic boundary conditions)
 !
 ! 75 2007-03-22 09:54:05Z raasch
Index: /palm/trunk/SOURCE/diffusion_u.f90
===================================================================
--- /palm/trunk/SOURCE/diffusion_u.f90	(revision 109)
+++ /palm/trunk/SOURCE/diffusion_u.f90	(revision 110)
@@ -4,10 +4,13 @@
 ! Actual revisions:
 ! -----------------
-! Momentumflux at top (uswst) included as boundary condition,
-! i loop is starting from nxlu (needed for non-cyclic boundary conditions)
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 106 2007-08-16 14:30:26Z raasch
+! Momentumflux at top (uswst) included as boundary condition,
+! i loop is starting from nxlu (needed for non-cyclic boundary conditions)
 !
 ! 75 2007-03-22 09:54:05Z raasch
Index: /palm/trunk/SOURCE/diffusion_v.f90
===================================================================
--- /palm/trunk/SOURCE/diffusion_v.f90	(revision 109)
+++ /palm/trunk/SOURCE/diffusion_v.f90	(revision 110)
@@ -4,10 +4,13 @@
 ! Actual revisions:
 ! -----------------
-! Momentumflux at top (vswst) included as boundary condition,
-! j loop is starting from nysv (needed for non-cyclic boundary conditions)
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 106 2007-08-16 14:30:26Z raasch
+! Momentumflux at top (vswst) included as boundary condition,
+! j loop is starting from nysv (needed for non-cyclic boundary conditions)
 !
 ! 75 2007-03-22 09:54:05Z raasch
Index: /palm/trunk/SOURCE/flow_statistics.f90
===================================================================
--- /palm/trunk/SOURCE/flow_statistics.f90	(revision 109)
+++ /palm/trunk/SOURCE/flow_statistics.f90	(revision 110)
@@ -4,10 +4,13 @@
 ! Actual revisions:
 ! -----------------
-! Prescribed momentum fluxes at the top surface are used,
-! profiles for w*p* and w"e are calculated
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 106 2007-08-16 14:30:26Z raasch
+! Prescribed momentum fluxes at the top surface are used,
+! profiles for w*p* and w"e are calculated
 !
 ! 97 2007-06-21 08:23:15Z raasch
Index: /palm/trunk/SOURCE/header.f90
===================================================================
--- /palm/trunk/SOURCE/header.f90	(revision 109)
+++ /palm/trunk/SOURCE/header.f90	(revision 110)
@@ -4,11 +4,14 @@
 ! Actual revisions:
 ! -----------------
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 108 2007-08-24 15:10:38Z letzel
 ! Output of informations for coupled model runs (boundary conditions etc.)
 ! + output of momentumfluxes at the top boundary
 ! Rayleigh damping for ocean, e_init
-!
-! Former revisions:
-! -----------------
-! $Id$
 !
 ! 97 2007-06-21 08:23:15Z raasch
Index: /palm/trunk/SOURCE/init_3d_model.f90
===================================================================
--- /palm/trunk/SOURCE/init_3d_model.f90	(revision 109)
+++ /palm/trunk/SOURCE/init_3d_model.f90	(revision 110)
@@ -7,13 +7,16 @@
 ! Actual revisions:
 ! -----------------
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 108 2007-08-24 15:10:38Z letzel
 ! Flux initialization in case of coupled runs, +momentum fluxes at top boundary,
 ! +arrays for phase speed c_u, c_v, c_w, indices for u|v|w_m_l|r changed
 ! +qswst_remote in case of atmosphere model with humidity coupled to ocean
-! Rayleigh damping for ocean
-! optionally calculate km and kh from initial TKE e_init
-!
-! Former revisions:
-! -----------------
-! $Id$
+! Rayleigh damping for ocean, optionally calculate km and kh from initial
+! TKE e_init
 !
 ! 97 2007-06-21 08:23:15Z raasch
Index: /palm/trunk/SOURCE/init_particles.f90
===================================================================
--- /palm/trunk/SOURCE/init_particles.f90	(revision 109)
+++ /palm/trunk/SOURCE/init_particles.f90	(revision 110)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! variable iran replaced by iran_part
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 106 2007-08-16 14:30:26Z raasch
+! variable iran replaced by iran_part
 !
 ! 82 2007-04-16 15:40:52Z raasch
Index: /palm/trunk/SOURCE/init_pegrid.f90
===================================================================
--- /palm/trunk/SOURCE/init_pegrid.f90	(revision 109)
+++ /palm/trunk/SOURCE/init_pegrid.f90	(revision 110)
@@ -5,12 +5,14 @@
 ! -----------------
 ! TEST OUTPUT (TO BE REMOVED) logging mpi2 ierr values
-! Intercommunicator (comm_inter) and derived data type (type_xy) for
-! coupled model runs created, assign coupling_mode_remote
-! indices nxlu and nysv are calculated (needed for non-cyclic boundary
-! conditions)
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 108 2007-08-24 15:10:38Z letzel
+! Intercommunicator (comm_inter) and derived data type (type_xy) for
+! coupled model runs created, assign coupling_mode_remote,
+! indices nxlu and nysv are calculated (needed for non-cyclic boundary
+! conditions)
 !
 ! 82 2007-04-16 15:40:52Z raasch
@@ -894,6 +896,6 @@
 #endif
 !
-!-- At the outflow, u or v, respectively, have to be calculated for one grid
-!-- point less.
+!-- At the outflow, u or v, respectively, have to be calculated for one more
+!-- grid point.
     IF ( outflow_l )  THEN
        nxlu = nxl + 1
Index: /palm/trunk/SOURCE/init_rankine.f90
===================================================================
--- /palm/trunk/SOURCE/init_rankine.f90	(revision 109)
+++ /palm/trunk/SOURCE/init_rankine.f90	(revision 110)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! Initial profiles are reset to constant profiles
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 107 2007-08-17 13:54:45Z raasch
+! Initial profiles are reset to constant profiles
 !
 ! 75 2007-03-22 09:54:05Z raasch
Index: /palm/trunk/SOURCE/local_stop.f90
===================================================================
--- /palm/trunk/SOURCE/local_stop.f90	(revision 109)
+++ /palm/trunk/SOURCE/local_stop.f90	(revision 110)
@@ -4,5 +4,5 @@
 ! Actual revisions:
 ! -----------------
-! modifications to terminate coupled runs
+! 
 ! 
 !
@@ -10,4 +10,8 @@
 ! -----------------
 ! $Id$
+!
+! 108 2007-08-24 15:10:38Z letzel
+! modifications to terminate coupled runs
+!
 ! RCS Log replace by Id keyword, revision history cleaned up
 !
Index: /palm/trunk/SOURCE/modules.f90
===================================================================
--- /palm/trunk/SOURCE/modules.f90	(revision 109)
+++ /palm/trunk/SOURCE/modules.f90	(revision 110)
@@ -5,4 +5,11 @@
 ! Actual revisions:
 ! -----------------
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 108 2007-08-24 15:10:38Z letzel
 ! +comm_inter, constant_top_momentumflux, coupling_char, coupling_mode,
 ! coupling_mode_remote, c_u, c_v, c_w, dt_coupling, e_init, humidity_remote, 
@@ -10,8 +17,4 @@
 ! terminate_coupled_remote, time_coupling, top_momentumflux_u|v, type_xy, 
 ! uswst*, vswst*
-!
-! Former revisions:
-! -----------------
-! $Id$
 !
 ! 97 2007-06-21 08:23:15Z raasch
Index: /palm/trunk/SOURCE/palm.f90
===================================================================
--- /palm/trunk/SOURCE/palm.f90	(revision 109)
+++ /palm/trunk/SOURCE/palm.f90	(revision 110)
@@ -4,10 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! Get coupling mode from environment variable
-! Change localtion of debug output
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 108 2007-08-24 15:10:38Z letzel
+! Get coupling mode from environment variable, change location of debug output
 !
 ! 75 2007-03-22 09:54:05Z raasch
@@ -64,5 +66,5 @@
     INTEGER           ::  i, run_description_header_i(80)
 
-    version = 'PALM 3.3a'
+    version = 'PALM 3.4'
 
 #if defined( __parallel )
Index: /palm/trunk/SOURCE/parin.f90
===================================================================
--- /palm/trunk/SOURCE/parin.f90	(revision 109)
+++ /palm/trunk/SOURCE/parin.f90	(revision 110)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! +e_init, top_momentumflux_u|v in inipar, +dt_coupling in d3par
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 108 2007-08-24 15:10:38Z letzel
+! +e_init, top_momentumflux_u|v in inipar, +dt_coupling in d3par
 !
 ! 95 2007-06-02 16:48:38Z raasch
Index: /palm/trunk/SOURCE/prandtl_fluxes.f90
===================================================================
--- /palm/trunk/SOURCE/prandtl_fluxes.f90	(revision 109)
+++ /palm/trunk/SOURCE/prandtl_fluxes.f90	(revision 110)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! assume saturation at k=nzb_s_inner(j,i) for atmosphere coupled to ocean
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 108 2007-08-24 15:10:38Z letzel
+! assume saturation at k=nzb_s_inner(j,i) for atmosphere coupled to ocean
 !
 ! 75 2007-03-22 09:54:05Z raasch
Index: /palm/trunk/SOURCE/pres.f90
===================================================================
--- /palm/trunk/SOURCE/pres.f90	(revision 109)
+++ /palm/trunk/SOURCE/pres.f90	(revision 110)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! Volume flow conservation added for the remaining three outflow boundaries
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 106 2007-08-16 14:30:26Z raasch
+! Volume flow conservation added for the remaining three outflow boundaries
 !
 ! 85 2007-05-11 09:35:14Z raasch
Index: /palm/trunk/SOURCE/production_e.f90
===================================================================
--- /palm/trunk/SOURCE/production_e.f90	(revision 109)
+++ /palm/trunk/SOURCE/production_e.f90	(revision 110)
@@ -4,4 +4,11 @@
 ! Actual revisions:
 ! -----------------
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 108 2007-08-24 15:10:38Z letzel
 ! Bugfix: wrong sign removed from the buoyancy production term in the case
 ! use_reference = .T.,
@@ -9,8 +16,4 @@
 ! not available in case of non-cyclic boundary conditions)
 ! Bugfix for ocean density flux at bottom
-!
-! Former revisions:
-! -----------------
-! $Id$
 !
 ! 97 2007-06-21 08:23:15Z raasch
Index: /palm/trunk/SOURCE/prognostic_equations.f90
===================================================================
--- /palm/trunk/SOURCE/prognostic_equations.f90	(revision 109)
+++ /palm/trunk/SOURCE/prognostic_equations.f90	(revision 110)
@@ -4,11 +4,14 @@
 ! Actual revisions:
 ! -----------------
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 106 2007-08-16 14:30:26Z raasch
 ! +uswst, vswst as arguments in calls of diffusion_u|v,
 ! loops for u and v are starting from index nxlu, nysv, respectively (needed
 ! for non-cyclic boundary conditions)
-!
-! Former revisions:
-! -----------------
-! $Id$
 !
 ! 97 2007-06-21 08:23:15Z raasch
Index: /palm/trunk/SOURCE/read_3d_binary.f90
===================================================================
--- /palm/trunk/SOURCE/read_3d_binary.f90	(revision 109)
+++ /palm/trunk/SOURCE/read_3d_binary.f90	(revision 110)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! +uswst, uswst_m, vswst, vswst_m
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 102 2007-07-27 09:09:17Z raasch
+! +uswst, uswst_m, vswst, vswst_m
 !
 ! 96 2007-06-04 08:07:41Z raasch
Index: /palm/trunk/SOURCE/read_var_list.f90
===================================================================
--- /palm/trunk/SOURCE/read_var_list.f90	(revision 109)
+++ /palm/trunk/SOURCE/read_var_list.f90	(revision 110)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! +time_coupling, top_momentumflux_u|v
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 102 2007-07-27 09:09:17Z raasch
+! +time_coupling, top_momentumflux_u|v
 !
 ! 95 2007-06-02 16:48:38Z raasch
Index: /palm/trunk/SOURCE/surface_coupler.f90
===================================================================
--- /palm/trunk/SOURCE/surface_coupler.f90	(revision 109)
+++ /palm/trunk/SOURCE/surface_coupler.f90	(revision 110)
@@ -4,6 +4,4 @@
 ! Actual revisions:
 ! -----------------
-! include latent heatflux and salinity flux for atmosphere runs with moisture
-! modifications to terminate coupled runs
 ! 
 !
@@ -12,4 +10,5 @@
 ! $Id$
 !
+! 109 2007-08-28 15:26:47Z letzel
 ! Initial revision
 !
Index: /palm/trunk/SOURCE/swap_timelevel.f90
===================================================================
--- /palm/trunk/SOURCE/swap_timelevel.f90	(revision 109)
+++ /palm/trunk/SOURCE/swap_timelevel.f90	(revision 110)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! swaping of uswst, vswst included
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 102 2007-07-27 09:09:17Z raasch
+! swaping of uswst, vswst included
 !
 ! 95 2007-06-02 16:48:38Z raasch
Index: /palm/trunk/SOURCE/time_integration.f90
===================================================================
--- /palm/trunk/SOURCE/time_integration.f90	(revision 109)
+++ /palm/trunk/SOURCE/time_integration.f90	(revision 110)
@@ -4,4 +4,11 @@
 ! Actual revisions:
 ! -----------------
+! 
+!
+! Former revisions:
+! -----------------
+! $Id$
+!
+! 108 2007-08-24 15:10:38Z letzel
 ! Call of new routine surface_coupler,
 ! presure solver is called after the first Runge-Kutta substep instead of the
@@ -9,10 +16,4 @@
 ! random perturbation has to be added to the velocity fields also after the
 ! first substep
-! modifications to terminate coupled runs
-! 
-!
-! Former revisions:
-! -----------------
-! $Id$
 !
 ! 97 2007-06-21 08:23:15Z raasch
Index: /palm/trunk/SOURCE/timestep.f90
===================================================================
--- /palm/trunk/SOURCE/timestep.f90	(revision 109)
+++ /palm/trunk/SOURCE/timestep.f90	(revision 110)
@@ -4,5 +4,4 @@
 ! Actual revisions:
 ! -----------------
-! modifications to terminate coupled runs
 ! 
 !
@@ -10,4 +9,8 @@
 ! -----------------
 ! $Id$
+!
+! 108 2007-08-24 15:10:38Z letzel
+! modifications to terminate coupled runs
+!
 ! RCS Log replace by Id keyword, revision history cleaned up
 !
Index: /palm/trunk/SOURCE/user_interface.f90
===================================================================
--- /palm/trunk/SOURCE/user_interface.f90	(revision 109)
+++ /palm/trunk/SOURCE/user_interface.f90	(revision 110)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! +dots_num_palm in module user, +module netcdf_control in user_init
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 105 2007-08-08 07:12:55Z raasch
+! +dots_num_palm in module user, +module netcdf_control in user_init
 !
 ! 95 2007-06-02 16:48:38Z raasch
Index: /palm/trunk/SOURCE/write_3d_binary.f90
===================================================================
--- /palm/trunk/SOURCE/write_3d_binary.f90	(revision 109)
+++ /palm/trunk/SOURCE/write_3d_binary.f90	(revision 110)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! +uswst, uswst_m, vswst, vswst_m
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 102 2007-07-27 09:09:17Z raasch
+! +uswst, uswst_m, vswst, vswst_m
 !
 ! 96 2007-06-04 08:07:41Z raasch
Index: /palm/trunk/SOURCE/write_var_list.f90
===================================================================
--- /palm/trunk/SOURCE/write_var_list.f90	(revision 109)
+++ /palm/trunk/SOURCE/write_var_list.f90	(revision 110)
@@ -4,9 +4,12 @@
 ! Actual revisions:
 ! -----------------
-! +top_momentumflux_u|v, time_coupling
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+!
+! 102 2007-07-27 09:09:17Z raasch
+! +top_momentumflux_u|v, time_coupling
 !
 ! 95 2007-06-02 16:48:38Z raasch
