import os
from opendrift.readers import reader_ROMS_native
from opendrift.models.oceandrift import OceanDrift


os.chdir('E:/')


o = OceanDrift(loglevel=0)

filename_nc = 'croco_var.nc';
mosa_native = reader_ROMS_native.Reader(filename_nc)
o.add_reader([mosa_native])



time_ini = mosa_native.start_time


lon = -109.9; lat = 25.6;  # P4

o.seed_elements(lon, lat, number=200, radius = 3000, time=time_ini)

lon = -110.0; lat = 25.7;  # P4

o.seed_elements(lon, lat, number=200, radius = 3000, time=time_ini)


o.run(time_step=60,outfile='GCalifornia.nc')

o.plot(filename='GCalifornia.png')