Changeset 3609

Timestamp:

Dec 7, 2018 1:37:59 PM (6 years ago)

Author:

suehring

Message:

Further correction in initialization of surfaces in cyclic-fill case

File:

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

palm/trunk/SOURCE/init_3d_model.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Furhter correction in initialization of surfaces in cyclic-fill case
+! 
+! 3608 2018-12-07 12:59:57Z suehring
 ! Bugfix in initialization of surfaces in cyclic-fill case
 ! 
@@ -1607 +1610 @@
                               .FALSE. )
 !
+!--    Initialize surface elements and its attributes, e.g. heat- and 
+!--    momentumfluxes, roughness, scaling parameters. As number of surface 
+!--    elements might be different between runs, e.g. in case of cyclic fill, 
+!--    and not all surface elements are read, surface elements need to be 
+!--    initialized before. 
+!--    Please note, in case of cyclic fill, surfaces should be initialized 
+!--    after restart data is read, else, individual settings of surface 
+!--    parameters will be overwritten from data of precursor run, hence, 
+!--    init_surfaces is called a second time after reading the restart data. 
+       CALL init_surfaces                        
+!
 !--    When reading data for cyclic fill of 3D prerun data files, read
 !--    some of the global variables from the restart file which are required
@@ -1634 +1648 @@
        ENDDO
 !
-!--    Initialize surface elements and its attributes, e.g. heat- and 
-!--    momentumfluxes, roughness, scaling parameters. As number of surface 
-!--    elements might be different between runs, e.g. in case of cyclic fill, 
-!--    and not all surface elements are read, surface elements need to be 
-!--    initialized before. 
-!--    Please note, surfaces should be initialized after restart data is read,
-!--    else, individual settings of surface parameters will be overwritten
-!--    from data of precursor run. 
-       CALL init_surfaces
+!--    In case of cyclic fill, call init_surfaces a second time, so that 
+!--    surface properties such as heat fluxes are initialized as prescribed. 
+       IF ( TRIM( initializing_actions ) == 'cyclic_fill' )                    &
+          CALL init_surfaces
 
 !