/[MITgcm]/MITgcm/pkg/exf/exf_set_obcs.F
ViewVC logotype

Diff of /MITgcm/pkg/exf/exf_set_obcs.F

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

revision 1.13 by dimitri, Fri May 9 10:00:40 2008 UTC revision 1.16 by jmc, Tue Sep 1 19:37:46 2009 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "EXF_OPTIONS.h"  #include "EXF_OPTIONS.h"
5    
6        subroutine exf_set_obcs_xz (  C--  File exf_set_obcs.F:
7       &       obcs_fld_xz, obcs_xz_0, obcs_xz_1  C--   Contents
8       I     , obcs_file, obcsmask  C--   o EXF_SET_OBCS_XZ
9       I     , fac, first, changed, useYearlyFields, obcs_period  C--   o EXF_SET_OBCS_YZ
10       I     , count0, count1, year0, year1  C--   o EXF_SET_OBCS_X
11       I     , mycurrenttime, mycurrentiter, mythid  C--   o EXF_SET_OBCS_Y
12       &                           )  
13    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
14    
15          SUBROUTINE EXF_SET_OBCS_XZ (
16         U       obcs_fld_xz, obcs_xz_0, obcs_xz_1,
17         I       obcs_file, obcsmask,
18         I       fac, first, changed, useYearlyFields, obcs_period,
19         I       count0, count1, year0, year1,
20         I       myTime, myIter, myThid )
21    
22  c     ==================================================================  c     ==================================================================
23  c     SUBROUTINE exf_set_obcs_xz  c     SUBROUTINE EXF_SET_OBCS_XZ
24  c     ==================================================================  c     ==================================================================
25  c  c
26  c     o set open boundary conditions  c     o set open boundary conditions
# Line 21  c     started: heimbach@mit.edu 01-May-2 Line 29  c     started: heimbach@mit.edu 01-May-2
29  c     mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002  c     mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002
30    
31  c     ==================================================================  c     ==================================================================
32  c     SUBROUTINE exf_set_obcs_xz  c     SUBROUTINE EXF_SET_OBCS_XZ
33  c     ==================================================================  c     ==================================================================
34    
35        implicit none        IMPLICIT NONE
36    
37  c     == global variables ==  c     == global variables ==
38    
# Line 36  c     == global variables == Line 44  c     == global variables ==
44    
45  c     == routine arguments ==  c     == routine arguments ==
46    
47        _RL obcs_fld_xz(1-olx:snx+olx,Nr,nsx,nsy)        _RL obcs_fld_xz(1-OLx:sNx+OLx,Nr,nSx,nSy)
48        _RL obcs_xz_0(1-olx:snx+olx,Nr,nsx,nsy)        _RL obcs_xz_0(1-OLx:sNx+OLx,Nr,nSx,nSy)
49        _RL obcs_xz_1(1-olx:snx+olx,Nr,nsx,nsy)        _RL obcs_xz_1(1-OLx:sNx+OLx,Nr,nSx,nSy)
50    
51        character*(128) obcs_file        CHARACTER*(128) obcs_file
52        character*1 obcsmask        CHARACTER*1 obcsmask
53        logical first, changed        LOGICAL first, changed
54        logical useYearlyFields        LOGICAL useYearlyFields
55        _RL     obcs_period        _RL     obcs_period
56        integer count0, count1, year0, year1        INTEGER count0, count1, year0, year1
57        _RL     fac        _RL     fac
58        _RL     mycurrenttime        _RL     myTime
59        integer mycurrentiter        INTEGER myIter
60        integer mythid        INTEGER myThid
61    
62  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
63    
64  c     == local variables ==  c     == local variables ==
65    
66        integer bi, bj        CHARACTER*(128) obcs_file0, obcs_file1
67        integer i, k        INTEGER bi, bj
68          INTEGER i, k
       integer il  
       character*(128) obcs_file0, obcs_file1  
   
 c     == external ==  
   
       integer  ilnblnk  
       external ilnblnk  
69    
70  c     == end of interface ==  c     == end of interface ==
71    
72        if ( obcs_file .NE. ' ' ) then        IF ( obcs_file .NE. ' ' ) THEN
73    
74           if ( first ) then           IF ( first ) THEN
75    
76              call exf_GetYearlyFieldName(              CALL exf_GetYearlyFieldName(
77       I         useYearlyFields, twoDigitYear, obcs_period, year0,       I         useYearlyFields, twoDigitYear, obcs_period, year0,
78       I         obcs_file,       I         obcs_file,
79       O         obcs_file0,       O         obcs_file0,
80       I         mycurrenttime, mycurrentiter, mythid )       I         myTime, myIter, myThid )
81    
82              call mdsreadfieldxz( obcs_file0, exf_iprec, exf_yftype, Nr              CALL READ_REC_XZ_RL( obcs_file0, exf_iprec_obcs, Nr,
83       &           , obcs_xz_1, count0, mythid       &                           obcs_xz_1, count0, myIter, myThid )
84       &           )           ENDIF
          endif  
85    
86           if (( first ) .or. ( changed )) then           IF (( first ) .OR. ( changed )) THEN
87              call exf_swapffields_xz( obcs_xz_0, obcs_xz_1, mythid )              CALL exf_swapffields_xz( obcs_xz_0, obcs_xz_1, myThid )
88    
89              call exf_GetYearlyFieldName(              CALL exf_GetYearlyFieldName(
90       I         useYearlyFields, twoDigitYear, obcs_period, year1,       I         useYearlyFields, twoDigitYear, obcs_period, year1,
91       I         obcs_file,       I         obcs_file,
92       O         obcs_file1,       O         obcs_file1,
93       I         mycurrenttime, mycurrentiter, mythid )       I         myTime, myIter, myThid )
94    
95              call mdsreadfieldxz( obcs_file1, exf_iprec, exf_yftype, Nr              CALL READ_REC_XZ_RL( obcs_file1, exf_iprec_obcs, Nr,
96       &           , obcs_xz_1, count1, mythid       &                           obcs_xz_1, count1, myIter, myThid )
97       &           )           ENDIF
98           endif  
99             DO bj = myByLo(myThid),myByHi(myThid)
100           do bj = mybylo(mythid),mybyhi(mythid)              DO bi = myBxLo(myThid),myBxHi(myThid)
101              do bi = mybxlo(mythid),mybxhi(mythid)                 DO k = 1,Nr
102                 do k = 1,Nr                    DO i = 1,sNx
103                    do i = 1,snx                       obcs_fld_xz(i,k,bi,bj) =
                      obcs_fld_xz(i,k,bi,bj) =  
104       &                    fac * obcs_xz_0(i,k,bi,bj) +       &                    fac * obcs_xz_0(i,k,bi,bj) +
105       &                    (exf_one - fac) * obcs_xz_1(i,k,bi,bj)       &                    (exf_one - fac) * obcs_xz_1(i,k,bi,bj)
106                    enddo                    ENDDO
107                 enddo                 ENDDO
108              enddo              ENDDO
109           enddo           ENDDO
110    
111        endif        ENDIF
112    
113  #endif /* ALLOW_OBCS */  #endif /* ALLOW_OBCS */
114    
115        end        RETURN
116          END
117    
118        subroutine exf_set_obcs_yz (  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
119       &       obcs_fld_yz, obcs_yz_0, obcs_yz_1  
120       I     , obcs_file, obcsmask        SUBROUTINE EXF_SET_OBCS_YZ (
121       I     , fac, first, changed, useYearlyFields, obcs_period       U       obcs_fld_yz, obcs_yz_0, obcs_yz_1,
122       I     , count0, count1, year0, year1       I       obcs_file, obcsmask,
123       I     , mycurrenttime, mycurrentiter, mythid       I       fac, first, changed, useYearlyFields, obcs_period,
124       &                           )       I       count0, count1, year0, year1,
125         I       myTime, myIter, myThid)
126    
127  c     ==================================================================  c     ==================================================================
128  c     SUBROUTINE exf_set_obcs_yz  c     SUBROUTINE EXF_SET_OBCS_YZ
129  c     ==================================================================  c     ==================================================================
130  c  c
131  c     o set open boundary conditions  c     o set open boundary conditions
# Line 132  c Line 133  c
133  c     started: heimbach@mit.edu 01-May-2001  c     started: heimbach@mit.edu 01-May-2001
134    
135  c     ==================================================================  c     ==================================================================
136  c     SUBROUTINE exf_set_obcs_yz  c     SUBROUTINE EXF_SET_OBCS_YZ
137  c     ==================================================================  c     ==================================================================
138    
139        implicit none        IMPLICIT NONE
140    
141  c     == global variables ==  c     == global variables ==
142    
# Line 147  c     == global variables == Line 148  c     == global variables ==
148    
149  c     == routine arguments ==  c     == routine arguments ==
150    
151        _RL obcs_fld_yz(1-oly:sny+oly,Nr,nsx,nsy)        _RL obcs_fld_yz(1-OLy:sNy+OLy,Nr,nSx,nSy)
152        _RL obcs_yz_0(1-oly:sny+oly,Nr,nsx,nsy)        _RL obcs_yz_0(1-OLy:sNy+OLy,Nr,nSx,nSy)
153        _RL obcs_yz_1(1-oly:sny+oly,Nr,nsx,nsy)        _RL obcs_yz_1(1-OLy:sNy+OLy,Nr,nSx,nSy)
154        character*(MAX_LEN_FNAM) obcs_file        CHARACTER*(MAX_LEN_FNAM) obcs_file
155        character*1 obcsmask        CHARACTER*1 obcsmask
156        logical first, changed        LOGICAL first, changed
157        logical useYearlyFields        LOGICAL useYearlyFields
158        _RL     obcs_period        _RL     obcs_period
159        integer count0, count1, year0, year1        INTEGER count0, count1, year0, year1
160        _RL     fac        _RL     fac
161        _RL     mycurrenttime        _RL     myTime
162        integer mycurrentiter        INTEGER myIter
163        integer mythid        INTEGER myThid
164    
165  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
166    
167  c     == local variables ==  c     == local variables ==
168    
169        integer bi, bj        CHARACTER*(128) obcs_file0, obcs_file1
170        integer j, k        INTEGER bi, bj
171        integer il        INTEGER j, k
       character*(128) obcs_file0, obcs_file1  
   
 c     == external ==  
   
       integer  ilnblnk  
       external ilnblnk  
