Changeset 2296 for palm/trunk/SOURCE/urban_surface_mod.f90

Timestamp:

Jun 28, 2017 7:53:56 AM (6 years ago)

Author:

maronga

Message:

added new spinup mechanism for surface/radiation models

File:

• palm/trunk/SOURCE/urban_surface_mod.f90 (modified) (6 diffs)

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Bugfix: adjust output of shf. Added support for spinups
+! 
+! 2287 2017-06-15 16:46:30Z suehring
 ! Bugfix in determination topography-top index
 ! 
@@ -139 +142 @@
     
     USE control_parameters,                                                    &
-        ONLY:  dz, topography, dt_3d, intermediate_timestep_count,             &
-               initializing_actions, intermediate_timestep_count_max,          &
-               simulated_time, end_time, timestep_scheme, tsc,                 &
-               coupling_char, io_blocks, io_group, message_string,             &
-               time_since_reference_point, surface_pressure,                   &
-               g, pt_surface, large_scale_forcing, lsf_surf,                   &
-               time_do3d, dt_do3d, average_count_3d, varnamelength,            &
-               urban_surface
+        ONLY:  coupling_start_time, dz, topography, dt_3d,                     &
+               intermediate_timestep_count, initializing_actions,              &
+               intermediate_timestep_count_max, simulated_time, end_time,      &
+               timestep_scheme, tsc, coupling_char, io_blocks, io_group,       &
+               message_string, time_since_reference_point, surface_pressure,   &
+               g, pt_surface, large_scale_forcing, lsf_surf, spinup,           &
+               spinup_pt_mean, spinup_time, time_do3d, dt_do3d,                &
+               average_count_3d, varnamelength, urban_surface
 
     USE cpulog,                                                                &
@@ -3489 +3492 @@
 !--    Set flag that indicates that the land surface model is switched on
        urban_surface = .TRUE.
+
+!
+!--    Activate spinup
+       IF ( spinup_time > 0.0_wp )  THEN
+          coupling_start_time = spinup_time
+          end_time = end_time + spinup_time
+          IF ( spinup_pt_mean == 9999999.9_wp )  THEN
+             spinup_pt_mean = pt_surface
+          ENDIF
+          spinup = .TRUE.
+       ENDIF
 
  10    CONTINUE
@@ -5161 +5175 @@
 !--        store kinematic surface heat fluxes for utilization in other processes
 !--        diffusion_s, surface_layer_fluxes,...
-           surf_usm_h%shf(m) = surf_usm_h%wshf_eb(m) / rho_cp
+           surf_usm_h%shf(m) = surf_usm_h%wshf_eb(m) / cp
 
        ENDDO
@@ -5282 +5296 @@
 !--           store kinematic surface heat fluxes for utilization in other processes
 !--           diffusion_s, surface_layer_fluxes,...
-              surf_usm_v(l)%shf(m) = surf_usm_v(l)%wshf_eb(m) / rho_cp
+              surf_usm_v(l)%shf(m) = surf_usm_v(l)%wshf_eb(m) / cp
 
            ENDDO
@@ -5316 +5330 @@
 !--              Instead of this, we can adjust shf in case AH only at surface
                  surf_usm_h%shf(m) = surf_usm_h%shf(m) +                       &
-                                   aheat(j,i) * acoef * ddx * ddy / rho_cp
+                                   aheat(j,i) * acoef * ddx * ddy / cp
               ENDIF
            ENDDO