Changeset 4573 for palm/trunk/SOURCE/time_integration.f90

Timestamp:

Jun 24, 2020 1:08:47 PM (4 years ago)

Author:

oliver.maas

Message:

pt_surface_heating_rate: calculate pt(0) incrementally by using dt_3d instead of calculating it absolutely by using time_since_reference_point, because time_since_reference_point is set to zero for initializing_actions = 'cyclic_fill', add statement for pt_surface_heating_rate in header.f90

File:

• palm/trunk/SOURCE/time_integration.f90 (modified) (2 diffs)

palm/trunk/SOURCE/time_integration.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! calculate pt(0) incrementally by using dt_3d instead of calculating it absolutely
+! by using time_since_reference_point, because time_since_reference_point is set
+! to zero for initializing_actions = 'cyclic_fill'
+! 
+! 4565 2020-06-15 08:30:38Z oliver.maas
 ! added new surface temperature forcing method for bc_pt_b = 'dirichlet':
 ! surface temperature pt(0) can be linearly increased by pt_surface_heating_rate (in K/h)
@@ -950 +955 @@
 !
 !--       Increase temperature pt(0) according to pt_surface_heating_rate (convert from K/h to K/s)
-          IF ( pt_surface_heating_rate /= 0.0_wp ) THEN
-             pt(0,:,:) = pt_surface + pt_surface_initial_change +                                  &
-                         time_since_reference_point * pt_surface_heating_rate / 3600.0_wp
+          IF ( pt_surface_heating_rate /= 0.0_wp .AND. intermediate_timestep_count == 1 ) THEN
+             pt(0,:,:) = pt(0,:,:) + dt_3d * pt_surface_heating_rate / 3600.0_wp
           ENDIF