[v4.2.x][v4.0.x] Wrong calculation of `zcoef` in *sbcssm.f90*

I have noticed a very minor bug in sbcssm.F90 L.225 :

zcoef = REAL( nn_fsbc - 1, wp ) / zf_sbc

should be :

zcoef = REAL( nn_fsbc - 1, wp ) / ( zf_sbc - 1)

Changing the frequency of call to SI³ is very rare and the error only affects the first time step, however, for the sake of robustness, better to fix it at some point.

Best regards
Virginie

1 Like

I agree. I will commit it.

fixed in r4.0-HEAD, branch_4.2 and main
https://forge.ipsl.jussieu.fr/nemo/changeset/15769