Changeset 1968 for palm/trunk/SOURCE/exchange_horiz_2d.f90

Timestamp:

Jul 18, 2016 12:01:49 PM (8 years ago)

Author:

suehring

Message:

PE-wise reading of topography file

File:

• palm/trunk/SOURCE/exchange_horiz_2d.f90 (modified) (8 diffs)

palm/trunk/SOURCE/exchange_horiz_2d.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! 2D-INTEGER exchange adopted for different multigrid level
 ! 
 ! Former revisions:
@@ -203 +203 @@
 !------------------------------------------------------------------------------!
  
- SUBROUTINE exchange_horiz_2d_int( ar )
+ SUBROUTINE exchange_horiz_2d_int( ar, nys_l, nyn_l, nxl_l, nxr_l, nbgp_local )
 
 
     USE control_parameters,                                                    &
-        ONLY:  bc_lr_cyc, bc_ns_cyc, nest_bound_l, nest_bound_n, nest_bound_r, &
-               nest_bound_s
+        ONLY:  bc_lr_cyc, bc_ns_cyc, grid_level, nest_bound_l, nest_bound_n,   &
+               nest_bound_r, nest_bound_s
         
     USE cpulog,                                                                &
         ONLY:  cpu_log, log_point_s
-        
-    USE indices,                                                               &
-        ONLY:  nbgp, nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg
-        
+               
     USE kinds
     
@@ -222 +219 @@
     IMPLICIT NONE
 
-    INTEGER(iwp) ::  i
-    INTEGER(iwp) ::  ar(nysg:nyng,nxlg:nxrg)  !<
+    INTEGER(iwp) ::  i           !< dummy index to zero-gradient conditions at in/outflow boundaries
+    INTEGER(iwp) ::  nxl_l       !< local index bound at current grid level, left side
+    INTEGER(iwp) ::  nxr_l       !< local index bound at current grid level, right side
+    INTEGER(iwp) ::  nyn_l       !< local index bound at current grid level, north side
+    INTEGER(iwp) ::  nys_l       !< local index bound at current grid level, south side
+    INTEGER(iwp) ::  nbgp_local  !< number of ghost layers to be exchanged
+
+    INTEGER(iwp), DIMENSION(nys_l-nbgp_local:nyn_l+nbgp_local,                 &
+                            nxl_l-nbgp_local:nxr_l+nbgp_local) ::  ar  !< treated array
 
     CALL cpu_log( log_point_s(13), 'exchange_horiz_2d', 'start' )
@@ -236 +240 @@
 !--    One-dimensional decomposition along y, boundary values can be exchanged
 !--    within the PE memory
-       ar(:,nxlg:nxl-1) = ar(:,nxr-nbgp+1:nxr)
-       ar(:,nxr+1:nxrg) = ar(:,nxl:nxl+nbgp-1)
-
+       ar(:,nxl_l-nbgp_local:nxl_l-1) = ar(:,nxr_l-nbgp_local+1:nxr_l)
+       ar(:,nxr_l+1:nxr_l+nbgp_local) = ar(:,nxl_l:nxl_l+nbgp_local-1)
 
     ELSE
 !
 !--    Send left boundary, receive right one
