File indexing completed on 2024-01-06 06:09:33 UTC
view on githubraw file Latest commit aef96ed3 on 2024-01-05 19:00:20 UTC
aea29c8517 Alis*0001 #include "CPP_OPTIONS.h"
0002
9366854e02 Chri*0003
0004
0005
aea29c8517 Alis*0006 SUBROUTINE INI_VEL( myThid )
0007
9366854e02 Chri*0008
0009
0010
f9ff59e0c6 Jean*0011
9366854e02 Chri*0012
0013
aea29c8517 Alis*0014
9366854e02 Chri*0015
0016 IMPLICIT NONE
aea29c8517 Alis*0017
0018 #include "SIZE.h"
0019 #include "EEPARAMS.h"
0020 #include "PARAMS.h"
0021 #include "GRID.h"
0022 #include "DYNVARS.h"
0023
9366854e02 Chri*0024
aef96ed361 Jean*0025
aea29c8517 Alis*0026 INTEGER myThid
0027
9366854e02 Chri*0028
b568a02ce1 Jean*0029
0030
aef96ed361 Jean*0031 INTEGER i, j, k
0032 INTEGER bi, bj
9366854e02 Chri*0033
aea29c8517 Alis*0034
0035
0036
0037
0038
0039
0040
0041
0042 DO bj = myByLo(myThid), myByHi(myThid)
0043 DO bi = myBxLo(myThid), myBxHi(myThid)
0044 DO k=1,Nr
aef96ed361 Jean*0045 DO j=1-OLy,sNy+OLy
0046 DO i=1-OLx,sNx+OLx
0047 uVel(i,j,k,bi,bj) = 0. _d 0
0048 vVel(i,j,k,bi,bj) = 0. _d 0
aea29c8517 Alis*0049 ENDDO
0050 ENDDO
0051 ENDDO
0052 ENDDO
0053 ENDDO
0054
aef96ed361 Jean*0055 IF ( uVelInitFile .NE. ' ' .OR. vVelInitFile .NE. ' ' ) THEN
aea29c8517 Alis*0056
aef96ed361 Jean*0057 IF ( uVelInitFile .NE. ' ' )
b568a02ce1 Jean*0058 & CALL READ_FLD_XYZ_RL( uVelInitFile, ' ', uVel, 0, myThid )
aea29c8517 Alis*0059
aef96ed361 Jean*0060 IF ( vVelInitFile .NE. ' ' )
b568a02ce1 Jean*0061 & CALL READ_FLD_XYZ_RL( vVelInitFile, ' ', vVel, 0, myThid )
aea29c8517 Alis*0062
48955d42d9 Alis*0063 CALL EXCH_UV_XYZ_RL(uVel,vVel,.TRUE.,myThid)
b568a02ce1 Jean*0064 ENDIF
48955d42d9 Alis*0065
aea29c8517 Alis*0066 DO bj = myByLo(myThid), myByHi(myThid)
0067 DO bi = myBxLo(myThid), myBxHi(myThid)
0068 DO k=1,Nr
aef96ed361 Jean*0069 DO j=1-OLy,sNy+OLy
0070 DO i=1-OLx,sNx+OLx
0071 uVel(i,j,k,bi,bj) = uVel(i,j,k,bi,bj)*_maskW(i,j,k,bi,bj)
0072 vVel(i,j,k,bi,bj) = vVel(i,j,k,bi,bj)*_maskS(i,j,k,bi,bj)
aea29c8517 Alis*0073 ENDDO
0074 ENDDO
0075 ENDDO
0076 ENDDO
0077 ENDDO
0078
0079 RETURN
0080 END