Changeset 4335 for palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90

Timestamp:

Dec 12, 2019 4:39:05 PM (4 years ago)

Author:

suehring

Message:

Plant canopy: implement fix for LAD on building edges also for the vectorl-optimized branch

File:

• palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90 (modified) (2 diffs)

palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Commentation of last commit
+! 
+! 4332 2019-12-10 19:44:12Z suehring
 ! Limit initial velocity seeds in restart runs, if not the seed calculation 
 ! may become unstable. Further, minor bugfix in initial velocity seed 
@@ -939 +942 @@
           DO  i = nxlg, nxrg
              DO  k = nzb, nzt+1
-
+!
+!--             In case the correlation coefficients are very small, the 
+!--             velocity seeds may become very large finally creating 
+!--             numerical instabilities in the synthetic turbulence generator.
+!--             Empirically, a value of 10E-8 seems to be sufficient.
                 IF  ( a11(k) > 10E-8_wp )  THEN
                    fu_xz(k,i) = ( u(k,j,i) - u_init(k) ) / a11(k)