Ignore:
Timestamp:
Mar 18, 2007 11:46:30 PM (15 years ago)
Author:
raasch
Message:

bugs fixed for particle code and bc-scheme

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/prognostic_equations.f90

    r63 r70  
    316316!
    317317!-- pt-tendency terms with communication
     318    sat = tsc(1)
     319    sbt = tsc(2)
    318320    IF ( scalar_advec == 'bc-scheme' )  THEN
    319 !
    320 !--    Bott-Chlond scheme always uses Euler time step. Thus:
    321        sat = 1.0
    322        sbt = 1.0
     321
     322       IF ( timestep_scheme(1:5) /= 'runge' )  THEN
     323!
     324!--       Bott-Chlond scheme always uses Euler time step when leapfrog is
     325!--       switched on. Thus:
     326          sat = 1.0
     327          sbt = 1.0
     328       ENDIF
    323329       tend = 0.0
    324330       CALL advec_s_bc( pt, 'pt' )
    325331    ELSE
    326        sat = tsc(1)
    327        sbt = tsc(2)
    328332       IF ( tsc(2) /= 2.0  .AND.  scalar_advec == 'ups-scheme' )  THEN
    329333          tend = 0.0
     
    411415!
    412416!--    Scalar/q-tendency terms with communication
     417       sat = tsc(1)
     418       sbt = tsc(2)
    413419       IF ( scalar_advec == 'bc-scheme' )  THEN
    414 !
    415 !--       Bott-Chlond scheme always uses Euler time step. Thus:
    416           sat = 1.0
    417           sbt = 1.0
     420
     421          IF ( timestep_scheme(1:5) /= 'runge' )  THEN
     422!
     423!--          Bott-Chlond scheme always uses Euler time step when leapfrog is
     424!--          switched on. Thus:
     425             sat = 1.0
     426             sbt = 1.0
     427          ENDIF
    418428          tend = 0.0
    419429          CALL advec_s_bc( q, 'q' )
    420430       ELSE
    421           sat = tsc(1)
    422           sbt = tsc(2)
    423431          IF ( tsc(2) /= 2.0 )  THEN
    424432             IF ( scalar_advec == 'ups-scheme' )  THEN
     
    507515!--    TKE-tendency terms with communication
    508516       CALL production_e_init
     517
     518       sat = tsc(1)
     519       sbt = tsc(2)
    509520       IF ( .NOT. use_upstream_for_tke )  THEN
    510521          IF ( scalar_advec == 'bc-scheme' )  THEN
    511 !
    512 !--          Bott-Chlond scheme always uses Euler time step. Thus:
    513              sat = 1.0
    514              sbt = 1.0
     522
     523             IF ( timestep_scheme(1:5) /= 'runge' )  THEN
     524!
     525!--             Bott-Chlond scheme always uses Euler time step when leapfrog is
     526!--             switched on. Thus:
     527                sat = 1.0
     528                sbt = 1.0
     529             ENDIF
    515530             tend = 0.0
    516531             CALL advec_s_bc( e, 'e' )
    517532          ELSE
    518              sat = tsc(1)
    519              sbt = tsc(2)
    520533             IF ( tsc(2) /= 2.0 )  THEN
    521534                IF ( scalar_advec == 'ups-scheme' )  THEN
Note: See TracChangeset for help on using the changeset viewer.