-       CALL MPI_SENDRECV( ar(nysg,nxl), 1, type_y_int, pleft,  0,             &
-                          ar(nysg,nxr+1), 1, type_y_int, pright, 0,           &
+       CALL MPI_SENDRECV( ar(nys_l-nbgp_local,nxl_l),   1,                     &
+                          type_y_int(grid_level), pleft,  0,                   &
+                          ar(nys_l-nbgp_local,nxr_l+1), 1,                     &
+                          type_y_int(grid_level), pright, 0,                   &
                           comm2d, status, ierr )
 !
 !--    Send right boundary, receive left one
-       CALL MPI_SENDRECV( ar(nysg,nxr+1-nbgp), 1, type_y_int, pright,  1,     &
-                          ar(nysg,nxlg), 1, type_y_int, pleft,   1,           &
-                          comm2d, status, ierr )
+       CALL MPI_SENDRECV( ar(nys_l-nbgp_local,nxr_l+1-nbgp_local), 1,          &
+                          type_y_int(grid_level), pright, 1,                   &
+                          ar(nys_l-nbgp_local,nxl_l-nbgp_local),   1,          & 
+                          type_y_int(grid_level), pleft,  1,                   &
+                          comm2d, status, ierr )                         
 
     ENDIF
@@ -258 +265 @@
 !--    One-dimensional decomposition along x, boundary values can be exchanged
 !--    within the PE memory
-       ar(nysg:nys-1,:) = ar(nyn+1-nbgp:nyn,:)
-       ar(nyn+1:nyng,:) = ar(nys:nys-1+nbgp,:)
+       ar(nys_l-nbgp_local:nys_l-1,:) = ar(nyn_l+1-nbgp_local:nyn_l,:)
+       ar(nyn_l+1:nyn_l+nbgp_local,:) = ar(nys_l:nys_l-1+nbgp_local,:)
 
 
@@ -265 +272 @@
 !
 !--    Send front boundary, receive rear one
-       CALL MPI_SENDRECV( ar(nys,nxlg), 1, type_x_int, psouth, 0,             &
-                          ar(nyn+1,nxlg), 1, type_x_int, pnorth, 0,           &
+       CALL MPI_SENDRECV( ar(nys_l,nxl_l-nbgp_local),   1,                    &
+                          type_x_int(grid_level), psouth, 0,                  &
+                          ar(nyn_l+1,nxl_l-nbgp_local), 1,                    &
+                          type_x_int(grid_level), pnorth, 0,                  &
                           comm2d, status, ierr )                          
 
 !
 !--    Send rear boundary, receive front one
-       CALL MPI_SENDRECV( ar(nyn+1-nbgp,nxlg), 1, type_x_int, pnorth, 1,      &
-                          ar(nysg,nxlg), 1, type_x_int, psouth, 1,            &
+       CALL MPI_SENDRECV( ar(nyn_l+1-nbgp_local,nxl_l-nbgp_local), 1,         &
+                          type_x_int(grid_level), pnorth, 1,                  &
+                          ar(nys_l-nbgp_local,nxl_l-nbgp_local),   1,         &
+                          type_x_int(grid_level), psouth, 1,                  &
                           comm2d, status, ierr )
 
@@ -282 +293 @@
 !-- Lateral boundary conditions in the non-parallel case
     IF ( bc_lr_cyc )  THEN
-       ar(:,nxlg:nxl-1) = ar(:,nxr-nbgp+1:nxr)
-       ar(:,nxr+1:nxrg) = ar(:,nxl:nxl+nbgp-1)
+       ar(:,nxl_l-nbgp_local:nxl_l-1) = ar(:,nxr_l-nbgp_local+1:nxr_l)
+       ar(:,nxr_l+1:nxr_l+nbgp_local) = ar(:,nxl_l:nxl_l+nbgp_local-1)
     ENDIF
 
     IF ( bc_ns_cyc )  THEN
-       ar(nysg:nys-1,:) = ar(nyn+1-nbgp:nyn,:)
-       ar(nyn+1:nyng,:) = ar(nys:nys-1+nbgp,:)
+       ar(nys_l-nbgp_local:nys_l-1,:) = ar(nyn_l+1-nbgp_local:nyn_l,:)
+       ar(nyn_l+1:nyn_l+nbgp_local,:) = ar(nys_l:nys_l-1+nbgp_local,:)
     ENDIF
 
@@ -295 +306 @@
 !-- Neumann-conditions at inflow/outflow/nested boundaries
     IF ( nest_bound_l )  THEN
-       DO  i = nbgp, 1, -1
-         ar(:,nxl-i) = ar(:,nxl)
+       DO  i = nbgp_local, 1, -1
+         ar(:,nxl_l-i) = ar(:,nxl_l)
        ENDDO
     ENDIF
     IF ( nest_bound_r )  THEN
-       DO  i = 1, nbgp
-          ar(:,nxr+i) = ar(:,nxr)
+       DO  i = 1, nbgp_local
+          ar(:,nxr_l+i) = ar(:,nxr_l)
        ENDDO
     ENDIF
     IF ( nest_bound_s )  THEN
-       DO  i = nbgp, 1, -1
-         ar(nys-i,:) = ar(nys,:)
+       DO  i = nbgp_local, 1, -1
+         ar(nys_l-i,:) = ar(nys_l,:)
        ENDDO
     ENDIF
     IF ( nest_bound_n )  THEN
-       DO  i = 1, nbgp
-         ar(nyn+i,:) = ar(nyn,:)
+       DO  i = 1, nbgp_local
+         ar(nyn_l+i,:) = ar(nyn_l,:)
        ENDDO
     ENDIF