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

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

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

revision 1.5 by cnh, Wed Sep 26 18:09:16 2001 UTC revision 1.5.6.2 by heimbach, Tue Jun 24 23:05:29 2003 UTC
# Line 11  C     !INTERFACE: Line 11  C     !INTERFACE:
11  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
12  C     *=================================================================  C     *=================================================================
13  C     | SUBROUTINE INI_VEL  C     | SUBROUTINE INI_VEL
14  C     | o Initialize 3D flow field (either to zero or from input files)  C     | o Initialize flow field (either to zero or from input files)
15  C     *=================================================================  C     *=================================================================
16  C     \ev  C     \ev
17    
# Line 26  C     === Global variables === Line 26  C     === Global variables ===
26    
27  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
28  C     == Routine arguments ==  C     == Routine arguments ==
29  C     myThid -  Number of this instance of INI_UVEL  C     myThid -  Number of this instance of INI_VEL
30        INTEGER myThid        INTEGER myThid
31    
32  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
# Line 49  C Line 49  C
49            DO i=1-Olx,sNx+Olx            DO i=1-Olx,sNx+Olx
50             uVel(i,j,k,bi,bj)=0. _d 0             uVel(i,j,k,bi,bj)=0. _d 0
51             vVel(i,j,k,bi,bj)=0. _d 0             vVel(i,j,k,bi,bj)=0. _d 0
            wVel(i,j,k,bi,bj)=0. _d 0  
52            ENDDO            ENDDO
53           ENDDO           ENDDO
54          ENDDO          ENDDO
# Line 73  c      _EXCH_XYZ_R8(vVel  , myThid ) Line 72  c      _EXCH_XYZ_R8(vVel  , myThid )
72    
73  c     IF (uVelInitFile .NE. ' ' .OR. vVelInitFile .NE. ' ') THEN  c     IF (uVelInitFile .NE. ' ' .OR. vVelInitFile .NE. ' ') THEN
74         CALL EXCH_UV_XYZ_RL(uVel,vVel,.TRUE.,myThid)         CALL EXCH_UV_XYZ_RL(uVel,vVel,.TRUE.,myThid)
        _EXCH_XYZ_R8(wVel,myThid)  
75  c     ENDIF  c     ENDIF
76    
77        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
# Line 85  c     ENDIF Line 83  c     ENDIF
83             vVel(i,j,k,bi,bj)=vVel(i,j,k,bi,bj)*_maskS(i,j,k,bi,bj)             vVel(i,j,k,bi,bj)=vVel(i,j,k,bi,bj)*_maskS(i,j,k,bi,bj)
84            ENDDO            ENDDO
85           ENDDO           ENDDO
          CALL INTEGRATE_FOR_W(  
      I                        bi, bj, k, uVel, vVel,  
      O                        wVel,  
      I                        myThid )  
86          ENDDO          ENDDO
87         ENDDO         ENDDO
88        ENDDO        ENDDO
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.5.6.2
  ViewVC Help
Powered by ViewVC 1.1.22