This routine will provide sea surface elevation - either by an analytical expression or read from file. The subroutine is called in the get_all_obs() subroutine as part of the main integration loop. The spatial interpolation is done via the reading routine and the temporal interpolation is done in this routine.
use time, only: time_diff,julian_day,fsecs use observations, only: pi,init_saved_vars,read_obs use observations, only: period_1,amp_1,phase_1,period_2,amp_2,phase_2 use observations, only: zeta,zeta_0 IMPLICIT NONEINPUT PARAMETERS:
integer, intent(in) :: method,unit,jul,secsREVISION HISTORY:
Original author(s): Karsten Bolding