| 72 | | Bugfix: 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). Hence, logical operator use_cmax is now used instead of bc_lr_dirneu/_neudir (boundary_conds.f90, check_parameters.f90, header.f90, init_3d_model.f90, init_pegrid.f90, modules.f90, parin.f90, poismg.f90). |
| 73 | | |
| 74 | | Bugfix: In case of use_cmax at the outflow, u, v, w are replaced by u_p, v_p, w_p (boundary_conds.f90). |
| | 72 | Bugfix: If radiation boundary conditions are used, the calculation method of the phase velocity is now set by the logical operator use_cmax. The case use_cmax = ''.T.'' was wrongly defined as Neumann boundary condition before. (boundary_conds, check_parameters, header, init_3d_model, init_pegrid, modules, parin, poismg). |
| | 73 | |
