Ignore:
Timestamp:
Mar 5, 2015 4:43:04 PM (7 years ago)
Author:
suehring
Message:

Enable monotone advection for scalars in combination with fifth-order scheme using monotonic limiter.

File:
1 edited

Legend:

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

    r1556 r1557  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Added checks for monotonic limiter
    2323!
    2424! Former revisions:
     
    605605    IF ( topography /= 'flat' )  THEN
    606606       action = ' '
    607        IF ( scalar_advec /= 'pw-scheme' .AND. scalar_advec /= 'ws-scheme')  THEN
     607       IF ( scalar_advec /= 'pw-scheme' .AND. scalar_advec /= 'ws-scheme'     &
     608      .AND. scalar_advec /= 'ws-scheme-mono' )  THEN
    608609          WRITE( action, '(A,A)' )  'scalar_advec = ', scalar_advec
    609610       ENDIF
     
    751752       CALL message( 'check_parameters', 'PA0022', 1, 2, 0, 6, 0 )
    752753    ENDIF
    753     IF ( ( momentum_advec == 'ws-scheme' .OR.  scalar_advec == 'ws-scheme' )   &
     754    IF ( ( momentum_advec == 'ws-scheme' .OR.  scalar_advec == 'ws-scheme'     &
     755           .OR. scalar_advec == 'ws-scheme-mono' )                             &
    754756           .AND. ( timestep_scheme == 'euler' .OR.                             &
    755757                   timestep_scheme == 'runge-kutta-2' ) )                      &
     
    761763    ENDIF
    762764    IF ( scalar_advec /= 'pw-scheme'  .AND.  scalar_advec /= 'ws-scheme' .AND. &
    763          scalar_advec /= 'bc-scheme' )                                         &
     765         scalar_advec /= 'ws-scheme-mono' .AND. scalar_advec /= 'bc-scheme' )  &
    764766    THEN
    765767       message_string = 'unknown advection scheme: scalar_advec = "' // &
     
    776778
    777779    IF ( use_sgs_for_particles  .AND.  .NOT. use_upstream_for_tke .AND.        &
    778          scalar_advec /= 'ws-scheme' )  THEN
     780         ( scalar_advec /= 'ws-scheme' .OR.                                    &
     781           scalar_advec /= 'ws-scheme-mono' )                                  &
     782       )  THEN
    779783       use_upstream_for_tke = .TRUE.
    780784       message_string = 'use_upstream_for_tke set .TRUE. because ' //          &
     
    792796!
    793797!-- Set LOGICAL switches to enhance performance
    794     IF ( momentum_advec == 'ws-scheme' )    ws_scheme_mom = .TRUE.
    795     IF ( scalar_advec   == 'ws-scheme'   )  ws_scheme_sca = .TRUE.
     798    IF ( momentum_advec == 'ws-scheme' )       ws_scheme_mom = .TRUE.
     799    IF ( scalar_advec   == 'ws-scheme' .OR.                                   &
     800         scalar_advec   == 'ws-scheme-mono' )  ws_scheme_sca = .TRUE.
     801    IF ( scalar_advec   == 'ws-scheme-mono' )  monotonic_adjustment = .TRUE.
     802
    796803
    797804!
     
    17071714          CALL message( 'check_parameters', 'PA0051', 1, 2, 0, 6, 0 )
    17081715       ENDIF
    1709        IF ( momentum_advec /= 'pw-scheme' .AND. &
    1710             momentum_advec /= 'ws-scheme')  THEN
     1716       IF ( scalar_advec /= 'pw-scheme' .AND.                                 &
     1717          ( scalar_advec /= 'ws-scheme' .OR.                                  &
     1718            scalar_advec /= 'ws-scheme-mono' )                                &
     1719          )  THEN
     1720
    17111721          message_string = 'non-cyclic lateral boundaries do not allow ' // &
    17121722                           'momentum_advec = "' // TRIM( momentum_advec ) // '"'
    17131723          CALL message( 'check_parameters', 'PA0052', 1, 2, 0, 6, 0 )
    17141724       ENDIF
    1715        IF ( scalar_advec /= 'pw-scheme' .AND. &
    1716             scalar_advec /= 'ws-scheme' )  THEN
     1725       IF ( scalar_advec /= 'pw-scheme' .AND.                                  &
     1726          ( scalar_advec /= 'ws-scheme' .OR.                                   &
     1727            scalar_advec /= 'ws-scheme-mono' )                                 &
     1728          )  THEN
    17171729          message_string = 'non-cyclic lateral boundaries do not allow ' // &
    17181730                           'scalar_advec = "' // TRIM( scalar_advec ) // '"'
Note: See TracChangeset for help on using the changeset viewer.