Blow-Up on the boundaries

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
antoinemarie

Blow-Up on the boundaries

#1 Unread post by antoinemarie »

Hi ROMS'users!

I try to run a realistic case. The aim of this case is to estimate currents produce by the wind along the Reunion Island. For this, I decided to use ANA_SMFLUX instead of ANA_WIND. I have to impose the sea level too.
The problem is I don't know which boundary to use. I try a lot of combination of boundaries, but I always have instability along the boundaries. Currently, I use the following boudaries:

Code: Select all

#define FSOBC_REDUCED
#define ANA_FSOBC

#define SOUTHERN_WALL

#define NORTH_FSRADIATION
#define NORTH_M2RADIATION
#define NORTH_M3RADIATION

#define EAST_FSREDUCED
#define EAST_M2REDUCED
#define EAST_M3RADIATION

#define WEST_FSRADIATION
#define WEST_M2RADIATION
#define WEST_M3RADIATION
I impose EAST_FSREDUCED and EAST_M2REDUCED to keep the sea level equal to 2.1 meters.

I don't think it comes from my time step, because I have dx=dy=300 meters, and my time step is equal to 15 seconds...

I don't manage to have a stable model with a Wind equal to 10m/s. And I would like to simulate an hurricane Wind... :?

Someone can help me?

Thanks in advance

Antoine

User avatar
kate
Posts: 4091
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: Blow-Up on the boundaries

#2 Unread post by kate »

It has been years since I've tried anything besides a combination of radiation and nudging. Which boundary is yours blowing up along? The eastern or the others? What about the tracer boundary conditions? I don't see them in your list.

antoinemarie

Re: Blow-Up on the boundaries

#3 Unread post by antoinemarie »

Thank you for your reply Kate,

For the tracers, I use radiation conditions for each boundary.

The Northern and Western boundaries are blowing up. But I have tried to put radiation conditions (FS, M2 and M3) for each boundary, and sea level doesn't stay equal to 2.1m, it increases rapidly.

I have also tried to put reduced conditions (FS, M2) + M3RADIATION for each boundary, but speed blows up...

Post Reply