164 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
165 |
strain(i,j) = 0. _d 0 |
strain(i,j) = 0. _d 0 |
166 |
tension(i,j) = 0. _d 0 |
tension(i,j) = 0. _d 0 |
167 |
|
hFacZ(i,j) = 0. _d 0 |
168 |
#endif |
#endif |
169 |
ENDDO |
ENDDO |
170 |
ENDDO |
ENDDO |
352 |
ENDDO |
ENDDO |
353 |
ENDDO |
ENDDO |
354 |
ENDIF |
ENDIF |
355 |
|
#ifdef ALLOW_SHELFICE |
356 |
|
IF (useShelfIce.AND.momViscosity.AND.bottomDragTerms) THEN |
357 |
|
CALL SHELFICE_U_DRAG(bi,bj,k,uFld,KE,KappaRU,vF,myThid) |
358 |
|
DO j=jMin,jMax |
359 |
|
DO i=iMin,iMax |
360 |
|
guDiss(i,j) = guDiss(i,j) + vF(i,j) |
361 |
|
ENDDO |
362 |
|
ENDDO |
363 |
|
ENDIF |
364 |
|
#endif /* ALLOW_SHELFICE */ |
365 |
|
|
366 |
|
|
367 |
C--- Other dissipation terms in Meridional momentum equation |
C--- Other dissipation terms in Meridional momentum equation |
368 |
|
|
417 |
ENDDO |
ENDDO |
418 |
ENDDO |
ENDDO |
419 |
ENDIF |
ENDIF |
420 |
|
#ifdef ALLOW_SHELFICE |
421 |
|
IF (useShelfIce.AND.momViscosity.AND.bottomDragTerms) THEN |
422 |
|
CALL SHELFICE_V_DRAG(bi,bj,k,vFld,KE,KappaRU,vF,myThid) |
423 |
|
DO j=jMin,jMax |
424 |
|
DO i=iMin,iMax |
425 |
|
gvDiss(i,j) = gvDiss(i,j) + vF(i,j) |
426 |
|
ENDDO |
427 |
|
ENDDO |
428 |
|
ENDIF |
429 |
|
#endif /* ALLOW_SHELFICE */ |
430 |
|
|
431 |
|
|
432 |
C- Vorticity diagnostics: |
C- Vorticity diagnostics: |
433 |
IF ( writeDiag ) THEN |
IF ( writeDiag ) THEN |