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 $Log: get_zeta.F90,v $ Revision 1.8 2010-09-17 12:53:51 jorn extensive code clean-up to ensure proper initialization and clean-up of all variables Revision 1.7 2006-11-27 09:25:18 kbk use logical var init_saved_vars to initialise saved variables Revision 1.6 2005/06/27 13:44:07 kbk modified + removed traling blanks Revision 1.5 2003/03/28 09:20:35 kbk added new copyright to files Revision 1.4 2003/03/28 09:02:09 kbk removed tabs Revision 1.3 2003/03/10 08:51:58 gotm Improved documentation and cleaned up code Revision 1.2 2001/11/18 16:06:31 gotm Avoid namelist member clashes by changing names in zetaspec Revision 18.104.22.168 2001/02/12 15:55:58 gotm initial import into CVS