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

Timestamp:

May 3, 2020 2:29:30 PM (5 years ago)

Author:

raasch

Message:

added restart with MPI-IO for reading local arrays

File:

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

palm/trunk/SOURCE/surface_data_output_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! added restart with MPI-IO for reading local arrays
+! 
+! 4502 2020-04-17 16:14:16Z schwenkel
 ! Implementation of ice microphysics
 ! 
@@ -147 +150 @@
 
    USE restart_data_mpi_io_mod,                                                                    &
-       ONLY:  rrd_mpi_io, wrd_mpi_io
+       ONLY:  rd_mpi_io_check_array, rrd_mpi_io, wrd_mpi_io
 
    USE surface_mod,                                                                                &
@@ -267 +270 @@
 
    INTERFACE  surface_data_output_rrd_local
-      MODULE PROCEDURE surface_data_output_rrd_local
+      MODULE PROCEDURE surface_data_output_rrd_local_ftn
+      MODULE PROCEDURE surface_data_output_rrd_local_mpi
    END INTERFACE  surface_data_output_rrd_local
 
@@ -4883 +4887 @@
 ! Description:
 ! ------------
-!> This routine reads the respective restart data.
+!> Read module-specific local restart data arrays (Fortran binary format).
 !------------------------------------------------------------------------------!
-    SUBROUTINE surface_data_output_rrd_local( found )
+    SUBROUTINE surface_data_output_rrd_local_ftn( found )
 
 
@@ -4895 +4899 @@
        LOGICAL, INTENT(OUT)  ::  found
 
-!
-!--    Here the reading of user-defined restart data follows:
-!--    Sample for user-defined output
+
        found = .TRUE.
 
@@ -4912 +4914 @@
 
 
-    END SUBROUTINE surface_data_output_rrd_local
+    END SUBROUTINE surface_data_output_rrd_local_ftn
+
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Read module-specific local restart data arrays (MPI-IO).
+!------------------------------------------------------------------------------!
+    SUBROUTINE surface_data_output_rrd_local_mpi
+
+       IMPLICIT NONE
+
+       LOGICAL ::  array_found  !<
+
+
+       CALL rd_mpi_io_check_array( 'surfaces%var_av' , found = array_found )
+
+!> does not work this way: surface%var_av has non-standard dimensions
+!       IF ( array_found )  THEN
+!          IF ( .NOT. ALLOCATED( surfaces%var_av ) )  ALLOCATE( ....... )
+!          CALL rrd_mpi_io( 'surfaces%var_av', surfaces%var_av )
+!       ENDIF
+
+    END SUBROUTINE surface_data_output_rrd_local_mpi
+
 
 !------------------------------------------------------------------------------!
@@ -4936 +4962 @@
     END SUBROUTINE surface_data_output_wrd_global
 
+
 !------------------------------------------------------------------------------!
 ! Description: