fluidsim-ipy-load doesn't work correctly

Some issues appeared with fluidsim-ipy-load:

• for usual simulations, plots are created only when trying to close with Ctrl+D
• for simulations from restart, can't load files, give this message:

ValueError                                Traceback (most recent call last)
Cell In[1], line 1
----> 1 import numpy as np; import matplotlib.pyplot as plt; import pandas as pd; from fluidsim import load; print("Loading simulation"); sim = load(); params = sim.params; print("`sim`, `params`, `np`, `plt` and `pd` variables are available")

File /.fsdyn_people/agoua5nn/project/24STRATURBKOLMO/fluidsim/fluidsim/util/util.py:246, in load_sim_for_plot(name_dir, merge_missing_params, hide_stdout)
    243     Simul = solver.Simul
    245 with stdout_redirected(hide_stdout):
--> 246     sim = Simul(params)
    247 return sim

File /.fsdyn_people/agoua5nn/project/24STRATURBKOLMO/fluidsim/fluidsim/base/solvers/base.py:129, in SimulBase.__init__(self, params)
    127 # initialization operators and grid
    128 Operators = dict_classes["Operators"]
--> 129 self.oper = Operators(params=params)
    131 # initialization output
    132 Output = dict_classes["Output"]

File /.fsdyn_people/agoua5nn/project/24STRATURBKOLMO/fluidsim/fluidsim/operators/operators3d.py:222, in OperatorsPseudoSpectral3D.__init__(self, params)
    219     ny = params.oper.ny
    220     nz = params.oper.nz
--> 222 super().__init__(
    223     nx,
    224     ny,
    225     nz,
    226     params.oper.Lx,
    227     params.oper.Ly,
    228     params.oper.Lz,
    229     fft=params.oper.type_fft,
    230     coef_dealiasing=params.oper.coef_dealiasing,
    231 )
    233 # compatibility for fluidfft <= 0.3.0
    234 if not hasattr(self, "oper_fft") and hasattr(self, "_op_fft"):

File /.fsdyn_people/agoua5nn/project/24STRATURBKOLMO/fluidsim/.venv/lib/python3.9/site-packages/fluidfft/fft3d/operators.py:219, in OperatorsPseudoSpectral3D.__init__(self, nx, ny, nz, lx, ly, lz, fft, coef_dealiasing)
    216 self._is_mpi_lib = self.shapeX_seq != self.shapeX_loc and mpi.nb_proc > 1
    218 Lx = self.Lx = float(lx)
--> 219 Ly = self.Ly = float(ly)
    220 Lz = self.Lz = float(lz)
    222 self.deltax = Lx / nx

ValueError: could not convert string to float: 'np.float64(3.0)'

• doesn't plot figures on gricad clusters.