/[MITgcm]/MITgcm/pkg/obcs/obcs_init_variables.F
ViewVC logotype

Diff of /MITgcm/pkg/obcs/obcs_init_variables.F

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

revision 1.9 by adcroft, Thu Jul 8 16:54:43 2004 UTC revision 1.16 by mlosch, Wed Dec 14 18:01:49 2005 UTC
# Line 18  C     === Global variables === Line 18  C     === Global variables ===
18  #include "PARAMS.h"  #include "PARAMS.h"
19  #include "DYNVARS.h"  #include "DYNVARS.h"
20  #include "OBCS.h"  #include "OBCS.h"
21    #ifdef ALLOW_PTRACERS.h
22    #include "PTRACERS_SIZE.h"
23    #include "PTRACERS.h"
24    #include "OBCS_PTRACERS.h"
25    #endif /* ALLOW_PTRACERS */
26    
27  C     == Routine arguments ==  C     == Routine arguments ==
28  C     myThid -  Number of this instance of INI_DEPTHS  C     myThid -  Number of this instance of INI_DEPTHS
# Line 30  C     == Local variables == Line 35  C     == Local variables ==
35        INTEGER I, J, K        INTEGER I, J, K
36        CHARACTER*(10) suff        CHARACTER*(10) suff
37        INTEGER prec        INTEGER prec
38    #ifdef ALLOW_PTRACERS
39          INTEGER iTracer
40    CML      _RL ptracerLoc(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
41    #endif /* ALLOW_PTRACERS */
42    
43    #ifdef ALLOW_DEBUG
44          IF (debugMode) CALL DEBUG_ENTER('OBCS_INIT_VARIABLES',myThid)
45    #endif
46    
47        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
48         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
49    
50          DO K=1,Nr          DO K=1,Nr
51           DO I=1-Olx,sNx+Olx           DO I=1-Olx,sNx+Olx
52            OBNu(I,K,bi,bj)=0.  #ifdef ALLOW_OBCS_NORTH
53            OBNv(I,K,bi,bj)=0.            OBNu(I,K,bi,bj)=0. _d 0
54            OBNt(I,K,bi,bj)=0.            OBNv(I,K,bi,bj)=0. _d 0
55            OBNs(I,K,bi,bj)=0.            OBNt(I,K,bi,bj)=0. _d 0
56            OBSu(I,K,bi,bj)=0.            OBNs(I,K,bi,bj)=0. _d 0
57            OBSv(I,K,bi,bj)=0.  # ifdef ALLOW_NONHYDROSTATIC
58            OBSt(I,K,bi,bj)=0.            OBNw(I,K,bi,bj)=0. _d 0
59            OBSs(I,K,bi,bj)=0.  # endif
60  #ifdef ALLOW_NONHYDROSTATIC  # ifdef ALLOW_OBCS_PRESCRIBE
61            OBNw(I,K,bi,bj)=0.            OBNu0(I,K,bi,bj)=0. _d 0
62            OBSw(I,K,bi,bj)=0.            OBNv0(I,K,bi,bj)=0. _d 0
63  #endif            OBNt0(I,K,bi,bj)=0. _d 0
64              OBNs0(I,K,bi,bj)=0. _d 0
65              OBNu1(I,K,bi,bj)=0. _d 0
66              OBNv1(I,K,bi,bj)=0. _d 0
67              OBNt1(I,K,bi,bj)=0. _d 0
68              OBNs1(I,K,bi,bj)=0. _d 0
69    # endif
70    #endif /* ALLOW_OBCS_NORTH */
71    
72    #ifdef ALLOW_OBCS_SOUTH
73              OBSu(I,K,bi,bj)=0. _d 0
74              OBSv(I,K,bi,bj)=0. _d 0
75              OBSt(I,K,bi,bj)=0. _d 0
76              OBSs(I,K,bi,bj)=0. _d 0
77    # ifdef ALLOW_NONHYDROSTATIC
78              OBSw(I,K,bi,bj)=0. _d 0
79    # endif
80    # ifdef ALLOW_OBCS_PRESCRIBE
81              OBSu0(I,K,bi,bj)=0. _d 0
82              OBSv0(I,K,bi,bj)=0. _d 0
83              OBSt0(I,K,bi,bj)=0. _d 0
84              OBSs0(I,K,bi,bj)=0. _d 0
85              OBSu1(I,K,bi,bj)=0. _d 0
86              OBSv1(I,K,bi,bj)=0. _d 0
87              OBSt1(I,K,bi,bj)=0. _d 0
88              OBSs1(I,K,bi,bj)=0. _d 0
89    # endif
90    #endif /* ALLOW_OBCS_SOUTH */
91           ENDDO           ENDDO
92    
93           DO J=1-Oly,sNy+Oly           DO J=1-Oly,sNy+Oly
94            OBEu(J,K,bi,bj)=0.  #ifdef ALLOW_OBCS_EAST
95            OBEv(J,K,bi,bj)=0.            OBEu(J,K,bi,bj)=0. _d 0
96            OBEt(J,K,bi,bj)=0.            OBEv(J,K,bi,bj)=0. _d 0
97            OBEs(J,K,bi,bj)=0.            OBEt(J,K,bi,bj)=0. _d 0
98            OBWu(J,K,bi,bj)=0.            OBEs(J,K,bi,bj)=0. _d 0
99            OBWv(J,K,bi,bj)=0.  # ifdef ALLOW_NONHYDROSTATIC
100            OBWt(J,K,bi,bj)=0.            OBEw(J,K,bi,bj)=0. _d 0
101            OBWs(J,K,bi,bj)=0.  # endif
102  #ifdef ALLOW_NONHYDROSTATIC  # ifdef ALLOW_OBCS_PRESCRIBE
103            OBEw(J,K,bi,bj)=0.            OBEu0(J,K,bi,bj)=0. _d 0
104            OBWw(J,K,bi,bj)=0.            OBEv0(J,K,bi,bj)=0. _d 0
105  #endif            OBEt0(J,K,bi,bj)=0. _d 0
106              OBEs0(J,K,bi,bj)=0. _d 0
107              OBEu1(J,K,bi,bj)=0. _d 0
108              OBEv1(J,K,bi,bj)=0. _d 0
109              OBEt1(J,K,bi,bj)=0. _d 0
110              OBEs1(J,K,bi,bj)=0. _d 0
111    # endif
112    #endif /* ALLOW_OBCS_EAST */
113    
114    #ifdef ALLOW_OBCS_WEST
115              OBWu(J,K,bi,bj)=0. _d 0
116              OBWv(J,K,bi,bj)=0. _d 0
117              OBWt(J,K,bi,bj)=0. _d 0
118              OBWs(J,K,bi,bj)=0. _d 0
119    # ifdef ALLOW_NONHYDROSTATIC
120              OBWw(J,K,bi,bj)=0. _d 0
121    # endif
122    # ifdef ALLOW_OBCS_PRESCRIBE
123              OBWu0(J,K,bi,bj)=0. _d 0
124              OBWv0(J,K,bi,bj)=0. _d 0
125              OBWt0(J,K,bi,bj)=0. _d 0
126              OBWs0(J,K,bi,bj)=0. _d 0
127              OBWu1(J,K,bi,bj)=0. _d 0
128              OBWv1(J,K,bi,bj)=0. _d 0
129              OBWt1(J,K,bi,bj)=0. _d 0
130              OBWs1(J,K,bi,bj)=0. _d 0
131    # endif
132    #endif /* ALLOW_OBCS_WEST */
133           ENDDO           ENDDO
134          ENDDO          ENDDO
135    
136    #ifdef ALLOW_PTRACERS
137            IF ( usePTRACERS ) THEN
138             DO iTracer=1,PTRACERS_numInUse
139              DO K=1,Nr
140               DO I=1-Olx,sNx+Olx
141    #ifdef ALLOW_OBCS_NORTH
142                OBNptr (I,K,bi,bj,iTracer)=0. _d 0
143    # ifdef ALLOW_OBCS_PRESCRIBE
144                OBNptr0(I,K,bi,bj,iTracer)=0. _d 0
145                OBNptr1(I,K,bi,bj,iTracer)=0. _d 0
146    # endif
147    #endif /* ALLOW_OBCS_NORTH */
148                
149    #ifdef ALLOW_OBCS_SOUTH
150                OBSptr (I,K,bi,bj,iTracer)=0. _d 0
151    # ifdef ALLOW_OBCS_PRESCRIBE
152                OBSptr0(I,K,bi,bj,iTracer)=0. _d 0
153                OBSptr1(I,K,bi,bj,iTracer)=0. _d 0
154    # endif
155    #endif /* ALLOW_OBCS_SOUTH */
156               ENDDO
157              
158               DO J=1-Oly,sNy+Oly
159    #ifdef ALLOW_OBCS_EAST
160                OBEptr (J,K,bi,bj,iTracer)=0. _d 0
161    # ifdef ALLOW_OBCS_PRESCRIBE
162                OBEptr0(J,K,bi,bj,iTracer)=0. _d 0
163                OBEptr1(J,K,bi,bj,iTracer)=0. _d 0
164    # endif
165    #endif /* ALLOW_OBCS_EAST */
166                
167    #ifdef ALLOW_OBCS_WEST
168                OBWptr (J,K,bi,bj,iTracer)=0. _d 0
169    # ifdef ALLOW_OBCS_PRESCRIBE
170                OBWptr0(J,K,bi,bj,iTracer)=0. _d 0
171                OBWptr1(J,K,bi,bj,iTracer)=0. _d 0
172    # endif
173    #endif /* ALLOW_OBCS_WEST */
174               ENDDO
175              ENDDO
176             ENDDO
177            ENDIF
178    #endif /* ALLOW_PTRACERS */
179    
180  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
181          DO I=1-Olx,sNx+Olx          DO I=1-Olx,sNx+Olx
182            OBNeta(I,bi,bj)=0.            OBNeta(I,bi,bj)=0.
# Line 79  C     == Local variables == Line 191  C     == Local variables ==
191  #ifdef ALLOW_ORLANSKI  #ifdef ALLOW_ORLANSKI
192          IF (useOrlanskiNorth.OR.useOrlanskiSouth.OR.          IF (useOrlanskiNorth.OR.useOrlanskiSouth.OR.
193       &      useOrlanskiEast.OR.useOrlanskiWest) THEN       &      useOrlanskiEast.OR.useOrlanskiWest) THEN
194    #ifdef ALLOW_DEBUG
195          IF (debugMode) CALL DEBUG_CALL('ORLANSKI_INIT',myThid)
196    #endif
197            CALL ORLANSKI_INIT(bi, bj, myThid)            CALL ORLANSKI_INIT(bi, bj, myThid)
198          ENDIF          ENDIF
199  #endif /* ALLOW_ORLANSKI */  #endif /* ALLOW_ORLANSKI */
# Line 102  c       CALL OBCS_READ_CHECKPOINT(prec, Line 217  c       CALL OBCS_READ_CHECKPOINT(prec,
217  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
218    
219  C--   Apply OBCS values to initial conditions for consistancy  C--   Apply OBCS values to initial conditions for consistancy
220    #ifdef ALLOW_DEBUG
221          IF (debugMode) CALL DEBUG_CALL('OBCS_CALC',myThid)
222    #endif
223        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
224         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
225          CALL OBCS_CALC( bi, bj, startTime, nIter0,          CALL OBCS_CALC( bi, bj, startTime, nIter0,
226       &        uVel, vVel, wVel, theta, salt, myThid )       &        uVel, vVel, wVel, theta, salt, 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  
227         ENDDO         ENDDO
228        ENDDO        ENDDO
229    
230          IF ( startTime .EQ. baseTime .AND.  nIter0 .EQ. 0
231         &     .AND. pickupSuff .EQ. ' ' ) THEN
232    #ifdef ALLOW_DEBUG
233           IF (debugMode)
234         &      CALL DEBUG_CALL('OBCS_APPLY_UV + OBCS_APPLY_TS',myThid)
235    #endif
236           DO bj = myByLo(myThid), myByHi(myThid)
237            DO bi = myBxLo(myThid), myBxHi(myThid)
238             DO K=1,Nr
239              CALL OBCS_APPLY_UV( bi, bj, k, uVel, vVel, myThid )
240              CALL OBCS_APPLY_TS( bi, bj, k, theta, salt, myThid )
241             ENDDO
242            ENDDO
243           ENDDO
244    #ifdef ALLOW_PTRACERS
245    C     repeat everything for passive tracers
246           IF ( usePTRACERS ) THEN
247    #ifdef ALLOW_DEBUG
248            IF (debugMode)
249         &       CALL DEBUG_CALL('OBCS_APPLY_PTRACER',myThid)
250    #endif
251            DO iTracer=1,PTRACERS_numInUse
252             DO bj = myByLo(myThid), myByHi(myThid)
253              DO bi = myBxLo(myThid), myBxHi(myThid)
254               DO K=1,Nr
255                CALL OBCS_APPLY_PTRACER(
256         I           bi, bj, K, iTracer,
257         U           ptracer(1-Olx,1-Oly,K,bi,bj,iTracer),
258         I           myThid )
259               ENDDO
260              ENDDO
261             ENDDO
262            ENDDO
263           ENDIF
264    #endif /* ALLOW_PTRACERS */
265           IF (useOBCSprescribe) THEN
266    C     After applying the boundary conditions exchange the 3D-fields.
267    C     This is only necessary of the boudnary values have been read
268    C     from a file.
269    #ifdef ALLOW_DEBUG
270            IF (debugMode)
271         &       CALL DEBUG_CALL('EXCHANGES in OBCS_INIT_VARIABLES',myThid)
272    #endif
273            CALL EXCH_UV_XYZ_RL(uVel,vVel,.TRUE.,myThid)
274            _EXCH_XYZ_R8( theta, myThid )
275            _EXCH_XYZ_R8( salt , myThid )
276           ENDIF
277    #ifdef ALLOW_PTRACERS
278           IF (usePTRACERS) THEN
279    C     Loop over passive tracers
280            DO iTracer=1,PTRACERS_numInUse
281             _EXCH_XYZ_R8( pTracer(1-Olx,1-Oly,1,1,1,iTracer),
282         &        myThid )
283            ENDDO
284           ENDIF
285    #endif /* ALLOW_PTRACERS */
286    C     endif start from rest
287          ENDIF
288    
289  #endif /* ALLOW_OBCS */  #endif /* ALLOW_OBCS */
290    
291    #ifdef ALLOW_DEBUG
292          IF (debugMode) CALL DEBUG_LEAVE('OBCS_INIT_VARIABLES',myThid)
293    #endif
294        RETURN        RETURN
295        END        END
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22