[MITgcm-devel] bloc and bloc in shelfice_thermodynamics

[Menemenlis, Dimitris (329D)]

Hi Martin, Yoshi and I wonder why you don’t multiply uLoc and vLoc by 0.5 in lines below:

          uLoc(I,J) = recip_hFacC(I,J,K,bi,bj) *
     &         ( uVel(I,  J,K,bi,bj) * _hFacW(I,  J,K,bi,bj)
     &         + uVel(I+1,J,K,bi,bj) * _hFacW(I+1,J,K,bi,bj) )
          vLoc(I,J) = recip_hFacC(I,J,K,bi,bj) *
     &         ( vVel(I,J,  K,bi,bj) * _hFacS(I,J,  K,bi,bj)
     &         + vVel(I,J+1,K,bi,bj) * _hFacS(I,J+1,K,bi,bj) )

            uLoc(I,J) = ( uLoc(I,J) * drF(K)*_hFacC(I,J,K,bi,bj)
     &           + drKp1 * recip_hFacC(I,J,Kp1,bi,bj) *
     &           ( uVel(I,  J,Kp1,bi,bj) * _hFacW(I,  J,Kp1,bi,bj)
     &           + uVel(I+1,J,Kp1,bi,bj) * _hFacW(I+1,J,Kp1,bi,bj) )
     &           ) * recip_drLoc
            vLoc(I,J) = ( vLoc(I,J) * drF(K)*_hFacC(I,J,K,bi,bj)
     &           + drKp1 * recip_hFacC(I,J,Kp1,bi,bj) *
     &           ( vVel(I,J,  Kp1,bi,bj) * _hFacS(I,J,  Kp1,bi,bj)
     &           + vVel(I,J+1,Kp1,bi,bj) * _hFacS(I,J+1,Kp1,bi,bj) )
     &           ) * recip_drLoc

You do it later on in line:

            ustarSq = shiCdrag * MAX( 1.D-6,
     &           0.25 _d 0 *(uLoc(I,J)*uLoc(I,J)+vLoc(I,J)*vLoc(I,J)) )

The problem with this formulation is that if one is trying to modify
drag formulation, it’s very easy to forget that uLoc = 2 * u
(Yoshi spent days tracking this missing factor :-(

Shall we change for clarity?

Cheers, Dimitris