| 374 |
I myCurrentTime,myThid) |
I myCurrentTime,myThid) |
| 375 |
|
|
| 376 |
C-- Metric terms for curvilinear grid systems |
C-- Metric terms for curvilinear grid systems |
| 377 |
IF (usingSphericalPolarMTerms) THEN |
IF (useNHMTerms) THEN |
| 378 |
C o Spherical polar grid metric terms |
C o Non-hydrosatic metric terms |
| 379 |
CALL MOM_U_METRIC_NH(bi,bj,k,uFld,wVel,mT,myThid) |
CALL MOM_U_METRIC_NH(bi,bj,k,uFld,wVel,mT,myThid) |
| 380 |
DO j=jMin,jMax |
DO j=jMin,jMax |
| 381 |
DO i=iMin,iMax |
DO i=iMin,iMax |
| 382 |
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) |
| 383 |
ENDDO |
ENDDO |
| 384 |
ENDDO |
ENDDO |
| 385 |
|
ENDIF |
| 386 |
|
IF (usingSphericalPolarMTerms) THEN |
| 387 |
CALL MOM_U_METRIC_SPHERE(bi,bj,k,uFld,vFld,mT,myThid) |
CALL MOM_U_METRIC_SPHERE(bi,bj,k,uFld,vFld,mT,myThid) |
| 388 |
DO j=jMin,jMax |
DO j=jMin,jMax |
| 389 |
DO i=iMin,iMax |
DO i=iMin,iMax |
| 524 |
I myCurrentTime,myThid) |
I myCurrentTime,myThid) |
| 525 |
|
|
| 526 |
C-- Metric terms for curvilinear grid systems |
C-- Metric terms for curvilinear grid systems |
| 527 |
IF (usingSphericalPolarMTerms) THEN |
IF (useNHMTerms) THEN |
| 528 |
C o Spherical polar grid metric terms |
C o Spherical polar grid metric terms |
| 529 |
CALL MOM_V_METRIC_NH(bi,bj,k,vFld,wVel,mT,myThid) |
CALL MOM_V_METRIC_NH(bi,bj,k,vFld,wVel,mT,myThid) |
| 530 |
DO j=jMin,jMax |
DO j=jMin,jMax |
| 532 |
gV(i,j,k,bi,bj) = gV(i,j,k,bi,bj)+mTFacV*mT(i,j) |
gV(i,j,k,bi,bj) = gV(i,j,k,bi,bj)+mTFacV*mT(i,j) |
| 533 |
ENDDO |
ENDDO |
| 534 |
ENDDO |
ENDDO |
| 535 |
|
ENDIF |
| 536 |
|
IF (usingSphericalPolarMTerms) THEN |
| 537 |
CALL MOM_V_METRIC_SPHERE(bi,bj,k,uFld,mT,myThid) |
CALL MOM_V_METRIC_SPHERE(bi,bj,k,uFld,mT,myThid) |
| 538 |
DO j=jMin,jMax |
DO j=jMin,jMax |
| 539 |
DO i=iMin,iMax |
DO i=iMin,iMax |