Bug with RK3, flux form momentum advection and linear free surface
I observed weird results with key_RK3
and key_linssh
defined with a flux form momentum advection. This comes from a wrong conditional statement in dynspg_ts.F90
here
Instead of
IF(.NOT.ln_dynadv_vec .AND. ln_bt_av ) THEN
!
! Blabla, transform time averaged fluxes into velocities
!
One should have:
IF( (.NOT.(ln_dynadv_vec .OR. ln_linssh)) .AND. ln_bt_av ) THEN
!
! Blabla, transform time averaged fluxes into velocities
!