Python run_assimilator Examples

Example #1

File: assimilate_csv_on_advector.py Project: gaelforget/Ocean-Plastic-Assimilator

from src.run_assimilator import run_assimilator

if __name__ == "__main__":
    run_assimilator(
        particles_dataset_path="data/advector_output_rivers_2012.nc",
        observations_type="from_csv",
        observations_source_path="data/observations.csv",
        assimilation_domain_coords=(195 - 360, 23, 235 - 360, 45),
        assimilation_grid_size=(int(40 / 0.5), int(22 / 0.5)),
        size_ensemble=10,
        initial_ensemble_spread_percent=0.05,
        t_start=0,
        t_end=300,
        radius_observation=50,
    )

Example #2

File: assimilate_advector.py Project: gaelforget/Ocean-Plastic-Assimilator

from src.run_assimilator import run_assimilator

if __name__ == "__main__":
    run_assimilator(
        particles_dataset_path="data/advector_output_rivers_2012.nc",
        observations_type="from_simulation",
        observations_source_path="data/advector_output_coastal_2012.nc",
        assimilation_domain_coords=(195 - 360, 23, 235 - 360, 45),
        assimilation_grid_size=(int(40 / 0.5), int(22 / 0.5)),
        size_ensemble=10,
        initial_ensemble_spread_percent=0.05,
        observations_error_percent=0.01,
        observation_locations=[(12, 4), (55, 27)],
        t_start=0,
        t_end=300,
        initial_mass_multiplicator=2,
        radius_observation=50,
    )

Example #3

File: assimilate_advector.py Project: TheOceanCleanup/OceanPlasticAssimilator

import numpy as np

cell_area = lambda dlon, dlat, lat: dlat * 111 * dlon * 111 * np.cos(lat * np.
                                                                     pi / 180)

GRID_COORDS = (-165, 23, -125, 45)
GRID_RESOLUTION = 0.5

cells_area = np.array([[
    cell_area(0.5, 0.5, lat_id * GRID_RESOLUTION + GRID_COORDS[1])
    for lat_id in range(44)
] for lon_id in range(80)])

if __name__ == "__main__":
    run_assimilator(
        particles_dataset_path="data/advector_output_rivers_2012.nc",
        observations_type="from_simulation",
        observations_source_path="data/advector_output_coastal_2012.nc",
        assimilation_domain_coords=GRID_COORDS,
        assimilation_grid_size=(int(40 / GRID_RESOLUTION),
                                int(22 / GRID_RESOLUTION)),
        size_ensemble=10,
        initial_ensemble_spread=50,
        observations_error_percent=0.01,
        observation_locations=[(12, 4), (55, 27)],
        t_start=0,
        t_end=365,
        radius_observation=np.inf,
        initial_mass_multiplicator=1000,  # Set to your initial particles mass
        cells_area=cells_area)

Example #4

File: assimilate_csv_on_advector.py Project: TheOceanCleanup/OceanPlasticAssimilator

from src.run_assimilator import run_assimilator
import numpy as np

cell_area = lambda dlon, dlat, lat: dlat * 111 * dlon * 111 * np.cos(lat * np.
                                                                     pi / 180)

GRID_COORDS = (-165, 23, -125, 45)
GRID_RESOLUTION = 0.5

cells_area = np.array([[
    cell_area(0.5, 0.5, lat_id * GRID_RESOLUTION + GRID_COORDS[1])
    for lat_id in range(44)
] for lon_id in range(80)])

if __name__ == "__main__":
    run_assimilator(
        particles_dataset_path="data/advector_output_rivers_2012.nc",
        observations_type="from_csv",
        observations_source_path="data/observations.csv",
        assimilation_domain_coords=GRID_COORDS,
        assimilation_grid_size=(int(40 / GRID_RESOLUTION),
                                int(22 / GRID_RESOLUTION)),
        size_ensemble=10,
        initial_ensemble_spread=
        50,  # Define, a tested, functionning value can be 0.5% of your initial particle mass
        t_start=0,
        t_end=365,
        radius_observation=np.inf,
        initial_mass_multiplicator=1000,  # Set to your initial particles mass
        cells_area=cells_area)

Example #5

File: assimilate_double_gyre.py Project: TheOceanCleanup/OceanPlasticAssimilator

from src.run_assimilator import run_assimilator

mu = 2
A = 0.1175
epsilon = 0.25
sigma_rel = 0.01

if __name__ == "__main__":
    run_assimilator(
        particles_dataset_path="data/parts_double_gyre_standard_as_3.nc",
        observations_type="from_simulation",
        observations_source_path=
        f"data/parts_double_gyre_ref_eps_{epsilon}_A_{A}_as_2.nc",
        simulation_name=f"mu_{mu}_A_{A}_eps_{epsilon}_sigma_rel_{sigma_rel}",
        assimilation_domain_coords=(195, 20, 225, 40),
        assimilation_grid_size=(60, 40),
        size_ensemble=10,
        initial_ensemble_spread=0.05,
        t_start=0,
        t_end=2000,
        observations_error_percent=sigma_rel,
        observation_locations=[(12, 4), (55, 27)],
        initial_mass_multiplicator=mu,
    )