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

Timestamp:

May 21, 2013 11:58:22 AM (8 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/modules.f90 (modified) (3 diffs)

palm/trunk/SOURCE/modules.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-!
+! -bc_lr_dirneu, bc_lr_neudir, bc_ns_dirneu, bc_ns_neudir
+! +use_cmax
 !
 ! Former revisions:
@@ -730 +731 @@
 
     LOGICAL ::  avs_output = .FALSE., &
-                bc_lr_cyc =.TRUE., bc_lr_dirneu = .FALSE., &
-                bc_lr_dirrad = .FALSE., bc_lr_neudir = .FALSE., &
+                bc_lr_cyc =.TRUE., bc_lr_dirrad = .FALSE., &
                 bc_lr_raddir = .FALSE., bc_ns_cyc = .TRUE., &
-                bc_ns_dirneu = .FALSE., bc_ns_dirrad = .FALSE., &
-                bc_ns_neudir = .FALSE., bc_ns_raddir = .FALSE., &
+                bc_ns_dirrad = .FALSE., bc_ns_raddir = .FALSE.,&
                 call_psolver_at_all_substeps = .TRUE., &
                 cloud_droplets = .FALSE., cloud_physics = .FALSE., &
@@ -767 +766 @@
                 stop_dt = .FALSE., synchronous_exchange = .FALSE., &
                 terminate_run = .FALSE., turbulence = .FALSE., &
-                turbulent_inflow = .FALSE., &
+                turbulent_inflow = .FALSE., use_cmax = .TRUE., &
                 use_prescribed_profile_data = .FALSE., use_reference = .FALSE.,&
                 use_surface_fluxes = .FALSE., use_top_fluxes = .FALSE., &