/[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.46 by heimbach, Wed Jun 7 01:55:13 2006 UTC revision 1.47 by jmc, Fri Oct 19 14:41:39 2007 UTC
# Line 26  C     == Global variables == Line 26  C     == Global variables ==
26  #include "EEPARAMS.h"  #include "EEPARAMS.h"
27  #include "PARAMS.h"  #include "PARAMS.h"
28  #include "GRID.h"  #include "GRID.h"
29    #include "RESTART.h"
30  #include "SURFACE.h"  #include "SURFACE.h"
31    
32  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
# Line 49  C     gvDissip :: dissipation tendency ( Line 50  C     gvDissip :: dissipation tendency (
50    
51  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
52  C     == Local variables ==  C     == Local variables ==
       LOGICAL momStartAB  
53        INTEGER i,j        INTEGER i,j
54        _RL phxFac,phyFac, psFac        _RL phxFac,phyFac, psFac
55        _RL     gUtmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL     gUtmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
# Line 68  C    Allow diagnosis of external forcing Line 68  C    Allow diagnosis of external forcing
68  #endif  #endif
69  CEOP  CEOP
70    
 C     Start AB with low-order timestepping weights  
       momStartAB = nIter0.EQ.0  
   
71  C-- explicit part of the surface potential gradient is added in this S/R  C-- explicit part of the surface potential gradient is added in this S/R
72        psFac = pfFacMom*(1. _d 0 - implicSurfPress)        psFac = pfFacMom*(1. _d 0 - implicSurfPress)
73    
# Line 181  C     and save gU,gV_[n] into guNm1,gvNm Line 178  C     and save gU,gV_[n] into guNm1,gvNm
178        CALL ADAMS_BASHFORTH3(        CALL ADAMS_BASHFORTH3(
179       I                        bi, bj, k,       I                        bi, bj, k,
180       U                        gU, guNm,       U                        gU, guNm,
181       I                        momStartAB, myIter, myThid )       I                        mom_StartAB, myIter, myThid )
182        CALL ADAMS_BASHFORTH3(        CALL ADAMS_BASHFORTH3(
183       I                        bi, bj, k,       I                        bi, bj, k,
184       U                        gV, gvNm,       U                        gV, gvNm,
185       I                        momStartAB, myIter, myThid )       I                        mom_StartAB, myIter, myThid )
186  #else /* ALLOW_ADAMSBASHFORTH_3 */  #else /* ALLOW_ADAMSBASHFORTH_3 */
187        CALL ADAMS_BASHFORTH2(        CALL ADAMS_BASHFORTH2(
188       I                        bi, bj, k,       I                        bi, bj, k,
189       U                        gU, guNm1,       U                        gU, guNm1,
190       I                        myIter, myThid )       I                        mom_StartAB, myIter, myThid )
191        CALL ADAMS_BASHFORTH2(        CALL ADAMS_BASHFORTH2(
192       I                        bi, bj, k,       I                        bi, bj, k,
193       U                        gV, gvNm1,       U                        gV, gvNm1,
194       I                        myIter, myThid )       I                        mom_StartAB, myIter, myThid )
195  #endif /* ALLOW_ADAMSBASHFORTH_3 */  #endif /* ALLOW_ADAMSBASHFORTH_3 */
196                
197  C--   Forcing term outside the Adams-Bashforth:  C--   Forcing term outside the Adams-Bashforth:
Legend:
Removed from v.1.46  
changed lines
  Added in v.1.47
  ViewVC Help
Powered by ViewVC 1.1.22