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

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! added restart with MPI-IO for reading local arrays
+! 
+! 4510 2020-04-29 14:19:18Z raasch
 ! Further re-formatting to follow the PALM coding standard
 !
@@ -368 +371 @@
     USE restart_data_mpi_io_mod,                                                                   &
         ONLY:  rd_mpi_io_surface_filetypes,                                                        &
+               rrd_mpi_io,                                                                         &
+               rrd_mpi_io_surface,                                                                 &
                wrd_mpi_io,                                                                         &
                wrd_mpi_io_surface
@@ -617 +622 @@
 
 
-
-
 !
 !-- Building facade/wall/green/window properties (partly according to PIDS).
@@ -840 +843 @@
 
     INTERFACE usm_rrd_local
-       MODULE PROCEDURE usm_rrd_local
+       MODULE PROCEDURE usm_rrd_local_ftn
+       MODULE PROCEDURE usm_rrd_local_mpi
     END INTERFACE usm_rrd_local
 
@@ -6110 +6114 @@
 ! Description:
 ! ------------
-!> Soubroutine reads t_surf and t_wall data from restart files
+!> Read module-specific local restart data arrays (Fortran binary format).
+!> Soubroutine reads t_surf and t_wall.
 !--------------------------------------------------------------------------------------------------!
- SUBROUTINE usm_rrd_local( k, nxlf, nxlc, nxl_on_file, nxrf, nxr_on_file, nynf, nyn_on_file,       &
-                           nysf, nysc, nys_on_file, found )
+ SUBROUTINE usm_rrd_local_ftn( k, nxlf, nxlc, nxl_on_file, nxrf, nxr_on_file, nynf, nyn_on_file,   &
+                               nysf, nysc, nys_on_file, found )
 
 
@@ -6654 +6659 @@
     END SELECT
 
