get_wave

INTERFACE:

    subroutine get_wave(unit,jul,secs)

DESCRIPTION:

This routine is responsible for providing sane values to `observed' wind generated waves. The observations consist of significant wave height (Hs), mean zero-crossing period (Tz) and mean direction (phiw). The variables can be set to constant values (wave_method=1) or read from file (wave_method=2). For wave_method=0 nothing is done. The subroutine is called in the get_all_obs() subroutine as part of the main integration loop. In case of observations from file the temporal interpolation is done in this routine.

USES:

    use time,         only: time_diff,julian_day
    use observations, only: init_saved_vars,read_obs
    use observations, only: Hs,Tz,phiw
    IMPLICIT NONE

INPUT PARAMETERS:

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

REVISION HISTORY:

    Original author(s): Karsten Bolding
    $Log: get_wave.F90,v $
    Revision 1.3  2010-09-17 12:53:51  jorn
    extensive code clean-up to ensure proper initialization and clean-up of all variables
    Revision 1.2  2007-01-04 12:19:09  kbk
    updated documentation
    Revision 1.1  2007-01-04 12:08:12  kbk
    adding surface waves

LOCAL VARIABLES:

    integer                   :: yy,mm,dd,hh,min,ss
    REALTYPE                  :: t
    REALTYPE, save            :: dt
    integer, save             :: jul1,secs1
    integer, save             :: jul2,secs2
    REALTYPE, save            :: alpha(3)
    REALTYPE, save            :: obs1(3),obs2(3)
    integer                   :: rc