14 |
c o set open boundary conditions |
c o set open boundary conditions |
15 |
c |
c |
16 |
c started: heimbach@mit.edu 01-May-2001 |
c started: heimbach@mit.edu 01-May-2001 |
17 |
|
c mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002 |
18 |
|
|
19 |
c ================================================================== |
c ================================================================== |
20 |
c SUBROUTINE exf_set_obcs_xz |
c SUBROUTINE exf_set_obcs_xz |
54 |
|
|
55 |
c == end of interface == |
c == end of interface == |
56 |
|
|
57 |
if ( first ) then |
if ( obcs_file .NE. ' ' ) then |
|
if ( obcs_file .NE. ' ' ) |
|
|
& call mdsreadfieldxz( obcs_file, exf_iprec, exf_yftype, Nr |
|
|
& , obcs_xz_1, count0, mythid |
|
|
& ) |
|
|
endif |
|
58 |
|
|
59 |
if (( first ) .or. ( changed )) then |
if ( first ) then |
60 |
call exf_swapffields_xz( obcs_xz_0, obcs_xz_1, mythid ) |
call mdsreadfieldxz( obcs_file, exf_iprec, exf_yftype, Nr |
61 |
|
& , obcs_xz_1, count0, mythid |
62 |
|
& ) |
63 |
|
endif |
64 |
|
|
65 |
|
if (( first ) .or. ( changed )) then |
66 |
|
call exf_swapffields_xz( obcs_xz_0, obcs_xz_1, mythid ) |
67 |
|
|
68 |
|
call mdsreadfieldxz( obcs_file, exf_iprec, exf_yftype, Nr |
69 |
|
& , obcs_xz_1, count1, mythid |
70 |
|
& ) |
71 |
|
endif |
72 |
|
|
73 |
|
do bj = mybylo(mythid),mybyhi(mythid) |
74 |
|
do bi = mybxlo(mythid),mybxhi(mythid) |
75 |
|
do k = 1,Nr |
76 |
|
do i = 1,snx |
77 |
|
obcs_fld_xz(i,k,bi,bj) = |
78 |
|
& fac * obcs_xz_0(i,k,bi,bj) + |
79 |
|
& (exf_one - fac) * obcs_xz_1(i,k,bi,bj) |
80 |
|
enddo |
81 |
|
enddo |
82 |
|
enddo |
83 |
|
enddo |
84 |
|
|
|
if ( obcs_file .NE. ' ' ) |
|
|
& call mdsreadfieldxz( obcs_file, exf_iprec, exf_yftype, Nr |
|
|
& , obcs_xz_1, count1, mythid |
|
|
& ) |
|
85 |
endif |
endif |
86 |
|
|
|
do bj = mybylo(mythid),mybyhi(mythid) |
|
|
do bi = mybxlo(mythid),mybxhi(mythid) |
|
|
do k = 1,Nr |
|
|
do i = 1,snx |
|
|
obcs_fld_xz(i,k,bi,bj) = |
|
|
& fac *obcs_xz_0(i,k,bi,bj) + |
|
|
& (exf_one - fac) *obcs_xz_1(i,k,bi,bj) |
|
|
enddo |
|
|
enddo |
|
|
enddo |
|
|
enddo |
|
|
|
|
87 |
#endif |
#endif |
88 |
|
|
89 |
end |
end |