172    
173  c     == end of interface ==  c     == end of interface ==
174    
175        if ( obcs_file .NE. ' ' ) then        IF ( obcs_file .NE. ' ' ) THEN
176    
177           if ( first ) then           IF ( first ) THEN
178    
179              call exf_GetYearlyFieldName(              CALL exf_GetYearlyFieldName(
180       I         useYearlyFields, twoDigitYear, obcs_period, year0,       I         useYearlyFields, twoDigitYear, obcs_period, year0,
181       I         obcs_file,       I         obcs_file,
182       O         obcs_file0,       O         obcs_file0,
183       I         mycurrenttime, mycurrentiter, mythid )       I         myTime, myIter, myThid )
184    
185              call mdsreadfieldyz( obcs_file0, exf_iprec, exf_yftype, Nr              CALL READ_REC_YZ_RL( obcs_file0, exf_iprec_obcs, Nr,
186       &           , obcs_yz_1, count0, mythid       &                           obcs_yz_1, count0, myIter, myThid )
187       &           )           ENDIF
          endif  
188    
189           if (( first ) .or. ( changed )) then           IF (( first ) .OR. ( changed )) THEN
190              call exf_swapffields_yz( obcs_yz_0, obcs_yz_1, mythid )              CALL exf_swapffields_yz( obcs_yz_0, obcs_yz_1, myThid )
191    
192              call exf_GetYearlyFieldName(              CALL exf_GetYearlyFieldName(
193       I         useYearlyFields, twoDigitYear, obcs_period, year1,       I         useYearlyFields, twoDigitYear, obcs_period, year1,
194       I         obcs_file,       I         obcs_file,
195       O         obcs_file1,       O         obcs_file1,
196       I         mycurrenttime, mycurrentiter, mythid )       I         myTime, myIter, myThid )
197    
198              call mdsreadfieldyz( obcs_file1, exf_iprec, exf_yftype, Nr              CALL READ_REC_YZ_RL( obcs_file1, exf_iprec_obcs, Nr,
199       &           , obcs_yz_1, count1, mythid       &                           obcs_yz_1, count1, myIter, myThid )
200       &           )           ENDIF
201           endif  
202             DO bj = myByLo(myThid),myByHi(myThid)
203           do bj = mybylo(mythid),mybyhi(mythid)              DO bi = myBxLo(myThid),myBxHi(myThid)
204              do bi = mybxlo(mythid),mybxhi(mythid)                 DO k = 1,Nr
205                 do k = 1,Nr                    DO j = 1,sNy
206                    do j = 1,sny                       obcs_fld_yz(j,k,bi,bj) =
                      obcs_fld_yz(j,k,bi,bj) =  
