[MITgcm-devel] kpp: ustarY

[Jean-Michel Campin]

Hi,

I found a strange line of code in kpp_calc.F, that is there for quiet
a long time (in version 1.2, from 2000/09/11, line 501 & 502):

499 heimbach 1.2                   ustarX = ( SurfaceTendencyU(i,  j,bi,bj) + 
500                   &                       SurfaceTendencyU(ip1,j,bi,bj) ) * p5
501                                ustarY = ( SurfaceTendencyV(i,j,  bi,bj) +
502                   &                       SurfaceTendencyU(i,jp1,bi,bj) ) * p5

looks like ustarY should be function of SurfaceTendencyV only, no ?

Any way, if I change line 502 to:
 &                       SurfaceTendencyV(i,jp1,bi,bj) ) * p5
it does not seem to have a big impact on the 2 exp. that are using KPP,
lab_sea & natl_box (get identical output.txt)

Jean-Michel