Caspian / lake

from datetime import datetime, timedelta
from opendrift.models.oceandrift import OceanDrift
from opendrift.readers import reader_shape
import cartopy

The default GSHHG landmask reader classifies lakes (e.g. Caspian Sea) as land and not water. For simulations in lakes, we can use a custom reader for lakes only.

level = 'h'  # using high resolution coastline
reader_lakes = reader_shape.Reader.from_shpfiles(
    f'{cartopy.config['data_dir']}/shapefiles/gshhs/{level}/GSHHS_{level}_L2.shp',
    invert=True)  # Inverting since inside of polygons is water (lake) and not land.

Disabling the default landmask reader, and using the above reader instead

o = OceanDrift(loglevel=20)
o.add_reader(reader_lakes)
o.set_config('general:use_auto_landmask', False)  # To use custom landmask instead
o.set_config('environment:constant:x_sea_water_velocity', 1)
o.set_config('drift:horizontal_diffusivity', 100)
o.seed_elements(lon=48.819, lat=44.959, radius=5000, number=100, time=datetime.now())

o.run(steps=10)
o.plot(fast=False, buffer=2)

09:34:52 INFO    opendrift:569: OpenDriftSimulation initialised (version 1.14.3 / v1.14.3-46-g60beecd)
09:34:52 INFO    opendrift.models.basemodel.environment:229: Fallback values will be used for the following variables which have no readers:
09:34:52 INFO    opendrift.models.basemodel.environment:232:    y_sea_water_velocity: 0.000000
09:34:52 INFO    opendrift.models.basemodel.environment:232:    sea_surface_height: 0.000000
09:34:52 INFO    opendrift.models.basemodel.environment:232:    x_wind: 0.000000
09:34:52 INFO    opendrift.models.basemodel.environment:232:    y_wind: 0.000000
09:34:52 INFO    opendrift.models.basemodel.environment:232:    upward_sea_water_velocity: 0.000000
09:34:52 INFO    opendrift.models.basemodel.environment:232:    ocean_vertical_diffusivity: 0.000000
09:34:52 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_significant_height: 0.000000
09:34:52 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
09:34:52 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_y_velocity: 0.000000
09:34:52 INFO    opendrift.models.basemodel.environment:232:    ocean_mixed_layer_thickness: 50.000000
09:34:52 INFO    opendrift.models.basemodel.environment:232:    sea_floor_depth_below_sea_level: 10000.000000
09:34:52 INFO    opendrift:1783: Skipping environment variable ocean_vertical_diffusivity because of condition ['drift:vertical_mixing', 'is', False]
09:34:52 INFO    opendrift:1783: Skipping environment variable ocean_mixed_layer_thickness because of condition ['drift:vertical_mixing', 'is', False]
09:34:52 INFO    opendrift:1794: Storing previous values of element property lon because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
09:34:52 INFO    opendrift:1794: Storing previous values of element property lat because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
09:34:52 INFO    opendrift:1802: Storing previous values of environment variable sea_surface_height because of condition ['drift:vertical_advection', 'is', True]
09:34:52 INFO    opendrift:945: Using existing reader for land_binary_mask to move elements to ocean
09:34:52 INFO    opendrift.readers.reader_shape:149: Building KDTree from 6601 geometries with buffer distance 0.000000e+00
09:34:57 INFO    opendrift:2090: 2025-10-22 09:34:52.257217 - step 1 of 10 - 100 active elements (0 deactivated)
09:34:57 INFO    opendrift:2090: 2025-10-22 10:34:52.257217 - step 2 of 10 - 100 active elements (0 deactivated)
09:34:57 INFO    opendrift:2090: 2025-10-22 11:34:52.257217 - step 3 of 10 - 100 active elements (0 deactivated)
09:34:57 INFO    opendrift:2090: 2025-10-22 12:34:52.257217 - step 4 of 10 - 100 active elements (0 deactivated)
09:34:57 INFO    opendrift:2090: 2025-10-22 13:34:52.257217 - step 5 of 10 - 100 active elements (0 deactivated)
09:34:57 INFO    opendrift:2090: 2025-10-22 14:34:52.257217 - step 6 of 10 - 100 active elements (0 deactivated)
09:34:57 INFO    opendrift:2090: 2025-10-22 15:34:52.257217 - step 7 of 10 - 100 active elements (0 deactivated)
09:34:57 INFO    opendrift:2090: 2025-10-22 16:34:52.257217 - step 8 of 10 - 100 active elements (0 deactivated)
09:34:57 INFO    opendrift:2090: 2025-10-22 17:34:52.257217 - step 9 of 10 - 100 active elements (0 deactivated)
09:34:57 INFO    opendrift:2090: 2025-10-22 18:34:52.257217 - step 10 of 10 - 100 active elements (0 deactivated)

(<GeoAxes: title={'center': 'OpenDrift - OceanDrift\n2025-10-22 09:34 to 2025-10-22 19:34 UTC (11 steps)'}>, <Figure size 1100x755.105 with 1 Axes>)