207       &                    fac             *obcs_yz_0(j,k,bi,bj) +       &                    fac             *obcs_yz_0(j,k,bi,bj) +
208       &                    (exf_one - fac) *obcs_yz_1(j,k,bi,bj)       &                    (exf_one - fac) *obcs_yz_1(j,k,bi,bj)
209                    enddo                    ENDDO
210                 enddo                 ENDDO
211              enddo              ENDDO
212           enddo           ENDDO
213    
214        endif        ENDIF
215    
216  #endif /* ALLOW_OBCS */  #endif /* ALLOW_OBCS */
217    
218        end        RETURN
219          END
220    
221        subroutine exf_set_obcs_x (  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
222       &       obcs_fld_x, obcs_x_0, obcs_x_1  
223       I     , obcs_file, obcsmask        SUBROUTINE EXF_SET_OBCS_X (
224       I     , fac, first, changed, useYearlyFields, obcs_period       U       obcs_fld_x, obcs_x_0, obcs_x_1,
225       I     , count0, count1, year0, year1       I       obcs_file, obcsmask,
226       I     , mycurrenttime, mycurrentiter, mythid       I       fac, first, changed, useYearlyFields, obcs_period,
227       &                           )       I       count0, count1, year0, year1,
228         I       myTime, myIter, myThid )
229    
230  c     ==================================================================  c     ==================================================================
231  c     SUBROUTINE exf_set_obcs_x  c     SUBROUTINE EXF_SET_OBCS_X
232  c     ==================================================================  c     ==================================================================
233  c  c
234  c     o set open boundary conditions  c     o set open boundary conditions
235  c       same as exf_set_obcs_xz but for NR=1  c       same as EXF_SET_OBCS_XZ but for Nr=1
236  c  c
237  c     ==================================================================  c     ==================================================================
238  c     SUBROUTINE exf_set_obcs_x  c     SUBROUTINE EXF_SET_OBCS_X
239  c     ==================================================================  c     ==================================================================
240    
241        implicit none        IMPLICIT NONE
242    
243  c     == global variables ==  c     == global variables ==
244    
# Line 255  c     == global variables == Line 250  c     == global variables ==
250    
251  c     == routine arguments ==  c     == routine arguments ==
252    
253        _RL obcs_fld_x(1-olx:snx+olx,nsx,nsy)        _RL obcs_fld_x(1-OLx:sNx+OLx,nSx,nSy)
254        _RL obcs_x_0(1-olx:snx+olx,nsx,nsy)        _RL obcs_x_0(1-OLx:sNx+OLx,nSx,nSy)
255        _RL obcs_x_1(1-olx:snx+olx,nsx,nsy)        _RL obcs_x_1(1-OLx:sNx+OLx,nSx,nSy)
256    
257        character*(128) obcs_file        CHARACTER*(128) obcs_file
258        character*1 obcsmask        CHARACTER*1 obcsmask
259        logical first, changed        LOGICAL first, changed
260        logical useYearlyFields        LOGICAL useYearlyFields
261        _RL     obcs_period        _RL     obcs_period
262        integer count0, count1, year0, year1        INTEGER count0, count1, year0, year1
263        _RL     fac        _RL     fac
264        _RL     mycurrenttime        _RL     myTime
265        integer mycurrentiter        INTEGER myIter
266        integer mythid        INTEGER myThid
267    
268  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
269    
270  c     == local variables ==  c     == local variables ==
271    
272        integer bi, bj, i        CHARACTER*(128) obcs_file0, obcs_file1
273          INTEGER bi, bj, i
       integer il  
       character*(128) obcs_file0, obcs_file1  
   
 c     == external ==  
   
       integer  ilnblnk  
       external ilnblnk  
274    
275  c     == end of interface ==  c     == end of interface ==
276    
277        if ( obcs_file .NE. ' ' ) then        IF ( obcs_file .NE. ' ' ) THEN
278    
279           if ( first ) then           IF ( first ) THEN
280    
281              call exf_GetYearlyFieldName(              CALL exf_GetYearlyFieldName(
282       I         useYearlyFields, twoDigitYear, obcs_period, year0,       I         useYearlyFields, twoDigitYear, obcs_period, year0,
283       I         obcs_file,       I         obcs_file,
284       O         obcs_file0,       O         obcs_file0,
285       I         mycurrenttime, mycurrentiter, mythid )       I         myTime, myIter, myThid )
286    
287              call mdsreadfieldxz( obcs_file0, exf_iprec, exf_yftype, 1              CALL READ_REC_XZ_RL( obcs_file0, exf_iprec_obcs, 1,
288       &           , obcs_x_1, count0, mythid       &                           obcs_x_1, count0, myIter, myThid )
289       &           )           ENDIF
          endif  
