Restarting a Tidal Flow

[thomas.roc]

Hi Guys,

I face troubles restarting a tidal flow simulation from an *_his.nc file.
The model restart perfectly using the *_his.nc but the tidal forcing doesn't follow and damp to zero after the first steps. Notice that the tidal forcing is specified through the ana_fobc file as follow:
fac=0.75_r8 ! zeta0
omega=2.0_r8*pi/(12.42_r8*3600.0_r8) ! M2 Tide period
# ifdef WEST_FSOBC
val=fac*SIN(omega*time(ng))
val=val*(1-exp(-time(ng)/2000.d0)) !Ramp up over ~3 hours
IF (WESTERN_EDGE) THEN
DO j=JstrR,JendR
BOUNDARY(ng)%zeta_west(j)=val
END DO
END IF
# endif
# ifdef EAST_FSOBC
phase=0.6266 ! radians, calculated by dist*freq/speed of wave
val=fac*SIN(omega*time(ng)+phase)
val=val*(1-exp(-time(ng)/2000.d0)) !Ramp up over ~3 hours
IF (EASTERN_EDGE) THEN
DO j=JstrR,JendR
BOUNDARY(ng)%zeta_east(j)=val
END DO
END IF
# endif

Any thoughts ?

Cheers,


-T-