31 |
CHARACTER*(10) suff |
CHARACTER*(10) suff |
32 |
INTEGER prec |
INTEGER prec |
33 |
|
|
34 |
|
#ifdef ALLOW_DEBUG |
35 |
|
IF (debugMode) CALL DEBUG_ENTER('OBCS_INIT_VARIABLES',myThid) |
36 |
|
#endif |
37 |
|
|
38 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
39 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
40 |
|
|
138 |
#ifdef ALLOW_ORLANSKI |
#ifdef ALLOW_ORLANSKI |
139 |
IF (useOrlanskiNorth.OR.useOrlanskiSouth.OR. |
IF (useOrlanskiNorth.OR.useOrlanskiSouth.OR. |
140 |
& useOrlanskiEast.OR.useOrlanskiWest) THEN |
& useOrlanskiEast.OR.useOrlanskiWest) THEN |
141 |
|
#ifdef ALLOW_DEBUG |
142 |
|
IF (debugMode) CALL DEBUG_CALL('ORLANSKI_INIT',myThid) |
143 |
|
#endif |
144 |
CALL ORLANSKI_INIT(bi, bj, myThid) |
CALL ORLANSKI_INIT(bi, bj, myThid) |
145 |
ENDIF |
ENDIF |
146 |
#endif /* ALLOW_ORLANSKI */ |
#endif /* ALLOW_ORLANSKI */ |
164 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
165 |
|
|
166 |
C-- Apply OBCS values to initial conditions for consistancy |
C-- Apply OBCS values to initial conditions for consistancy |
167 |
|
#ifdef ALLOW_DEBUG |
168 |
|
IF (debugMode) CALL DEBUG_CALL('OBCS_CALC',myThid) |
169 |
|
#endif |
170 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
171 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
172 |
CALL OBCS_CALC( bi, bj, startTime, nIter0, |
CALL OBCS_CALC( bi, bj, startTime, nIter0, |
174 |
ENDDO |
ENDDO |
175 |
ENDDO |
ENDDO |
176 |
|
|
177 |
|
#ifdef ALLOW_DEBUG |
178 |
|
IF (debugMode) CALL DEBUG_CALL('OBCS_APPLY_UV + OBCS_APPLY_TS',myThid) |
179 |
|
#endif |
180 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
181 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
182 |
DO K=1,Nr |
DO K=1,Nr |
190 |
C After applying the boundary conditions exchange the 3D-fields. |
C After applying the boundary conditions exchange the 3D-fields. |
191 |
C This is only necessary of the boudnary values have been read |
C This is only necessary of the boudnary values have been read |
192 |
C from a file. |
C from a file. |
193 |
|
#ifdef ALLOW_DEBUG |
194 |
|
IF (debugMode) CALL DEBUG_CALL('EXCHANGES in OBCS_INIT_VARIABLES',myThid) |
195 |
|
#endif |
196 |
CALL EXCH_UV_XYZ_RL(uVel,vVel,.TRUE.,myThid) |
CALL EXCH_UV_XYZ_RL(uVel,vVel,.TRUE.,myThid) |
197 |
_EXCH_XYZ_R8( theta, myThid ) |
_EXCH_XYZ_R8( theta, myThid ) |
198 |
_EXCH_XYZ_R8( salt , myThid ) |
_EXCH_XYZ_R8( salt , myThid ) |
200 |
|
|
201 |
#endif /* ALLOW_OBCS */ |
#endif /* ALLOW_OBCS */ |
202 |
|
|
203 |
|
#ifdef ALLOW_DEBUG |
204 |
|
IF (debugMode) CALL DEBUG_LEAVE('OBCS_INIT_VARIABLES',myThid) |
205 |
|
#endif |
206 |
RETURN |
RETURN |
207 |
END |
END |