Changeset 70 for palm/trunk/SOURCE/prognostic_equations.f90

Timestamp:

Mar 18, 2007 11:46:30 PM (17 years ago)

Author:

raasch

Message:

bugs fixed for particle code and bc-scheme

File:

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

palm/trunk/SOURCE/prognostic_equations.f90

@@ -316 +316 @@
 !
 !-- pt-tendency terms with communication
+    sat = tsc(1)
+    sbt = tsc(2)
     IF ( scalar_advec == 'bc-scheme' )  THEN
-!
-!--    Bott-Chlond scheme always uses Euler time step. Thus:
-       sat = 1.0
-       sbt = 1.0
+
+       IF ( timestep_scheme(1:5) /= 'runge' )  THEN
+!
+!--       Bott-Chlond scheme always uses Euler time step when leapfrog is
+!--       switched on. Thus:
+          sat = 1.0
+          sbt = 1.0
+       ENDIF
        tend = 0.0
        CALL advec_s_bc( pt, 'pt' )
     ELSE
-       sat = tsc(1)
-       sbt = tsc(2)
        IF ( tsc(2) /= 2.0  .AND.  scalar_advec == 'ups-scheme' )  THEN
           tend = 0.0
@@ -411 +415 @@
 !
 !--    Scalar/q-tendency terms with communication
+       sat = tsc(1)
+       sbt = tsc(2)
        IF ( scalar_advec == 'bc-scheme' )  THEN
-!
-!--       Bott-Chlond scheme always uses Euler time step. Thus:
-          sat = 1.0
-          sbt = 1.0
+
+          IF ( timestep_scheme(1:5) /= 'runge' )  THEN
+!
+!--          Bott-Chlond scheme always uses Euler time step when leapfrog is
+!--          switched on. Thus:
+             sat = 1.0
+             sbt = 1.0
+          ENDIF
           tend = 0.0
           CALL advec_s_bc( q, 'q' )
        ELSE
-          sat = tsc(1)
-          sbt = tsc(2)
           IF ( tsc(2) /= 2.0 )  THEN
              IF ( scalar_advec == 'ups-scheme' )  THEN
@@ -507 +515 @@
 !--    TKE-tendency terms with communication
        CALL production_e_init
+
+       sat = tsc(1)
+       sbt = tsc(2)
        IF ( .NOT. use_upstream_for_tke )  THEN
           IF ( scalar_advec == 'bc-scheme' )  THEN
-!
-!--          Bott-Chlond scheme always uses Euler time step. Thus:
-             sat = 1.0
-             sbt = 1.0
+
+             IF ( timestep_scheme(1:5) /= 'runge' )  THEN
+!
+!--             Bott-Chlond scheme always uses Euler time step when leapfrog is
+!--             switched on. Thus:
+                sat = 1.0
+                sbt = 1.0
+             ENDIF
              tend = 0.0
              CALL advec_s_bc( e, 'e' )
           ELSE
-             sat = tsc(1)
-             sbt = tsc(2)
              IF ( tsc(2) /= 2.0 )  THEN
                 IF ( scalar_advec == 'ups-scheme' )  THEN