/[MITgcm]/MITgcm/pkg/mom_fluxform/mom_fluxform.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_fluxform/mom_fluxform.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.3 by adcroft, Wed Sep 26 19:05:21 2001 UTC revision 1.6 by adcroft, Tue Nov 5 19:58:21 2002 UTC
# Line 4  C $Name$ Line 4  C $Name$
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:
# Line 375  C--   Forcing term Line 374  C--   Forcing term
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
# Line 523  C--   Forcing term Line 524  C--   Forcing term
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
# Line 531  C      o Spherical polar grid metric ter Line 532  C      o Spherical polar grid metric ter
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
# Line 564  C     Note. As coded here, coriolis will Line 567  C     Note. As coded here, coriolis will
567         ENDDO         ENDDO
568        ENDDO        ENDDO
569  #endif /* INCLUDE_CD_CODE */  #endif /* INCLUDE_CD_CODE */
570          IF (nonHydrostatic) THEN
571           CALL MOM_U_CORIOLIS_NH(bi,bj,k,wVel,cf,myThid)
572           DO j=jMin,jMax
573            DO i=iMin,iMax
574             gU(i,j,k,bi,bj) = gU(i,j,k,bi,bj)+fuFac*cf(i,j)
575            ENDDO
576           ENDDO
577          ENDIF
578    
579        RETURN        RETURN
580        END        END
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22