pkg/obcs/obcs_init_variables.F

C $Header: /u/gcmpack/MITgcm/pkg/obcs/obcs_init_variables.F,v 1.18 2005/12/15 14:30:45 mlosch Exp $
C $Name:  $

#include "OBCS_OPTIONS.h"

      SUBROUTINE OBCS_INIT_VARIABLES( myThid )
C     /==========================================================\
C     | SUBROUTINE OBCS_INIT_VARIABLES                           |
C     | o Initialise OBCs variable data                          |
C     |==========================================================|
C     |                                                          |
C     \==========================================================/
      IMPLICIT NONE

C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "DYNVARS.h"
#include "OBCS.h"
#ifdef ALLOW_PTRACERS
#include "PTRACERS_SIZE.h"
#include "PTRACERS.h"
#include "OBCS_PTRACERS.h"
#endif /* ALLOW_PTRACERS */

C     == Routine arguments ==
C     myThid -  Number of this instance of INI_DEPTHS
      INTEGER myThid

#ifdef ALLOW_OBCS

C     == Local variables ==
      INTEGER bi, bj
      INTEGER I, J, K
      CHARACTER*(10) suff
      INTEGER prec
#ifdef ALLOW_PTRACERS
      INTEGER iTracer
#endif /* ALLOW_PTRACERS */

#ifdef ALLOW_DEBUG
      IF (debugMode) CALL DEBUG_ENTER('OBCS_INIT_VARIABLES',myThid)
#endif

      DO bj = myByLo(myThid), myByHi(myThid)
       DO bi = myBxLo(myThid), myBxHi(myThid)

        DO K=1,Nr
         DO I=1-Olx,sNx+Olx
#ifdef ALLOW_OBCS_NORTH
          OBNu(I,K,bi,bj)=0. _d 0
          OBNv(I,K,bi,bj)=0. _d 0
          OBNt(I,K,bi,bj)=0. _d 0
          OBNs(I,K,bi,bj)=0. _d 0
# ifdef ALLOW_NONHYDROSTATIC
          OBNw(I,K,bi,bj)=0. _d 0
# endif
# ifdef ALLOW_OBCS_PRESCRIBE
          OBNu0(I,K,bi,bj)=0. _d 0
          OBNv0(I,K,bi,bj)=0. _d 0
          OBNt0(I,K,bi,bj)=0. _d 0
          OBNs0(I,K,bi,bj)=0. _d 0
          OBNu1(I,K,bi,bj)=0. _d 0
          OBNv1(I,K,bi,bj)=0. _d 0
          OBNt1(I,K,bi,bj)=0. _d 0
          OBNs1(I,K,bi,bj)=0. _d 0
# endif
#endif /* ALLOW_OBCS_NORTH */

#ifdef ALLOW_OBCS_SOUTH
          OBSu(I,K,bi,bj)=0. _d 0
          OBSv(I,K,bi,bj)=0. _d 0
          OBSt(I,K,bi,bj)=0. _d 0
          OBSs(I,K,bi,bj)=0. _d 0
# ifdef ALLOW_NONHYDROSTATIC
          OBSw(I,K,bi,bj)=0. _d 0
# endif
# ifdef ALLOW_OBCS_PRESCRIBE
          OBSu0(I,K,bi,bj)=0. _d 0
          OBSv0(I,K,bi,bj)=0. _d 0
          OBSt0(I,K,bi,bj)=0. _d 0
          OBSs0(I,K,bi,bj)=0. _d 0
          OBSu1(I,K,bi,bj)=0. _d 0
          OBSv1(I,K,bi,bj)=0. _d 0
          OBSt1(I,K,bi,bj)=0. _d 0
          OBSs1(I,K,bi,bj)=0. _d 0
# endif
#endif /* ALLOW_OBCS_SOUTH */
         ENDDO

         DO J=1-Oly,sNy+Oly
#ifdef ALLOW_OBCS_EAST
          OBEu(J,K,bi,bj)=0. _d 0
          OBEv(J,K,bi,bj)=0. _d 0
          OBEt(J,K,bi,bj)=0. _d 0
          OBEs(J,K,bi,bj)=0. _d 0
# ifdef ALLOW_NONHYDROSTATIC
          OBEw(J,K,bi,bj)=0. _d 0
