Changeset 1159 for palm/trunk/SOURCE/header.f90

Timestamp:

May 21, 2013 11:58:22 AM (11 years ago)

Author:

fricke

Message:

Bugfix: In case of non-cyclic lateral boundary conditions, Neumann boundary conditions for the velocity components at the outflow are in fact radiation boundary conditions using the maximum phase velocity that ensures numerical stability (CFL-condition).
Logical operator use_cmax is now used instead of bc_lr_dirneu/_neudir.

File:

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

palm/trunk/SOURCE/header.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
+! +use_cmax
 !
 !
@@ -742 +743 @@
     WRITE ( io, 317 )  bc_lr, bc_ns
     IF ( .NOT. bc_lr_cyc  .OR.  .NOT. bc_ns_cyc )  THEN
-       WRITE ( io, 318 )  pt_damping_width, pt_damping_factor       
+       WRITE ( io, 318 )  use_cmax, pt_damping_width, pt_damping_factor       
        IF ( turbulent_inflow )  THEN
           WRITE ( io, 319 )  recycling_width, recycling_plane, &
@@ -1781 +1782 @@
             '       left/right:  ',A/    &
             '       north/south: ',A)
-318 FORMAT (/'       pt damping layer width = ',F7.2,' m, pt ', &
+318 FORMAT (/'       use_cmax: ',L1 / &
+            '       pt damping layer width = ',F8.2,' m, pt ', &
                     'damping factor = ',F6.4)
 319 FORMAT ('       turbulence recycling at inflow switched on'/ &