Changeset 3608

Timestamp:

Dec 7, 2018 12:59:57 PM (5 years ago)

Author:

suehring

Message:

Bugfix in initialization of surfaces in cyclic-fill case + bugfix in radiation output

Location:

palm/trunk/SOURCE

Files:

• init_3d_model.f90 (modified) (3 diffs, 1 prop)
• radiation_model_mod.f90 (modified) (2 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix in initialization of surfaces in cyclic-fill case
+! 
+! 3589 2018-11-30 15:09:51Z suehring
 ! Move the control parameter "salsa" from salsa_mod to control_parameters 
 ! (M. Kurppa)
@@ -1604 +1607 @@
                               .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.     
-       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
@@ -1637 +1633 @@
 #endif
        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
 
 !

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Bugfix radiation output
+! 
+! 3607 2018-12-07 11:56:58Z suehring
 ! Output of radiation-related quantities migrated to radiation_model_mod.
 ! 
@@ -10372 +10375 @@
 
             IF ( surf(ix,isurfs) == is  .AND.  surf(iy,isurfs) == js  .AND. surf(iz,isurfs) == ks  .AND. &
-                 (surf(id,isurfs) == idsint_u .OR. surfl(id,isurf) == idsint_l ) ) THEN
+                 (surf(id,isurfs) == idsint_u .OR. surfl(id,isurfs) == idsint_l ) ) THEN
 !--            correct source surface
                local_pf(surfl(ix,isurflt),surfl(iy,isurflt),surfl(iz,isurflt)) = svf(k,isvf)