Changeset 2618 for palm/trunk/SOURCE/init_3d_model.f90

Timestamp:

Nov 16, 2017 3:37:30 PM (6 years ago)

Author:

suehring

Message:

Reorder calls for initializing surface quantities and give example of user-defined initialization of surface quantities

File:

• palm/trunk/SOURCE/init_3d_model.f90 (modified) (5 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Reorder calls of init_surfaces.
+! 
+! 2564 2017-10-19 15:56:56Z Giersch
 ! Variable wind_turbine was added to control_parameters.
 ! 
@@ -1135 +1138 @@
 
           ENDIF
+!
+!--       Initialize surface variables, e.g. friction velocity, momentum 
+!--       fluxes, etc. 
+          CALL init_surfaces
 
           CALL location_message( 'finished', .TRUE. )
@@ -1247 +1254 @@
 !--       of a sloping surface
           IF ( sloping_surface )  CALL init_slope
+!
+!--       Initialize surface variables, e.g. friction velocity, momentum 
+!--       fluxes, etc. 
+          CALL init_surfaces
 
           CALL location_message( 'finished', .TRUE. )
@@ -1254 +1265 @@
 
           CALL location_message( 'initializing by user', .FALSE. )
+!
+!--       Pre-initialize surface variables, i.e. setting start- and end-indices
+!--       at each (j,i)-location. Please note, this does not supersede 
+!--       user-defined initialization of surface quantities. 
+          CALL init_surfaces
 !
 !--       Initialization will completely be done by the user
@@ -1924 +1940 @@
     ENDIF
 !
-!-- Initialize surface elements and its attributes, e.g. heat- and 
-!-- momentumfluxes, roughness, scaling parameters. 
-!-- This is already done in case of restart data.  
+!-- Finally, if random_heatflux is set, disturb shf at horizontal
+!-- surfaces. Actually, this should be done in surface_mod, where all other 
+!-- initializations of surface quantities are done. However, this
+!-- would create a ring dependency, hence, it is done here. Maybe delete
+!-- disturb_heatflux and tranfer the respective code directly into the 
+!-- initialization in surface_mod.         
     IF ( TRIM( initializing_actions ) /= 'read_restart_data'  .AND.            &
          TRIM( initializing_actions ) /= 'cyclic_fill' )  THEN
-       CALL init_surfaces
-!
-!--    Finally, if random_heatflux is set, disturb shf at horizontal
-!--    surfaces. Actually, this should be done in surface_mod, where all other 
-!--    initializations of surface quantities are done. However, this
-!--    would create a ring dependency, hence, it is done here. Maybe delete
-!--    disturb_heatflux and tranfer the respective code directly into the 
-!--    initialization in surface_mod.         
+ 
        IF ( use_surface_fluxes  .AND.  constant_heatflux  .AND.                &
             random_heatflux )  THEN