Changeset 1411

Timestamp:

May 26, 2014 4:01:51 PM (10 years ago)

Author:

suehring

Message:

Initial horizontal velocity profiles were not set to zero at the first vertical

File:

• palm/trunk/SOURCE/init_3d_model.f90 (modified) (2 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Initial horizontal velocity profiles were not set to zero at the first vertical 
+! grid level in case of non-cyclic lateral boundary conditions. 
 ! 
 ! Former revisions:
@@ -786 +787 @@
 !--       in the limiting formula!). The original values are stored to be later
 !--       used for volume flow control.
-          IF ( ibc_uv_b /= 1 )  THEN
+          IF ( ibc_uv_b /= 1 .AND. ( .NOT. bc_lr_dirrad .AND.                 &
+                                     .NOT. bc_lr_raddir .AND.                 &
+                                     .NOT. bc_ns_dirrad .AND.                 &
+                                     .NOT. bc_ns_raddir ) )  THEN    
              DO  i = nxlg, nxrg
                 DO  j = nysg, nyng