179 |
fMer(i,j) = 0. |
fMer(i,j) = 0. |
180 |
rTransU(i,j) = 0. |
rTransU(i,j) = 0. |
181 |
rTransV(i,j) = 0. |
rTransV(i,j) = 0. |
|
#ifdef ALLOW_AUTODIFF_TAMC |
|
|
C- jmc: this is wrong, but at least with #ifdef/endif TAMC, it does not break |
|
|
C the forward code ; (same thing in mom_vectinv) |
|
|
fVerU(i,j,1) = 0. _d 0 |
|
|
fVerU(i,j,2) = 0. _d 0 |
|
|
fVerV(i,j,1) = 0. _d 0 |
|
|
fVerV(i,j,2) = 0. _d 0 |
|
|
#endif |
|
182 |
ENDDO |
ENDDO |
183 |
ENDDO |
ENDDO |
184 |
|
|
253 |
ENDDO |
ENDDO |
254 |
ENDDO |
ENDDO |
255 |
|
|
256 |
CALL MOM_CALC_KE(bi,bj,k,uFld,vFld,KE,myThid) |
CALL MOM_CALC_KE(bi,bj,k,3,uFld,vFld,KE,myThid) |
257 |
|
|
258 |
C--- First call (k=1): compute vertical adv. flux fVerU(kUp) & fVerV(kUp) |
C--- First call (k=1): compute vertical adv. flux fVerU(kUp) & fVerV(kUp) |
259 |
IF (momAdvection.AND.k.EQ.1) THEN |
IF (momAdvection.AND.k.EQ.1) THEN |