/[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.4 by adcroft, Thu Apr 30 14:03:29 1998 UTC revision 1.5 by cnh, Mon May 25 20:05:55 1998 UTC
# Line 26  C     == Routine Arguments == Line 26  C     == Routine Arguments ==
26  C     == Local variables ==  C     == Local variables ==
27        INTEGER i,j        INTEGER i,j
28        REAL ab15,ab05        REAL ab15,ab05
29          _RL  hxFac, hyFac
30    
31  C     Adams-Bashforth timestepping weights  C     Adams-Bashforth timestepping weights
32        ab15=1.5+abeps        ab15=1.5+abeps
33        ab05=-0.5-abeps        ab05=-0.5-abeps
34    
35    C     On/off scaling paramters
36          hxFac = pfFacMom
37          hyFac = pfFacMom
38    
39  C     Step forward zonal velocity  C     Step forward zonal velocity
40         DO j=jMin,jMax         DO j=jMin,jMax
41          DO i=iMin,iMax          DO i=iMin,iMax
42           uVel(i,j,k,bi,bj)=uVel(i,j,k,bi,bj)           uVel(i,j,k,bi,bj)=uVel(i,j,k,bi,bj)
43       &    +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)
44       &                -pSurfX(i,j)/rhonil       &                -hxFac*pSurfX(i,j)/rhonil
45       &    )*maskW(i,j,k,bi,bj)       &    )*maskW(i,j,k,bi,bj)
46           gUNm1(i,j,k,bi,bj)=gU(i,j,k,bi,bj)           gUNm1(i,j,k,bi,bj)=gU(i,j,k,bi,bj)
47          ENDDO          ENDDO
# Line 46  C     Step forward meridional velocity Line 51  C     Step forward meridional velocity
51          DO i=iMin,iMax          DO i=iMin,iMax
52           vVel(i,j,k,bi,bj)=vVel(i,j,k,bi,bj)           vVel(i,j,k,bi,bj)=vVel(i,j,k,bi,bj)
53       &    +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)
54       &                -pSurfY(i,j)/rhonil       &                -hyFac*pSurfY(i,j)/rhonil
55       &               )*maskS(i,j,k,bi,bj)       &               )*maskS(i,j,k,bi,bj)
56           gVNm1(i,j,k,bi,bj)=gV(i,j,k,bi,bj)           gVNm1(i,j,k,bi,bj)=gV(i,j,k,bi,bj)
57          ENDDO          ENDDO
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22