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

palm/trunk/SOURCE/plant_canopy_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! added restart with MPI-IO for reading local arrays
+! 
+! 4515 2020-04-30 16:37:18Z suehring
 ! Rename error number again since this was recently given in -r 4511
 ! 
@@ -237 +240 @@
 
     USE restart_data_mpi_io_mod,                                                                   &
-        ONLY:  rrd_mpi_io,                                                                         &
+        ONLY:  rd_mpi_io_check_array,                                                              &
+               rrd_mpi_io,                                                                         &
                wrd_mpi_io
 
@@ -361 +365 @@
 
     INTERFACE pcm_rrd_local
-       MODULE PROCEDURE pcm_rrd_local
+       MODULE PROCEDURE pcm_rrd_local_ftn
+       MODULE PROCEDURE pcm_rrd_local_mpi
     END INTERFACE pcm_rrd_local
 
@@ -1685 +1690 @@
     END SUBROUTINE pcm_rrd_global_mpi
 
+
 !------------------------------------------------------------------------------!
 ! Description:
 ! ------------
-!> Subroutine reads local (subdomain) restart data
-!------------------------------------------------------------------------------!
-    SUBROUTINE pcm_rrd_local( k, nxlf, nxlc, nxl_on_file, nxrf, nxrc,          &
-                              nxr_on_file, nynf, nync, nyn_on_file, nysf,      &
-                              nysc, nys_on_file, found )
+!> Read module-specific local restart data arrays (Fortran binary format).
+!------------------------------------------------------------------------------!
+    SUBROUTINE pcm_rrd_local_ftn( k, nxlf, nxlc, nxl_on_file, nxrf, nxrc,          &
+                                  nxr_on_file, nynf, nync, nyn_on_file, nysf,      &
+                                  nysc, nys_on_file, found )
 
        INTEGER(iwp) ::  k               !< 
@@ -1751 +1757 @@
        END SELECT
 
-    END SUBROUTINE pcm_rrd_local
+    END SUBROUTINE pcm_rrd_local_ftn
+
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Read module-specific local restart data arrays (MPI-IO).
+!------------------------------------------------------------------------------!
+    SUBROUTINE pcm_rrd_local_mpi
+
+       IMPLICIT NONE
+
+       LOGICAL ::  array_found  !<
+
+
+!>     TODO: following code does not work because arrays are sized 0:pch_index along k!!!!!
+       CALL rd_mpi_io_check_array( 'pcm_heatrate_av' , found = array_found )
+       IF ( array_found )  THEN
+          IF ( .NOT. ALLOCATED( pcm_heatrate_av ) )  ALLOCATE( pcm_heatrate_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+          CALL rrd_mpi_io( 'pcm_heatrate_av', pcm_heatrate_av )
+       ENDIF
+
+       CALL rd_mpi_io_check_array( 'pcm_latentrate_av' , found = array_found )
+       IF ( array_found )  THEN
+          IF ( .NOT. ALLOCATED( pcm_latentrate_av ) )  ALLOCATE( pcm_latentrate_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+          CALL rrd_mpi_io( 'pcm_latentrate_av', pcm_latentrate_av )
+       ENDIF
+
+       CALL rd_mpi_io_check_array( 'pcm_transpirationrate_av' , found = array_found )
+       IF ( array_found )  THEN
+          IF ( .NOT. ALLOCATED( pcm_transpirationrate_av ) )  ALLOCATE( pcm_transpirationrate_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+          CALL rrd_mpi_io( 'pcm_transpirationrate_av', pcm_transpirationrate_av )
+       ENDIF
+
+    END SUBROUTINE pcm_rrd_local_mpi
+
 
 !------------------------------------------------------------------------------!