Changeset 4517 for palm/trunk/SOURCE/ocean_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/ocean_mod.f90 (modified) (5 diffs)

palm/trunk/SOURCE/ocean_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! added restart with MPI-IO for reading local arrays
+! 
+! 4495 2020-04-13 20:11:20Z raasch
 ! restart data handling with MPI-IO added
 ! 
@@ -102 +105 @@
 
     USE restart_data_mpi_io_mod,                                                                   &
-        ONLY:  rrd_mpi_io, rrd_mpi_io_global_array, wrd_mpi_io, wrd_mpi_io_global_array
+        ONLY:  rd_mpi_io_check_array, rrd_mpi_io, rrd_mpi_io_global_array, wrd_mpi_io,             &
+               wrd_mpi_io_global_array
 
     USE surface_mod,                                                           &
@@ -232 +236 @@
 
     INTERFACE ocean_rrd_local
-       MODULE PROCEDURE ocean_rrd_local
+       MODULE PROCEDURE ocean_rrd_local_ftn
+       MODULE PROCEDURE ocean_rrd_local_mpi
     END INTERFACE ocean_rrd_local
 
@@ -2039 +2044 @@
 ! Description:
 ! ------------
-!> This routine reads the respective restart data for the ocean module.
-!------------------------------------------------------------------------------!
- SUBROUTINE ocean_rrd_local( k, nxlf, nxlc, nxl_on_file, nxrf, nxrc,           &
-                             nxr_on_file, nynf, nync, nyn_on_file, nysf,       &
-                             nysc, nys_on_file, tmp_3d, found )
+!> Read module-specific local restart data arrays (Fortran binary format).
+!------------------------------------------------------------------------------!
+ SUBROUTINE ocean_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_3d, found )
 
     USE averaging,                                                             &
@@ -2108 +2113 @@
     END SELECT
 
- END SUBROUTINE ocean_rrd_local
+ END SUBROUTINE ocean_rrd_local_ftn
+
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Read module-specific local restart data arrays (MPI-IO).
+!------------------------------------------------------------------------------!
+ SUBROUTINE ocean_rrd_local_mpi
+
+    USE averaging,                                                             &
+        ONLY:  rho_ocean_av, sa_av
+
+    USE indices,                                                               &
+        ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzt
+
+    IMPLICIT NONE
+
+    LOGICAL ::  array_found  !<
+
+
+    CALL rd_mpi_io_check_array( 'rho_ocean_av' , found = array_found )
+    IF ( array_found )  THEN
+       IF ( .NOT. ALLOCATED( rho_ocean_av ) )  ALLOCATE( rho_ocean_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+       CALL rrd_mpi_io( 'rho_ocean_av', rho_ocean_av )
+    ENDIF
+
+    CALL rrd_mpi_io( 'sa', sa )
+
+    CALL rd_mpi_io_check_array( 'sa_av' , found = array_found )
+    IF ( array_found )  THEN
+       IF ( .NOT. ALLOCATED( sa_av ) )  ALLOCATE( sa_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+       CALL rrd_mpi_io( 'sa_av', sa_av )
+    ENDIF
+
+ END SUBROUTINE ocean_rrd_local_mpi