Ignore:
Timestamp:
Aug 16, 2007 2:30:26 PM (14 years ago)
Author:
raasch
Message:

preliminary update of bugfixes and extensions for non-cyclic BCs

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/time_integration.f90

    r102 r106  
    44! Actual revisions:
    55! -----------------
    6 ! Call of new routine surface_coupler
     6! Call of new routine surface_coupler,
     7! presure solver is called after the first Runge-Kutta substep instead of the
     8! last in case that call_psolver_at_all_substeps = .F.; for this case, the
     9! random perturbation has to be added to the velocity fields also after the
     10! first substep
    711!
    812! Former revisions:
     
    194198!
    195199!--       Impose a random perturbation on the horizontal velocity field
    196           IF ( create_disturbances  .AND.  &
     200          IF ( create_disturbances  .AND.                                      &
     201               ( call_psolver_at_all_substeps  .AND.                           &
    197202               intermediate_timestep_count == intermediate_timestep_count_max )&
     203          .OR. ( .NOT. call_psolver_at_all_substeps  .AND.                     &
     204               intermediate_timestep_count == 1 ) )                            &
    198205          THEN
    199206             time_disturb = time_disturb + dt_3d
     
    218225!--       Reduce the velocity divergence via the equation for perturbation
    219226!--       pressure.
    220           IF ( intermediate_timestep_count == intermediate_timestep_count_max &
    221                .OR. call_psolver_at_all_substeps )  THEN
     227          IF ( intermediate_timestep_count == 1  .OR. &
     228                call_psolver_at_all_substeps )  THEN
    222229             CALL pres
    223230          ENDIF
    224 
    225 !
    226 !--       In case of a non-cyclic lateral wall, set the boundary conditions for
    227 !--       the velocities at the outflow
    228 !          IF ( bc_lr /= 'cyclic'  .OR.  bc_ns /= 'cyclic' )  THEN
    229 !             CALL boundary_conds( 'outflow_uvw' )
    230 !          ENDIF
    231231
    232232!
Note: See TracChangeset for help on using the changeset viewer.