/[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.13 by edhill, Wed Oct 20 14:16:23 2004 UTC revision 1.14 by mlosch, Mon Oct 10 05:53:48 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  #ifdef ALLOW_DEBUG
44        IF (debugMode) CALL DEBUG_ENTER('OBCS_INIT_VARIABLES',myThid)        IF (debugMode) CALL DEBUG_ENTER('OBCS_INIT_VARIABLES',myThid)
# Line 124  C     == Local variables == Line 133  C     == Local variables ==
133           ENDDO           ENDDO
134          ENDDO          ENDDO
135    
136    #ifdef ALLOW_PTRACERS
137            DO iTracer=1,PTRACERS_numInUse
138             DO K=1,Nr
139              DO I=1-Olx,sNx+Olx
140    #ifdef ALLOW_OBCS_NORTH
141               OBNptr (I,K,bi,bj,iTracer)=0. _d 0
142    # ifdef ALLOW_OBCS_PRESCRIBE
143               OBNptr0(I,K,bi,bj,iTracer)=0. _d 0
144               OBNptr1(I,K,bi,bj,iTracer)=0. _d 0
145    # endif
146    #endif /* ALLOW_OBCS_NORTH */
147    
148    #ifdef ALLOW_OBCS_SOUTH
149               OBSptr (I,K,bi,bj,iTracer)=0. _d 0
150    # ifdef ALLOW_OBCS_PRESCRIBE
151               OBSptr0(I,K,bi,bj,iTracer)=0. _d 0
152               OBSptr1(I,K,bi,bj,iTracer)=0. _d 0
153    # endif
154    #endif /* ALLOW_OBCS_SOUTH */
155              ENDDO
156    
157              DO J=1-Oly,sNy+Oly
158    #ifdef ALLOW_OBCS_EAST
159               OBEptr (J,K,bi,bj,iTracer)=0. _d 0
160    # ifdef ALLOW_OBCS_PRESCRIBE
161               OBEptr0(J,K,bi,bj,iTracer)=0. _d 0
162               OBEptr1(J,K,bi,bj,iTracer)=0. _d 0
163    # endif
164    #endif /* ALLOW_OBCS_EAST */
165    
166    #ifdef ALLOW_OBCS_WEST
167               OBWptr (J,K,bi,bj,iTracer)=0. _d 0
168    # ifdef ALLOW_OBCS_PRESCRIBE
169               OBWptr0(J,K,bi,bj,iTracer)=0. _d 0
170               OBWptr1(J,K,bi,bj,iTracer)=0. _d 0
171    # endif
172    #endif /* ALLOW_OBCS_WEST */
173              ENDDO
174             ENDDO
175            ENDDO
176    #endif /* ALLOW_PTRACERS */
177    
178  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
179          DO I=1-Olx,sNx+Olx          DO I=1-Olx,sNx+Olx
180            OBNeta(I,bi,bj)=0.            OBNeta(I,bi,bj)=0.
# Line 183  C--   Apply OBCS values to initial condi Line 234  C--   Apply OBCS values to initial condi
234          DO K=1,Nr          DO K=1,Nr
235           CALL OBCS_APPLY_UV( bi, bj, k, uVel, vVel, myThid )           CALL OBCS_APPLY_UV( bi, bj, k, uVel, vVel, myThid )
236           CALL OBCS_APPLY_TS( bi, bj, k, theta, salt, myThid )           CALL OBCS_APPLY_TS( bi, bj, k, theta, salt, myThid )
237    #ifdef ALLOW_PTRACERS
238             DO iTracer=1,PTRACERS_numInUse
239              CALL OBCS_APPLY_PTRACER(
240         I         bi, bj, K, iTracer,
241         U         ptracer(1-Olx,1-Oly,K,bi,bj,iTracer),
242         I         myThid )
243             ENDDO
244    #endif /* ALLOW_PTRACERS */
245          ENDDO          ENDDO
246         ENDDO         ENDDO
247        ENDDO        ENDDO
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22