Tide open boundry?

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
yang
Posts: 56
Joined: Mon Sep 25, 2006 2:37 pm
Location: Institue of oceanology ,Chinese acedemy of scinece

Tide open boundry?

#1 Unread post by yang »

I am runing a basin scale climatology simulation in west Pacific ocean with grid 5'x5'.
It will be ok for runing 5 years with the following configuration without tide.
#define EAST_FSCHAPMAN
#define EAST_M2FLATHER
#define EAST_M3RADIATION
#define EAST_TNUDGING
If included tide,(that is to say , include the configuration as following),the model will always blow up nomatter how small the dt is (even dt=60 second;ndtfast=90)?


#define SSH_TIDES
#defing UV_TIDES
#define RAMP_TIDES
#ifdef SSH_TIDES
#define ADD_FSOBC
#endif

#ifdef UV_TIDES
#undef FSOBC_REDUCED
#define ADD_M2OBC
#endif
Can i use #define FSFLATHER; #define M2NUDGING to simulate the tide ??

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

Re: Tide open boundry?

#2 Unread post by kate »

I don't see anything wrong with your tides setup - can you run in 2-D with just tides? How about without the ADD_xxx? Can you plot the ROMS output to see what it is doing?

For the 3-D, what I use is:

Code: Select all

# define NORTH_FSCHAPMAN
# define NORTH_M2FLATHER
# ifdef SOLVE3D
#  define NORTH_M3RADIATION
#  define NORTH_M3NUDGING
#  define NORTH_TRADIATION
#  define NORTH_TNUDGING
# endif
In other words, both nudging and radiation for both M3 and T. You would then need some external source of information for both, plus the nudging coefficients and the OBCFAC.

yang
Posts: 56
Joined: Mon Sep 25, 2006 2:37 pm
Location: Institue of oceanology ,Chinese acedemy of scinece

Re: Tide open boundry?

#3 Unread post by yang »

Thank you a lot. Kate

In the above runing , I have changed OBCFAC from 0 to 2.0d+1.It always blows up in the second day.
OBCFAC == 2.0d+1 ! nondimensional Ydz changed from 0.0d0
In addition, i have make a configurating for nudging as following.
TNUDG == 2*3.0d+1 ! days Ydz changed from 0.0d0
ZNUDG == 3.0d+1 ! days Ydz changed from 0.0d0
M2NUDG == 3.0d+1 ! days Ydz changed from 0.0d0
M3NUDG == 3.0d+1 ! days Ydz changed from 0.0d0

Now, i am runing the tide again with the combination of boudary CPP(#define M3Radiation and #define M3Nudging;#define TRadiation and #define TNudging; ) for T and M3.

yang
Posts: 56
Joined: Mon Sep 25, 2006 2:37 pm
Location: Institue of oceanology ,Chinese acedemy of scinece

Re: Tide open boundry?

#4 Unread post by yang »

Kate,please, accept my very best wishes for the holiday season.

According your advice, it is ok now that the tide simulation can be run coupled with current for a whole year.

Another question, are there some other boundary CPPs which can be used in the combination form(for exmple: for M3 ,use both #define M3RADIATION and #define M3NUDGING)?

Where can i find the information about the CPP combination ?

I have read a lot of information about ROMS, but i have not find the introduction about the CPP combination. In addition, i didn't find the CPP combination in the *.h in ROMS/Include

User avatar
susonic
Posts: 171
Joined: Tue Aug 21, 2007 5:44 pm
Location: UST21 / Korea
Contact:

Re: Tide open boundry?

#5 Unread post by susonic »

Hi, yang
please check out below address and download nobuhitomori's pdf file.
There is a boundary code information in that file.

viewtopic.php?f=17&t=930&p=4200#p4200

Hope that it helps to you.

-Peter
Joonho Lee

yang
Posts: 56
Joined: Mon Sep 25, 2006 2:37 pm
Location: Institue of oceanology ,Chinese acedemy of scinece

Re: Tide open boundry?

#6 Unread post by yang »

Thanks.
I have downded that 3 days ago and it is helpful for me. But there are not information about combination of open boundary CCP.

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

Re: Tide open boundry?

#7 Unread post by kate »

You have all I know about the open boundaries. Someone told me to try that many years ago and it's what I've been using ever since. It would be great if someone more informed would write to the wiki on boundary conditions, but don't hold your breath.

Post Reply