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

Timestamp:

Jul 19, 2018 1:26:52 PM (6 years ago)

Author:

suehring

Message:

Further adjustments for surface structure

File:

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

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Initialize q_surface, which might be used in surface_layer_fluxes
+! 
+! 3151 2018-07-19 08:45:38Z raasch
 ! remaining preprocessor define strings __check removed
 ! 
@@ -313 +316 @@
     
     USE control_parameters,                                                    &
-        ONLY:  coupling_start_time, topography, dt_3d,                         &
+        ONLY:  coupling_start_time, topography, dt_3d, humidity,               &
                intermediate_timestep_count, initializing_actions,              &
                intermediate_timestep_count_max, simulated_time, end_time,      &
@@ -4706 +4709 @@
             ENDDO
 #endif
+!
+!--         For the sake of correct initialization, set also q_surface. 
+!--         Note, at urban surfaces q_surface is initialized with 0.
+            IF ( humidity )  THEN
+               DO  m = 1, surf_usm_h%ns
+                  surf_usm_h%q_surface(m) = 0.0_wp
+               ENDDO
+               DO  l = 0, 3
+                  DO  m = 1, surf_usm_v(l)%ns
+                     surf_usm_v(l)%q_surface(m) = 0.0_wp
+                  ENDDO
+               ENDDO
+            ENDIF
       
 !--         initial values for t_wall