roms-tools 1.2.0__py3-none-any.whl → 1.4.0__py3-none-any.whl

{roms_tools-1.2.0.dist-info → roms_tools-1.4.0.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (73.0.1)
+Generator: setuptools (75.1.0)
 Root-Is-Purelib: true
 Tag: py3-none-any
 

roms_tools/tests/test_boundary_forcing.py

@@ -1,658 +0,0 @@
-import pytest
-from datetime import datetime
-from roms_tools import BoundaryForcing, Grid
-import numpy as np
-import tempfile
-import os
-import textwrap
-from roms_tools.setup.download import download_test_data
-
-
-@pytest.fixture
-def example_grid():
-    """
-    Fixture for creating a Grid object.
-    """
-    grid = Grid(
-        nx=2,
-        ny=2,
-        size_x=500,
-        size_y=1000,
-        center_lon=0,
-        center_lat=55,
-        rot=10,
-        N=3,  # number of vertical levels
-        theta_s=5.0,  # surface control parameter
-        theta_b=2.0,  # bottom control parameter
-        hc=250.0,  # critical depth
-    )
-
-    return grid
-
-
-@pytest.fixture
-def boundary_forcing(example_grid):
-    """
-    Fixture for creating a BoundaryForcing object.
-    """
-
-    fname = download_test_data("GLORYS_test_data.nc")
-
-    return BoundaryForcing(
-        grid=example_grid,
-        start_time=datetime(2021, 6, 29),
-        end_time=datetime(2021, 6, 30),
-        physics_source={"name": "GLORYS", "path": fname},
-    )
-
-
-@pytest.fixture
-def boundary_forcing_with_bgc_from_climatology(example_grid):
-    """
-    Fixture for creating a BoundaryForcing object.
-    """
-
-    fname = download_test_data("GLORYS_test_data.nc")
-    fname_bgc = download_test_data("CESM_regional_test_data_climatology.nc")
-
-    return BoundaryForcing(
-        grid=example_grid,
-        start_time=datetime(2021, 6, 29),
-        end_time=datetime(2021, 6, 30),
-        physics_source={"name": "GLORYS", "path": fname},
-        bgc_source={"path": fname_bgc, "name": "CESM_REGRIDDED", "climatology": True},
-    )
-
-
-@pytest.mark.parametrize(
-    "bdry_forcing_fixture",
-    [
-        "boundary_forcing",
-        "boundary_forcing_with_bgc_from_climatology",
-    ],
-)
-def test_boundary_forcing_creation(bdry_forcing_fixture, request):
-    """
-    Test the creation of the BoundaryForcing object.
-    """
-
-    bdry_forcing = request.getfixturevalue(bdry_forcing_fixture)
-
-    assert bdry_forcing.start_time == datetime(2021, 6, 29)
-    assert bdry_forcing.end_time == datetime(2021, 6, 30)
-
-    assert bdry_forcing.ds["physics"].physics_source == "GLORYS"
-    for direction in ["south", "east", "north", "west"]:
-        assert f"temp_{direction}" in bdry_forcing.ds["physics"]
-        assert f"salt_{direction}" in bdry_forcing.ds["physics"]
-        assert f"u_{direction}" in bdry_forcing.ds["physics"]
-        assert f"v_{direction}" in bdry_forcing.ds["physics"]
-        assert f"zeta_{direction}" in bdry_forcing.ds["physics"]
-    assert len(bdry_forcing.ds["physics"].bry_time) == 1
-
-
-def test_boundary_forcing_creation_with_bgc(boundary_forcing_with_bgc_from_climatology):
-    """
-    Test the creation of the BoundaryForcing object.
-    """
-
-    assert (
-        boundary_forcing_with_bgc_from_climatology.ds["bgc"].bgc_source
-        == "CESM_REGRIDDED"
-    )
-    for direction in ["south", "east", "north", "west"]:
-        for var in ["ALK", "PO4"]:
-            assert (
-                f"{var}_{direction}"
-                in boundary_forcing_with_bgc_from_climatology.ds["bgc"]
-            )
-
-    assert len(boundary_forcing_with_bgc_from_climatology.ds["bgc"].bry_time) == 12
-
-
-def test_boundary_forcing_data_consistency_plot_save(
-    boundary_forcing_with_bgc_from_climatology,
-):
-    """
-    Test that the data within the BoundaryForcing object remains consistent.
-    Also test plot and save methods in the same test since we dask arrays are already computed.
-    """
-    # Define the expected data
-    expected_zeta_south = np.array(
-        [[-0.30468762, -0.29416865, -0.30391693, -0.32985148]], dtype=np.float32
-    )
-    expected_zeta_east = np.array(
-        [[-0.32985148, -0.36176518, -0.40663475, -0.40699923]], dtype=np.float32
-    )
-    expected_zeta_north = np.array(
-        [[-0.5534979, -0.5270749, -0.45107934, -0.40699923]], dtype=np.float32
-    )
-    expected_zeta_west = np.array(
-        [[-0.30468762, -0.34336275, -0.3699948, -0.5534979]], dtype=np.float32
-    )
-
-    expected_temp_south = np.array(
-        [
-            [
-                [16.84414, 18.088203, 18.431192, 19.294329],
-                [16.905312, 18.121834, 18.496748, 19.30358],
-                [16.967817, 18.315424, 18.718002, 19.439777],
-            ]
-        ],
-        dtype=np.float32,
-    )
-    expected_temp_east = np.array(
-        [
-            [
-                [19.294329, 18.633307, 8.710737, 9.20282],
-                [19.30358, 18.637077, 11.25943, 10.667074],
-                [19.439777, 18.667465, 13.111585, 11.752404],
-            ]
-        ],
-        dtype=np.float32,
-    )
-    expected_temp_north = np.array(
-        [
-            [
-                [10.233599, 10.147332, 10.458557, 9.20282],
-                [10.546486, 10.502733, 11.209945, 10.667074],
-                [10.671082, 10.68275, 11.377164, 11.752404],
-            ]
-        ],
-        dtype=np.float32,
-    )
-    expected_temp_west = np.array(
-        [
-            [
-                [16.84414, 12.639833, 11.027701, 10.233599],
-                [16.905312, 13.479691, 11.650267, 10.546486],
-                [16.967817, 14.426711, 12.200586, 10.671082],
-            ]
-        ],
-        dtype=np.float32,
-    )
-
-    expected_u_south = np.array(
-        [[[-0.0, -0.0, -0.0], [-0.0, -0.0, -0.0], [0.0, -0.0, -0.0]]], dtype=np.float32
-    )
-    expected_u_east = np.array(
-        [
-            [
-                [0.0, -0.0, 0.06979556, 0.0211786],
-                [-0.0, -0.0, 0.06167743, 0.03679834],
-                [-0.0, -0.0, -0.02247071, 0.0274788],
-            ]
-        ],
-        dtype=np.float32,
-    )
-
-    expected_u_north = np.array(
-        [
-            [
-                [0.04268532, 0.04645353, 0.0211786],
-                [0.03889201, 0.04914769, 0.03679834],
-                [0.03351666, 0.03673013, 0.0274788],
-            ]
-        ],
-        dtype=np.float32,
-    )
-
-    expected_u_west = np.array(
-        [
-            [
-                [-0.0, 0.0, 0.0, 0.04268532],
-                [-0.0, -0.0, 0.0, 0.03889201],
-                [-0.0, -0.0, -0.0, 0.03351666],
-            ]
-        ],
-        dtype=np.float32,
-    )
-
-    expected_v_south = np.array(
-        [[[0.0, 0.0, -0.0, -0.0], [0.0, 0.0, -0.0, -0.0], [0.0, -0.0, -0.0, -0.0]]],
-        dtype=np.float32,
-    )
-
-    expected_v_east = np.array(
-        [
-            [
-                [-0.0, -0.0, -0.06720348],
-                [-0.0, -0.0, -0.08354441],
-                [-0.0, -0.0, -0.13835917],
-            ]
-        ],
-        dtype=np.float32,
-    )
-
-    expected_v_north = np.array(
-        [
-            [
-                [-0.00951457, -0.0, 0.01915873, -0.06720348],
-                [-0.00576979, -0.0, 0.02625698, -0.08354441],
-                [-0.02147919, -0.0, 0.01757628, -0.13835917],
-            ]
-        ],
-        dtype=np.float32,
-    )
-    expected_v_west = np.array(
-        [
-            [
-                [0.0, -0.0, -0.00951457],
-                [0.0, -0.0, -0.00576979],
-                [0.0, -0.0, -0.02147919],
-            ]
-        ],
-        dtype=np.float32,
-    )
-
-    expected_ubar_south = np.array([[0.0, 0.0, 0.0]], dtype=np.float32)
-    expected_ubar_east = np.array(
-        [[0.0, 0.0, 0.04028399, 0.02812303]], dtype=np.float32
-    )
-    expected_ubar_north = np.array(
-        [[0.03866891, 0.04446249, 0.02812303]], dtype=np.float32
-    )
-    expected_ubar_west = np.array([[0.0, 0.0, 0.0, 0.03866891]], dtype=np.float32)
-
-    expected_vbar_south = np.array([[0.0, 0.0, 0.0, 0.0]], dtype=np.float32)
-    expected_vbar_east = np.array([[0.0, 0.0, -0.09326097]], dtype=np.float32)
-    expected_vbar_north = np.array(
-        [[-0.01189703, 0.0, 0.02102064, -0.09326097]], dtype=np.float32
-    )
-    expected_vbar_west = np.array([[0.0, 0.0, -0.01189703]], dtype=np.float32)
-
-    expected_alk_south = np.array(
-        [
-            [
-                [2352.1636, 2333.1094, 2308.618, 2286.2327],
-                [2352.128, 2332.9932, 2308.4705, 2286.0442],
-                [2352.091, 2332.8738, 2308.3032, 2285.8623],
-            ],
-            [
-                [2354.7297, 2337.391, 2314.0232, 2294.0876],
-                [2354.6863, 2337.259, 2313.8638, 2293.8918],
-                [2354.6414, 2337.1223, 2313.6616, 2293.6626],
-            ],
-            [
-                [2355.02, 2336.321, 2312.7905, 2292.2637],
-                [2354.9158, 2336.117, 2312.4849, 2291.8752],
-                [2354.817, 2335.9214, 2312.2444, 2291.6702],
-            ],
-            [
-                [2355.2583, 2334.3098, 2309.5273, 2287.4226],
-                [2355.052, 2333.913, 2308.292, 2285.5466],
-                [2354.9163, 2333.651, 2307.8772, 2285.0479],
-            ],
-            [
-                [2354.1216, 2330.6604, 2304.88, 2280.6],
-                [2353.7737, 2329.948, 2302.105, 2276.278],
-                [2353.6316, 2329.6023, 2301.4136, 2275.2761],
-            ],
-            [
-                [2346.881, 2322.218, 2301.2046, 2279.729],
-                [2345.7993, 2320.3225, 2291.6438, 2265.0488],
-                [2345.4675, 2319.691, 2290.308, 2263.0654],
-            ],
-            [
-                [2336.2632, 2312.938, 2293.6711, 2272.6418],
-                [2335.273, 2310.7412, 2279.314, 2250.8848],
-                [2334.9324, 2310.055, 2277.593, 2248.0718],
-            ],
-            [
-                [2333.8801, 2311.3232, 2289.1548, 2266.9172],
-                [2333.265, 2309.5337, 2276.1438, 2246.1714],
-                [2333.0667, 2309.034, 2274.885, 2244.1096],
-            ],
-            [
-                [2336.8193, 2314.2544, 2288.887, 2263.7944],
-                [2336.5115, 2313.1343, 2280.0654, 2249.8926],
-                [2336.4402, 2312.9004, 2279.546, 2249.0842],
-            ],
-            [
-                [2341.18, 2317.8801, 2289.4714, 2262.6213],
-                [2341.0679, 2317.4556, 2286.2214, 2257.4502],
-                [2341.025, 2317.3325, 2285.9263, 2256.9663],
-            ],
-            [
-                [2344.9026, 2323.6404, 2297.1284, 2273.7996],
-                [2344.86, 2323.5337, 2296.8164, 2273.2507],
-                [2344.8218, 2323.442, 2296.6846, 2273.0938],
-            ],
-            [
-                [2349.1636, 2329.3665, 2304.1536, 2281.6301],
-                [2349.1338, 2329.2893, 2304.0178, 2281.4062],
-                [2349.1035, 2329.2102, 2303.9006, 2281.269],
-            ],
-        ],
-        dtype=np.float32,
-    )
-
-    expected_alk_east = np.array(
-        [
-            [
-                [2286.2327, 2268.2373, 2297.2825, 2349.247],
-                [2286.0442, 2268.162, 2296.5747, 2349.1697],
-                [2285.8623, 2268.0366, 2296.2695, 2349.0938],
-            ],
-            [
-                [2294.0876, 2272.6245, 2301.072, 2352.814],
-                [2293.8918, 2272.5269, 2300.6504, 2352.756],
-                [2293.6626, 2272.3384, 2300.4255, 2352.6907],
-            ],
-            [
-                [2292.2637, 2272.901, 2304.492, 2352.8616],
-                [2291.8752, 2272.7153, 2303.8516, 2352.795],
-                [2291.6702, 2272.5413, 2303.6313, 2352.7444],
-            ],
-            [
-                [2287.4226, 2269.8926, 2310.2466, 2350.9084],
-                [2285.5466, 2269.2112, 2307.2576, 2350.3499],
-                [2285.0479, 2268.8284, 2306.0806, 2350.0247],
-            ],
-            [
-                [2280.6, 2266.4937, 2309.5537, 2350.1707],
-                [2276.278, 2264.7148, 2297.2288, 2348.5422],
-                [2275.2761, 2263.8596, 2293.6697, 2347.7415],
-            ],
-            [
-                [2279.729, 2263.0278, 2311.0305, 2350.3872],
-                [2265.0488, 2254.806, 2290.098, 2345.0693],
-                [2263.0654, 2253.1594, 2284.9631, 2342.4902],
-            ],
-            [
-                [2272.6418, 2256.4866, 2309.7383, 2349.7207],
-                [2250.8848, 2239.6753, 2284.9133, 2339.6501],
-                [2248.0718, 2237.9412, 2277.729, 2335.4924],
-            ],
-            [
-                [2266.9172, 2249.1013, 2307.639, 2347.7605],
-                [2246.1714, 2234.224, 2275.4429, 2330.1702],
-                [2244.1096, 2232.722, 2266.299, 2325.427],
-            ],
-            [
-                [2263.7944, 2246.2002, 2301.0461, 2338.5127],
-                [2249.8926, 2236.0933, 2269.4448, 2321.7393],
-                [2249.0842, 2235.4578, 2263.889, 2316.971],
-            ],
-            [
-                [2262.6213, 2244.3662, 2292.1316, 2331.162],
-                [2257.4502, 2241.8052, 2274.763, 2322.5522],
-                [2256.9663, 2241.589, 2272.2124, 2320.8901],
-            ],
-            [
-                [2273.7996, 2254.6848, 2290.5637, 2333.583],
-                [2273.2507, 2254.4407, 2286.885, 2332.7646],
-                [2273.0938, 2254.3308, 2286.088, 2332.3955],
-            ],
-            [
-                [2281.6301, 2262.4211, 2296.095, 2340.852],
-                [2281.4062, 2262.3313, 2293.6157, 2340.6072],
-                [2281.269, 2262.2532, 2293.1213, 2340.4978],
-            ],
-        ],
-        dtype=np.float32,
-    )
-
-    expected_alk_north = np.array(
-        [
-            [
-                [2376.7993, 2375.0688, 2372.2307, 2349.247],
-                [2376.7961, 2375.0647, 2372.2188, 2349.1697],
-                [2376.7932, 2375.0596, 2372.201, 2349.0938],
-            ],
-            [
-                [2376.757, 2374.8325, 2371.8342, 2352.814],
-                [2376.7537, 2374.8286, 2371.8257, 2352.756],
-                [2376.7478, 2374.825, 2371.814, 2352.6907],
-            ],
-            [
-                [2377.0188, 2374.9753, 2371.7104, 2352.8616],
-                [2377.0137, 2374.9707, 2371.7085, 2352.795],
-                [2377.0073, 2374.9666, 2371.707, 2352.7444],
-            ],
-            [
-                [2377.3914, 2375.5757, 2371.869, 2350.9084],
-                [2377.383, 2375.5718, 2371.8706, 2350.3499],
-                [2377.3792, 2375.57, 2371.8743, 2350.0247],
-            ],
-            [
-                [2378.4722, 2377.151, 2373.7805, 2350.1707],
-                [2378.5, 2377.2292, 2373.9, 2348.5422],
-                [2378.5017, 2377.2666, 2374.0066, 2347.7415],
-            ],
-            [
-                [2381.3555, 2380.8108, 2378.7131, 2350.3872],
-                [2383.3838, 2382.6125, 2379.7717, 2345.0693],
-                [2383.9692, 2383.3374, 2380.4094, 2342.4902],
-            ],
-            [
-                [2388.7502, 2387.0815, 2381.9631, 2349.7207],
-                [2389.4453, 2388.0613, 2383.4429, 2339.6501],
-                [2389.6758, 2388.3323, 2383.521, 2335.4924],
-            ],
-            [
-                [2386.29, 2383.9194, 2380.237, 2347.7605],
-                [2385.9524, 2383.5654, 2377.2883, 2330.1702],
-                [2385.863, 2383.4429, 2376.6057, 2325.427],
-            ],
-            [
-                [2375.1409, 2370.7205, 2371.8071, 2338.5127],
-                [2373.7124, 2368.2014, 2365.7, 2321.7393],
-                [2373.2695, 2367.7415, 2365.1226, 2316.971],
-            ],
-            [
-                [2371.951, 2365.2358, 2365.8464, 2331.162],
-                [2371.8308, 2364.598, 2362.7102, 2322.5522],
-                [2371.786, 2364.4639, 2362.5464, 2320.8901],
-            ],
-            [
-                [2373.2239, 2370.007, 2367.3286, 2333.583],
-                [2373.2004, 2369.9714, 2367.2654, 2332.7646],
-                [2373.1711, 2369.9507, 2367.2217, 2332.3955],
-            ],
-            [
-                [2374.99, 2373.9058, 2370.9688, 2340.852],
-                [2374.9817, 2373.8901, 2370.9421, 2340.6072],
-                [2374.9717, 2373.8748, 2370.907, 2340.4978],
-            ],
-        ],
-        dtype=np.float32,
-    )
-
-    expected_alk_west = np.array(
-        [
-            [
-                [2352.1636, 2335.311, 2362.5132, 2376.7993],
-                [2352.128, 2335.1162, 2362.5088, 2376.7961],
-                [2352.091, 2334.924, 2362.4663, 2376.7932],
-            ],
-            [
-                [2354.7297, 2335.165, 2363.018, 2376.757],
-                [2354.6863, 2334.9695, 2363.0068, 2376.7537],
-                [2354.6414, 2334.78, 2362.965, 2376.7478],
-            ],
-            [
-                [2355.02, 2334.8015, 2365.7913, 2377.0188],
-                [2354.9158, 2334.5276, 2365.8025, 2377.0137],
-                [2354.817, 2334.272, 2365.7432, 2377.0073],
-            ],
-            [
-                [2355.2583, 2335.5789, 2370.7808, 2377.3914],
-                [2355.052, 2335.267, 2370.754, 2377.383],
-                [2354.9163, 2335.0227, 2370.6836, 2377.3792],
-            ],
-            [
-                [2354.1216, 2335.7122, 2375.5793, 2378.4722],
-                [2353.7737, 2335.1372, 2375.3132, 2378.5],
-                [2353.6316, 2334.8125, 2375.1387, 2378.5017],
-            ],
-            [
-                [2346.881, 2332.2158, 2378.4954, 2381.3555],
-                [2345.7993, 2331.5295, 2376.4001, 2383.3838],
-                [2345.4675, 2331.2346, 2375.9187, 2383.9692],
-            ],
-            [
-                [2336.2632, 2328.7188, 2374.763, 2388.7502],
-                [2335.273, 2327.8726, 2370.0576, 2389.4453],
-                [2334.9324, 2327.576, 2369.491, 2389.6758],
-            ],
-            [
-                [2333.8801, 2326.699, 2356.8801, 2386.29],
-                [2333.265, 2326.0037, 2352.0618, 2385.9524],
-                [2333.0667, 2325.771, 2351.4526, 2385.863],
-            ],
-            [
-                [2336.8193, 2325.9893, 2346.735, 2375.1409],
-                [2336.5115, 2325.512, 2344.9722, 2373.7124],
-                [2336.4402, 2325.2942, 2344.6653, 2373.2695],
-            ],
-            [
-                [2341.18, 2328.2231, 2352.159, 2371.951],
-                [2341.0679, 2327.9817, 2351.6226, 2371.8308],
-                [2341.025, 2327.8281, 2351.4314, 2371.786],
-            ],
-            [
-                [2344.9026, 2331.0703, 2356.0095, 2373.2239],
-                [2344.86, 2330.9097, 2355.9424, 2373.2004],
-                [2344.8218, 2330.7537, 2355.8538, 2373.1711],
-            ],
-            [
-                [2349.1636, 2334.3328, 2359.3823, 2374.99],
-                [2349.1338, 2334.043, 2359.3276, 2374.9817],
-                [2349.1035, 2333.749, 2359.2634, 2374.9717],
-            ],
-        ],
-        dtype=np.float32,
-    )
-
-    ds = boundary_forcing_with_bgc_from_climatology.ds["physics"].compute()
-    ds_bgc = boundary_forcing_with_bgc_from_climatology.ds["bgc"]
-
-    # Check the values in the dataset
-    assert np.allclose(ds["zeta_south"].values, expected_zeta_south)
-    assert np.allclose(ds["zeta_east"].values, expected_zeta_east)
-    assert np.allclose(ds["zeta_north"].values, expected_zeta_north)
-    assert np.allclose(ds["zeta_west"].values, expected_zeta_west)
-    assert np.allclose(ds["temp_south"].values, expected_temp_south)
-    assert np.allclose(ds["temp_east"].values, expected_temp_east)
-    assert np.allclose(ds["temp_north"].values, expected_temp_north)
-    assert np.allclose(ds["temp_west"].values, expected_temp_west)
-    assert np.allclose(ds["u_south"].values, expected_u_south)
-    assert np.allclose(ds["u_east"].values, expected_u_east)
-    assert np.allclose(ds["u_north"].values, expected_u_north)
-    assert np.allclose(ds["u_west"].values, expected_u_west)
-    assert np.allclose(ds["v_south"].values, expected_v_south)
-    assert np.allclose(ds["v_east"].values, expected_v_east)
-    assert np.allclose(ds["v_north"].values, expected_v_north)
-    assert np.allclose(ds["v_west"].values, expected_v_west)
-    assert np.allclose(ds["ubar_south"].values, expected_ubar_south)
-    assert np.allclose(ds["ubar_east"].values, expected_ubar_east)
-    assert np.allclose(ds["ubar_north"].values, expected_ubar_north)
-    assert np.allclose(ds["ubar_west"].values, expected_ubar_west)
-    assert np.allclose(ds["vbar_south"].values, expected_vbar_south)
-    assert np.allclose(ds["vbar_east"].values, expected_vbar_east)
-    assert np.allclose(ds["vbar_north"].values, expected_vbar_north)
-    assert np.allclose(ds["vbar_west"].values, expected_vbar_west)
-    assert np.allclose(ds_bgc["ALK_south"].values, expected_alk_south)
-    assert np.allclose(ds_bgc["ALK_east"].values, expected_alk_east)
-    assert np.allclose(ds_bgc["ALK_north"].values, expected_alk_north)
-    assert np.allclose(ds_bgc["ALK_west"].values, expected_alk_west)
-
-    boundary_forcing_with_bgc_from_climatology.plot(
-        varname="temp_south", layer_contours=True
-    )
-    boundary_forcing_with_bgc_from_climatology.plot(
-        varname="temp_east", layer_contours=True
-    )
-    boundary_forcing_with_bgc_from_climatology.plot(
-        varname="temp_north", layer_contours=True
-    )
-    boundary_forcing_with_bgc_from_climatology.plot(
-        varname="temp_west", layer_contours=True
-    )
-    boundary_forcing_with_bgc_from_climatology.plot(varname="zeta_south")
-    boundary_forcing_with_bgc_from_climatology.plot(varname="zeta_east")
-    boundary_forcing_with_bgc_from_climatology.plot(varname="zeta_north")
-    boundary_forcing_with_bgc_from_climatology.plot(varname="zeta_west")
-    boundary_forcing_with_bgc_from_climatology.plot(varname="vbar_north")
-    boundary_forcing_with_bgc_from_climatology.plot(varname="ubar_west")
-    boundary_forcing_with_bgc_from_climatology.plot(varname="ALK_south")
-    boundary_forcing_with_bgc_from_climatology.plot(varname="ALK_east")
-    boundary_forcing_with_bgc_from_climatology.plot(varname="ALK_north")
-    boundary_forcing_with_bgc_from_climatology.plot(varname="ALK_west")
-
-    # Create a temporary file
-    with tempfile.NamedTemporaryFile(delete=True) as tmpfile:
-        filepath = tmpfile.name
-
-    boundary_forcing_with_bgc_from_climatology.save(filepath)
-    physics_filepath = filepath + "_physics_20210629-29.nc"
-    bgc_filepath = filepath + "_bgc_clim.nc"
-
-    try:
-        assert os.path.exists(physics_filepath)
-        assert os.path.exists(bgc_filepath)
-    finally:
-        os.remove(physics_filepath)
-        os.remove(bgc_filepath)
-
-
-@pytest.mark.parametrize(
-    "bdry_forcing_fixture",
-    [
-        "boundary_forcing",
-        "boundary_forcing_with_bgc_from_climatology",
-    ],
-)
-def test_roundtrip_yaml(bdry_forcing_fixture, request):
-    """Test that creating a BoundaryForcing object, saving its parameters to yaml file, and re-opening yaml file creates the same object."""
-
-    bdry_forcing = request.getfixturevalue(bdry_forcing_fixture)
-
-    # Create a temporary file
-    with tempfile.NamedTemporaryFile(delete=False) as tmpfile:
-        filepath = tmpfile.name
-
-    try:
-        bdry_forcing.to_yaml(filepath)
-
-        boundary_forcing_from_file = BoundaryForcing.from_yaml(filepath)
-
-        assert bdry_forcing == boundary_forcing_from_file
-
-    finally:
-        os.remove(filepath)
-
-
-def test_from_yaml_missing_boundary_forcing():
-    yaml_content = textwrap.dedent(
-        """\
-    ---
-    roms_tools_version: 0.0.0
-    ---
-    Grid:
-      nx: 100
-      ny: 100
-      size_x: 1800
-      size_y: 2400
-      center_lon: -10
-      center_lat: 61
-      rot: -20
-      topography_source: ETOPO5
-      smooth_factor: 8
-      hmin: 5.0
-      rmax: 0.2
-    """
-    )
-
-    with tempfile.NamedTemporaryFile(delete=False) as tmp_file:
-        yaml_filepath = tmp_file.name
-        tmp_file.write(yaml_content.encode())
-
-    try:
-        with pytest.raises(
-            ValueError, match="No BoundaryForcing configuration found in the YAML file."
-        ):
-            BoundaryForcing.from_yaml(yaml_filepath)
-    finally:
-        os.remove(yaml_filepath)