- END SUBROUTINE usm_rrd_local
+ END SUBROUTINE usm_rrd_local_ftn
+
+
+!--------------------------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Read module-specific local restart data arrays (MPI-IO).
+!> Soubroutine reads t_surf and t_wall.
+!>
+!> This read routine is a counterpart of usm_wrd_local.
+!> In usm_wrd_local, all array are unconditionally written, therefore all arrays are read here.
+!> This is a preliminary version of reading usm data. The final version has to be discussed with
+!> the developers.
+!>
+!> If it is possible to call usm_allocate_surface before reading the restart file, this reading
+!> routine would become much simpler, because no checking for allocation will be necessary any more.
+!--------------------------------------------------------------------------------------------------!
+ SUBROUTINE usm_rrd_local_mpi
+
+
+    CHARACTER(LEN=1) ::  dum  !< dummy string to create input-variable name
+
+    INTEGER(iwp) ::  l  !< loop index for surface types
+
+    INTEGER(iwp), DIMENSION(nys:nyn,nxl:nxr) ::  global_start
+
+    LOGICAL ::  ldum  !< dummy variable
+
+
+    CALL rrd_mpi_io( 'usm_start_index_h',  surf_usm_h%start_index )
+    CALL rrd_mpi_io( 'usm_end_index_h', surf_usm_h%end_index )
+    CALL rrd_mpi_io( 'usm_global_start_h', global_start )
+
+    CALL rd_mpi_io_surface_filetypes( surf_usm_h%start_index, surf_usm_h%end_index, ldum,          &
+                                      global_start )
+
+    IF ( .NOT.  ALLOCATED( t_surf_wall_h_1 ) )  ALLOCATE( t_surf_wall_h_1(1:surf_usm_h%ns) )
+    CALL rrd_mpi_io_surface( 't_surf_wall_h', t_surf_wall_h_1 )
+
+    IF ( .NOT.  ALLOCATED( t_surf_window_h_1 ) )  ALLOCATE( t_surf_window_h_1(1:surf_usm_h%ns) )
+    CALL rrd_mpi_io_surface( 't_surf_window_h', t_surf_window_h_1 )
+
+    IF ( .NOT.  ALLOCATED( t_surf_green_h_1 ) )  ALLOCATE( t_surf_green_h_1(1:surf_usm_h%ns) )
+    CALL rrd_mpi_io_surface( 't_surf_green_h', t_surf_green_h_1 )
+
+    DO  l = 0, 3
+
+       WRITE( dum, '(I1)' )  l
+
+       CALL rrd_mpi_io( 'usm_start_index_v_' //dum, surf_usm_v(l)%start_index )
+       CALL rrd_mpi_io( 'usm_end_index_v_' // dum, surf_usm_v(l)%end_index )
+       CALL rrd_mpi_io( 'usm_global_start_v_' // dum, global_start )
+
+       CALL rd_mpi_io_surface_filetypes( surf_usm_v(l)%start_index, surf_usm_v(l)%end_index, ldum, &
+                                         global_start )
+
+       IF ( .NOT.  ALLOCATED( t_surf_wall_v_1(l)%t ) )                                             &
+          ALLOCATE( t_surf_wall_v_1(l)%t(1:surf_usm_v(l)%ns) )
+       CALL rrd_mpi_io_surface( 't_surf_wall_v(' // dum // ')', t_surf_wall_v_1(l)%t )
+
+       IF ( .NOT.  ALLOCATED( t_surf_window_v_1(l)%t ) )                                           &
+          ALLOCATE( t_surf_window_v_1(l)%t(1:surf_usm_v(l)%ns) )
+       CALL rrd_mpi_io_surface( 't_surf_window_v(' // dum // ')', t_surf_window_v_1(l)%t )
+
+       IF ( .NOT.  ALLOCATED( t_surf_green_v_1(l)%t ) )                                            &
+          ALLOCATE( t_surf_green_v_1(l)%t(1:surf_usm_v(l)%ns) )
+       CALL rrd_mpi_io_surface( 't_surf_green_v(' // dum // ')', t_surf_green_v_1(l)%t)
+
+    ENDDO
+
+    CALL rrd_mpi_io( 'usm_start_index_h_2',  surf_usm_h%start_index )
+    CALL rrd_mpi_io( 'usm_end_index_h_2', surf_usm_h%end_index )
+    CALL rrd_mpi_io( 'usm_global_start_h_2', global_start )
+
+    CALL rd_mpi_io_surface_filetypes( surf_usm_h%start_index, surf_usm_h%end_index, ldum,          &
+                                      global_start )
+
+    IF ( .NOT.  ALLOCATED( t_wall_h_1 ) )                                                          &
+       ALLOCATE( t_wall_h_1(nzb_wall:nzt_wall+1,1:surf_usm_h%ns) )
+    CALL rrd_mpi_io_surface( 't_wall_h', t_wall_h_1 )
+
+    IF ( .NOT.  ALLOCATED( t_window_h_1 ) )                                                        &
+       ALLOCATE( t_window_h_1(nzb_wall:nzt_wall+1,1:surf_usm_h%ns) )
+    CALL rrd_mpi_io_surface( 't_window_h', t_window_h_1 )
+
+    IF ( .NOT.  ALLOCATED( t_green_h_1 ) )                                                         &
+       ALLOCATE( t_green_h_1(nzb_wall:nzt_wall+1,1:surf_usm_h%ns) )
+    CALL rrd_mpi_io_surface( 't_green_h', t_green_h_1 )
+
+    DO  l = 0, 3
+
+       WRITE( dum, '(I1)' )  l
+
+       CALL rrd_mpi_io( 'usm_start_index_v_2_' //dum, surf_usm_v(l)%start_index )
+       CALL rrd_mpi_io( 'usm_end_index_v_2_' // dum, surf_usm_v(l)%end_index )
+       CALL rrd_mpi_io( 'usm_global_start_v_2_' // dum, global_start )
+
+       CALL rd_mpi_io_surface_filetypes( surf_usm_v(l)%start_index, surf_usm_v(l)%end_index, ldum, &
+                                         global_start )
+
+       IF ( .NOT. ALLOCATED( t_wall_v_1(l)%t ) )                                                   &
+          ALLOCATE ( t_wall_v_1(l)%t(nzb_wall:nzt_wall+1,1:surf_usm_v(l)%ns) )
+       CALL rrd_mpi_io_surface( 't_wall_v(' // dum // ')', t_wall_v_1(l)%t )
+
+       IF ( .NOT. ALLOCATED( t_window_v_1(l)%t ) )                                                 &
+          ALLOCATE ( t_window_v_1(l)%t(nzb_wall:nzt_wall+1,1:surf_usm_v(l)%ns) )
+       CALL rrd_mpi_io_surface( 't_window_v(' // dum // ')', t_window_v_1(l)%t )
+
+       IF ( .NOT. ALLOCATED( t_green_v_1(l)%t ) )                                                  &
+          ALLOCATE ( t_green_v_1(l)%t(nzb_wall:nzt_wall+1,1:surf_usm_v(l)%ns) )
+       CALL rrd_mpi_io_surface( 't_green_v(' // dum // ')', t_green_v_1(l)%t )
+
+    ENDDO
+
+ END SUBROUTINE usm_rrd_local_mpi
+