[MITgcm-support] OBCS for ocean, but periodic for sea ice?
dngoldberg at gmail.com
Wed Feb 14 14:03:44 EST 2018
Thanks martin!! Will give it a shot..
Sent from my iPhone. Apologies for brevity.
> On Feb 14, 2018, at 6:57 PM, Martin Losch <Martin.Losch at awi.de> wrote:
> Hi Dan,
> The model is periodic by default. That means it should be enough to make sure the the routines obcs_*ice* don’t do anything (there are 4: obcs_adjust_uvice.F obcs_apply_seaice.F obcs_apply_uvice.F obcs_seaice_sponge.F), i.e. put a EXIT as the first executable statement the code, or a CPP-FLAG or something like that.
>> On 13. Feb 2018, at 12:02, Daniel Goldberg <dngoldberg at gmail.com> wrote:
>> I am helping to set up an idealised continental shelf circulation model -- and OBCS is being used for east and west boundaries. For the sea ice model, though, we would like to impose periodic boundary conditions. Has that ever been done before? From what I can see, defining OBCS boundaries will then necessarily apply to both codes.
>> If this is not doable with the present code, a "hack" i thought of was to modify lines 88-120 of obcs_apply_seaice.F such that the appropriate values from the opposite boundary are used instead of OBEh/a/sn, etc, as appropriate (depending on the velocity direction). This would be very crude, but might get the job done.. would anyone be able to comment on this (vague) strategy?
>> Many thanks
>> Daniel Goldberg, PhD
>> Lecturer in Glaciology
>> School of Geosciences, University of Edinburgh
>> Geography Building, Drummond Street, Edinburgh EH8 9XP
>> em: Dan.Goldberg at ed.ac.uk
>> web: https://www.geos.ed.ac.uk/homes/dgoldber
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
More information about the MITgcm-support