290    
291           if (( first ) .or. ( changed )) then           IF (( first ) .OR. ( changed )) THEN
292              call exf_swapffields_x( obcs_x_0, obcs_x_1, mythid )              CALL exf_swapffields_x( obcs_x_0, obcs_x_1, myThid )
293    
294              call exf_GetYearlyFieldName(              CALL exf_GetYearlyFieldName(
295       I         useYearlyFields, twoDigitYear, obcs_period, year1,       I         useYearlyFields, twoDigitYear, obcs_period, year1,
296       I         obcs_file,       I         obcs_file,
297       O         obcs_file1,       O         obcs_file1,
298       I         mycurrenttime, mycurrentiter, mythid )       I         myTime, myIter, myThid )
299    
300              call mdsreadfieldxz( obcs_file1, exf_iprec, exf_yftype, 1              CALL READ_REC_XZ_RL( obcs_file1, exf_iprec_obcs, 1,
301       &           , obcs_x_1, count1, mythid       &                           obcs_x_1, count1, myIter, myThid )
302       &           )           ENDIF
303           endif  
304             DO bj = myByLo(myThid),myByHi(myThid)
305           do bj = mybylo(mythid),mybyhi(mythid)              DO bi = myBxLo(myThid),myBxHi(myThid)
306              do bi = mybxlo(mythid),mybxhi(mythid)                 DO i = 1,sNx
307                 do i = 1,snx                    obcs_fld_x(i,bi,bj) =
                   obcs_fld_x(i,bi,bj) =  
308       &                 fac * obcs_x_0(i,bi,bj) +       &                 fac * obcs_x_0(i,bi,bj) +
309       &                 (exf_one - fac) * obcs_x_1(i,bi,bj)       &                 (exf_one - fac) * obcs_x_1(i,bi,bj)
310                 enddo                 ENDDO
311              enddo              ENDDO
312           enddo           ENDDO
313    
314        endif        ENDIF
315    
316  #endif /* ALLOW_OBCS */  #endif /* ALLOW_OBCS */
317    
318        end        RETURN
319          END
320    
321    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
322    
323        subroutine exf_set_obcs_y (        SUBROUTINE EXF_SET_OBCS_Y (
324       &       obcs_fld_y, obcs_y_0, obcs_y_1       U       obcs_fld_y, obcs_y_0, obcs_y_1,
325       I     , obcs_file, obcsmask       I       obcs_file, obcsmask,
326       I     , fac, first, changed, useYearlyFields, obcs_period       I       fac, first, changed, useYearlyFields, obcs_period,
327       I     , count0, count1, year0, year1       I       count0, count1, year0, year1,
328       I     , mycurrenttime, mycurrentiter, mythid       I       myTime, myIter, myThid )
      &                           )  
