Changeset 1406 for palm/trunk/SOURCE/init_3d_model.f90

Timestamp:

May 16, 2014 1:47:01 PM (10 years ago)

Author:

raasch

Message:

 

File:

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

palm/trunk/SOURCE/init_3d_model.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! bugfix: setting of initial velocities at k=1 to zero not in case of a
+! no-slip boundary condition for uv
 ! 
 ! Former revisions:
@@ -782 +783 @@
 !--       in the limiting formula!). The original values are stored to be later
 !--       used for volume flow control.
-          DO  i = nxlg, nxrg
-             DO  j = nysg, nyng
-                u(nzb:nzb_u_inner(j,i)+1,j,i) = 0.0_wp
-                v(nzb:nzb_v_inner(j,i)+1,j,i) = 0.0_wp
-             ENDDO
-          ENDDO
+          IF ( ibc_uv_b /= 1 )  THEN
+             DO  i = nxlg, nxrg
+                DO  j = nysg, nyng
+                   u(nzb:nzb_u_inner(j,i)+1,j,i) = 0.0_wp
+                   v(nzb:nzb_v_inner(j,i)+1,j,i) = 0.0_wp
+                ENDDO
+             ENDDO
+          ENDIF
 
           IF ( humidity  .OR.  passive_scalar )  THEN