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

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

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

revision 1.6 by adcroft, Tue Oct 19 17:52:03 2004 UTC revision 1.13 by dimitri, Tue Oct 23 08:25:28 2007 UTC
# Line 25  c     == global variables == Line 25  c     == global variables ==
25  #include "GRID.h"  #include "GRID.h"
26  #include "OBCS.h"  #include "OBCS.h"
27  #ifdef ALLOW_EXF  #ifdef ALLOW_EXF
28  # include "exf_param.h"  # include "EXF_PARAM.h"
29  #endif  #endif
30    #if (defined (ALLOW_SEAICE) || defined (ALLOW_PTRACERS))
31    # include "PARAMS.h"
32    #endif
33    #ifdef ALLOW_PTRACERS
34    # include "PTRACERS_SIZE.h"
35    # include "PTRACERS.h"
36    # include "OBCS_PTRACERS.h"
37    #endif /* ALLOW_PTRACERS */
38    
39  c     == routine arguments ==  c     == routine arguments ==
40    
# Line 38  c     == routine arguments == Line 46  c     == routine arguments ==
46    
47  c     == local variables ==  c     == local variables ==
48    
49    #ifdef ALLOW_EXF
50        logical first, changed        logical first, changed
51        integer count0, count1        integer count0, count1
52        integer year0, year1        integer year0, year1
53        _RL     fac        _RL     fac
54    #ifdef ALLOW_PTRACERS
55          integer iTracer, i,j,k
56    #endif /* ALLOW_PTRACERS */
57    #endif /* ALLOW_EXF */
58    
59  c     == end of interface ==  c     == end of interface ==
60    
# Line 68  c     == end of interface == Line 81  c     == end of interface ==
81        call exf_set_obcs_xz(  OBNs, OBNs0, OBNs1, OBNsfile, 's'        call exf_set_obcs_xz(  OBNs, OBNs0, OBNs1, OBNsfile, 's'
82       I                     , fac, first, changed, count0, count1       I                     , fac, first, changed, count0, count1
83       I                     , mycurrenttime, mycurrentiter, mythid )       I                     , mycurrenttime, mycurrentiter, mythid )
84  #endif  #ifdef ALLOW_SEAICE
85          IF (useSEAICE) THEN
86           call exf_set_obcs_x (  OBNa, OBNa0, OBNa1, OBNafile, 's'
87         I                     , fac, first, changed, count0, count1
88         I                     , mycurrenttime, mycurrentiter, mythid )
89           call exf_set_obcs_x (  OBNh, OBNh0, OBNh1, OBNhfile, 's'
90         I                     , fac, first, changed, count0, count1
91         I                     , mycurrenttime, mycurrentiter, mythid )
92           call exf_set_obcs_x (  OBNsl, OBNsl0, OBNsl1, OBNslfile, 's'
93         I                     , fac, first, changed, count0, count1
94         I                     , mycurrenttime, mycurrentiter, mythid )
95           call exf_set_obcs_x (  OBNsn, OBNsn0, OBNsn1, OBNsnfile, 's'
96         I                     , fac, first, changed, count0, count1
97         I                     , mycurrenttime, mycurrentiter, mythid )
98          ENDIF
99    #endif /* ALLOW_SEAICE */
100    #ifdef ALLOW_PTRACERS
101          if ( usePTRACERS ) then
102           do itracer = 1, PTRACERS_numInUse
103            call exf_set_obcs_xz(  OBNptr (1-Olx,1,1,1,iTracer)
104         I                       , OBNptr0(1-Olx,1,1,1,iTracer)
105         I                       , OBNptr1(1-Olx,1,1,1,iTracer)
106         I                       , OBNptrFile(iTracer), 's'
107         I                       , fac, first, changed, count0, count1
108         I                       , mycurrenttime, mycurrentiter, mythid )
109           enddo
110          endif
111    #endif /* ALLOW_PTRACERS */
112    #endif /* ALLOW_OBCS_NORTH */
113    
114  #ifdef ALLOW_OBCS_SOUTH  #ifdef ALLOW_OBCS_SOUTH
115        call exf_getffieldrec(        call exf_getffieldrec(
# Line 92  c     == end of interface == Line 133  c     == end of interface ==
133        call exf_set_obcs_xz(  OBSs, OBSs0, OBSs1, OBSsfile, 's'        call exf_set_obcs_xz(  OBSs, OBSs0, OBSs1, OBSsfile, 's'
134       I                     , fac, first, changed, count0, count1       I                     , fac, first, changed, count0, count1
135       I                     , mycurrenttime, mycurrentiter, mythid )       I                     , mycurrenttime, mycurrentiter, mythid )
136    #ifdef ALLOW_SEAICE
137  #endif        IF (useSEAICE) THEN
138           call exf_set_obcs_x (  OBSa, OBSa0, OBSa1, OBSafile, 's'
139         I                     , fac, first, changed, count0, count1
140         I                     , mycurrenttime, mycurrentiter, mythid )
141           call exf_set_obcs_x (  OBSh, OBSh0, OBSh1, OBShfile, 's'
142         I                     , fac, first, changed, count0, count1
143         I                     , mycurrenttime, mycurrentiter, mythid )
144           call exf_set_obcs_x (  OBSsl, OBSsl0, OBSsl1, OBSslfile, 's'
145         I                     , fac, first, changed, count0, count1
146         I                     , mycurrenttime, mycurrentiter, mythid )
147           call exf_set_obcs_x (  OBSsn, OBSsn0, OBSsn1, OBSsnfile, 's'
148         I                     , fac, first, changed, count0, count1
149         I                     , mycurrenttime, mycurrentiter, mythid )
150          ENDIF
151    #endif /* ALLOW_SEAICE */
152    #ifdef ALLOW_PTRACERS
153          if ( usePTRACERS ) then
154           do itracer = 1, PTRACERS_numInUse
155            call exf_set_obcs_xz(  OBSptr (1-Olx,1,1,1,iTracer)
156         I                       , OBSptr0(1-Olx,1,1,1,iTracer)
157         I                       , OBSptr1(1-Olx,1,1,1,iTracer)
158         I                       , OBSptrFile(iTracer), 's'
159         I                       , fac, first, changed, count0, count1
160         I                       , mycurrenttime, mycurrentiter, mythid )
161           enddo
162          endif
163    #endif /* ALLOW_PTRACERS */
164    #endif /* ALLOW_OBCS_SOUTH */
165    
166  #ifdef ALLOW_OBCS_EAST  #ifdef ALLOW_OBCS_EAST
167        call exf_getffieldrec(        call exf_getffieldrec(
# Line 117  c     == end of interface == Line 185  c     == end of interface ==
185        call exf_set_obcs_yz(  OBEs, OBEs0, OBEs1, OBEsfile, 's'        call exf_set_obcs_yz(  OBEs, OBEs0, OBEs1, OBEsfile, 's'
186       I                     , fac, first, changed, count0, count1       I                     , fac, first, changed, count0, count1
187       I                     , mycurrenttime, mycurrentiter, mythid )       I                     , mycurrenttime, mycurrentiter, mythid )
188    #ifdef ALLOW_SEAICE
189  #endif        IF (useSEAICE) THEN
190           call exf_set_obcs_y (  OBEa, OBEa0, OBEa1, OBEafile, 's'
191         I                     , fac, first, changed, count0, count1
192         I                     , mycurrenttime, mycurrentiter, mythid )
193           call exf_set_obcs_y (  OBEh, OBEh0, OBEh1, OBEhfile, 's'
194         I                     , fac, first, changed, count0, count1
195         I                     , mycurrenttime, mycurrentiter, mythid )
196           call exf_set_obcs_y (  OBEsl, OBEsl0, OBEsl1, OBEslfile, 's'
197         I                     , fac, first, changed, count0, count1
198         I                     , mycurrenttime, mycurrentiter, mythid )
199           call exf_set_obcs_y (  OBEsn, OBEsn0, OBEsn1, OBEsnfile, 's'
200         I                     , fac, first, changed, count0, count1
201         I                     , mycurrenttime, mycurrentiter, mythid )
202          ENDIF
203    #endif /* ALLOW_SEAICE */
204    #ifdef ALLOW_PTRACERS
205          if ( usePTRACERS ) then
206           do itracer = 1, PTRACERS_numInUse
207            call exf_set_obcs_yz(  OBEptr (1-Oly,1,1,1,iTracer)
208         I                       , OBEptr0(1-Oly,1,1,1,iTracer)
209         I                       , OBEptr1(1-Oly,1,1,1,iTracer)
210         I                       , OBEptrFile(iTracer), 's'
211         I                       , fac, first, changed, count0, count1
212         I                       , mycurrenttime, mycurrentiter, mythid )
213           enddo
214          endif
215    #endif /* ALLOW_PTRACERS */
216    #endif /* ALLOW_OBCS_EAST */
217    
218  #ifdef ALLOW_OBCS_WEST  #ifdef ALLOW_OBCS_WEST
219        call exf_getffieldrec(        call exf_getffieldrec(
# Line 142  c     == end of interface == Line 237  c     == end of interface ==
237        call exf_set_obcs_yz(  OBWs, OBWs0, OBWs1, OBWsfile, 's'        call exf_set_obcs_yz(  OBWs, OBWs0, OBWs1, OBWsfile, 's'
238       I                     , fac, first, changed, count0, count1       I                     , fac, first, changed, count0, count1
239       I                     , mycurrenttime, mycurrentiter, mythid )       I                     , mycurrenttime, mycurrentiter, mythid )
240    #ifdef ALLOW_SEAICE
241  #endif        IF (useSEAICE) THEN
242           call exf_set_obcs_y (  OBWa, OBWa0, OBWa1, OBWafile, 's'
243         I                     , fac, first, changed, count0, count1
244         I                     , mycurrenttime, mycurrentiter, mythid )
245           call exf_set_obcs_y (  OBWh, OBWh0, OBWh1, OBWhfile, 's'
246         I                     , fac, first, changed, count0, count1
247         I                     , mycurrenttime, mycurrentiter, mythid )
248           call exf_set_obcs_y (  OBWsl, OBWsl0, OBWsl1, OBWslfile, 's'
249         I                     , fac, first, changed, count0, count1
250         I                     , mycurrenttime, mycurrentiter, mythid )
251           call exf_set_obcs_y (  OBWsn, OBWsn0, OBWsn1, OBWsnfile, 's'
252         I                     , fac, first, changed, count0, count1
253         I                     , mycurrenttime, mycurrentiter, mythid )
254          ENDIF
255    #endif /* ALLOW_SEAICE */
256    #ifdef ALLOW_PTRACERS
257          if ( usePTRACERS ) then
258           do itracer = 1, PTRACERS_numInUse
259            call exf_set_obcs_yz(  OBWptr (1-Oly,1,1,1,iTracer)
260         I                       , OBWptr0(1-Oly,1,1,1,iTracer)
261         I                       , OBWptr1(1-Oly,1,1,1,iTracer)
262         I                       , OBWptrFile(iTracer), 's'
263         I                       , fac, first, changed, count0, count1
264         I                       , mycurrenttime, mycurrentiter, mythid )
265           enddo
266          endif
267    #endif /* ALLOW_PTRACERS */
268    #endif /* ALLOW_OBCS_WEST */
269    
270  #ifdef ALLOW_OBCS_CONTROL  #ifdef ALLOW_OBCS_CONTROL
271  cgg   WARNING: Assuming North Open Boundary exists and has same  cgg   WARNING: Assuming North Open Boundary exists and has same
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22