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.
use time, only: time_diff,julian_day use observations, only: init_saved_vars,read_obs use observations, only: Hs,Tz,phiw IMPLICIT NONEINPUT PARAMETERS:
integer, intent(in) :: unit,jul,secsREVISION HISTORY:
Original author(s): Karsten BoldingLOCAL 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