Issues with Simulating Salinity Plume in SWASH: Model Crashes When Approaching the Coast

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
maddayout
Posts: 1
Joined: Mon Mar 10, 2025 2:57 am
Location: utfsm

Issues with Simulating Salinity Plume in SWASH: Model Crashes When Approaching the Coast

#1 Unread post by maddayout »

Hello ROMS community,
I am trying to simulate the interaction between the ocean and a river to analyze the salinity plume on SWASH MODEL. Currently, I am working only with a river discharge flowing into the sea (no waves or tides involved). However, when the plume starts approaching the coast, the model crashes. Honestly, I don't fully understand why this is happening.
My domain is 1200x999 meters with cells of 3x3 meters each. The simulation time is 12 hours, but it crashes at 6 hours. According to the SWASH manual, when dealing with salinity and sediment analysis, it is suggested to use TVD as the numerical scheme. I applied it, but the model still crashes.
I will attach the .sws, .bot, and .sal files in case anyone can help me.
Could it be the bathymetry? In my opinion, it shouldn't cause any issues.
The .sws file contains the following:

$*************HEADING***********************************
$
PROJ 'tester' '01'
$
$***********MODEL INPUT*********************************
$
MODE NONSTationary TWOD
COORDinates CARTesian
$
CGRID 0. 0. 0. 999 1200 333 400
VERT 2
$
INPGRID BOTTOM 0. 0. 0. 333 400 3.0 3.0
READINP BOTTOM 1. 'bathymetry.bot' 1 0 FREE
$
INITial ZERO
$
INPTRANS SALINITY 0. 0. 0. 333 400 3.0 3.0 EXCEPTION -999
READTRANS SALINITY 1. 'salinity.sal' 1 0 FREE
$
BOU SEGM XY 999 321 999 303 BTYPE DISCH SMOO 10 MIN CON FOUR -1
BOU SIDE W CCW BTYPE RADIATION
BOU SIDE N CCW BTYPE RADIATION
BOU SIDE S CCW BTYPE RADIATION
$
VISC H CON 0.001
FRIC MANN 0.019
NONHYDROSTATIC STAndard 1
$
DISCRET UPW MOM
DISCRET UPW WMOM
DISCRET CORRDEP
DISCRET MINMOD
$
TIMEI METH EXPL 0.3 0.6
$
$************ OUTPUT REQUESTS *************************
$
QUANTITY VELK
QUANTITY SALK
$
GROUP 'OUTPUT' SUBGRID 1 334 1 401
BLOCK 'OUTPUT' NOHEAD 'salinidad_limiter.mat' LAY-OUT 4 SALK OUTPUT 000000.0 60 SEC
BLOCK 'OUTPUT' NOHEAD 'velocidad_limiter.mat' LAY-OUT 4 VELK OUTPUT 000000.0 60 SEC
$
TEST 1,0
COMPUTE 000000.0 0.09 SEC 043200.0
STOP

I hope someone can help me.
Best regards,
Leonora Pereira
+56976494317
leonora.pereira@usm.cl
Civil Engineering Student, UTFSM

jcwarner
Posts: 1215
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: Issues with Simulating Salinity Plume in SWASH: Model Crashes When Approaching the Coast

#2 Unread post by jcwarner »

i think this is the wrong site for your question.
Does SWASH have their own user help site?
-j

Post Reply