Dev

[brey]

No description provided.

[pmav99]

Not a biggie, but would it make sense to move this log message inside execute_schism_mpirun_script()?

I guess the answer depends on whether the function is being used elsewhere

[brey]

I guess, if we put it in tools the logging reference will also be different.

pyposeidon.schism run:671 executing model
pyposeidon.tools :xxx model finished successfully

which might be misleading as well, but I am open to it if you think it's better.

[pmav99]

With the first 6 days of May 2023 we get an exception:

  5%|█████████▋                                                                                                                                                                                       | 48/961 [00:00<00:04, 187.58it/s]
---------------------------------------------------------------------------
IndexError                                Traceback (most recent call last)
File /mnt/to_thalassa/run.py:50
     47 end_time = pd.Timestamp("2023-05-06T00:00:00")
     48 stations = gpd.read_file("base_model/stations.json")
---> 50 ods = compute_obs(stations=stations, start_time=start_time, end_time=end_time)
     51 save_ods(ods)
     52 1/0

File /mnt/micromamba/envs/pos_dev/lib/python3.10/site-packages/pyposeidon/utils/post.py:55, in compute_obs(stations, start_time, end_time)
     52     if oi[var].isnull().all().values == True:
     53         oi = oi.drop(var)
---> 55 var = [k for k, v in oi.data_vars.items() if v.dims == ("time",)][0]
     57 obs = oi[var].to_dataframe().drop(["location"], axis=1)  # Get observational data
     59 # de-tide obs
     60 #        if not obs[var].dropna().empty | (obs[var].dropna() == obs[var].dropna()[0]).all():
     61 #            obs = get_ss(obs, oi.lat.values)
     62 #            oi[var].values = obs.elev.values

IndexError: list index out of range