[MITgcm-support] [EXTERNAL] topography with OBCS

[Martin Losch]

Hi Andreas,

check the masks. 
You should be able to specify salinity everywhere and the masks should take care of this, i.e. if you specify salinity over land points, they will be masked.
In your case it looks like your mask and obcs file disagree and a land point is not masked but interpreted as zero salinity (but you knew that). I usually extrapolate my boundary field into land to avoid this.

Martin

> On 19. Jun 2020, at 05:12, Dimitris Menemenlis <dmenemenlis at gmail.com> wrote:
> 
> Hi Andreas, are you sure that you are extracting your boundary conditions at the correct location?
> You could be off by one grid cell, which could cause this problem.  Dimitris
> 
> 
>> On Jun 18, 2020, at 6:21 PM, Andreas Klocker <andreas.klocker at utas.edu.au> wrote:
>> 
>> Hi,
>> 
>> I am trying to set up a regional model configuration with OBCS. For some reason the model seems to introduce salinities where I actually have bathymetry (see attached figure; bathymetry left, salinities right). Salinities this low definitely do not exist in the files containing the bounday conditions. Is it necessary to set the open boundary exactly to where I have ocean, or should the bathymetry file take care of this by telling the model there should be nothing happening in case of zero water depth?
>> 
>> Here's my data.obcs file:
>> 
>> # Open-boundaries
>> &OBCS_PARM01
>> 
>> OB_Iwest  = 600*1,
>> OB_Ieast  = 600*-1,
>> OB_Jsouth = 2600*1,
>> OB_Jnorth = 2600*-1,
>> 
>> useOBCSsponge    = .TRUE.,
>> useOBCSbalance   = .TRUE.,
>> useOBCSprescribe = .TRUE.,
>> 
>> useOrlanskiWest  = .FALSE.,
>> useOrlanskiEast  = .FALSE.,
>> useOrlanskiSouth = .FALSE.,
>> useOrlanskiNorth = .FALSE.,
>> 
>> OBWuFile = 'OWU',
>> OBWvFile = 'OWV',
>> OBWtFile = 'OWT',
>> OBWsFile = 'OWS',
>> 
>> OBEuFile = 'OEU',
>> OBEvFile = 'OEV',
>> OBEtFile = 'OET',
>> OBEsFile = 'OES',
>> 
>> OBSuFile = 'OSU',
>> OBSvFile = 'OSV',
>> OBStFile = 'OST',
>> OBSsFile = 'OSS',
>> 
>> OBNuFile = 'ONU',
>> OBNvFile = 'ONV',
>> OBNtFile = 'ONT',
>> OBNsFile = 'ONS',
>> 
>> # OBCS_monitorFreq= 86400.,
>> # OBCS_monSelect = 1,
>> &
>> 
>> # Orlanski parameters
>> &OBCS_PARM02
>> # Cmax=0.45,
>> # cVelTimeScale=1000.,
>> &
>> 
>> &OBCS_PARM03
>> Urelaxobcsinner=864000.,
>> Urelaxobcsbound=86400.,
>> Vrelaxobcsinner=864000.,
>> Vrelaxobcsbound=86400.,
>> spongeThickness=40.,
>> &
>> 
>> Thanks in advance!
>> Andreas
>> 
>> 
>> University of Tasmania Electronic Communications Policy (December, 2014). 
>> This email is confidential, and is for the intended recipient only. Access, disclosure, copying, distribution, or reliance on any of it by anyone outside the intended recipient organisation is prohibited and may be a criminal offence. Please delete if obtained in error and email confirmation to the sender. The views expressed in this email are not necessarily the views of the University of Tasmania, unless clearly intended otherwise.
>> <obcs.png>_______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support