/[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.3 by heimbach, Mon Oct 11 16:41:01 2004 UTC revision 1.12 by dimitri, Thu Oct 11 10:30:34 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
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 49  c     == end of interface == Line 63  c     == end of interface ==
63        call exf_getffieldrec(        call exf_getffieldrec(
64       I                       obcsNstartdate, obcsNperiod       I                       obcsNstartdate, obcsNperiod
65       I                     , obcsNstartdate1, obcsNstartdate2       I                     , obcsNstartdate1, obcsNstartdate2
66       I                     , false.       I                     , .false.
67       O                     , fac, first, changed,       O                     , fac, first, changed
68       O                     , count0, count1, year0, year1       O                     , count0, count1, year0, year1
69       I                     , mycurrenttime, mycurrentiter, mythid       I                     , mycurrenttime, mycurrentiter, mythid
70       &                     )       &                     )
# Line 67  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          ENDIF
93    #endif /* ALLOW_SEAICE */
94    #ifdef ALLOW_PTRACERS
95          if ( usePTRACERS ) then
96           do itracer = 1, PTRACERS_numInUse
97            call exf_set_obcs_xz(  OBNptr (1-Olx,1,1,1,iTracer)
98         I                       , OBNptr0(1-Olx,1,1,1,iTracer)
99         I                       , OBNptr1(1-Olx,1,1,1,iTracer)
100         I                       , OBNptrFile(iTracer), 's'
101         I                       , fac, first, changed, count0, count1
102         I                       , mycurrenttime, mycurrentiter, mythid )
103           enddo
104          endif
105    #endif /* ALLOW_PTRACERS */
106    #endif /* ALLOW_OBCS_NORTH */
107    
108  #ifdef ALLOW_OBCS_SOUTH  #ifdef ALLOW_OBCS_SOUTH
109        call exf_getffieldrec(        call exf_getffieldrec(
110       I                       obcsSstartdate, obcsSperiod       I                       obcsSstartdate, obcsSperiod
111       I                     , obcsSstartdate1, obcsSstartdate2       I                     , obcsSstartdate1, obcsSstartdate2
112       I                     , false.       I                     , .false.
113       O                     , fac, first, changed, count0, count1       O                     , fac, first, changed
114         O                     , count0, count1, year0, year1
115       I                     , mycurrenttime, mycurrentiter, mythid       I                     , mycurrenttime, mycurrentiter, mythid
116       &                     )       &                     )
117    
# Line 90  c     == end of interface == Line 127  c     == end of interface ==
127        call exf_set_obcs_xz(  OBSs, OBSs0, OBSs1, OBSsfile, 's'        call exf_set_obcs_xz(  OBSs, OBSs0, OBSs1, OBSsfile, 's'
128       I                     , fac, first, changed, count0, count1       I                     , fac, first, changed, count0, count1
129       I                     , mycurrenttime, mycurrentiter, mythid )       I                     , mycurrenttime, mycurrentiter, mythid )
130    #ifdef ALLOW_SEAICE
131  #endif        IF (useSEAICE) THEN
132           call exf_set_obcs_x (  OBSa, OBSa0, OBSa1, OBSafile, 's'
133         I                     , fac, first, changed, count0, count1
134         I                     , mycurrenttime, mycurrentiter, mythid )
135           call exf_set_obcs_x (  OBSh, OBSh0, OBSh1, OBShfile, 's'
136         I                     , fac, first, changed, count0, count1
137         I                     , mycurrenttime, mycurrentiter, mythid )
138          ENDIF
139    #endif /* ALLOW_SEAICE */
140    #ifdef ALLOW_PTRACERS
141          if ( usePTRACERS ) then
142           do itracer = 1, PTRACERS_numInUse
143            call exf_set_obcs_xz(  OBSptr (1-Olx,1,1,1,iTracer)
144         I                       , OBSptr0(1-Olx,1,1,1,iTracer)
145         I                       , OBSptr1(1-Olx,1,1,1,iTracer)
146         I                       , OBSptrFile(iTracer), 's'
147         I                       , fac, first, changed, count0, count1
148         I                       , mycurrenttime, mycurrentiter, mythid )
149           enddo
150          endif
151    #endif /* ALLOW_PTRACERS */
152    #endif /* ALLOW_OBCS_SOUTH */
153    
154  #ifdef ALLOW_OBCS_EAST  #ifdef ALLOW_OBCS_EAST
155        call exf_getffieldrec(        call exf_getffieldrec(
156       I                       obcsEstartdate, obcsEperiod       I                       obcsEstartdate, obcsEperiod
157       I                     , obcsEstartdate1, obcsEstartdate2       I                     , obcsEstartdate1, obcsEstartdate2
158       I                     , false.       I                     , .false.
159       O                     , fac, first, changed, count0, count1       O                     , fac, first, changed
160         O                     , count0, count1, year0, year1
161       I                     , mycurrenttime, mycurrentiter, mythid       I                     , mycurrenttime, mycurrentiter, mythid
162       &                     )       &                     )
163    
# Line 114  c     == end of interface == Line 173  c     == end of interface ==
173        call exf_set_obcs_yz(  OBEs, OBEs0, OBEs1, OBEsfile, 's'        call exf_set_obcs_yz(  OBEs, OBEs0, OBEs1, OBEsfile, 's'
174       I                     , fac, first, changed, count0, count1       I                     , fac, first, changed, count0, count1
175       I                     , mycurrenttime, mycurrentiter, mythid )       I                     , mycurrenttime, mycurrentiter, mythid )
176    #ifdef ALLOW_SEAICE
177  #endif        IF (useSEAICE) THEN
178           call exf_set_obcs_y (  OBEa, OBEa0, OBEa1, OBEafile, 's'
179         I                     , fac, first, changed, count0, count1
180         I                     , mycurrenttime, mycurrentiter, mythid )
181           call exf_set_obcs_y (  OBEh, OBEh0, OBEh1, OBEhfile, 's'
182         I                     , fac, first, changed, count0, count1
183         I                     , mycurrenttime, mycurrentiter, mythid )
184          ENDIF
185    #endif /* ALLOW_SEAICE */
186    #ifdef ALLOW_PTRACERS
187          if ( usePTRACERS ) then
188           do itracer = 1, PTRACERS_numInUse
189            call exf_set_obcs_yz(  OBEptr (1-Oly,1,1,1,iTracer)
190         I                       , OBEptr0(1-Oly,1,1,1,iTracer)
191         I                       , OBEptr1(1-Oly,1,1,1,iTracer)
192         I                       , OBEptrFile(iTracer), 's'
193         I                       , fac, first, changed, count0, count1
194         I                       , mycurrenttime, mycurrentiter, mythid )
195           enddo
196          endif
197    #endif /* ALLOW_PTRACERS */
198    #endif /* ALLOW_OBCS_EAST */
199    
200  #ifdef ALLOW_OBCS_WEST  #ifdef ALLOW_OBCS_WEST
201        call exf_getffieldrec(        call exf_getffieldrec(
202       I                       obcsWstartdate, obcsWperiod       I                       obcsWstartdate, obcsWperiod
203       I                     , obcsWstartdate1, obcsWstartdate2       I                     , obcsWstartdate1, obcsWstartdate2
204       I                     , false.       I                     , .false.
205       O                     , fac, first, changed, count0, count1       O                     , fac, first, changed
206         O                     , count0, count1, year0, year1
207       I                     , mycurrenttime, mycurrentiter, mythid       I                     , mycurrenttime, mycurrentiter, mythid
208       &                     )       &                     )
209    
# Line 138  c     == end of interface == Line 219  c     == end of interface ==
219        call exf_set_obcs_yz(  OBWs, OBWs0, OBWs1, OBWsfile, 's'        call exf_set_obcs_yz(  OBWs, OBWs0, OBWs1, OBWsfile, 's'
220       I                     , fac, first, changed, count0, count1       I                     , fac, first, changed, count0, count1
221       I                     , mycurrenttime, mycurrentiter, mythid )       I                     , mycurrenttime, mycurrentiter, mythid )
222    #ifdef ALLOW_SEAICE
223  #endif        IF (useSEAICE) THEN
224           call exf_set_obcs_y (  OBWa, OBWa0, OBWa1, OBWafile, 's'
225         I                     , fac, first, changed, count0, count1
226         I                     , mycurrenttime, mycurrentiter, mythid )
227           call exf_set_obcs_y (  OBWh, OBWh0, OBWh1, OBWhfile, 's'
228         I                     , fac, first, changed, count0, count1
229         I                     , mycurrenttime, mycurrentiter, mythid )
230          ENDIF
231    #endif /* ALLOW_SEAICE */
232    #ifdef ALLOW_PTRACERS
233          if ( usePTRACERS ) then
234           do itracer = 1, PTRACERS_numInUse
235            call exf_set_obcs_yz(  OBWptr (1-Oly,1,1,1,iTracer)
236         I                       , OBWptr0(1-Oly,1,1,1,iTracer)
237         I                       , OBWptr1(1-Oly,1,1,1,iTracer)
238         I                       , OBWptrFile(iTracer), 's'
239         I                       , fac, first, changed, count0, count1
240         I                       , mycurrenttime, mycurrentiter, mythid )
241           enddo
242          endif
243    #endif /* ALLOW_PTRACERS */
244    #endif /* ALLOW_OBCS_WEST */
245    
246  #ifdef ALLOW_OBCS_CONTROL  #ifdef ALLOW_OBCS_CONTROL
247  cgg   WARNING: Assuming North Open Boundary exists and has same  cgg   WARNING: Assuming North Open Boundary exists and has same
# Line 170  cgg    calendar information as other bou Line 272  cgg    calendar information as other bou
272    
273  #endif /* ALLOW_OBCS */  #endif /* ALLOW_OBCS */
274    
275        end        RETURN
276          END
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22