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

Timestamp:

Feb 25, 2016 12:31:13 PM (8 years ago)

Author:

hellstea

Message:

Introduction of nested domain system

File:

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

palm/trunk/SOURCE/init_3d_model.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-!
+! Introduction of nested domain feature
 !
 ! Former revisions:
@@ -669 +669 @@
 !-- is done). Further below the values needed within the timestep scheme
 !-- will be set.
+!-- TO_DO: zeroth element added to weight_pres because in case of nesting
+!--        pres may need to be called  outside the RK-substeps. Antti will
+!--        check if this is really required.
     ALLOCATE( weight_substep(1:intermediate_timestep_count_max), &
-              weight_pres(1:intermediate_timestep_count_max) )
+              weight_pres(0:intermediate_timestep_count_max) )
     weight_substep = 1.0_wp
     weight_pres    = 1.0_wp
@@ -774 +777 @@
 
 !
+!--       Inside buildings set velocities back to zero
           IF ( topography /= 'flat' )  THEN
-             DO  i = nxl-1, nxr+1
-                DO  j = nys-1, nyn+1
+             DO  i = nxlg, nxrg
+                DO  j = nysg, nyng
                    u(nzb:nzb_u_inner(j,i),j,i) = 0.0_wp
                    v(nzb:nzb_v_inner(j,i),j,i) = 0.0_wp
@@ -841 +845 @@
 !--       in the limiting formula!). The original values are stored to be later
 !--       used for volume flow control.
-          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
+          ! 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
 
           IF ( humidity  .OR.  passive_scalar )  THEN
@@ -1669 +1674 @@
 !
 !-- Setting weighting factors for calculation of perturbation pressure 
-!-- and turbulent quantities from the RK substeps                
+!-- and turbulent quantities from the RK substeps
+!-- TO_DO: zeroth element is added to weight_pres because in nesting pres
+!--        may need to be called outside the RK-substeps
+    weight_pres(0) = 1.0_wp
     IF ( TRIM(timestep_scheme) == 'runge-kutta-3' )  THEN      ! for RK3-method