Changeset 1729 for palm

Timestamp:

Nov 20, 2015 11:01:48 AM (8 years ago)

Author:

gronemeier

Message:

bugfix for subsidence in combination with Rayleigh damping

File:

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

palm/trunk/SOURCE/subsidence.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: shifting of initial profiles only once each time step
 ! 
 ! Former revisions:
@@ -235 +235 @@
 !--    Shifting of the initial profile is especially necessary with Rayleigh 
 !--    damping switched on
-       IF ( scalar_rayleigh_damping ) THEN
+       IF ( scalar_rayleigh_damping .AND.                                      &
+            intermediate_timestep_count == 1 )  THEN
           DO  k = nzb, nzt
              IF ( w_subs(k) < 0.0_wp )  THEN      ! large-scale subsidence
@@ -334 +335 @@
 !--    Shifting of the initial profile is especially necessary with Rayleigh 
 !--    damping switched on
-       IF ( scalar_rayleigh_damping ) THEN
+       IF ( scalar_rayleigh_damping .AND.                                      &
+            intermediate_timestep_count == 1 )  THEN
           IF ( i == nxl .AND. j == nys )  THEN ! shifting only once per PE