/[MITgcm]/MITgcm/model/src/forward_step.F
ViewVC logotype

Diff of /MITgcm/model/src/forward_step.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.20 by adcroft, Thu Sep 13 17:43:55 2001 UTC revision 1.21 by jmc, Wed Sep 19 13:58:08 2001 UTC
# Line 50  C     mythid - thread number for this in Line 50  C     mythid - thread number for this in
50        integer mythid        integer mythid
51        integer myiter        integer myiter
52        _RL     mytime        _RL     mytime
53          INTEGER bi,bj
54    
55  C     == local variables ==  C     == local variables ==
56    
# Line 70  C       (DO_GTERM_BLOCKING_EXCHANGES). Line 71  C       (DO_GTERM_BLOCKING_EXCHANGES).
71          CALL DUMMY_IN_STEPPING( myTime, myIter, myThid )          CALL DUMMY_IN_STEPPING( myTime, myIter, myThid )
72  #endif  #endif
73    
74    #ifdef EXACT_CONSERV
75          IF (exactConserv) THEN
76    C--   Update etaH(n+1) :
77            DO bj=myByLo(myThid),myByHi(myThid)
78             DO bi=myBxLo(myThid),myBxHi(myThid)
79               CALL CALC_EXACT_ETA( .FALSE., bi,bj, uVel,vVel,
80         I                          startTime, nIter0, myThid )
81             ENDDO
82            ENDDO
83            IF (implicDiv2Dflow .NE. 1. _d 0 )
84         &     _EXCH_XY_R8(etaH, myThid )
85          ENDIF
86    #endif /* EXACT_CONSERV */
87    
88  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
89  C--   compute the future surface level thickness  C--   compute the future surface level thickness
90  C      according to the current Eta field  C      according to etaH(n+1)
91          IF ( nonlinFreeSurf.GT.0) THEN          IF ( nonlinFreeSurf.GT.0) THEN
92            CALL CALC_SURF_DR(etaN, myTime, myIter, myThid )            CALL CALC_SURF_DR(etaH, myTime, myIter, myThid )
93          ENDIF          ENDIF
94  #endif  #endif /* NONLIN_FRSURF */
95    
96  C--     Load forcing/external data fields.  C--     Load forcing/external data fields.
97  #ifdef INCLUDE_EXTERNAL_FORCING_PACKAGE  #ifdef INCLUDE_EXTERNAL_FORCING_PACKAGE
# Line 114  C--   Step forward W field in N-H algori Line 129  C--   Step forward W field in N-H algori
129  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
130    
131  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
132  C--   update hfacC,W,S and recip_hFac according to the current Eta field  C--   update hfacC,W,S and recip_hFac according to etaH(n+1) :
133        IF ( momStepping ) THEN        IF ( momStepping ) THEN
134        IF ( nonlinFreeSurf.GT.0) THEN        IF ( nonlinFreeSurf.GT.0) THEN
135          CALL UPDATE_SURF_DR( myTime, myIter, myThid )          CALL UPDATE_SURF_DR( myTime, myIter, myThid )
Legend:
Removed from v.1.20  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.22