[MITgcm-support] OBCS for deep sewage water source
Martin Losch
mlosch at awi-bremerhaven.de
Mon Dec 20 11:45:47 EST 2004
Andrea,
sorry for the delay, I didn't get around to doing anything last week.
And now I only have a short answer: I don't think that you can use OBCS
in the way you want to. You have a source in the middle of the domain
at the bottom, and I don't see why or how OBCS should work in this
case. If I were you, I would prescribe a source at one grid point in
ptracers_forcing.F (in tracer units/time), --- and if you want to
simulate outflow through a sewer at the same point, some velocity
source in external_forcing.F (S/R external_forcing_u), may some
restoring--- and NOT use OBCS for the sewer problem. OBCS was not made
for this type of situation, I guess. Any other opinions, anyone?
Martin
On Dec 13, 2004, at 4:14 PM, cimatori at libero.it wrote:
> Hi Martin, some time ago we wrote about this case: we are trying to
> model the diffusion of a deep
> source of wastewater in a semi-enclosed basin (gulf of Trieste,
> northern Adriatic) using obcs
> prescribe to attain the effect of lighter water, ptracers (for
> modeling the polluttant) and exf. We
> use a western obcs at the open side of the gulf (we tried both
> Orlanski and sponge (nudging) layer)
> and tried to use an eastern boundary for the sewage source. This
> source is a diffuser offshore on
> the seabed, so we need to prescribe the boundary values only at a
> particular level (15) leaving the
> rest of the column without any constraint.
> This works only if we use a closed domain (without any kind of western
> open boundary): we see a
> plume of lighter water wich follows the wind-driven currents.
> When we open the western boundary we obtain irrealistic velocity
> fields on the column above the
> diffuser.
> We corrected the routines (see attachments) making some improvements
> but the main problem
> remains.
> The problem lies in the use of another open boundary together with the
> one in the middle of the
> domain (we are going to try using northern or southern obcs instead of
> eastern for the source:
> maybe it could be a C-grid problem?).
> We will check again the source code (expecially our modifications),
> but, of course, any further
> suggestion is welcome!
> Thanks!
> Andrea
>
> PS: the passive tracer with obcs seems to give sensible results, at
> least qualitatively, without
> using the deep source of light water (using only ptracer forcing). See
> image attached.
>
>
>
> ____________________________________________________________
> Libero ADSL: 3 mesi gratis e navighi a 1.2 Mega, senza costi di
> attivazione.
> Abbonati subito su http://www.libero.it
> <U_hor_sect_log(ptrac).png><obcs_apply_uv.F><obcs_apply_w.F>___________
> ____________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://dev.mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list