/[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.23 by jmc, Thu Mar 8 20:27:33 2001 UTC revision 1.25 by adcroft, Tue Jul 31 15:01:33 2001 UTC
# Line 37  C     == Local variables == Line 37  C     == Local variables ==
37        _RL phxFac,phyFac, psFac        _RL phxFac,phyFac, psFac
38    
39  C     Adams-Bashforth timestepping weights  C     Adams-Bashforth timestepping weights
40  Caja  IF (myIter .EQ. 0) THEN        IF (myIter .EQ. 0) THEN
41  Caja   ab15=1.0         ab15=1.0
42  Caja   ab05=0.0         ab05=0.0
43  Caja  ELSE        ELSE
44         ab15=1.5+abeps         ab15=1.5+abeps
45         ab05=-0.5-abeps         ab05=-0.5-abeps
46  Caja  ENDIF        ENDIF
47    
48  C     Step forward zonal velocity (store in Gu)  C     Step forward zonal velocity (store in Gu)
49        psFac = pfFacMom*(1. _d 0 - implicSurfPress)        psFac = pfFacMom*(1. _d 0 - implicSurfPress)
# Line 93  C--   -grad Phi_Hyd has not been incorpo Line 93  C--   -grad Phi_Hyd has not been incorpo
93            DO i=iMin,iMax            DO i=iMin,iMax
94              gVNm1(i,j,k,bi,bj)=gVNm1(i,j,k,bi,bj)              gVNm1(i,j,k,bi,bj)=gVNm1(i,j,k,bi,bj)
95       &       - _recip_dyC(i,j,bi,bj)       &       - _recip_dyC(i,j,bi,bj)
96       &         *(phiHyd(i,j,k)-phiHyd(i,j-1,k))*phyFac                             &         *(phiHyd(i,j,k)-phiHyd(i,j-1,k))*phyFac
97       &         *_maskS(i,j,k,bi,bj)       &         *_maskS(i,j,k,bi,bj)
98            ENDDO            ENDDO
99          ENDDO          ENDDO
Legend:
Removed from v.1.23  
changed lines
  Added in v.1.25
  ViewVC Help
Powered by ViewVC 1.1.22