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

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

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

revision 1.4 by jmc, Sat Feb 8 02:13:02 2003 UTC revision 1.6 by jmc, Tue Feb 18 15:36:45 2003 UTC
# Line 8  C !ROUTINE: MOM_CDSCHEME Line 8  C !ROUTINE: MOM_CDSCHEME
8    
9  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
10        SUBROUTINE MOM_CDSCHEME(        SUBROUTINE MOM_CDSCHEME(
11       I        bi,bj,k,phi_hyd,dPhiHydX,dPhiHydY,       I        bi,bj,k,dPhiHydX,dPhiHydY,
12       I        myThid)       I        myThid)
13    
14  C !DESCRIPTION:  C !DESCRIPTION:
# Line 27  C     == Global variables == Line 27  C     == Global variables ==
27  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
28  C  bi,bj                :: tile indices  C  bi,bj                :: tile indices
29  C  k                    :: vertical level  C  k                    :: vertical level
 C  phi_hyd              :: hydrostatic pressure  
30  C     dPhiHydX,Y        :: Gradient (X & Y dir.) of Hydrostatic Potential  C     dPhiHydX,Y        :: Gradient (X & Y dir.) of Hydrostatic Potential
31  C  myThid               :: thread number  C  myThid               :: thread number
32        INTEGER bi,bj,k        INTEGER bi,bj,k
       _RL phi_hyd(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)  
33        _RL dPhiHydX(1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL dPhiHydX(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
34        _RL dPhiHydY(1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL dPhiHydY(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
35        INTEGER myThid        INTEGER myThid
# Line 62  C     Adams-Bashforth weighting factors Line 60  C     Adams-Bashforth weighting factors
60        ab05   = -0.5 - abEps        ab05   = -0.5 - abEps
61    
62  C-- stagger time stepping: grad Phi_Hyp is not in gU,gV and needs to be added:  C-- stagger time stepping: grad Phi_Hyp is not in gU,gV and needs to be added:
63        IF (.FALSE.) THEN        IF (staggerTimeStep) THEN
 c     IF (staggerTimeStep) THEN  
64          phxFac = pfFacMom          phxFac = pfFacMom
65          phyFac = pfFacMom          phyFac = pfFacMom
66        ELSE        ELSE
# Line 79  C     Pressure extrapolated forward in t Line 76  C     Pressure extrapolated forward in t
76       &  +ab05*(etaNm1(i,j,bi,bj)*Bo_surf(i,j,bi,bj) )       &  +ab05*(etaNm1(i,j,bi,bj)*Bo_surf(i,j,bi,bj) )
77         ENDDO         ENDDO
78        ENDDO        ENDDO
       IF (staggerTimeStep) THEN  
        DO j=jMin,jMax  
         DO i=iMin,iMax  
          pf(i,j) = pf(i,j)+phi_hyd(i,j,k)  
         ENDDO  
        ENDDO  
       ENDIF  
79    
80  C--   Zonal velocity coriolis term  C--   Zonal velocity coriolis term
81  C     Note. As coded here, coriolis will not work with "thin walls"  C     Note. As coded here, coriolis will not work with "thin walls"
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22