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 |
|
|
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: |
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 |
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) |
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 |