get_zeta

INTERFACE:

    subroutine get_zeta(method,unit,jul,secs)

DESCRIPTION:

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.

USES:

    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 NONE

INPUT PARAMETERS:

    integer, intent(in)                 :: method,unit,jul,secs

REVISION 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 1.1.1.1  2001/02/12 15:55:58  gotm
    initial import into CVS