# endif
# ifdef ALLOW_OBCS_PRESCRIBE
          OBEu0(J,K,bi,bj)=0. _d 0
          OBEv0(J,K,bi,bj)=0. _d 0
          OBEt0(J,K,bi,bj)=0. _d 0
          OBEs0(J,K,bi,bj)=0. _d 0
          OBEu1(J,K,bi,bj)=0. _d 0
          OBEv1(J,K,bi,bj)=0. _d 0
          OBEt1(J,K,bi,bj)=0. _d 0
          OBEs1(J,K,bi,bj)=0. _d 0
# endif
#endif /* ALLOW_OBCS_EAST */

#ifdef ALLOW_OBCS_WEST
          OBWu(J,K,bi,bj)=0. _d 0
          OBWv(J,K,bi,bj)=0. _d 0
          OBWt(J,K,bi,bj)=0. _d 0
          OBWs(J,K,bi,bj)=0. _d 0
# ifdef ALLOW_NONHYDROSTATIC
          OBWw(J,K,bi,bj)=0. _d 0
# endif
# ifdef ALLOW_OBCS_PRESCRIBE
          OBWu0(J,K,bi,bj)=0. _d 0
          OBWv0(J,K,bi,bj)=0. _d 0
          OBWt0(J,K,bi,bj)=0. _d 0
          OBWs0(J,K,bi,bj)=0. _d 0
          OBWu1(J,K,bi,bj)=0. _d 0
          OBWv1(J,K,bi,bj)=0. _d 0
          OBWt1(J,K,bi,bj)=0. _d 0
          OBWs1(J,K,bi,bj)=0. _d 0
# endif
#endif /* ALLOW_OBCS_WEST */
         ENDDO
        ENDDO

#ifdef ALLOW_PTRACERS
        IF ( usePTRACERS ) THEN
         DO iTracer=1,PTRACERS_numInUse
          DO K=1,Nr
           DO I=1-Olx,sNx+Olx
#ifdef ALLOW_OBCS_NORTH
            OBNptr (I,K,bi,bj,iTracer)=0. _d 0
# ifdef ALLOW_OBCS_PRESCRIBE
            OBNptr0(I,K,bi,bj,iTracer)=0. _d 0
            OBNptr1(I,K,bi,bj,iTracer)=0. _d 0
# endif
#endif /* ALLOW_OBCS_NORTH */
            
#ifdef ALLOW_OBCS_SOUTH
            OBSptr (I,K,bi,bj,iTracer)=0. _d 0
# ifdef ALLOW_OBCS_PRESCRIBE
            OBSptr0(I,K,bi,bj,iTracer)=0. _d 0
            OBSptr1(I,K,bi,bj,iTracer)=0. _d 0
# endif
#endif /* ALLOW_OBCS_SOUTH */
           ENDDO
           
           DO J=1-Oly,sNy+Oly
#ifdef ALLOW_OBCS_EAST
            OBEptr (J,K,bi,bj,iTracer)=0. _d 0
# ifdef ALLOW_OBCS_PRESCRIBE
            OBEptr0(J,K,bi,bj,iTracer)=0. _d 0
            OBEptr1(J,K,bi,bj,iTracer)=0. _d 0
# endif
#endif /* ALLOW_OBCS_EAST */
            
#ifdef ALLOW_OBCS_WEST
            OBWptr (J,K,bi,bj,iTracer)=0. _d 0
# ifdef ALLOW_OBCS_PRESCRIBE
            OBWptr0(J,K,bi,bj,iTracer)=0. _d 0
            OBWptr1(J,K,bi,bj,iTracer)=0. _d 0
# endif
#endif /* ALLOW_OBCS_WEST */
           ENDDO
          ENDDO
         ENDDO
        ENDIF
#endif /* ALLOW_PTRACERS */

#ifdef NONLIN_FRSURF
        DO I=1-Olx,sNx+Olx
          OBNeta(I,bi,bj)=0.
          OBSeta(I,bi,bj)=0.
        ENDDO
        DO J=1-Oly,sNy+Oly
          OBEeta(J,bi,bj)=0.
          OBWeta(J,bi,bj)=0.
        ENDDO
#endif /* NONLIN_FRSURF */

