455 |
gU(i,j,k,bi,bj) = gU(i,j,k,bi,bj)+mTFacU*mT(i,j) |
gU(i,j,k,bi,bj) = gU(i,j,k,bi,bj)+mTFacU*mT(i,j) |
456 |
ENDDO |
ENDDO |
457 |
ENDDO |
ENDDO |
458 |
|
|
459 |
|
IF (usingCylindricalGrid) THEN |
460 |
|
CALL MOM_U_METRIC_CYLINDER(bi,bj,k,uFld,vFld,mT,myThid) |
461 |
|
DO j=jMin,jMax |
462 |
|
DO i=iMin,iMax |
463 |
|
gU(i,j,k,bi,bj) = gU(i,j,k,bi,bj)+mTFacU*mT(i,j) |
464 |
|
ENDDO |
465 |
|
ENDDO |
466 |
|
ENDIF |
467 |
|
|
468 |
ENDIF |
ENDIF |
|
|
|
469 |
C-- Set du/dt on boundaries to zero |
C-- Set du/dt on boundaries to zero |
470 |
DO j=jMin,jMax |
DO j=jMin,jMax |
471 |
DO i=iMin,iMax |
DO i=iMin,iMax |
617 |
ENDDO |
ENDDO |
618 |
ENDDO |
ENDDO |
619 |
ENDIF |
ENDIF |
620 |
|
IF (usingCylindricalGrid) THEN |
621 |
|
CALL MOM_V_METRIC_CYLINDER(bi,bj,k,uFld,vFld,mT,myThid) |
622 |
|
DO j=jMin,jMax |
623 |
|
DO i=iMin,iMax |
624 |
|
gV(i,j,k,bi,bj) = gV(i,j,k,bi,bj)+mTFacV*mT(i,j) |
625 |
|
ENDDO |
626 |
|
ENDDO |
627 |
|
ENDIF |
628 |
|
|
629 |
C-- Set dv/dt on boundaries to zero |
C-- Set dv/dt on boundaries to zero |
630 |
DO j=jMin,jMax |
DO j=jMin,jMax |