4 |
CBOI |
CBOI |
5 |
C !TITLE: pkg/mom\_advdiff |
C !TITLE: pkg/mom\_advdiff |
6 |
C !AUTHORS: adcroft@mit.edu |
C !AUTHORS: adcroft@mit.edu |
7 |
C !INTRODUCTION: |
C !INTRODUCTION: Flux-form Momentum Equations Package |
|
C \section{Flux-form Momentum Equations Package} |
|
8 |
C |
C |
9 |
C Package "mom\_fluxform" provides methods for calculating explicit terms |
C Package "mom\_fluxform" provides methods for calculating explicit terms |
10 |
C in the momentum equation cast in flux-form: |
C in the momentum equation cast in flux-form: |
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 |