Changeset 3822 for palm/trunk/SOURCE

Timestamp:

Mar 27, 2019 1:10:23 PM (6 years ago)

Author:

hellstea

Message:

Temporary increase of the vertical dimension of the parent-grid arrays and workarrc_t is cancelled as unnecessary

File:

• palm/trunk/SOURCE/pmc_interface_mod.f90 (modified) (6 diffs)

palm/trunk/SOURCE/pmc_interface_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Temporary increase of the vertical dimension of the parent-grid arrays and 
+! workarrc_t is cancelled as unnecessary.
+! 
+! 3819 2019-03-27 11:01:36Z hellstea
 ! Adjustable anterpolation buffer introduced on all nest boundaries, it is controlled
 ! by the new nesting_parameters parameter anterpolation_buffer_width.
@@ -879 +883 @@
 !--       transfer
           DO  k = 1, nz                 
-!AH  Let's try to make the parent-grid arrays higher by one parent dz 
-!AH             IF ( zw(k) > fval(1) )  THEN
-             IF ( zw(k) > fval(1)+dz(1) )  THEN
+             IF ( zw(k) > fval(1) )  THEN
                 nz_cl = k
                 EXIT
              ENDIF
           ENDDO
-
-!AH  Let's make the parent-grid arrays higher by one parent dz 
-!AH          zmax_coarse = fval(1:2)
-!AH          cl_height   = fval(1)
-          zmax_coarse(1) = fval(1) + dz(1)
-          zmax_coarse(2) = fval(2) + dz(1)
-          cl_height   = fval(1) + dz(1)
-!AH
-
+          zmax_coarse = fval(1:2)
+          cl_height   = fval(1)
 !    
 !--       Get absolute coordinates from the child
@@ -1903 +1898 @@
 !
 !--    Top boundary.
-!AH       ALLOCATE( workarrc_t(0:2,jcsw:jcnw,iclw:icrw) )
-       ALLOCATE( workarrc_t(-2:3,jcsw:jcnw,iclw:icrw) )
+       ALLOCATE( workarrc_t(0:2,jcsw:jcnw,iclw:icrw) )
 
     END SUBROUTINE pmci_allocate_workarrays
@@ -1928 +1922 @@
 !
 !--    For the top-boundary x-slices
-!AH       CALL MPI_TYPE_VECTOR( icrw-iclw+1, 3, 3*(jcnw-jcsw+1), MPI_REAL,         &
-!AH            workarrc_t_exchange_type_x, ierr )
-       CALL MPI_TYPE_VECTOR( icrw-iclw+1, 6, 6*(jcnw-jcsw+1), MPI_REAL,         &
+       CALL MPI_TYPE_VECTOR( icrw-iclw+1, 3, 3*(jcnw-jcsw+1), MPI_REAL,         &
             workarrc_t_exchange_type_x, ierr )
        CALL MPI_TYPE_COMMIT( workarrc_t_exchange_type_x, ierr )
 !
 !--    For the top-boundary y-slices
-!AH       CALL MPI_TYPE_VECTOR( 1, 3*(jcnw-jcsw+1), 3*(jcnw-jcsw+1), MPI_REAL,     &
-!AH            workarrc_t_exchange_type_y, ierr )
-       CALL MPI_TYPE_VECTOR( 1, 6*(jcnw-jcsw+1), 6*(jcnw-jcsw+1), MPI_REAL,     &
+       CALL MPI_TYPE_VECTOR( 1, 3*(jcnw-jcsw+1), 3*(jcnw-jcsw+1), MPI_REAL,     &
             workarrc_t_exchange_type_y, ierr )
        CALL MPI_TYPE_COMMIT( workarrc_t_exchange_type_y, ierr )
@@ -3968 +3958 @@
       ENDIF
       workarrc_t = 0.0_wp
-      workarrc_t(-2:3,jcsc:jcnc,iclc:icrc) = fc(kct-2:kct+3,jcsc:jcnc,iclc:icrc)
+!      workarrc_t(-2:3,jcsc:jcnc,iclc:icrc) = fc(kct-2:kct+3,jcsc:jcnc,iclc:icrc)
+      workarrc_t(0:2,jcsc:jcnc,iclc:icrc) = fc(kct:kct+2,jcsc:jcnc,iclc:icrc)
 !
 !--   Left-right exchange if more than one PE subdomain in the x-direction.
@@ -4208 +4199 @@
       kcb  = 0
       IF ( nesting_mode /= 'vertical' )  THEN
-!         IF ( bc_dirichlet_l )  THEN
-!!            iclant = icl + 3
-!            iclant = icl + 3 + anterpolation_buffer_width
-!         ENDIF
-!         IF ( bc_dirichlet_r )  THEN
-!!            icrant = icr - 3
-!            icrant = icr - 3 - anterpolation_buffer_width
-!         ENDIF
-!
-!         IF ( bc_dirichlet_s )  THEN
-!!            jcsant = jcs + 3
-!            jcsant = jcs + 3 + anterpolation_buffer_width
-!         ENDIF
-!         IF ( bc_dirichlet_n )  THEN
-!!            jcnant = jcn - 3
-!            jcnant = jcn - 3 - anterpolation_buffer_width
-!         ENDIF
-
-!
-!--      New method: 
+!
+!--      Set the anterpolation buffers on the lateral boundaries
          iclant = MAX( icl, iplg + 3 + anterpolation_buffer_width )
          icrant = MIN( icr, iprg - 3 - anterpolation_buffer_width )