Changeset 4694 for palm/trunk/SOURCE/urban_surface_mod.f90

Timestamp:

Sep 23, 2020 3:09:19 PM (4 years ago)

Author:

pavelkrc

Message:

Fix writing and reading of surface data to/from MPI restart file + small fixes

File:

• palm/trunk/SOURCE/urban_surface_mod.f90 (modified) (10 diffs)

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Fix writing and reading of surface data to/from MPI restart file
+!
+! 4693 2020-09-22 19:47:04Z maronga
 ! Bugfix for last commit
 !
@@ -6498 +6501 @@
                                          global_start )
 
+       IF ( MAXVAL( surf_usm_h(l)%end_index ) <= 0 )  CYCLE
+
        IF ( .NOT.  ALLOCATED( t_surf_wall_h_1(l)%val ) )                                             &
           ALLOCATE( t_surf_wall_h_1(l)%val(1:surf_usm_h(l)%ns) )
@@ -6510 +6515 @@
        CALL rrd_mpi_io_surface( 't_surf_green_h(' // dum // ')', t_surf_green_h_1(l)%val )
 
+       IF ( .NOT.  ALLOCATED( m_liq_usm_h_1(l)%val ) )                                            &
+          ALLOCATE( m_liq_usm_h_1(l)%val(1:surf_usm_h(l)%ns) )
+       CALL rrd_mpi_io_surface( 'm_liq_usm_h(' // dum // ')', m_liq_usm_h_1(l)%val )
+
+       IF ( indoor_model )  THEN
+          IF ( .NOT.  ALLOCATED( surf_usm_h(l)%waste_heat ) )                                     &
+             ALLOCATE( surf_usm_h(l)%waste_heat(1:surf_usm_h(l)%ns) )
+          CALL rrd_mpi_io_surface( 'waste_heat_h(' // dum // ')', surf_usm_h(l)%waste_heat )
+       ENDIF
+
     ENDDO
     DO  l = 0, 3
@@ -6522 +6537 @@
                                          global_start )
 
+       IF ( MAXVAL( surf_usm_v(l)%end_index ) <= 0 )  CYCLE
+
        IF ( .NOT.  ALLOCATED( t_surf_wall_v_1(l)%val ) )                                             &
           ALLOCATE( t_surf_wall_v_1(l)%val(1:surf_usm_v(l)%ns) )
@@ -6547 +6564 @@
                                          global_start )
 
+       IF ( MAXVAL( surf_usm_h(l)%end_index ) <= 0 )  CYCLE
+
        IF ( .NOT.  ALLOCATED( t_wall_h_1(l)%val ) )                                                          &
           ALLOCATE( t_wall_h_1(l)%val(nzb_wall:nzt_wall+1,1:surf_usm_h(l)%ns) )
@@ -6558 +6577 @@
           ALLOCATE( t_green_h_1(l)%val(nzb_wall:nzt_wall+1,1:surf_usm_h(l)%ns) )
        CALL rrd_mpi_io_surface( 't_green_h(' // dum // ')', t_green_h_1(l)%val )
+
     ENDDO
 
@@ -6570 +6590 @@
        CALL rd_mpi_io_surface_filetypes( surf_usm_v(l)%start_index, surf_usm_v(l)%end_index, ldum, &
                                          global_start )
+
+       IF ( MAXVAL( surf_usm_v(l)%end_index ) <= 0 )  CYCLE
 
        IF ( .NOT. ALLOCATED( t_wall_v_1(l)%val ) )                                                   &
@@ -7392 +7414 @@
           WRITE( dum, '(I1)')  l
 
-          END DO
           CALL rd_mpi_io_surface_filetypes( surf_usm_h(l)%start_index, surf_usm_h(l)%end_index,             &
                                             surface_data_to_write, global_start_index )
@@ -7400 +7421 @@
           CALL wrd_mpi_io( 'usm_global_start_h_' // dum, global_start_index )
 
-          CALL wrd_mpi_io_surface( 't_surf_wall_h_' // dum,  t_surf_wall_h(l)%val )
-          CALL wrd_mpi_io_surface( 't_surf_window_h_' // dum, t_surf_window_h(l)%val )
-          CALL wrd_mpi_io_surface( 't_surf_green_h_' // dum, t_surf_green_h(l)%val )
-
-          CALL wrd_mpi_io_surface( 'm_liq_usm_h_' // dum, m_liq_usm_h(l)%val )
+          IF ( .NOT. surface_data_to_write )  CYCLE
+
+          CALL wrd_mpi_io_surface( 't_surf_wall_h(' // dum // ')',  t_surf_wall_h(l)%val )
+          CALL wrd_mpi_io_surface( 't_surf_window_h(' // dum // ')', t_surf_window_h(l)%val )
+          CALL wrd_mpi_io_surface( 't_surf_green_h(' // dum // ')', t_surf_green_h(l)%val )
+
+          CALL wrd_mpi_io_surface( 'm_liq_usm_h(' // dum // ')', m_liq_usm_h(l)%val )
           IF ( indoor_model )  THEN
-             CALL wrd_mpi_io_surface( 'waste_heat_h_' // dum, surf_usm_h(l)%waste_heat ) ! NEED TO BE CHECKED!!!!!
-          ENDIF
+             CALL wrd_mpi_io_surface( 'waste_heat_h(' // dum // ')', surf_usm_h(l)%waste_heat ) ! NEED TO BE CHECKED!!!!!
+          ENDIF
+
+       ENDDO
 
        DO  l = 0, 3
@@ -7439 +7464 @@
           CALL wrd_mpi_io( 'usm_global_start_h_2_' // dum, global_start_index )
 
-          CALL wrd_mpi_io_surface( 't_wall_h_' // dum, t_wall_h(l)%val )
-          CALL wrd_mpi_io_surface( 't_window_h_' // dum, t_window_h(l)%val )
-          CALL wrd_mpi_io_surface( 't_green_h_' // dum, t_green_h(l)%val )
+          IF ( .NOT. surface_data_to_write )  CYCLE
+
+          CALL wrd_mpi_io_surface( 't_wall_h(' // dum // ')', t_wall_h(l)%val )
+          CALL wrd_mpi_io_surface( 't_window_h(' // dum // ')', t_window_h(l)%val )
+          CALL wrd_mpi_io_surface( 't_green_h(' // dum // ')', t_green_h(l)%val )
 
        ENDDO