#ifdef ALLOW_ORLANSKI
        IF (useOrlanskiNorth.OR.useOrlanskiSouth.OR.
     &      useOrlanskiEast.OR.useOrlanskiWest) THEN
#ifdef ALLOW_DEBUG
      IF (debugMode) CALL DEBUG_CALL('ORLANSKI_INIT',myThid)
#endif
          CALL ORLANSKI_INIT(bi, bj, myThid)
        ENDIF
#endif /* ALLOW_ORLANSKI */

       ENDDO
      ENDDO     

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
C jmc: here is the logical place to read OBCS-pickup files
C      but a) without Orlanski: pass the test 1+1=2 without reading pickup.
C          b) with Orlanski: 1+1=2 fail even with this bit of code
      IF ( nIter0.NE.0 ) THEN
        prec = precFloat64
        IF (pickupSuff.EQ.' ') THEN
         WRITE(suff,'(I10.10)') nIter0
        ELSE
         WRITE(suff,'(A10)') pickupSuff
        ENDIF
c       CALL OBCS_READ_CHECKPOINT(prec, nIter0, suff, myThid)
      ENDIF
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

C--   Load/compute OBCS values in all cases, although these values are
C--   only used for initialisation and never for restart:
#ifdef ALLOW_DEBUG
      IF (debugMode) CALL DEBUG_CALL('OBCS_CALC',myThid)
#endif
      DO bj = myByLo(myThid), myByHi(myThid)
       DO bi = myBxLo(myThid), myBxHi(myThid)
        CALL OBCS_CALC( bi, bj, startTime, nIter0,
     &        uVel, vVel, wVel, theta, salt, myThid )
       ENDDO
      ENDDO

      IF ( startTime .EQ. baseTime .AND.  nIter0 .EQ. 0
     &     .AND. pickupSuff .EQ. ' ' ) THEN
C--   Apply OBCS values to initial conditions for consistency
C      (but initial conditions only)
#ifdef ALLOW_DEBUG
       IF (debugMode) 
     &      CALL DEBUG_CALL('OBCS_APPLY_UV + OBCS_APPLY_TS',myThid)
#endif
       DO bj = myByLo(myThid), myByHi(myThid)
        DO bi = myBxLo(myThid), myBxHi(myThid)
         DO K=1,Nr
          CALL OBCS_APPLY_UV( bi, bj, k, uVel, vVel, myThid )
          CALL OBCS_APPLY_TS( bi, bj, k, theta, salt, myThid )
         ENDDO
        ENDDO
       ENDDO
       IF (useOBCSprescribe) THEN
C     After applying the boundary conditions exchange the 3D-fields.
C     This is only necessary of the boudnary values have been read
C     from a file.
#ifdef ALLOW_DEBUG
        IF (debugMode) 
     &       CALL DEBUG_CALL('EXCHANGES in OBCS_INIT_VARIABLES',myThid)
#endif
        CALL EXCH_UV_XYZ_RL(uVel,vVel,.TRUE.,myThid)
        _EXCH_XYZ_R8( theta, myThid )
        _EXCH_XYZ_R8( salt , myThid )
       ENDIF
C     endif start from rest
      ENDIF
#ifdef ALLOW_PTRACERS
C     repeat everything for passive tracers
      IF ( usePTRACERS ) THEN
C     catch the case when we do start from a pickup for dynamics variables
C     but initialise ptracers differently
       IF ( nIter0 .EQ. PTRACERS_Iter0 ) THEN
#ifdef ALLOW_DEBUG
        IF (debugMode) 
     &       CALL DEBUG_CALL('OBCS_APPLY_PTRACER',myThid)
#endif
        DO iTracer=1,PTRACERS_numInUse
         DO bj = myByLo(myThid), myByHi(myThid)
          DO bi = myBxLo(myThid), myBxHi(myThid)
           DO K=1,Nr
            CALL OBCS_APPLY_PTRACER( 
     I           bi, bj, K, iTracer, 
     U           ptracer(1-Olx,1-Oly,K,bi,bj,iTracer),
     I           myThid )
           ENDDO
          ENDDO
         ENDDO
        ENDDO
        IF (useOBCSprescribe) THEN
C     After applying the boundary conditions exchange the 3D-fields.
C     This is only necessary of the boudnary values have been read
C     from a file.
#ifdef ALLOW_DEBUG
         IF (debugMode) CALL DEBUG_CALL( 
     &        'PTRACERS EXCHANGES in OBCS_INIT_VARIABLES',myThid)
