/[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.22 by jmc, Tue Mar 6 17:10:29 2001 UTC revision 1.23 by jmc, Thu Mar 8 20:27:33 2001 UTC
# Line 20  C     == Global variables == Line 20  C     == Global variables ==
20  #include "GRID.h"  #include "GRID.h"
21    
22  C     == Routine Arguments ==  C     == Routine Arguments ==
23  C     phiHyd    - Hydrostatic pressure (ocean) or geopotentiel (atmos)  C     phiHyd    - Hydrostatic Potential (ocean: pressure/rho)
24    C                                       (atmos: geopotentiel)
25  C     phiSurfX, - gradient of Surface potentiel (Pressure/rho, ocean)  C     phiSurfX, - gradient of Surface potentiel (Pressure/rho, ocean)
26  C     phiSurfY             or geopotentiel (atmos) in X and Y direction  C     phiSurfY             or geopotentiel (atmos) in X and Y direction
27        INTEGER bi,bj,iMin,iMax,jMin,jMax        INTEGER bi,bj,iMin,iMax,jMin,jMax
# Line 60  C     Step forward zonal velocity (store Line 61  C     Step forward zonal velocity (store
61    
62        IF (staggerTimeStep) THEN        IF (staggerTimeStep) THEN
63  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:
64          phxFac = pfFacMom*deltaTmom*recip_rhoConst          phxFac = pfFacMom*deltaTmom
65          DO j=jMin,jMax          DO j=jMin,jMax
66            DO i=iMin,iMax            DO i=iMin,iMax
67              gUNm1(i,j,k,bi,bj)=gUNm1(i,j,k,bi,bj)              gUNm1(i,j,k,bi,bj)=gUNm1(i,j,k,bi,bj)
# Line 87  C     Step forward meridional velocity ( Line 88  C     Step forward meridional velocity (
88    
89        IF (staggerTimeStep) THEN        IF (staggerTimeStep) THEN
90  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:
91          phyFac = pfFacMom*deltaTmom*recip_rhoConst          phyFac = pfFacMom*deltaTmom
92          DO j=jMin,jMax          DO j=jMin,jMax
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)
Legend:
Removed from v.1.22  
changed lines
  Added in v.1.23
  ViewVC Help
Powered by ViewVC 1.1.22