Changeset 4573 for palm/trunk/SOURCE/header.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/header.f90 (modified) (3 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', &