#endif
         CALL PTRACERS_FIELDS_BLOCKING_EXCH( myThid )
        ENDIF
C     endif start from rest
       ENDIF
C     endif usePTRACERS
      ENDIF
#endif /* ALLOW_PTRACERS */

#endif /* ALLOW_OBCS */

#ifdef ALLOW_DEBUG
      IF (debugMode) CALL DEBUG_LEAVE('OBCS_INIT_VARIABLES',myThid)
#endif
      RETURN
      END
1	jmc	1.19	C $Header: /u/gcmpack/MITgcm/pkg/obcs/obcs_init_variables.F,v 1.18 2005/12/15 14:30:45 mlosch Exp $
2	adcroft	1.3	C $Name: $
3	adcroft	1.2
4			#include "OBCS_OPTIONS.h"
5
6			SUBROUTINE OBCS_INIT_VARIABLES( myThid )
7			C /==========================================================\
8			C \| SUBROUTINE OBCS_INIT_VARIABLES \|
9			C \| o Initialise OBCs variable data \|
10			C \|==========================================================\|
11			C \| \|
12			C \==========================================================/
13			IMPLICIT NONE
14
15			C === Global variables ===
16			#include "SIZE.h"
17			#include "EEPARAMS.h"
18			#include "PARAMS.h"
19			#include "DYNVARS.h"
20			#include "OBCS.h"
21	jmc	1.17	#ifdef ALLOW_PTRACERS
22	mlosch	1.14	#include "PTRACERS_SIZE.h"
23			#include "PTRACERS.h"
24			#include "OBCS_PTRACERS.h"
25			#endif /* ALLOW_PTRACERS */
26	adcroft	1.2
27			C == Routine arguments ==
28			C myThid - Number of this instance of INI_DEPTHS
29			INTEGER myThid
30
31			#ifdef ALLOW_OBCS
32
33			C == Local variables ==
34			INTEGER bi, bj
35			INTEGER I, J, K
36	adcroft	1.9	CHARACTER*(10) suff
37			INTEGER prec
38	mlosch	1.14	#ifdef ALLOW_PTRACERS
39			INTEGER iTracer
40			#endif /* ALLOW_PTRACERS */
41	adcroft	1.2
42	adcroft	1.12	#ifdef ALLOW_DEBUG
43			IF (debugMode) CALL DEBUG_ENTER('OBCS_INIT_VARIABLES',myThid)
44			#endif
45
46	adcroft	1.2	DO bj = myByLo(myThid), myByHi(myThid)
47			DO bi = myBxLo(myThid), myBxHi(myThid)
48
49			DO K=1,Nr
50			DO I=1-Olx,sNx+Olx
51	heimbach	1.10	#ifdef ALLOW_OBCS_NORTH
52			OBNu(I,K,bi,bj)=0. _d 0
53			OBNv(I,K,bi,bj)=0. _d 0
54			OBNt(I,K,bi,bj)=0. _d 0
55			OBNs(I,K,bi,bj)=0. _d 0
56			# ifdef ALLOW_NONHYDROSTATIC
57			OBNw(I,K,bi,bj)=0. _d 0
58			# endif
59			# ifdef ALLOW_OBCS_PRESCRIBE
60			OBNu0(I,K,bi,bj)=0. _d 0
61			OBNv0(I,K,bi,bj)=0. _d 0
62			OBNt0(I,K,bi,bj)=0. _d 0
63			OBNs0(I,K,bi,bj)=0. _d 0
64			OBNu1(I,K,bi,bj)=0. _d 0
65			OBNv1(I,K,bi,bj)=0. _d 0
66			OBNt1(I,K,bi,bj)=0. _d 0
67			OBNs1(I,K,bi,bj)=0. _d 0
68			# endif
69			#endif /* ALLOW_OBCS_NORTH */
70
71			#ifdef ALLOW_OBCS_SOUTH
72			OBSu(I,K,bi,bj)=0. _d 0
73			OBSv(I,K,bi,bj)=0. _d 0
74			OBSt(I,K,bi,bj)=0. _d 0
75			OBSs(I,K,bi,bj)=0. _d 0
76			# ifdef ALLOW_NONHYDROSTATIC
77			OBSw(I,K,bi,bj)=0. _d 0
78			# endif
79			# ifdef ALLOW_OBCS_PRESCRIBE
80			OBSu0(I,K,bi,bj)=0. _d 0
81			OBSv0(I,K,bi,bj)=0. _d 0
82			OBSt0(I,K,bi,bj)=0. _d 0
83			OBSs0(I,K,bi,bj)=0. _d 0
84			OBSu1(I,K,bi,bj)=0. _d 0
85			OBSv1(I,K,bi,bj)=0. _d 0
86			OBSt1(I,K,bi,bj)=0. _d 0
87			OBSs1(I,K,bi,bj)=0. _d 0
88			# endif
89			#endif /* ALLOW_OBCS_SOUTH */
90	adcroft	1.2	ENDDO
91	heimbach	1.10
92	adcroft	1.2	DO J=1-Oly,sNy+Oly
93	heimbach	1.10	#ifdef ALLOW_OBCS_EAST
94			OBEu(J,K,bi,bj)=0. _d 0
95			OBEv(J,K,bi,bj)=0. _d 0
96			OBEt(J,K,bi,bj)=0. _d 0
97			OBEs(J,K,bi,bj)=0. _d 0
98			# ifdef ALLOW_NONHYDROSTATIC
99			OBEw(J,K,bi,bj)=0. _d 0
100			# endif
101			# ifdef ALLOW_OBCS_PRESCRIBE
102			OBEu0(J,K,bi,bj)=0. _d 0
103			OBEv0(J,K,bi,bj)=0. _d 0
104			OBEt0(J,K,bi,bj)=0. _d 0
105			OBEs0(J,K,bi,bj)=0. _d 0
106			OBEu1(J,K,bi,bj)=0. _d 0
107			OBEv1(J,K,bi,bj)=0. _d 0
108			OBEt1(J,K,bi,bj)=0. _d 0
109			OBEs1(J,K,bi,bj)=0. _d 0
110			# endif
111			#endif /* ALLOW_OBCS_EAST */
112
113			#ifdef ALLOW_OBCS_WEST
114			OBWu(J,K,bi,bj)=0. _d 0
115			OBWv(J,K,bi,bj)=0. _d 0
116			OBWt(J,K,bi,bj)=0. _d 0
117			OBWs(J,K,bi,bj)=0. _d 0
118			# ifdef ALLOW_NONHYDROSTATIC
119			OBWw(J,K,bi,bj)=0. _d 0
120			# endif
121			# ifdef ALLOW_OBCS_PRESCRIBE
122			OBWu0(J,K,bi,bj)=0. _d 0
123			OBWv0(J,K,bi,bj)=0. _d 0
124			OBWt0(J,K,bi,bj)=0. _d 0
125			OBWs0(J,K,bi,bj)=0. _d 0
126			OBWu1(J,K,bi,bj)=0. _d 0
127			OBWv1(J,K,bi,bj)=0. _d 0
128			OBWt1(J,K,bi,bj)=0. _d 0
129			OBWs1(J,K,bi,bj)=0. _d 0
130			# endif
131			#endif /* ALLOW_OBCS_WEST */
132	adcroft	1.2	ENDDO
133			ENDDO
134	jmc	1.5
135	mlosch	1.14	#ifdef ALLOW_PTRACERS
136	mlosch	1.16	IF ( usePTRACERS ) THEN
137			DO iTracer=1,PTRACERS_numInUse
138			DO K=1,Nr
139			DO I=1-Olx,sNx+Olx
140	mlosch	1.14	#ifdef ALLOW_OBCS_NORTH
141	mlosch	1.16	OBNptr (I,K,bi,bj,iTracer)=0. _d 0
142	mlosch	1.14	# ifdef ALLOW_OBCS_PRESCRIBE
143	mlosch	1.16	OBNptr0(I,K,bi,bj,iTracer)=0. _d 0
144			OBNptr1(I,K,bi,bj,iTracer)=0. _d 0
145	mlosch	1.14	# endif
146			#endif /* ALLOW_OBCS_NORTH */
147	mlosch	1.16
148	mlosch	1.14	#ifdef ALLOW_OBCS_SOUTH
149	mlosch	1.16	OBSptr (I,K,bi,bj,iTracer)=0. _d 0
150	mlosch	1.14	# ifdef ALLOW_OBCS_PRESCRIBE
151	mlosch	1.16	OBSptr0(I,K,bi,bj,iTracer)=0. _d 0
152			OBSptr1(I,K,bi,bj,iTracer)=0. _d 0
153	mlosch	1.14	# endif
154			#endif /* ALLOW_OBCS_SOUTH */
155	mlosch	1.16	ENDDO
156
157			DO J=1-Oly,sNy+Oly
158	mlosch	1.14	#ifdef ALLOW_OBCS_EAST
159	mlosch	1.16	OBEptr (J,K,bi,bj,iTracer)=0. _d 0
160	mlosch	1.14	# ifdef ALLOW_OBCS_PRESCRIBE
161	mlosch	1.16	OBEptr0(J,K,bi,bj,iTracer)=0. _d 0
162			OBEptr1(J,K,bi,bj,iTracer)=0. _d 0
163	mlosch	1.14	# endif
164			#endif /* ALLOW_OBCS_EAST */
165	mlosch	1.16
166	mlosch	1.14	#ifdef ALLOW_OBCS_WEST
167	mlosch	1.16	OBWptr (J,K,bi,bj,iTracer)=0. _d 0
168	mlosch	1.14	# ifdef ALLOW_OBCS_PRESCRIBE
169	mlosch	1.16	OBWptr0(J,K,bi,bj,iTracer)=0. _d 0
170			OBWptr1(J,K,bi,bj,iTracer)=0. _d 0
171	mlosch	1.14	# endif
172			#endif /* ALLOW_OBCS_WEST */
173	mlosch	1.16	ENDDO
174	mlosch	1.14	ENDDO
175			ENDDO
176	mlosch	1.16	ENDIF
177	mlosch	1.14	#endif /* ALLOW_PTRACERS */
178
179	jmc	1.5	#ifdef NONLIN_FRSURF
180			DO I=1-Olx,sNx+Olx
181			OBNeta(I,bi,bj)=0.
182			OBSeta(I,bi,bj)=0.
183			ENDDO
184			DO J=1-Oly,sNy+Oly
185			OBEeta(J,bi,bj)=0.
186			OBWeta(J,bi,bj)=0.
187			ENDDO
188			#endif /* NONLIN_FRSURF */
189	adcroft	1.2
190			#ifdef ALLOW_ORLANSKI
191			IF (useOrlanskiNorth.OR.useOrlanskiSouth.OR.
192			& useOrlanskiEast.OR.useOrlanskiWest) THEN
193	adcroft	1.12	#ifdef ALLOW_DEBUG
194			IF (debugMode) CALL DEBUG_CALL('ORLANSKI_INIT',myThid)
195			#endif
196	adcroft	1.2	CALL ORLANSKI_INIT(bi, bj, myThid)
197			ENDIF
198			#endif /* ALLOW_ORLANSKI */
199	adcroft	1.3
200			ENDDO
201			ENDDO
202	jmc	1.7
203	adcroft	1.9	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
204			C jmc: here is the logical place to read OBCS-pickup files
205			C but a) without Orlanski: pass the test 1+1=2 without reading pickup.
206			C b) with Orlanski: 1+1=2 fail even with this bit of code
207			IF ( nIter0.NE.0 ) THEN
208			prec = precFloat64
209			IF (pickupSuff.EQ.' ') THEN
210			WRITE(suff,'(I10.10)') nIter0
211			ELSE
212			WRITE(suff,'(A10)') pickupSuff
213			ENDIF
214			c CALL OBCS_READ_CHECKPOINT(prec, nIter0, suff, myThid)
215			ENDIF
216			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
217
218	mlosch	1.18	C-- Load/compute OBCS values in all cases, although these values are
219			C-- only used for initialisation and never for restart:
220	adcroft	1.12	#ifdef ALLOW_DEBUG
221			IF (debugMode) CALL DEBUG_CALL('OBCS_CALC',myThid)
222			#endif
223	adcroft	1.2	DO bj = myByLo(myThid), myByHi(myThid)
224			DO bi = myBxLo(myThid), myBxHi(myThid)
225	jmc	1.6	CALL OBCS_CALC( bi, bj, startTime, nIter0,
226	adcroft	1.2	& uVel, vVel, wVel, theta, salt, myThid )
227	heimbach	1.10	ENDDO
228			ENDDO
229
230	mlosch	1.15	IF ( startTime .EQ. baseTime .AND. nIter0 .EQ. 0
231			& .AND. pickupSuff .EQ. ' ' ) THEN
232	jmc	1.17	C-- Apply OBCS values to initial conditions for consistency
233			C (but initial conditions only)
234	adcroft	1.12	#ifdef ALLOW_DEBUG
235	mlosch	1.15	IF (debugMode)
236			& CALL DEBUG_CALL('OBCS_APPLY_UV + OBCS_APPLY_TS',myThid)
237	adcroft	1.12	#endif
238	mlosch	1.15	DO bj = myByLo(myThid), myByHi(myThid)
239			DO bi = myBxLo(myThid), myBxHi(myThid)
240			DO K=1,Nr
241			CALL OBCS_APPLY_UV( bi, bj, k, uVel, vVel, myThid )
242			CALL OBCS_APPLY_TS( bi, bj, k, theta, salt, myThid )
243			ENDDO
244			ENDDO
245			ENDDO
246	mlosch	1.18	IF (useOBCSprescribe) THEN
247			C After applying the boundary conditions exchange the 3D-fields.
248			C This is only necessary of the boudnary values have been read
249			C from a file.
250			#ifdef ALLOW_DEBUG
251			IF (debugMode)
252			& CALL DEBUG_CALL('EXCHANGES in OBCS_INIT_VARIABLES',myThid)
253			#endif
254			CALL EXCH_UV_XYZ_RL(uVel,vVel,.TRUE.,myThid)
255			_EXCH_XYZ_R8( theta, myThid )
256			_EXCH_XYZ_R8( salt , myThid )
257			ENDIF
258			C endif start from rest
259			ENDIF
260	mlosch	1.14	#ifdef ALLOW_PTRACERS
261	mlosch	1.15	C repeat everything for passive tracers
262	mlosch	1.18	IF ( usePTRACERS ) THEN
263			C catch the case when we do start from a pickup for dynamics variables
264			C but initialise ptracers differently
265			IF ( nIter0 .EQ. PTRACERS_Iter0 ) THEN
266	mlosch	1.15	#ifdef ALLOW_DEBUG
267			IF (debugMode)
268			& CALL DEBUG_CALL('OBCS_APPLY_PTRACER',myThid)
269			#endif
270			DO iTracer=1,PTRACERS_numInUse
271			DO bj = myByLo(myThid), myByHi(myThid)
272			DO bi = myBxLo(myThid), myBxHi(myThid)
273			DO K=1,Nr
274			CALL OBCS_APPLY_PTRACER(
275			I bi, bj, K, iTracer,
276	mlosch	1.18	U ptracer(1-Olx,1-Oly,K,bi,bj,iTracer),
277	mlosch	1.15	I myThid )
278			ENDDO
279			ENDDO
280	mlosch	1.14	ENDDO
281	mlosch	1.15	ENDDO
282	mlosch	1.18	IF (useOBCSprescribe) THEN
283	mlosch	1.11	C After applying the boundary conditions exchange the 3D-fields.
284			C This is only necessary of the boudnary values have been read
285			C from a file.
286	adcroft	1.12	#ifdef ALLOW_DEBUG
287	mlosch	1.18	IF (debugMode) CALL DEBUG_CALL(
288			& 'PTRACERS EXCHANGES in OBCS_INIT_VARIABLES',myThid)
289	adcroft	1.12	#endif
290	jmc	1.19	CALL PTRACERS_FIELDS_BLOCKING_EXCH( myThid )
291	mlosch	1.18	ENDIF
292			C endif start from rest
293	mlosch	1.15	ENDIF
294	mlosch	1.18	C endif usePTRACERS
295			ENDIF
296	mlosch	1.15	#endif /* ALLOW_PTRACERS */
297	mlosch	1.11
298	adcroft	1.2	#endif /* ALLOW_OBCS */
299	heimbach	1.10
300	adcroft	1.12	#ifdef ALLOW_DEBUG
301			IF (debugMode) CALL DEBUG_LEAVE('OBCS_INIT_VARIABLES',myThid)
302			#endif
303	adcroft	1.2	RETURN
304			END