Changeset 4850 for palm/trunk/SOURCE/indoor_model_mod.f90

Timestamp:

Jan 21, 2021 5:59:25 PM (3 years ago)

Author:

suehring

Message:

Bugfix in indoor model: consider previous indoor temperature during restarts; Further bugfix in mpi-io restart mechanism for the waste-heat flux from buildings

File:

• palm/trunk/SOURCE/indoor_model_mod.f90 (modified) (7 diffs)

palm/trunk/SOURCE/indoor_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Enable restart mechanism for previous indoor temperature
+!
+! 4843 2021-01-15 15:22:11Z raasch
 ! local namelist parameter added to switch off the module although the respective module namelist
 ! appears in the namelist file
@@ -590 +593 @@
     INTEGER(iwp) ::  bt          !< local building type
     INTEGER(iwp) ::  day_of_year !< day of the year
+    INTEGER(iwp) ::  fa          !< running index for facade elements of each building
     INTEGER(iwp) ::  i           !< running index along x-direction
     INTEGER(iwp) ::  j           !< running index along y-direction
@@ -1185 +1189 @@
 !
 !-- Initialize indoor temperature. Actually only for output at initial state.
-    DO  nb = 1, num_build
-       IF ( buildings(nb)%on_pe )  THEN
-          buildings(nb)%t_in(:) = initial_indoor_temperature
-
-!
-!--       (after first loop, use theta_m_t as theta_m_t_prev)
-          buildings(nb)%theta_m_t_prev_h(:) = initial_indoor_temperature
-          buildings(nb)%theta_m_t_prev_v(:) = initial_indoor_temperature
-
-       ENDIF
-    ENDDO
+    IF ( TRIM( initializing_actions ) /= 'read_restart_data' )  THEN
+       DO  nb = 1, num_build
+          IF ( buildings(nb)%on_pe )  THEN
+             buildings(nb)%t_in(:) = initial_indoor_temperature
+
+!
+!--          (after first loop, use theta_m_t as theta_m_t_prev)
+             buildings(nb)%theta_m_t_prev_h(:) = initial_indoor_temperature
+             buildings(nb)%theta_m_t_prev_v(:) = initial_indoor_temperature
+
+          ENDIF
+       ENDDO
+!
+!-- Initialize indoor temperature at previous timestep.
+    ELSE
+       DO  nb = 1, num_build
+          IF ( buildings(nb)%on_pe )  THEN
+!
+!--          Mean indoor temperature can be initialized with initial value. This is just
+!--          used for output.
+             buildings(nb)%t_in(:) = initial_indoor_temperature
+!
+!--          Initialize theta_m_t_prev arrays. The respective data during the restart mechanism
+!--          is stored on the surface-data array.
+             DO  fa = 1, buildings(nb)%num_facades_per_building_h_l
+!
+!--             Determine indices where corresponding surface-type information is stored.
+                l = buildings(nb)%l_h(fa)
+                m = buildings(nb)%m_h(fa)
+                buildings(nb)%theta_m_t_prev_h(fa) = surf_usm_h(l)%t_prev(m)
+             ENDDO
+             DO  fa = 1, buildings(nb)%num_facades_per_building_v_l
+!
+!--             Determine indices where corresponding surface-type information is stored.
+                l = buildings(nb)%l_v(fa)
+                m = buildings(nb)%m_v(fa)
+                buildings(nb)%theta_m_t_prev_v(fa) = surf_usm_v(l)%t_prev(m)
+             ENDDO
+          ENDIF
+       ENDDO
+    ENDIF
 
     CALL location_message( 'initializing indoor model', 'finished' )
@@ -1303 +1337 @@
           ENDIF
 !
-!--       Initialize/reset indoor temperature
+!--       Initialize/reset indoor temperature - note, this is only for output
           buildings(nb)%t_in_l = 0.0_wp
 !
@@ -1548 +1582 @@
                               buildings(nb)%params_waste_heat_h * heating_on +                     &
                               buildings(nb)%params_waste_heat_c * cooling_on )                     &
-                            ) / facade_element_area                                             !< [W/m2] , observe the directional
-                                                                                                !< convention in PALM!
+                            ) / facade_element_area  !< [W/m2] , observe the directional convention in PALM!
+!
+!--          Store waste heat and previous previous indoor temperature on surface-data type.
+!--          These will be used in the urban-surface model.
+             surf_usm_h(l)%t_prev(m) = buildings(nb)%theta_m_t_prev_h(fa)
              surf_usm_h(l)%waste_heat(m) = q_waste_heat
           ENDDO !< Horizontal surfaces loop
@@ -1559 +1596 @@
              l = buildings(nb)%l_v(fa)
              m = buildings(nb)%m_v(fa)
+
 !
 !--          During spinup set window fraction to zero and add these to wall fraction.
@@ -1805 +1843 @@
              q_waste_heat = ( phi_hc_nd * ( buildings(nb)%params_waste_heat_h * heating_on +       &
                                             buildings(nb)%params_waste_heat_c * cooling_on )       &
-                                                    ) / facade_element_area  !< [W/m2] , observe the directional convention in
-                                                                             !< PALM!
+                            ) / facade_element_area  !< [W/m2] , observe the directional convention in PALM!
+!
+!--          Store waste heat and previous previous indoor temperature on surface-data type.
+!--          These will be used in the urban-surface model.
+             surf_usm_v(l)%t_prev(m)     = buildings(nb)%theta_m_t_prev_v(fa)
              surf_usm_v(l)%waste_heat(m) = q_waste_heat
           ENDDO !< Vertical surfaces loop