/[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.6 by cnh, Sat May 30 02:10:16 1998 UTC revision 1.7 by adcroft, Mon Jun 1 20:36:13 1998 UTC
# Line 7  C     | S/R TIMESTEP Line 7  C     | S/R TIMESTEP
7  C     | o Step model fields forward in time                      |  C     | o Step model fields forward in time                      |
8  C     \==========================================================/  C     \==========================================================/
9        SUBROUTINE TIMESTEP( bi, bj, iMin, iMax, jMin, jMax,        SUBROUTINE TIMESTEP( bi, bj, iMin, iMax, jMin, jMax,
10       I                     K, pSurfX, pSurfY,       I                     K,
11       I                     myThid )       I                     myThid )
12        implicit none        implicit none
13  ! Common  ! Common
# Line 15  C     \================================= Line 15  C     \=================================
15  #include "DYNVARS.h"  #include "DYNVARS.h"
16  #include "PARAMS.h"  #include "PARAMS.h"
17  #include "GRID.h"  #include "GRID.h"
 #include "EEPARAMS.h"  
 #include "CG2D.h"  
18  C     == Routine Arguments ==  C     == Routine Arguments ==
19        INTEGER bi,bj,iMin,iMax,jMin,jMax        INTEGER bi,bj,iMin,iMax,jMin,jMax
20        INTEGER K        INTEGER K
       _RL  pSurfX(1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
       _RL  pSurfY(1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
21        INTEGER myThid        INTEGER myThid
22  C     == Local variables ==  C     == Local variables ==
23        INTEGER i,j        INTEGER i,j
24        REAL ab15,ab05        _RL ab15,ab05
       _RL  hxFac, hyFac  
25    
26  C     Adams-Bashforth timestepping weights  C     Adams-Bashforth timestepping weights
27        ab15=1.5+abeps        ab15=1.5+abeps
28        ab05=-0.5-abeps        ab05=-0.5-abeps
29    
30  C     On/off scaling paramters  C     Step forward zonal velocity (store in Gu)
       hxFac = pfFacMom  
       hyFac = pfFacMom  
   
 C     Step forward zonal velocity  
31         DO j=jMin,jMax         DO j=jMin,jMax
32          DO i=iMin,iMax          DO i=iMin,iMax
33           uVel(i,j,k,bi,bj)=uVel(i,j,k,bi,bj)           gUNm1(i,j,k,bi,bj)=uVel(i,j,k,bi,bj)
34       &    +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)
35       &                -hxFac*pSurfX(i,j)/rhonil       &        )*_maskW(i,j,k,bi,bj)
      &    )*_maskW(i,j,k,bi,bj)  
          gUNm1(i,j,k,bi,bj)=gU(i,j,k,bi,bj)  
36          ENDDO          ENDDO
37         ENDDO         ENDDO
38  C     Step forward meridional velocity  C     Step forward meridional velocity (store in Gv)
39         DO j=jMin,jMax         DO j=jMin,jMax
40          DO i=iMin,iMax          DO i=iMin,iMax
41           vVel(i,j,k,bi,bj)=vVel(i,j,k,bi,bj)           gVNm1(i,j,k,bi,bj)=vVel(i,j,k,bi,bj)
42       &    +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)
43       &                -hyFac*pSurfY(i,j)/rhonil       &        )*_maskS(i,j,k,bi,bj)
      &               )*_maskS(i,j,k,bi,bj)  
          gVNm1(i,j,k,bi,bj)=gV(i,j,k,bi,bj)  
44          ENDDO          ENDDO
45         ENDDO         ENDDO
46  C     Step forward temperature  C     Step forward temperature
47         DO j=jMin,jMax         DO j=jMin,jMax
48          DO i=iMin,iMax          DO i=iMin,iMax
49           theta(i,j,k,bi,bj)=theta(i,j,k,bi,bj)           gTNm1(i,j,k,bi,bj)=theta(i,j,k,bi,bj)
50       &    +deltaTtracer*(ab15*gT(i,j,k,bi,bj)+ab05*gTNm1(i,j,k,bi,bj))       &    +deltaTtracer*(ab15*gT(i,j,k,bi,bj)+ab05*gTNm1(i,j,k,bi,bj))
          gTNm1(i,j,k,bi,bj)=gT(i,j,k,bi,bj)  
51          ENDDO          ENDDO
52         ENDDO         ENDDO
53    
       _BARRIER  
 C     CALL PLOT_FIELD_XYZR8( uVel, 'TIEMSTEP.1 uVel',Nz,1,myThid)  
 C     CALL PLOT_FIELD_XYZR8( vVel, 'TIEMSTEP.1 vVel',Nz,1,myThid)  
 C     CALL PLOT_FIELD_XYZR8( theta, 'TIEMSTEP.1 theta',Nz,1,myThid)  
   
   
54        RETURN        RETURN
55        END        END
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22