#%% # Setting some parameters duration = timedelta(seconds=12) # T time_step=timedelta(seconds=.5) time_step_output=timedelta(seconds=.5) delta=.02 # spatial resolution steps = int(duration.total_seconds()/ time_step_output.total_seconds() + 1) o = OceanDrift(loglevel=20) #%% # Note that Runge-Kutta here makes a difference to Euler scheme o.set_config('drift:advection_scheme', 'runge-kutta4') o.disable_vertical_motion() o.set_config('environment:fallback:land_binary_mask', 0) double_gyre = reader_double_gyre.Reader(epsilon=.25, omega=0.628, A=0.1) print(double_gyre) o.add_reader(double_gyre) #%% # Calculate Lyapunov exponents times = [double_gyre.initial_time + n*time_step_output for n in range(steps)] lcs = o.calculate_ftle(time=times, time_step=time_step, duration=duration, delta=delta, RLCS=False) #%% # Make run with particles for the same period
o.seed_elements(lon=174.046669, lat=-40.928116, radius=20, number=100, z=np.linspace(0, -10, 100), time=schism_native.start_time) o.seed_elements(lon=173.8839, lat=-40.9160, radius=20, number=100, z=np.linspace(0, -10, 100), time=schism_native.start_time) o.seed_elements(lon=174.2940, lat=-41.0888, radius=20, number=100, z=np.linspace(0, -10, 100), time=schism_native.start_time) o.disable_vertical_motion() #Deactivate any vertical processes/advection""" #%% # Running model o.run(time_step=900, end_time=schism_native.start_time + timedelta(days=1.0)) # outfile='schism_native_output.nc') # Print and plot results print(o) o.plot(fast=True) o.animation() o.animation_profile()