import pandas as pd
import numpy as np
import os
from datetime import datetime, timedelta
from opendrift.readers import reader_shape
from opendrift.readers import reader_ROMS_native
from opendrift.models.pelagicegg import PelagicEggDrift


os.chdir('C:/…')


siembraJunio = pd.read_excel('C:/....')
Lon = siembraMayo.iloc[:,3] 
Lat = siembraMayo.iloc[:,2]


o = PelagicEggDrift(loglevel=20)  
Landmask = reader_shape.Reader.from_shpfiles('C:/.../SHP/xx.shp')
uv = reader_ROMS_native.Reader('C:/…./croco06.nc')
o.add_reader([Landmask, uv])


o.set_config('drift:max_age_seconds', 1296000) 
start_time = uv.start_time

for i in range(0, 3, 2):  
    o.seed_elements(Lon, Lat, origin_marker=i, time=start_time + timedelta(days=i)) 

o.run(duration=timedelta(days=46), outfile='CROCOdrift.nc')  


o.plot(fast=True)
o.animation(fast=True, color='z', filename='animationCROCO.gif')