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

Timestamp:

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

Author:

raasch

Message:

added restart with MPI-IO for reading local arrays

File:

• palm/trunk/SOURCE/dynamics_mod.f90 (modified) (4 diffs)

palm/trunk/SOURCE/dynamics_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! added restart with MPI-IO for reading local arrays
+! 
+! 4505 2020-04-20 15:37:15Z schwenkel
 ! Add flag for saturation check
 ! 
@@ -279 +282 @@
 
     INTERFACE dynamics_rrd_local
-       MODULE PROCEDURE dynamics_rrd_local
+       MODULE PROCEDURE dynamics_rrd_local_ftn
+       MODULE PROCEDURE dynamics_rrd_local_mpi
     END INTERFACE dynamics_rrd_local
 
@@ -1871 +1875 @@
 ! Description:
 ! ------------
-!> Read module-specific processor specific restart data from file(s).
+!> Read module-specific local restart data arrays (Fortran binary format).
 !> Subdomain index limits on file are given by nxl_on_file, etc.
 !> Indices nxlc, etc. indicate the range of gridpoints to be mapped from the subdomain on file (f)
 !> to the subdomain of the current PE (c). They have been calculated in routine rrd_local.
 !--------------------------------------------------------------------------------------------------!
- SUBROUTINE dynamics_rrd_local( k, nxlf, nxlc, nxl_on_file, nxrf, nxrc, nxr_on_file, nynf, nync,   &
-                                nyn_on_file, nysf, nysc, nys_on_file, tmp_2d, tmp_3d, found )
+ SUBROUTINE dynamics_rrd_local_ftn( k, nxlf, nxlc, nxl_on_file, nxrf, nxrc, nxr_on_file, nynf, nync,   &
+                                    nyn_on_file, nysf, nysc, nys_on_file, tmp_2d, tmp_3d, found )
 
 
@@ -1920 +1924 @@
     END SELECT
 
- END SUBROUTINE dynamics_rrd_local
+ END SUBROUTINE dynamics_rrd_local_ftn
+
+
+!--------------------------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Read module-specific local restart data arrays (MPI-IO).
+!--------------------------------------------------------------------------------------------------!
+ SUBROUTINE dynamics_rrd_local_mpi
+
+    IMPLICIT NONE
+
+!    LOGICAL ::  array_found  !<
+
+
+!    CALL rd_mpi_io_check_array( 'u2_av' , found = array_found )
+!    IF ( array_found )  THEN
+!       IF ( .NOT. ALLOCATED( u2_av ) )  ALLOCATE( u2_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+!       CALL rrd_mpi_io( 'u2_av', u2_av )
+!    ENDIF
+
+    CONTINUE
+
+ END SUBROUTINE dynamics_rrd_local_mpi
+