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

Timestamp:

Apr 6, 2016 1:49:59 PM (8 years ago)

Author:

raasch

Message:

cpp-switches removed + cpp-bugfixes + zero-settings for velocities inside topography re-activated

File:

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

palm/trunk/SOURCE/init_3d_model.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-! 
+! zero-settings for velocities inside topography re-activated (was deactivated
+! in r1762)
 ! 
 ! Former revisions:
@@ -857 +858 @@
 !--       levels to zero in order to avoid too small time steps caused by the 
 !--       diffusion limit in the initial phase of a run (at k=1, dz/2 occurs
-!--       in the limiting formula!). The original values are stored to be later
-!--       used for volume flow control.
-          ! TO_DO:  Antti will check if this is really required
-          !AH IF ( ibc_uv_b /= 1 )  THEN
-          !AH    DO  i = nxlg, nxrg
-          !AH       DO  j = nysg, nyng
-          !AH          u(nzb:nzb_u_inner(j,i)+1,j,i) = 0.0_wp
-          !AH          v(nzb:nzb_v_inner(j,i)+1,j,i) = 0.0_wp
-          !AH       ENDDO
-          !AH    ENDDO
-          !AH ENDIF
+!--       in the limiting formula!).
+          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