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

I have 3 open boundary except north.when I run the model I get weird zeta vlaue in my result.#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

What seems be the problem?

Any tips would be great for me.

Best,

-Peter