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

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

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

revision 1.2 by jmc, Thu Feb 23 20:55:49 2006 UTC revision 1.3 by jmc, Fri Jul 7 20:10:36 2006 UTC
# Line 8  CBOP Line 8  CBOP
8  C     !ROUTINE: TIMESTEP_WVEL  C     !ROUTINE: TIMESTEP_WVEL
9  C     !INTERFACE:  C     !INTERFACE:
10        SUBROUTINE TIMESTEP_WVEL(        SUBROUTINE TIMESTEP_WVEL(
11       I                          myTime, myIter, myThid )       I                          bi,bj, myTime, myIter, myThid )
12  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
13  C     *==========================================================*  C     *==========================================================*
14  C     | S/R TIMESTEP_WVEL  C     | S/R TIMESTEP_WVEL
# Line 29  c #include "SURFACE.h" Line 29  c #include "SURFACE.h"
29    
30  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
31  C     == Routine Arguments ==  C     == Routine Arguments ==
32          INTEGER bi,bj
33        _RL     myTime        _RL     myTime
34        INTEGER myIter, myThid        INTEGER myIter, myThid
35    
36  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
37  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
38  C     == Local variables ==  C     == Local variables ==
       INTEGER k,bi,bj  
39        INTEGER iMin,iMax,jMin,jMax        INTEGER iMin,iMax,jMin,jMax
40        INTEGER i,j        INTEGER i,j,k
41        _RL     gWtmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL     gWtmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
42        _RL     tmpFac, nh_Fac, igwFac        _RL     tmpFac, nh_Fac, igwFac
43  CEOP  CEOP
# Line 52  CEOP Line 52  CEOP
52    
53  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
54    
55  C     tile loops :        IF ( nonHydrostatic ) THEN
       DO bj=myByLo(myThid),myByHi(myThid)  
        DO bi=myBxLo(myThid),myBxHi(myThid)  
   
        IF ( nonHydrostatic ) THEN  
56          nh_Fac = 0.          nh_Fac = 0.
57          IF ( nh_Am2.NE.0. ) nh_Fac = 1. _d 0 / nh_Am2          IF ( nh_Am2.NE.0. ) nh_Fac = 1. _d 0 / nh_Am2
58    
# Line 103  C     impact. This is purely for diagnos Line 99  C     impact. This is purely for diagnos
99          ENDIF          ENDIF
100  #endif /* ALLOW_OBCS */  #endif /* ALLOW_OBCS */
101    
102         ELSEIF ( implicitIntGravWave ) THEN        ELSEIF ( implicitIntGravWave ) THEN
103  C     keep a copy of wVel in gW:  C     keep a copy of wVel in gW:
104          DO k=1,Nr          DO k=1,Nr
105           DO j=1-Oly,sNy+Oly           DO j=1-Oly,sNy+Oly
# Line 114  C     keep a copy of wVel in gW: Line 110  C     keep a copy of wVel in gW:
110          ENDDO          ENDDO
111    
112  C-    End if nonHydrostatic / elseif implicitIntGravWave  C-    End if nonHydrostatic / elseif implicitIntGravWave
113         ENDIF        ENDIF
   
 C-    End of bi,bj loops  
        ENDDO  
       ENDDO  
114    
115  #endif /* ALLOW_NONHYDROSTATIC */  #endif /* ALLOW_NONHYDROSTATIC */
116    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22