program main_grid
    use mod_harm2grid
    implicit none
    integer :: itime

    call read_input
    allocate(theta(1:180), phi(1:360))
    call generate_theta(size(theta))
    call generate_phi(size(phi))
    do itime = 1, 10
        write(*,*) 'Processing time step:', itime
        call ocean_velocity_xyz(itime)
        write(*,*) 'Time step:', itime, 'done'
    end do
end program main_grid