/[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.11 by jmc, Tue Feb 18 15:36:45 2003 UTC revision 1.12 by jmc, Thu Apr 17 13:44:10 2003 UTC
# Line 409  C-    No-slip BCs impose a drag at botto Line 409  C-    No-slip BCs impose a drag at botto
409         ENDDO         ENDDO
410        ENDIF        ENDIF
411    
412  C--   Forcing term  C--   Forcing term (moved to timestep.F)
413        IF (momForcing)  c     IF (momForcing)
414       &  CALL EXTERNAL_FORCING_U(  c    &  CALL EXTERNAL_FORCING_U(
415       I     iMin,iMax,jMin,jMax,bi,bj,k,  c    I     iMin,iMax,jMin,jMax,bi,bj,k,
416       I     myTime,myThid)  c    I     myTime,myThid)
417    
418  C--   Metric terms for curvilinear grid systems  C--   Metric terms for curvilinear grid systems
419        IF (useNHMTerms) THEN        IF (useNHMTerms) THEN
# Line 561  C-    No-slip BCs impose a drag at botto Line 561  C-    No-slip BCs impose a drag at botto
561         ENDDO         ENDDO
562        ENDIF        ENDIF
563    
564  C--   Forcing term  C--   Forcing term (moved to timestep.F)
565        IF (momForcing)  c     IF (momForcing)
566       & CALL EXTERNAL_FORCING_V(  c    & CALL EXTERNAL_FORCING_V(
567       I     iMin,iMax,jMin,jMax,bi,bj,k,  c    I     iMin,iMax,jMin,jMax,bi,bj,k,
568       I     myTime,myThid)  c    I     myTime,myThid)
569    
570  C--   Metric terms for curvilinear grid systems  C--   Metric terms for curvilinear grid systems
571        IF (useNHMTerms) THEN        IF (useNHMTerms) THEN
# Line 595  C--   Set dv/dt on boundaries to zero Line 595  C--   Set dv/dt on boundaries to zero
595    
596  C--   Coriolis term  C--   Coriolis term
597  C     Note. As coded here, coriolis will not work with "thin walls"  C     Note. As coded here, coriolis will not work with "thin walls"
598  #ifdef INCLUDE_CD_CODE  c     IF (useCDscheme) THEN
599        CALL MOM_CDSCHEME(bi,bj,k,dPhiHydX,dPhiHydY,myThid)  c       CALL MOM_CDSCHEME(bi,bj,k,dPhiHydX,dPhiHydY,myThid)
600  #else  c     ELSE
601        CALL MOM_U_CORIOLIS(bi,bj,k,vFld,cf,myThid)        IF (.NOT.useCDscheme) THEN
602        DO j=jMin,jMax          CALL MOM_U_CORIOLIS(bi,bj,k,vFld,cf,myThid)
603         DO i=iMin,iMax          DO j=jMin,jMax
604          gU(i,j,k,bi,bj) = gU(i,j,k,bi,bj)+fuFac*cf(i,j)           DO i=iMin,iMax
605         ENDDO            gU(i,j,k,bi,bj) = gU(i,j,k,bi,bj)+fuFac*cf(i,j)
606        ENDDO           ENDDO
607        CALL MOM_V_CORIOLIS(bi,bj,k,uFld,cf,myThid)          ENDDO
608        DO j=jMin,jMax          CALL MOM_V_CORIOLIS(bi,bj,k,uFld,cf,myThid)
609         DO i=iMin,iMax          DO j=jMin,jMax
610          gV(i,j,k,bi,bj) = gV(i,j,k,bi,bj)+fvFac*cf(i,j)           DO i=iMin,iMax
611         ENDDO            gV(i,j,k,bi,bj) = gV(i,j,k,bi,bj)+fvFac*cf(i,j)
612        ENDDO           ENDDO
613  #endif /* INCLUDE_CD_CODE */          ENDDO
614          ENDIF
615    
616        IF (nonHydrostatic.OR.quasiHydrostatic) THEN        IF (nonHydrostatic.OR.quasiHydrostatic) THEN
617         CALL MOM_U_CORIOLIS_NH(bi,bj,k,wVel,cf,myThid)         CALL MOM_U_CORIOLIS_NH(bi,bj,k,wVel,cf,myThid)
618         DO j=jMin,jMax         DO j=jMin,jMax
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22