/[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.7 by mlosch, Mon Oct 10 05:53:48 2005 UTC revision 1.14 by dimitri, Fri Oct 26 02:00:47 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  #ifdef ALLOW_PTRACERS.h  #if (defined (ALLOW_SEAICE) || defined (ALLOW_PTRACERS))
31    # include "PARAMS.h"
32    #endif
33    #ifdef ALLOW_PTRACERS
34  # include "PTRACERS_SIZE.h"  # include "PTRACERS_SIZE.h"
35  # include "PTRACERS.h"  # include "PTRACERS.h"
36  # include "OBCS_PTRACERS.h"  # include "OBCS_PTRACERS.h"
# Line 43  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
# Line 50  c     == local variables == Line 54  c     == local variables ==
54  #ifdef ALLOW_PTRACERS  #ifdef ALLOW_PTRACERS
55        integer iTracer, i,j,k        integer iTracer, i,j,k
56  #endif /* ALLOW_PTRACERS */  #endif /* ALLOW_PTRACERS */
57    #endif /* ALLOW_EXF */
58    
59  c     == end of interface ==  c     == end of interface ==
60    
# Line 76  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    #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           call exf_set_obcs_x ( OBNuice,OBNuice0,OBNuice1,OBNuicefile,'s'
99         I                     , fac, first, changed, count0, count1
100         I                     , mycurrenttime, mycurrentiter, mythid )
101           call exf_set_obcs_x ( OBNvice,OBNvice0,OBNvice1,OBNvicefile,'s'
102         I                     , fac, first, changed, count0, count1
103         I                     , mycurrenttime, mycurrentiter, mythid )
104          ENDIF
105    #endif /* ALLOW_SEAICE */
106  #ifdef ALLOW_PTRACERS  #ifdef ALLOW_PTRACERS
107        if ( usePTRACERS ) then        if ( usePTRACERS ) then
108         do itracer = 1, PTRACERS_numInUse         do itracer = 1, PTRACERS_numInUse
# Line 88  c     == end of interface == Line 115  c     == end of interface ==
115         enddo         enddo
116        endif        endif
117  #endif /* ALLOW_PTRACERS */  #endif /* ALLOW_PTRACERS */
118  #endif  #endif /* ALLOW_OBCS_NORTH */
119    
120  #ifdef ALLOW_OBCS_SOUTH  #ifdef ALLOW_OBCS_SOUTH
121        call exf_getffieldrec(        call exf_getffieldrec(
# Line 112  c     == end of interface == Line 139  c     == end of interface ==
139        call exf_set_obcs_xz(  OBSs, OBSs0, OBSs1, OBSsfile, 's'        call exf_set_obcs_xz(  OBSs, OBSs0, OBSs1, OBSsfile, 's'
140       I                     , fac, first, changed, count0, count1       I                     , fac, first, changed, count0, count1
141       I                     , mycurrenttime, mycurrentiter, mythid )       I                     , mycurrenttime, mycurrentiter, mythid )
142    #ifdef ALLOW_SEAICE
143          IF (useSEAICE) THEN
144           call exf_set_obcs_x (  OBSa, OBSa0, OBSa1, OBSafile, 's'
145         I                     , fac, first, changed, count0, count1
146         I                     , mycurrenttime, mycurrentiter, mythid )
147           call exf_set_obcs_x (  OBSh, OBSh0, OBSh1, OBShfile, 's'
148         I                     , fac, first, changed, count0, count1
149         I                     , mycurrenttime, mycurrentiter, mythid )
150           call exf_set_obcs_x (  OBSsl, OBSsl0, OBSsl1, OBSslfile, 's'
151         I                     , fac, first, changed, count0, count1
152         I                     , mycurrenttime, mycurrentiter, mythid )
153           call exf_set_obcs_x (  OBSsn, OBSsn0, OBSsn1, OBSsnfile, 's'
154         I                     , fac, first, changed, count0, count1
155         I                     , mycurrenttime, mycurrentiter, mythid )
156           call exf_set_obcs_x ( OBSuice,OBSuice0,OBSuice1,OBSuicefile,'s'
157         I                     , fac, first, changed, count0, count1
158         I                     , mycurrenttime, mycurrentiter, mythid )
159           call exf_set_obcs_x ( OBSvice,OBSvice0,OBSvice1,OBSvicefile,'s'
160         I                     , fac, first, changed, count0, count1
161         I                     , mycurrenttime, mycurrentiter, mythid )
162          ENDIF
163    #endif /* ALLOW_SEAICE */
164  #ifdef ALLOW_PTRACERS  #ifdef ALLOW_PTRACERS
165        if ( usePTRACERS ) then        if ( usePTRACERS ) then
166         do itracer = 1, PTRACERS_numInUse         do itracer = 1, PTRACERS_numInUse
# Line 125  c     == end of interface == Line 173  c     == end of interface ==
173         enddo         enddo
174        endif        endif
175  #endif /* ALLOW_PTRACERS */  #endif /* ALLOW_PTRACERS */
176  #endif  #endif /* ALLOW_OBCS_SOUTH */
177    
178  #ifdef ALLOW_OBCS_EAST  #ifdef ALLOW_OBCS_EAST
179        call exf_getffieldrec(        call exf_getffieldrec(
# Line 149  c     == end of interface == Line 197  c     == end of interface ==
197        call exf_set_obcs_yz(  OBEs, OBEs0, OBEs1, OBEsfile, 's'        call exf_set_obcs_yz(  OBEs, OBEs0, OBEs1, OBEsfile, 's'
198       I                     , fac, first, changed, count0, count1       I                     , fac, first, changed, count0, count1
199       I                     , mycurrenttime, mycurrentiter, mythid )       I                     , mycurrenttime, mycurrentiter, mythid )
200    #ifdef ALLOW_SEAICE
201          IF (useSEAICE) THEN
202           call exf_set_obcs_y (  OBEa, OBEa0, OBEa1, OBEafile, 's'
203         I                     , fac, first, changed, count0, count1
204         I                     , mycurrenttime, mycurrentiter, mythid )
205           call exf_set_obcs_y (  OBEh, OBEh0, OBEh1, OBEhfile, 's'
206         I                     , fac, first, changed, count0, count1
207         I                     , mycurrenttime, mycurrentiter, mythid )
208           call exf_set_obcs_y (  OBEsl, OBEsl0, OBEsl1, OBEslfile, 's'
209         I                     , fac, first, changed, count0, count1
210         I                     , mycurrenttime, mycurrentiter, mythid )
211           call exf_set_obcs_y (  OBEsn, OBEsn0, OBEsn1, OBEsnfile, 's'
212         I                     , fac, first, changed, count0, count1
213         I                     , mycurrenttime, mycurrentiter, mythid )
214           call exf_set_obcs_y ( OBEuice,OBEuice0,OBEuice1,OBEuicefile,'s'
215         I                     , fac, first, changed, count0, count1
216         I                     , mycurrenttime, mycurrentiter, mythid )
217           call exf_set_obcs_y ( OBEvice,OBEvice0,OBEvice1,OBEvicefile,'s'
218         I                     , fac, first, changed, count0, count1
219         I                     , mycurrenttime, mycurrentiter, mythid )
220          ENDIF
221    #endif /* ALLOW_SEAICE */
222  #ifdef ALLOW_PTRACERS  #ifdef ALLOW_PTRACERS
223        if ( usePTRACERS ) then        if ( usePTRACERS ) then
224         do itracer = 1, PTRACERS_numInUse         do itracer = 1, PTRACERS_numInUse
# Line 161  c     == end of interface == Line 231  c     == end of interface ==
231         enddo         enddo
232        endif        endif
233  #endif /* ALLOW_PTRACERS */  #endif /* ALLOW_PTRACERS */
234  #endif  #endif /* ALLOW_OBCS_EAST */
235    
236  #ifdef ALLOW_OBCS_WEST  #ifdef ALLOW_OBCS_WEST
237        call exf_getffieldrec(        call exf_getffieldrec(
# Line 185  c     == end of interface == Line 255  c     == end of interface ==
255        call exf_set_obcs_yz(  OBWs, OBWs0, OBWs1, OBWsfile, 's'        call exf_set_obcs_yz(  OBWs, OBWs0, OBWs1, OBWsfile, 's'
256       I                     , fac, first, changed, count0, count1       I                     , fac, first, changed, count0, count1
257       I                     , mycurrenttime, mycurrentiter, mythid )       I                     , mycurrenttime, mycurrentiter, mythid )
258    #ifdef ALLOW_SEAICE
259          IF (useSEAICE) THEN
260           call exf_set_obcs_y (  OBWa, OBWa0, OBWa1, OBWafile, 's'
261         I                     , fac, first, changed, count0, count1
262         I                     , mycurrenttime, mycurrentiter, mythid )
263           call exf_set_obcs_y (  OBWh, OBWh0, OBWh1, OBWhfile, 's'
264         I                     , fac, first, changed, count0, count1
265         I                     , mycurrenttime, mycurrentiter, mythid )
266           call exf_set_obcs_y (  OBWsl, OBWsl0, OBWsl1, OBWslfile, 's'
267         I                     , fac, first, changed, count0, count1
268         I                     , mycurrenttime, mycurrentiter, mythid )
269           call exf_set_obcs_y (  OBWsn, OBWsn0, OBWsn1, OBWsnfile, 's'
270         I                     , fac, first, changed, count0, count1
271         I                     , mycurrenttime, mycurrentiter, mythid )
272           call exf_set_obcs_y ( OBWuice,OBWuice0,OBWuice1,OBWuicefile,'s'
273         I                     , fac, first, changed, count0, count1
274         I                     , mycurrenttime, mycurrentiter, mythid )
275           call exf_set_obcs_y ( OBWvice,OBWvice0,OBWvice1,OBWvicefile,'s'
276         I                     , fac, first, changed, count0, count1
277         I                     , mycurrenttime, mycurrentiter, mythid )
278          ENDIF
279    #endif /* ALLOW_SEAICE */
280  #ifdef ALLOW_PTRACERS  #ifdef ALLOW_PTRACERS
281        if ( usePTRACERS ) then        if ( usePTRACERS ) then
282         do itracer = 1, PTRACERS_numInUse         do itracer = 1, PTRACERS_numInUse
# Line 197  c     == end of interface == Line 289  c     == end of interface ==
289         enddo         enddo
290        endif        endif
291  #endif /* ALLOW_PTRACERS */  #endif /* ALLOW_PTRACERS */
292  #endif  #endif /* ALLOW_OBCS_WEST */
293    
294  #ifdef ALLOW_OBCS_CONTROL  #ifdef ALLOW_OBCS_CONTROL
295  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.7  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22