## Tidal effection with analytical expression

General scientific issues regarding ROMS

Moderators: arango, robertson

Message
Author
susonic
Posts: 160
Joined: Tue Aug 21, 2007 5:44 pm
Location: Jeju National University
Contact:

### Tidal effection with analytical expression

Hi all,
I want to use tidal effect using analytical expression in my domain. Since my domain is too small to use large scale data such as TPXO6. So, I think it's better to use the ana_fsobc.h to make tadal effect. Hence I put this
m2=0.757_r8
s2=0.337_r8
k1=0.245_r8
o1=0.177_r8
ttime=time(ng)+162.0_r8*3600.0_r8
m2t=2.0_r8*pi*ttime/(12.42_r8*3600.0_r8)
s2t=2.0_r8*pi*ttime/(12.00_r8*3600.0_r8)
k1t=2.0_r8*pi*ttime/(23.93_r8*3600.0_r8)
o1t=2.0_r8*pi*ttime/(25.82_r8*3600.0_r8)
IF (WESTERN_EDGE) THEN
DO j=JstrR,JendR
BOUNDARY(ng)%zeta_west(j)=m2*sin(m2t+274.6*pi/180.0)+ &
& s2*sin(s2t+294.4*pi/180.0)+ &
& k1*sin(k1t+201.3*pi/180.0)+ &
& o1*sin(o1t+177.6*pi/180.0)
END DO
END IF
IF (EASTERN_EDGE) THEN
DO j=JstrR,JendR
BOUNDARY(ng)%zeta_east(j)=m2*sin(m2t+274.6*pi/180.0)+ &
& s2*sin(s2t+294.4*pi/180.0)+ &
& k1*sin(k1t+201.3*pi/180.0)+ &
& o1*sin(o1t+177.6*pi/180.0)
IF (SOUTHERN_EDGE) THEN
DO i=IstrR,IendR
BOUNDARY(ng)%zeta_east(i)=m2*sin(m2t+274.6*pi/180.0)+ &
& s2*sin(s2t+294.4*pi/180.0)+ &
& k1*sin(k1t+201.3*pi/180.0)+ &
& o1*sin(o1t+177.6*pi/180.0)

END DO
END IF
#define NORTHERN_WALL
#define EAST_M2FLATHER
#define EAST_FSCHAPMAN
#define EAST_M3CLAMPED
#define EAST_TCLAMPED
#define SOUTH_M2FLATHER
#define SOUTH_FSCHAPMAN
#define SOUTH_M3CLAMPED
#define SOUTH_TCLAMPED
#define WEST_M2FLATHER
#define WEST_FSCHAPMAN
#define WEST_M3CLAMPED
#define WEST_TCLAMPED
#define ANA_M2OBC
#define ANA_FSOBC
I have 3 open boundary except north.when I run the model I get weird zeta vlaue in my result.
What seems be the problem?
Any tips would be great for me.
Best,
-Peter

wilkin
Posts: 526
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:

### Re: Tidal effection with analytical expression

I think you would find it much easier to build your "analytical" tidal forcing conditions into an external tides forcing netcdf file, and then give this as an input in ocean.in in the same way that most users apply the TPXO and ADCIRC (and other) tidal model forcings.
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu

johnluick
Posts: 97
Joined: Mon Jun 20, 2005 3:46 pm
Location: South Australian Research and Development Institute

### Re: Tidal effection with analytical expression

I agree with Wilkin, but if you insist on trying this (i.e. you don't care about accurately reproducing the tides, just want some tidal dispersion) just drive one boundary the way the kelvin test case does. Your tides are in phase all around the boundary! You are driving it like an eskimo blanket toss. Note that in the kelvin test case the tidal velocities (out of phase with sea levels) are also specified. But in the end, you'll never know how realistic your results are.

susonic
Posts: 160
Joined: Tue Aug 21, 2007 5:44 pm
Location: Jeju National University
Contact: