Changeset 4573 for palm

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

Location:

palm/trunk/SOURCE

Files:

• header.f90 (modified) (3 diffs)
• time_integration.f90 (modified) (2 diffs)

palm/trunk/SOURCE/header.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! added statement for pt_surface_heating_rate
+! 
+! 4536 2020-05-17 17:24:13Z raasch
 ! output of restart data format added
 ! 
@@ -1579 +1582 @@
        WRITE ( io, 475 )  pt_surface_initial_change
     ENDIF
+    IF ( pt_surface_heating_rate /= 0.0_wp )  THEN
+       WRITE ( io, 476 )  pt_surface_heating_rate
+    ENDIF
     IF ( humidity  .AND.  q_surface_initial_change /= 0.0_wp )  THEN
-       WRITE ( io, 476 )  q_surface_initial_change       
+       WRITE ( io, 477 )  q_surface_initial_change       
     ENDIF
     IF ( passive_scalar  .AND.  q_surface_initial_change /= 0.0_wp )  THEN
-       WRITE ( io, 477 )  q_surface_initial_change       
+       WRITE ( io, 478 )  q_surface_initial_change       
     ENDIF
 
@@ -1931 +1937 @@
             '    the value is negative) by ',F5.2,' K at the beginning of the',&
                  ' 3D-simulation'/)
-476 FORMAT ('    The surface humidity is increased (or decreased, ',&
+476 FORMAT ('    The surface temperature increases (or decreases, ', &
+                 'respectively, if'/ &
+            '    the value is negative) by ',F8.4,' K/h during the', &
+                 ' 3D-simulation'/)
+477 FORMAT ('    The surface humidity is increased (or decreased, ',&
                  'respectively, if the'/ &
             '    value is negative) by ',E8.1,' kg/kg at the beginning of', &
                  ' the 3D-simulation'/)
-477 FORMAT ('    The scalar value is increased at the surface (or decreased, ',&
+478 FORMAT ('    The scalar value is increased at the surface (or decreased, ',&
                  'respectively, if the'/ &
             '    value is negative) by ',E8.1,' kg/m**3 at the beginning of', &

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