Changeset 4517 for palm/trunk/SOURCE/bulk_cloud_model_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/bulk_cloud_model_mod.f90 (modified) (6 diffs)

palm/trunk/SOURCE/bulk_cloud_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! added restart with MPI-IO for reading local arrays
+! 
+! 4506 2020-04-21 10:57:45Z schwenkel
 ! Use correct magnus formula for liquid water temperature
 ! 
@@ -183 +186 @@
 
     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,                                                            &
@@ -404 +407 @@
 
     INTERFACE bcm_rrd_local
-       MODULE PROCEDURE bcm_rrd_local
+       MODULE PROCEDURE bcm_rrd_local_ftn
+       MODULE PROCEDURE bcm_rrd_local_mpi
     END INTERFACE bcm_rrd_local
 
@@ -3397 +3401 @@
 ! Description:
 ! ------------
-!> This routine reads the respective restart data for the bulk cloud module.
-!------------------------------------------------------------------------------!
-    SUBROUTINE bcm_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 )
+!> Read module-specific local restart data arrays (Fortran binary format).
+!------------------------------------------------------------------------------!
+    SUBROUTINE bcm_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 )
 
 
@@ -3432 +3436 @@
        REAL(wp), DIMENSION(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp,nxl_on_file-nbgp:nxr_on_file+nbgp) :: tmp_3d   !<
 
-!
-!-- Here the reading of user-defined restart data follows:
-!-- Sample for user-defined output
-
 
        found = .TRUE.
@@ -3561 +3561 @@
 
 
-    END SUBROUTINE bcm_rrd_local
+    END SUBROUTINE bcm_rrd_local_ftn
+
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Read module-specific local restart data arrays (MPI-IO).
+!------------------------------------------------------------------------------!
+    SUBROUTINE bcm_rrd_local_mpi
+
+       LOGICAL ::  array_found  !<
+
+
+       CALL rd_mpi_io_check_array( 'prr' , found = array_found )
+       IF ( array_found )  THEN
+          IF ( .NOT. ALLOCATED( prr ) )  ALLOCATE( prr(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+          CALL rrd_mpi_io( 'prr', prr )
+       ENDIF
+
+       CALL rd_mpi_io_check_array( 'prr_av' , found = array_found )
+       IF ( array_found )  THEN
+          IF ( .NOT. ALLOCATED( prr_av ) )  ALLOCATE( prr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+          CALL rrd_mpi_io( 'prr_av', prr_av )
+       ENDIF
+
+       CALL rd_mpi_io_check_array( 'precipitation_amount' , found = array_found )
+       IF ( array_found )  THEN
+          IF ( .NOT. ALLOCATED( precipitation_amount ) )  THEN
+             ALLOCATE( precipitation_amount(nysg:nyng,nxlg:nxrg) )
+          ENDIF
+          CALL rrd_mpi_io( 'precipitation_amount', precipitation_amount )
+       ENDIF
+
+       CALL rrd_mpi_io( 'ql', ql )
+
+       CALL rd_mpi_io_check_array( 'ql_av' , found = array_found )
+       IF ( array_found )  THEN
+          IF ( .NOT. ALLOCATED( ql_av ) )  ALLOCATE( ql_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+          CALL rrd_mpi_io( 'ql_av', ql_av )
+       ENDIF
+
+       CALL rrd_mpi_io( 'qc', qc )
+
+       CALL rd_mpi_io_check_array( 'qc_av' , found = array_found )
+       IF ( array_found )  THEN
+          IF ( .NOT. ALLOCATED( qc_av ) )  ALLOCATE( qc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+          CALL rrd_mpi_io( 'qc_av', qc_av )
+       ENDIF
+
+       IF ( microphysics_morrison )  THEN
+
+          CALL rrd_mpi_io( 'nc', nc )
+
+          CALL rd_mpi_io_check_array( 'nc_av' , found = array_found )
+          IF ( array_found )  THEN
+             IF ( .NOT. ALLOCATED( nc_av ) )  ALLOCATE( nc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+             CALL rrd_mpi_io( 'nc_av', nc_av )
+          ENDIF
+
+       ENDIF
+
+       IF ( microphysics_seifert )  THEN
+
+          CALL rrd_mpi_io( 'nr', nr )
+
+          CALL rd_mpi_io_check_array( 'nr_av' , found = array_found )
+          IF ( array_found )  THEN
+             IF ( .NOT. ALLOCATED( nr_av ) )  ALLOCATE( nr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+             CALL rrd_mpi_io( 'nr_av', nr_av )
+          ENDIF
+
+          CALL rrd_mpi_io( 'qr', qr )
+
+          CALL rd_mpi_io_check_array( 'qr_av' , found = array_found )
+          IF ( array_found )  THEN
+             IF ( .NOT. ALLOCATED( qr_av ) )  ALLOCATE( qr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+             CALL rrd_mpi_io( 'qr_av', qr_av )
+          ENDIF
+
+       ENDIF
+
+       IF ( microphysics_ice_extension )  THEN
+
+          CALL rrd_mpi_io( 'ni', ni )
+
+          CALL rd_mpi_io_check_array( 'ni_av' , found = array_found )
+          IF ( array_found )  THEN
+             IF ( .NOT. ALLOCATED( ni_av ) )  ALLOCATE( ni_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+             CALL rrd_mpi_io( 'ni_av', ni_av )
+          ENDIF
+
+          CALL rrd_mpi_io( 'qi', qi )
+
+          CALL rd_mpi_io_check_array( 'qi_av' , found = array_found )
+          IF ( array_found )  THEN
+             IF ( .NOT. ALLOCATED( qi_av ) )  ALLOCATE( qi_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+             CALL rrd_mpi_io( 'qi_av', qi_av )
+          ENDIF
+
+       ENDIF
+
+
+    END SUBROUTINE bcm_rrd_local_mpi