329    
330  c     ==================================================================  c     ==================================================================
331  c     SUBROUTINE exf_set_obcs_y  c     SUBROUTINE EXF_SET_OBCS_Y
332  c     ==================================================================  c     ==================================================================
333  c  c
334  c     o set open boundary conditions  c     o set open boundary conditions
335  c       same as exf_set_obcs_yz but for NR=1  c       same as EXF_SET_OBCS_YZ but for Nr=1
336  c  c
337  c     ==================================================================  c     ==================================================================
338  c     SUBROUTINE exf_set_obcs_y  c     SUBROUTINE EXF_SET_OBCS_Y
339  c     ==================================================================  c     ==================================================================
340    
341        implicit none        IMPLICIT NONE
342    
343  c     == global variables ==  c     == global variables ==
344    
# Line 362  c     == global variables == Line 350  c     == global variables ==
350    
351  c     == routine arguments ==  c     == routine arguments ==
352    
353        _RL obcs_fld_y(1-oly:sny+oly,nsx,nsy)        _RL obcs_fld_y(1-OLy:sNy+OLy,nSx,nSy)
354        _RL obcs_y_0(1-oly:sny+oly,nsx,nsy)        _RL obcs_y_0(1-OLy:sNy+OLy,nSx,nSy)
355        _RL obcs_y_1(1-oly:sny+oly,nsx,nsy)        _RL obcs_y_1(1-OLy:sNy+OLy,nSx,nSy)
356        character*(MAX_LEN_FNAM) obcs_file        CHARACTER*(MAX_LEN_FNAM) obcs_file
357        character*1 obcsmask        CHARACTER*1 obcsmask
358        logical first, changed        LOGICAL first, changed
359        logical useYearlyFields        LOGICAL useYearlyFields
360        _RL     obcs_period        _RL     obcs_period
361        integer count0, count1, year0, year1        INTEGER count0, count1, year0, year1
362        _RL     fac        _RL     fac
363        _RL     mycurrenttime        _RL     myTime
364        integer mycurrentiter        INTEGER myIter
365        integer mythid        INTEGER myThid
366    
367  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
368    
369  c     == local variables ==  c     == local variables ==
370    
371        integer bi, bj, j        CHARACTER*(128) obcs_file0, obcs_file1
372          INTEGER bi, bj, j
       integer il  
       character*(128) obcs_file0, obcs_file1  
   
 c     == external ==  
   
       integer  ilnblnk  
       external ilnblnk  
373    
374  c     == end of interface ==  c     == end of interface ==
375    
376        if ( obcs_file .NE. ' ' ) then        IF ( obcs_file .NE. ' ' ) THEN
377    
378           if ( first ) then           IF ( first ) THEN
379    
380              call exf_GetYearlyFieldName(              CALL exf_GetYearlyFieldName(
381       I         useYearlyFields, twoDigitYear, obcs_period, year0,       I         useYearlyFields, twoDigitYear, obcs_period, year0,
382       I         obcs_file,       I         obcs_file,
383       O         obcs_file0,       O         obcs_file0,
384       I         mycurrenttime, mycurrentiter, mythid )       I         myTime, myIter, myThid )
385    
386              call mdsreadfieldyz( obcs_file0, exf_iprec, exf_yftype, 1              CALL READ_REC_YZ_RL( obcs_file0, exf_iprec_obcs, 1,
387       &           , obcs_y_1, count0, mythid       &                           obcs_y_1, count0, myIter, myThid )
388       &           )           ENDIF
          endif  
389    
390           if (( first ) .or. ( changed )) then           IF (( first ) .OR. ( changed )) THEN
391              call exf_swapffields_y( obcs_y_0, obcs_y_1, mythid )              CALL exf_swapffields_y( obcs_y_0, obcs_y_1, myThid )
392    
393              call exf_GetYearlyFieldName(              CALL exf_GetYearlyFieldName(
394       I         useYearlyFields, twoDigitYear, obcs_period, year1,       I         useYearlyFields, twoDigitYear, obcs_period, year1,
395       I         obcs_file,       I         obcs_file,
396       O         obcs_file1,       O         obcs_file1,
397       I         mycurrenttime, mycurrentiter, mythid )       I         myTime, myIter, myThid )
398    
399              call mdsreadfieldyz( obcs_file1, exf_iprec, exf_yftype, 1              CALL READ_REC_YZ_RL( obcs_file1, exf_iprec_obcs, 1,
400       &           , obcs_y_1, count1, mythid       &                           obcs_y_1, count1, myIter, myThid )
401       &           )           ENDIF
402           endif  
403             DO bj = myByLo(myThid),myByHi(myThid)
404           do bj = mybylo(mythid),mybyhi(mythid)              DO bi = myBxLo(myThid),myBxHi(myThid)
405              do bi = mybxlo(mythid),mybxhi(mythid)                 DO j = 1,sNy
406                 do j = 1,sny                    obcs_fld_y(j,bi,bj) =
                   obcs_fld_y(j,bi,bj) =  
407       &                 fac             *obcs_y_0(j,bi,bj) +       &                 fac             *obcs_y_0(j,bi,bj) +
408       &                 (exf_one - fac) *obcs_y_1(j,bi,bj)       &                 (exf_one - fac) *obcs_y_1(j,bi,bj)
409                 enddo                 ENDDO
410              enddo              ENDDO
411           enddo           ENDDO
412    
413        endif        ENDIF
414    
415  #endif /* ALLOW_OBCS */  #endif /* ALLOW_OBCS */
416    
417        end        RETURN
418          END
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22