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

palm/trunk/SOURCE/lagrangian_particle_model_mod.f90

@@ -210 +210 @@
 
     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 statistics,                                                            &
@@ -413 +413 @@
 
     INTERFACE lpm_rrd_local
-       MODULE PROCEDURE lpm_rrd_local
+       MODULE PROCEDURE lpm_rrd_local_ftn
+       MODULE PROCEDURE lpm_rrd_local_mpi
     END INTERFACE lpm_rrd_local
 
@@ -3021 +3022 @@
  
  
- SUBROUTINE lpm_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 )
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Read module-specific local restart data arrays (Fortran binary format).
+!------------------------------------------------------------------------------!
+ SUBROUTINE lpm_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 )
 
 
@@ -3101 +3107 @@
        END SELECT
 
-
- END SUBROUTINE lpm_rrd_local
+ END SUBROUTINE lpm_rrd_local_ftn
+
  
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Read module-specific local restart data arrays (MPI-IO).
+!------------------------------------------------------------------------------!
+ SUBROUTINE lpm_rrd_local_mpi
+
+    IMPLICIT NONE
+
+    LOGICAL ::  array_found  !<
+
+    CALL rrd_mpi_io( 'iran', iran ) ! matching random numbers is still unresolved issue
+    CALL rrd_mpi_io( 'iran_part', iran_part )
+
+    CALL rd_mpi_io_check_array( 'pc_av' , found = array_found )
+    IF ( array_found )  THEN
+       IF ( .NOT. ALLOCATED( pc_av ) )  ALLOCATE( pc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+       CALL rrd_mpi_io( 'pc_av', pc_av )
+    ENDIF
+
+    CALL rd_mpi_io_check_array( 'pr_av' , found = array_found )
+    IF ( array_found )  THEN
+       IF ( .NOT. ALLOCATED( pr_av ) )  ALLOCATE( pr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+       CALL rrd_mpi_io( 'pr_av', pr_av )
+    ENDIF
+
+    CALL rd_mpi_io_check_array( 'ql_c_av' , found = array_found )
+    IF ( array_found )  THEN
+       IF ( .NOT. ALLOCATED( ql_c_av ) )  ALLOCATE( ql_c_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+       CALL rrd_mpi_io( 'ql_c_av', ql_c_av )
+    ENDIF
+
+    CALL rd_mpi_io_check_array( 'ql_v_av' , found = array_found )
+    IF ( array_found )  THEN
+       IF ( .NOT. ALLOCATED( ql_v_av ) )  ALLOCATE( ql_v_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+       CALL rrd_mpi_io( 'ql_v_av', ql_v_av )
+    ENDIF
+
+    CALL rd_mpi_io_check_array( 'ql_vp_av' , found = array_found )
+    IF ( array_found )  THEN
+       IF ( .NOT. ALLOCATED( ql_vp_av ) )  ALLOCATE( ql_vp_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+       CALL rrd_mpi_io( 'ql_vp_av', ql_vp_av )
+    ENDIF
+
+ END SUBROUTINE lpm_rrd_local_mpi
+
+
 !------------------------------------------------------------------------------!
 ! Description:
@@ -3168 +3221 @@
 #endif
 
-    CALL wrd_write_string( 'iran' ) 
-    WRITE ( 14 )  iran, iran_part 
-
+    IF ( TRIM( restart_data_format_output ) == 'fortran_binary' )  THEN
+
+       CALL wrd_write_string( 'iran' )  ! matching random numbers is still unresolved issue
+       WRITE ( 14 )  iran, iran_part
+
+    ELSEIF ( TRIM( restart_data_format_output ) == 'mpi' )  THEN
+
+       CALL wrd_mpi_io( 'iran', iran )  ! matching random numbers is still unresolved issue
+       CALL wrd_mpi_io( 'iran_part', iran_part )
+
+    ENDIF
 
  END SUBROUTINE lpm_wrd_local