Changeset 4893 for palm/trunk/SOURCE/surface_data_output_mod.f90

Timestamp:

Mar 2, 2021 4:39:14 PM (3 years ago)

Author:

raasch

Message:

revised output of surface data via MPI-IO for better performance

File:

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

palm/trunk/SOURCE/surface_data_output_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! revised output of surface data via MPI-IO for better performance
+!
+! 4892 2021-03-02 11:53:58Z suehring
 ! Remove outdated error message.
 !
@@ -4578 +4581 @@
 
     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  end_index           !< end index of surface data at (j,i)
-    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  global_start_index  !< index array for surface data (MPI-IO)
+    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  global_end_index    !< end index array for surface data (MPI-IO)
+    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  global_start_index  !< start index array for surface data (MPI-IO)
     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  num_surf            !< number of surface data at (j,i)
     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  start_index         !< start index of surface data at (j,i)
@@ -4600 +4604 @@
     ALLOCATE( surf_in(1:surfaces%ns) )
 
-    CALL rd_mpi_io_check_array( 'surfaces%start_index', found = array_found )
-    IF ( array_found )  CALL rrd_mpi_io( 'surfaces%start_index', start_index )
-
-    CALL rd_mpi_io_check_array( 'surfaces%end_index', found = array_found )
-    IF ( array_found )  CALL rrd_mpi_io( 'surfaces%end_index', end_index )
-
     CALL rd_mpi_io_check_array( 'surfaces%global_start_index', found = array_found )
     IF ( array_found )  CALL rrd_mpi_io( 'surfaces%global_start_index', global_start_index )
 
-    CALL rd_mpi_io_surface_filetypes( start_index, end_index, ldum, global_start_index )
+    CALL rd_mpi_io_check_array( 'surfaces%global_end_index', found = array_found )
+    IF ( array_found )  CALL rrd_mpi_io( 'surfaces%global_end_index', global_end_index )
+
+    CALL rd_mpi_io_surface_filetypes( start_index, end_index, ldum, global_start_index,            &
+                                      global_end_index )
 
     DO  nv = 1, dosurf_no(1)
@@ -4728 +4730 @@
 
     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  end_index           !< end index of surface data at (j,i)
-    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  global_start_index  !< index array for surface data (MPI-IO)
+    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  global_end_index    !< end index array for surface data (MPI-IO)
+    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  global_start_index  !< start index array for surface data (MPI-IO)
     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  num_surf            !< number of surface data at (j,i)
     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  start_index         !< start index of surface data at (j,i)
@@ -4811 +4814 @@
 
           CALL rd_mpi_io_surface_filetypes( start_index, end_index, surface_data_to_write,         &
-                                            global_start_index )
-          CALL wrd_mpi_io( 'surfaces%start_index',        start_index        )
-          CALL wrd_mpi_io( 'surfaces%end_index',          end_index          )
+                                            global_start_index, global_end_index )
+
           CALL wrd_mpi_io( 'surfaces%global_start_index', global_start_index )
+          CALL wrd_mpi_io( 'surfaces%global_end_index', global_end_index )
 
           DO  nv = 1, dosurf_no(1)