Ignore:
Timestamp:
May 21, 2013 11:58:22 AM (8 years ago)
Author:
fricke
Message:

Bugfix: In case of non-cyclic lateral boundary conditions, Neumann boundary conditions for the velocity components at the outflow are in fact radiation boundary conditions using the maximum phase velocity that ensures numerical stability (CFL-condition).
Logical operator use_cmax is now used instead of bc_lr_dirneu/_neudir.

File:
1 edited

Legend:

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

    r1116 r1159  
    2020! Current revisions:
    2121! -----------------
     22! +use_cmax
    2223!
    2324!
     
    742743    WRITE ( io, 317 )  bc_lr, bc_ns
    743744    IF ( .NOT. bc_lr_cyc  .OR.  .NOT. bc_ns_cyc )  THEN
    744        WRITE ( io, 318 )  pt_damping_width, pt_damping_factor       
     745       WRITE ( io, 318 )  use_cmax, pt_damping_width, pt_damping_factor       
    745746       IF ( turbulent_inflow )  THEN
    746747          WRITE ( io, 319 )  recycling_width, recycling_plane, &
     
    17811782            '       left/right:  ',A/    &
    17821783            '       north/south: ',A)
    1783 318 FORMAT (/'       pt damping layer width = ',F7.2,' m, pt ', &
     1784318 FORMAT (/'       use_cmax: ',L1 / &
     1785            '       pt damping layer width = ',F8.2,' m, pt ', &
    17841786                    'damping factor = ',F6.4)
    17851787319 FORMAT ('       turbulence recycling at inflow switched on'/ &
Note: See TracChangeset for help on using the changeset viewer.