Error in ROMS model

[yernisrinivasnekkali]

Dear all,
I am trying to run the COAWST model by coupling ROMS and WRF. However, ROMS encounters an error during the initial integration.

Based on my research, this issue is often caused by NaN values in the initial data or bathymetry.
1. I generated the initial and boundary conditions using HYCOM data ('https://tds.hycom.org/thredds/dodsC/GLBy0.08/expt_93.0') with MATLAB scripts provided in roms_clim. I also cross-checked the data, and there are no NaN values.
2. For bathymetry, I used ETOPO 2022 data and generated a grid using MATLAB GridBuilder.

For reference, I have attached the restart file, ocean.in file, and output log file in the attachment. In addition, pltooing for basic variables from the restart file.

I am new to ocean modeling, so I would appreciate any suggestions to resolve this issue.
Thank you,
Sree

https://drive.google.com/drive/folders/ ... sp=sharing

[jcwarner]

3 2023-06-07 00:01:30.00 2.974222E-03 1.787333E+04 1.787333E+04 7.224398E+15
(038,107,01) 8.201589E-03 5.863093E-01 6.355423E-04 1.714799E+02
DIAG speed trouble 171.479911641918
Found Error: 01 Line: 355 Source: ROMS/Nonlinear/main3d.F
Found Error: 01 Line: 330 Source: ROMS/Drivers/nl_roms.h, ROMS_run

Blowing-up: Saving latest model state into RESTART file
REASON: MaxSpeed = 1.7E+02

So you need to look at the ocean_rst file to see where this large velocity is occurring.
Is it at the open boundary? at the land sea boundary? at a strong bathy gradient?

[yernisrinivasnekkali]

"Thank you, Sir, for your reply.

The image is plotted from the RST file. I observed the gradient (small patch) and highlighted it with a circle in the image at the attached link.

1. There is a large gradient in Vbar and salt at the top of the East boundary (closed boundary).
2. Meanwhile, zeta values appear at the open boundary (West and South).
3. In addition, as per my knowledge, the bathymetry is not very steep at that (For your reference, I have attached an image)
What is the best solution to overcome this?

Thank you,
Sree

https://drive.google.com/file/d/1lMwZzC ... sp=sharing

[jcwarner]

lets start with that spot of high salt and vbar near the middle top.
what is the h value there? i see your min h is 0.4 m. that might be too small for the ~8km grid spacing.