/[MITgcm]/MITgcm/model/src/timestep.F
ViewVC logotype

Diff of /MITgcm/model/src/timestep.F

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

revision 1.29 by jmc, Sat Feb 8 02:09:20 2003 UTC revision 1.30 by jmc, Sun Feb 9 02:00:50 2003 UTC
# Line 61  C     Adams-Bashforth timestepping weigh Line 61  C     Adams-Bashforth timestepping weigh
61        ENDIF        ENDIF
62    
63  C-- stagger time step: grad Phi_Hyp is not in gU,gV => add it in this S/R  C-- stagger time step: grad Phi_Hyp is not in gU,gV => add it in this S/R
64  c     IF (staggerTimeStep) THEN        IF (staggerTimeStep) THEN
65        IF (.FALSE.) THEN  c     IF (.FALSE.) THEN
66          phxFac = pfFacMom          phxFac = pfFacMom
67          phyFac = pfFacMom          phyFac = pfFacMom
68        ELSE        ELSE
# Line 119  C     Step forward zonal velocity (store Line 119  C     Step forward zonal velocity (store
119          ENDDO          ENDDO
120        ENDDO        ENDDO
121    
122  c     IF (.FALSE.) THEN        IF (.FALSE.) THEN
123        IF (staggerTimeStep) THEN  c     IF (staggerTimeStep) THEN
124  C--   -grad Phi_Hyd has not been incorporated to gU and is added here:  C--   -grad Phi_Hyd has not been incorporated to gU and is added here:
125          phxFac = pfFacMom*deltaTmom          phxFac = pfFacMom*deltaTmom
126          DO j=jMin,jMax          DO j=jMin,jMax
127            DO i=iMin,iMax            DO i=iMin,iMax
128              gUNm1(i,j,k,bi,bj)=gUNm1(i,j,k,bi,bj)              gUNm1(i,j,k,bi,bj)=gUNm1(i,j,k,bi,bj)
 c    &       - phxFac*dPhiHydX(i,j)*_maskW(i,j,k,bi,bj)  
129       &       - _recip_dxC(i,j,bi,bj)       &       - _recip_dxC(i,j,bi,bj)
130       &         *(phiHyd(i,j,k)-phiHyd(i-1,j,k))*phxFac       &         *(phiHyd(i,j,k)-phiHyd(i-1,j,k))*phxFac
131       &         *_maskW(i,j,k,bi,bj)       &         *_maskW(i,j,k,bi,bj)
# Line 180  C     Step forward meridional velocity ( Line 179  C     Step forward meridional velocity (
179          ENDDO          ENDDO
180        ENDDO        ENDDO
181    
182  c     IF (.FALSE.) THEN        IF (.FALSE.) THEN
183        IF (staggerTimeStep) THEN  c     IF (staggerTimeStep) THEN
184  C--   -grad Phi_Hyd has not been incorporated to gV and is added here:  C--   -grad Phi_Hyd has not been incorporated to gV and is added here:
185          phyFac = pfFacMom*deltaTmom          phyFac = pfFacMom*deltaTmom
186          DO j=jMin,jMax          DO j=jMin,jMax
187            DO i=iMin,iMax            DO i=iMin,iMax
188              gVNm1(i,j,k,bi,bj)=gVNm1(i,j,k,bi,bj)              gVNm1(i,j,k,bi,bj)=gVNm1(i,j,k,bi,bj)
 c    &       - phyFac*dPhiHydY(i,j)*_maskS(i,j,k,bi,bj)  
189       &       - _recip_dyC(i,j,bi,bj)       &       - _recip_dyC(i,j,bi,bj)
190       &         *(phiHyd(i,j,k)-phiHyd(i,j-1,k))*phyFac       &         *(phiHyd(i,j,k)-phiHyd(i,j-1,k))*phyFac
191       &         *_maskS(i,j,k,bi,bj)       &         *_maskS(i,j,k,bi,bj)
Legend:
Removed from v.1.29  
changed lines
  Added in v.1.30
  ViewVC Help
Powered by ViewVC 1.1.22