/[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.2 by cnh, Fri Apr 24 02:05:42 1998 UTC revision 1.3 by cnh, Fri Apr 24 03:07:12 1998 UTC
# Line 45  C     Step forward zonal velocity Line 45  C     Step forward zonal velocity
45           uVel(i,j,k,bi,bj)=uVel(i,j,k,bi,bj)           uVel(i,j,k,bi,bj)=uVel(i,j,k,bi,bj)
46       &    +deltaTmom*(ab15*gU(i,j,k,bi,bj)+ab05*gUNm1(i,j,k,bi,bj)       &    +deltaTmom*(ab15*gU(i,j,k,bi,bj)+ab05*gUNm1(i,j,k,bi,bj)
47       &                -pg(i,j)/rhonil       &                -pg(i,j)/rhonil
 CcnhDebugStarts  
 Cdbg &                -rdxC(i,j,bi,bj)*  
 Cdbg &                 (pHSave(i,j,k,bi,bj)-pHSave(i-1,j,k,bi,bj))/rhonil  
 CcnhDebugEnds  
48       &    )*maskW(i,j,k,bi,bj)       &    )*maskW(i,j,k,bi,bj)
49           gUNm1(i,j,k,bi,bj)=gU(i,j,k,bi,bj)           gUNm1(i,j,k,bi,bj)=gU(i,j,k,bi,bj)
50          ENDDO          ENDDO
# Line 69  C     Step forward meridional velocity Line 65  C     Step forward meridional velocity
65           vVel(i,j,k,bi,bj)=vVel(i,j,k,bi,bj)           vVel(i,j,k,bi,bj)=vVel(i,j,k,bi,bj)
66       &    +deltaTmom*(ab15*gV(i,j,k,bi,bj)+ab05*gVNm1(i,j,k,bi,bj)       &    +deltaTmom*(ab15*gV(i,j,k,bi,bj)+ab05*gVNm1(i,j,k,bi,bj)
67       &                -pg(i,j)/rhonil       &                -pg(i,j)/rhonil
 CcnhDebugStarts  
 Cdbg &                -rdyC(i,j,bi,bj)*  
 Cdbg &                 (pHSave(i,j,k,bi,bj)-pHSave(i,j-1,k,bi,bj))/rhonil  
 CcnhDebugEnds  
68       &               )*maskS(i,j,k,bi,bj)       &               )*maskS(i,j,k,bi,bj)
69           gVNm1(i,j,k,bi,bj)=gV(i,j,k,bi,bj)           gVNm1(i,j,k,bi,bj)=gV(i,j,k,bi,bj)
70          ENDDO          ENDDO
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22