<div dir="ltr">Hi Martin,<div><br></div><div>Sorry for the unclear words. I set the western and northwestern as Stevens BC. I just checked the model's last step output, the model crashed at the western boundary, the third grid point from the western boundary with very large zonal and meridional velocity. These abnormal values didn't occupy all the western boundary but within a very small region.</div><div><br></div><div>Thanks,</div><div>Lambert</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, May 13, 2021 at 8:53 AM Martin Losch <<a href="mailto:Martin.Losch@awi.de">Martin.Losch@awi.de</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi Lambert,<br>
<br>
this is way to unspecific to help. You need to find out where the model crashes (near the boundaries?) and then possibly understand if your prescribed files are not correct.<br>
<br>
Martin<br>
<br>
> On 13. May 2021, at 14:37, 赤羽蓝 <<a href="mailto:lanchiyu12@gmail.com" target="_blank">lanchiyu12@gmail.com</a>> wrote:<br>
> <br>
> Hi all,<br>
> <br>
> I am running a regional ocean model with the OBCS package. When I use the Stevens boundary conditions, the model always crashes with "stops due to EXTREME Pot.Temp". The model is OK if I turn off the Stevens BC. I have tried different inflows(outflow) at the Stevens boundary, which encounter the same crashing. The sponge boundary layer didn't help either. I also tried different relaxation time, no help.<br>
> <br>
> Attaching my data.obcs file.<br>
>   <br>
> Any ideas will be grateful.<br>
> <br>
> Lambert<br>
> <br>
> # ***************<br>
> # Open boundaries<br>
> # ***************<br>
>  &OBCS_PARM01<br>
> # This flag turns off checking and fixing problematic topography across<br>
> # open boundaries.<br>
>  OB_Jnorth = 360*-1,<br>
>  OB_Jsouth = 360*1,<br>
>  OB_Iwest = 528*1,<br>
>  OB_Ieast = 528*-1,<br>
>  useOBCSprescribe=.TRUE.,<br>
>  useStevensWest = .TRUE.,<br>
>  useStevensNorth = .TRUE.,<br>
> # useOBCSsponge=.TRUE.,<br>
>  useOBCSbalance=.FALSE.,<br>
> # OBCS_balanceFacE = 1.,<br>
> # OBCS_balanceFacW = 1.,<br>
> # OBCS_balanceFacN = 1.,<br>
> # OBCS_balanceFacS = 1.,<br>
> <br>
>  OBNsFile='sn.scs',<br>
>  OBNtFile='tn.scs',<br>
>  OBNvFile='vn.scs',<br>
>  OBNuFile='un.scs',<br>
> <br>
>  OBSsFile='ss.scs',<br>
>  OBStFile='ts.scs',<br>
>  OBSvFile='vs.scs',<br>
>  OBSuFile='us.scs',<br>
> <br>
>  OBEsFile='se.scs',<br>
>  OBEtFile='te.scs',<br>
>  OBEvFile='ve.scs',<br>
>  OBEuFile='ue.scs',<br>
> <br>
>  OBWsFile='sw.scs',<br>
>  OBWtFile='tw.scs',<br>
>  OBWvFile='vw.scs',<br>
>  OBWuFile='uw.scs',<br>
>  &<br>
> <br>
> # &OBCS_PARM03<br>
> # Urelaxobcsinner=864000,<br>
> # Urelaxobcsbound=864000,<br>
> # Vrelaxobcsinner=864000,<br>
> # Vrelaxobcsbound=864000,<br>
> # spongeThickness=5,<br>
> # &<br>
> <br>
>  &OBCS_PARM04<br>
>  TrelaxStevens = 86400.,<br>
>  SrelaxStevens = 86400.,<br>
>  &<br>
> _______________________________________________<br>
> MITgcm-support mailing list<br>
> <a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
<br>
_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</blockquote></div>