/[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.12 by dimitri, Fri May 9 09:20:30 2008 UTC revision 1.13 by dimitri, Fri May 9 10:00:40 2008 UTC
# Line 39  c     == routine arguments == Line 39  c     == routine arguments ==
39        _RL obcs_fld_xz(1-olx:snx+olx,Nr,nsx,nsy)        _RL obcs_fld_xz(1-olx:snx+olx,Nr,nsx,nsy)
40        _RL obcs_xz_0(1-olx:snx+olx,Nr,nsx,nsy)        _RL obcs_xz_0(1-olx:snx+olx,Nr,nsx,nsy)
41        _RL obcs_xz_1(1-olx:snx+olx,Nr,nsx,nsy)        _RL obcs_xz_1(1-olx:snx+olx,Nr,nsx,nsy)
42        character*(MAX_LEN_FNAM) obcs_file  
43          character*(128) obcs_file
44        character*1 obcsmask        character*1 obcsmask
45        logical first, changed        logical first, changed
46        logical useYearlyFields        logical useYearlyFields
# Line 56  c     == local variables == Line 57  c     == local variables ==
57    
58        integer bi, bj        integer bi, bj
59        integer i, k        integer i, k
60    
61        integer il        integer il
62        character*(MAX_LEN_FNAM) obcs_file0, obcs_file1        character*(128) obcs_file0, obcs_file1
63    
64  c     == external ==  c     == external ==
65    
# Line 100  c     == end of interface == Line 102  c     == end of interface ==
102                 do k = 1,Nr                 do k = 1,Nr
103                    do i = 1,snx                    do i = 1,snx
104                       obcs_fld_xz(i,k,bi,bj) =                       obcs_fld_xz(i,k,bi,bj) =
105       &                    fac             * obcs_xz_0(i,k,bi,bj) +       &                    fac * obcs_xz_0(i,k,bi,bj) +
106       &                    (exf_one - fac) * obcs_xz_1(i,k,bi,bj)       &                    (exf_one - fac) * obcs_xz_1(i,k,bi,bj)
107                    enddo                    enddo
108                 enddo                 enddo
# Line 128  c Line 130  c
130  c     o set open boundary conditions  c     o set open boundary conditions
131  c  c
132  c     started: heimbach@mit.edu 01-May-2001  c     started: heimbach@mit.edu 01-May-2001
 c     mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002  
133    
134  c     ==================================================================  c     ==================================================================
135  c     SUBROUTINE exf_set_obcs_yz  c     SUBROUTINE exf_set_obcs_yz
# Line 167  c     == local variables == Line 168  c     == local variables ==
168        integer bi, bj        integer bi, bj
169        integer j, k        integer j, k
170        integer il        integer il
171        character*(MAX_LEN_FNAM) obcs_file0, obcs_file1        character*(128) obcs_file0, obcs_file1
172    
173  c     == external ==  c     == external ==
174    
# Line 181  c     == end of interface == Line 182  c     == end of interface ==
182           if ( first ) then           if ( first ) then
183    
184              call exf_GetYearlyFieldName(              call exf_GetYearlyFieldName(
185       I         useYearlyFields, twoDigitYear, obcs_period, year0,       I         useYearlyFields, twoDigitYear, obcs_period, year0,
186       I         obcs_file,       I         obcs_file,
187       O         obcs_file0,       O         obcs_file0,
188       I         mycurrenttime, mycurrentiter, mythid )       I         mycurrenttime, mycurrentiter, mythid )
# Line 210  c     == end of interface == Line 211  c     == end of interface ==
211                 do k = 1,Nr                 do k = 1,Nr
212                    do j = 1,sny                    do j = 1,sny
213                       obcs_fld_yz(j,k,bi,bj) =                       obcs_fld_yz(j,k,bi,bj) =
214       &                    fac             * obcs_yz_0(j,k,bi,bj) +       &                    fac             *obcs_yz_0(j,k,bi,bj) +
215       &                    (exf_one - fac) * obcs_yz_1(j,k,bi,bj)       &                    (exf_one - fac) *obcs_yz_1(j,k,bi,bj)
216                    enddo                    enddo
217                 enddo                 enddo
218              enddo              enddo
# Line 257  c     == routine arguments == Line 258  c     == routine arguments ==
258        _RL obcs_fld_x(1-olx:snx+olx,nsx,nsy)        _RL obcs_fld_x(1-olx:snx+olx,nsx,nsy)
259        _RL obcs_x_0(1-olx:snx+olx,nsx,nsy)        _RL obcs_x_0(1-olx:snx+olx,nsx,nsy)
260        _RL obcs_x_1(1-olx:snx+olx,nsx,nsy)        _RL obcs_x_1(1-olx:snx+olx,nsx,nsy)
261        character*(MAX_LEN_FNAM) obcs_file  
262          character*(128) obcs_file
263        character*1 obcsmask        character*1 obcsmask
264        logical first, changed        logical first, changed
265        logical useYearlyFields        logical useYearlyFields
# Line 273  c     == routine arguments == Line 275  c     == routine arguments ==
275  c     == local variables ==  c     == local variables ==
276    
277        integer bi, bj, i        integer bi, bj, i
278    
279        integer il        integer il
280        character*(MAX_LEN_FNAM) obcs_file0, obcs_file1        character*(128) obcs_file0, obcs_file1
281    
282  c     == external ==  c     == external ==
283    
# Line 316  c     == end of interface == Line 319  c     == end of interface ==
319              do bi = mybxlo(mythid),mybxhi(mythid)              do bi = mybxlo(mythid),mybxhi(mythid)
320                 do i = 1,snx                 do i = 1,snx
321                    obcs_fld_x(i,bi,bj) =                    obcs_fld_x(i,bi,bj) =
322       &                 fac             * obcs_x_0(i,bi,bj) +       &                 fac * obcs_x_0(i,bi,bj) +
323       &                 (exf_one - fac) * obcs_x_1(i,bi,bj)       &                 (exf_one - fac) * obcs_x_1(i,bi,bj)
324                 enddo                 enddo
325              enddo              enddo
# Line 378  c     == routine arguments == Line 381  c     == routine arguments ==
381  c     == local variables ==  c     == local variables ==
382    
383        integer bi, bj, j        integer bi, bj, j
384    
385        integer il        integer il
386        character*(MAX_LEN_FNAM) obcs_file0, obcs_file1        character*(128) obcs_file0, obcs_file1
387    
388  c     == external ==  c     == external ==
389    
# Line 421  c     == end of interface == Line 425  c     == end of interface ==
425              do bi = mybxlo(mythid),mybxhi(mythid)              do bi = mybxlo(mythid),mybxhi(mythid)
426                 do j = 1,sny                 do j = 1,sny
427                    obcs_fld_y(j,bi,bj) =                    obcs_fld_y(j,bi,bj) =
428       &                 fac             * obcs_y_0(j,bi,bj) +       &                 fac             *obcs_y_0(j,bi,bj) +
429       &                 (exf_one - fac) * obcs_y_1(j,bi,bj)       &                 (exf_one - fac) *obcs_y_1(j,bi,bj)
430                 enddo                 enddo
431              enddo              enddo
432           enddo           enddo
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22