Changeset 3055 for palm/trunk/SOURCE/surface_mod.f90

Timestamp:

Jun 1, 2018 5:15:08 PM (6 years ago)

Author:

suehring

Message:

Bugfix, initialization of surface elements also in case of restart runs and cyclic fill

File:

• palm/trunk/SOURCE/surface_mod.f90 (modified) (8 diffs)

palm/trunk/SOURCE/surface_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Initialize Obukhov length, friction velocity and momentum fluxes also 
+! in case of restarts and cyclic fill
+! 
+! 3026 2018-05-22 10:30:53Z schwenkel
 ! Changed the name specific humidity to mixing ratio, since we are computing
 ! mixing ratios.
@@ -2053 +2057 @@
                 ENDIF
 !
-!--          Initialization in case of constant profiles
-             ELSEIF ( INDEX(initializing_actions, 'set_constant_profiles') /= 0&
-                .OR.  INDEX(initializing_actions, 'inifor' ) /= 0 )  THEN
+!--          Initialization in all other cases
+             ELSE
 
                 surf%ol(num_h)   = surf%z_mo(num_h) / zeta_min
@@ -2067 +2070 @@
              ENDIF
 
-             surf%rib(num_h)   = 0.0_wp 
+             surf%rib(num_h)     = 0.0_wp 
              surf%uvw_abs(num_h) = 0.0_wp
 
@@ -3321 +3324 @@
 
        found              = .TRUE.
-
 
        SELECT CASE ( restart_string(1:length) )
@@ -3946 +3948 @@
 
        END SELECT
-
 !
 !--    Redistribute surface elements on its respective type. 
@@ -3952 +3953 @@
             .NOT. INDEX( restart_string(1:length), '%start_index' ) /= 0 )     &
        THEN
+       
           ic = nxlc
           DO  i = nxlf, nxrf
@@ -3965 +3967 @@
                                     surf_usm_h%start_index(jc,ic) )            &
                             .AND.  l == 0 
-
+                           
                 IF ( surf_match_def )  THEN
                    mm = surf_def_h(l)%start_index(jc,ic)
@@ -4017 +4019 @@
                 surf_match_usm  = surf_usm_v(l)%end_index(jc,ic) >=            &
                                   surf_usm_v(l)%start_index(jc,ic)
-
+                                  
                 IF ( surf_match_def )  THEN
                    mm = surf_def_v(l)%start_index(jc,ic)