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

Timestamp:

Jun 15, 2020 8:30:38 AM (4 years ago)

Author:

oliver.maas

Message:

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)

File:

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

palm/trunk/SOURCE/time_integration.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! 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)
+! 
+! 4564 2020-06-12 14:03:36Z raasch
 ! Vertical nesting method of Huq et al. (2019) removed
 ! 
@@ -269 +273 @@
                multi_agent_system_end, multi_agent_system_start, nesting_offline, neutral,         &
                nr_timesteps_this_run, nudging, ocean_mode, passive_scalar, pt_reference,           &
-               pt_slope_offset, random_heatflux, rans_tke_e, run_coupled, salsa,                   &
+               pt_slope_offset, pt_surface, pt_surface_heating_rate, pt_surface_initial_change,    &
+               random_heatflux, rans_tke_e, run_coupled, salsa,                                    &
                simulated_time, simulated_time_chr, skip_time_do2d_xy, skip_time_do2d_xz,           &
                skip_time_do2d_yz, skip_time_do3d, skip_time_domask, skip_time_dopr,                &
@@ -944 +949 @@
 
 !
+!--       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
+          ENDIF
+
+!
 !--       Impose a turbulent inflow using the recycling method
           IF ( turbulent_inflow )  CALL inflow_turbulence