program xbeach

use params
use spaceparams


IMPLICIT NONE

type(parameters)        ::par
type(spacepars)         ::s

integer*4				::nt
integer					::it
real*8					::degrad,thetamin,thetamax,tbegin,tend

call cpu_time(tbegin)

write(*,*) 'Welcome to Xbeach'
write(*,*) 'General Input Module'
	
! General input per module
call wave_input(par)
call flow_input(par)
call sed_input(par)


write(*,*) 'Building Grid and Bathymetry and....'
write(*,*) 'Distributing wave energy across the directional space ....'		
! Grid and bathymetry
call grid_bathy(s,par)	

! Jump into subroutine readtide
if(par%tideloc>=1)then 
   call readtide (s,par)  !Ap 15/10
end if

par%t=0.
it=0

write(*,*) 'Initializing .....'	
! Initialisations
call wave_init (s,par)
call flow_init (s,par)
call sed_init (s,par)
call init_output(s,par,it)

call readkey('params.txt','checkparams','')




write(*,*) 'Stepping into the time loop ....'	

do while (par%t<par%tstop)
    ! Calculate timestep
	call timestep(s,par,it)
	! Wave boundary conditions
    call wave_bc (s,par);
    ! Flow boundary conditions
    call flow_bc (s,par);
    ! Wave timestep
    if (par%instat==0) then
       if (abs(mod(par%t,real(par%wavint)))<1.e-4) call wave_stationary(s,par)
	else
       call wave_timestep(s,par)
    endif
    ! Flow timestep
    call flow_timestep (s,par)
    ! Suspended transport
    call transus(s,par)
    ! Bed level update
    call bed_update(s,par)
    ! Output
    call var_output(it,s,par)

enddo

call cpu_time(tend)
write(*,*)'Total cpu time: ',tend-tbegin,' seconds'


end program