Changeset 106 for palm/trunk/SOURCE/time_integration.f90

Timestamp:

Aug 16, 2007 2:30:26 PM (17 years ago)

Author:

raasch

Message:

preliminary update of bugfixes and extensions for non-cyclic BCs

File:

• palm/trunk/SOURCE/time_integration.f90 (modified) (3 diffs)

palm/trunk/SOURCE/time_integration.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! Call of new routine surface_coupler
+! Call of new routine surface_coupler,
+! presure solver is called after the first Runge-Kutta substep instead of the
+! last in case that call_psolver_at_all_substeps = .F.; for this case, the
+! random perturbation has to be added to the velocity fields also after the
+! first substep
 !
 ! Former revisions:
@@ -194 +198 @@
 !
 !--       Impose a random perturbation on the horizontal velocity field
-          IF ( create_disturbances  .AND.  &
+          IF ( create_disturbances  .AND.                                      &
+               ( call_psolver_at_all_substeps  .AND.                           &
                intermediate_timestep_count == intermediate_timestep_count_max )&
+          .OR. ( .NOT. call_psolver_at_all_substeps  .AND.                     &
+               intermediate_timestep_count == 1 ) )                            &
           THEN
              time_disturb = time_disturb + dt_3d
@@ -218 +225 @@
 !--       Reduce the velocity divergence via the equation for perturbation
 !--       pressure.
-          IF ( intermediate_timestep_count == intermediate_timestep_count_max &
-               .OR.  call_psolver_at_all_substeps )  THEN
+          IF ( intermediate_timestep_count == 1  .OR. &
+                call_psolver_at_all_substeps )  THEN
              CALL pres
           ENDIF
-
-!
-!--       In case of a non-cyclic lateral wall, set the boundary conditions for
-!--       the velocities at the outflow
-!          IF ( bc_lr /= 'cyclic'  .OR.  bc_ns /= 'cyclic' )  THEN
-!             CALL boundary_conds( 'outflow_uvw' )
-!          ENDIF
 
 !