2 SUBROUTINE ana_fsobc (ng, tile, model)
21 integer,
intent(in) :: ng, tile, model
25 character (len=*),
parameter :: MyFile = &
30 CALL ana_fsobc_tile (ng, tile, model, &
31 & lbi, ubi, lbj, ubj, &
32 & imins, imaxs, jmins, jmaxs)
48 SUBROUTINE ana_fsobc_tile (ng, tile, model, &
49 & LBi, UBi, LBj, UBj, &
50 & IminS, ImaxS, JminS, JmaxS)
61 integer,
intent(in) :: ng, tile, model
62 integer,
intent(in) :: LBi, UBi, LBj, UBj
63 integer,
intent(in) :: IminS, ImaxS, JminS, JmaxS
69 real(r8) :: cff, fac, omega, phase, val
71#include "set_bounds.h"
79 &
domain(ng)%Northern_Edge(tile))
THEN
80 cff=-1.0_r8*sin(2.0_r8*
pi*
time(ng)/(12.0_r8*3600.0_r8))
87 omega=2.0_r8*
pi/(12.42_r8*3600.0_r8)
89 &
domain(ng)%Western_Edge(tile))
THEN
91 val=fac*exp(-
grid(ng)%f(istr-1,j)*
grid(ng)%yp(istr-1,j)/ &
92 & sqrt(
g*
grid(ng)%h(istr-1,j)))
98 &
domain(ng)%Eastern_Edge(tile))
THEN
100 cff=1.0_r8/sqrt(
g*
grid(ng)%h(istr-1,j))
101 val=fac*exp(-
grid(ng)%f(istr-1,j)*
grid(ng)%yp(iend,j)*cff)
102 boundary(ng)%zeta_east(j)=val*cos(omega*
grid(ng)%xp(iend,j)* &
103 & cff-omega*
time(ng))
106#elif defined ESTUARY_TEST
108 &
domain(ng)%Western_Edge(tile))
THEN
109 cff=1.0_r8*sin(2.0_r8*
pi*
time(ng)/(12.0_r8*3600.0_r8))
114#elif defined SED_TEST1
116 &
domain(ng)%Western_Edge(tile))
THEN
119 boundary(ng)%zeta_west(j)=9.0e-06_r8*fac
124 &
domain(ng)%Eastern_Edge(tile))
THEN
127 boundary(ng)%zeta_east(j)=9.0e-06_r8*real(iend+1,r8)*fac
130#elif defined SHOREFACE
132 &
domain(ng)%Western_Edge(tile))
THEN
139#elif defined TEST_CHAN
141 &
domain(ng)%Western_Edge(tile))
THEN
149 &
domain(ng)%Eastern_Edge(tile))
THEN
150 cff=-0.4040_r8*min(
time(ng)/150000.0_r8,1.0_r8)
157 &
domain(ng)%Western_Edge(tile))
THEN
159 omega=2.0_r8*
pi*
time(ng)/(12.42_r8*3600.0_r8)
160 val=0.53_r8+(0.53_r8-0.48_r8)/real(iend+1,r8)
161 phase=(277.0_r8+(277.0_r8-240.0_r8)/real(iend+1,r8))*
deg2rad
163 boundary(ng)%zeta_west(j)=fac*val*cos(omega-phase)
168 &
domain(ng)%Eastern_Edge(tile))
THEN
170 omega=2.0_r8*
pi*
time(ng)/(12.42_r8*3600.0_r8)
171 val=0.53_r8+(0.53_r8-0.48_r8)
172 phase=(277.0_r8+(277.0_r8-240.0_r8))*
deg2rad
174 boundary(ng)%zeta_east(j)=fac*val*cos(omega-phase)
179 &
domain(ng)%Eastern_Edge(tile))
THEN
186 &
domain(ng)%Western_Edge(tile))
THEN
193 &
domain(ng)%Southern_Edge(tile))
THEN
200 &
domain(ng)%Northern_Edge(tile))
THEN
subroutine ana_fsobc_tile(ng, tile, model, lbi, ubi, lbj, ubj, imins, imaxs, jmins, jmaxs)