New test problem LSS_SHR

[vivi235711]

This pull request adds a new test problem named LSS_SHR. It builds upon the existing LSS_Hybrid problem and provides Python scripts for analyzing the soliton-halo relationship.

• Analysis Scripts: The included scripts reside in the plot_script subdirectory with README.md.
• Paper Reproduction: The provided initial condition for LSS_SHR can reproduce the halo in this paper.
• Scripts related to this paper's figures are located at plot_script/script_for_paper_fig.
• SHR Prediction Tool Integration: LSS_SHR utilizes the shr prediction tool managed through Git modules.

[hyschive]

@vivi235711 Thanks for the contribution! Please check my comment regarding a potential incompatibility issue with the latest SHR module.

[hyschive]

• Why is redshift fixed to 0?
• This line seems incompatible with the latest SHR modul. Note that revised_theo_c_FDM_Ms returns 8 quantities (ms, rs, peak_dens, Rh, NFW_Rs, c_theo, c_CDM, beta), whereas here we only need the core mass.
  • I think this issue is caused by the fact that git submodule update --init points to the commit hash 57832b41de1f4c0284ca929b52ee1712d7af68f3 rather than the latest commit d6de09ca7ca2449c3614d135fb8adb96ec0bbb46. Please update this PR to use the latest commit.
  • If this is indeed an incompatibility issue, please update all scripts in script_for_paper_fig/ as well.

Note that the second point above leads to the following error:

Command-line arguments:
-------------------------------------------------------------------
plot__all.py -s 30 -e 37 -halo 1
-------------------------------------------------------------------
Traceback (most recent call last):
  File "/projectV/fish/gamer_test/t21-pr-lss-shr/gamer.new/bin/run/plot_script/plot__all.py", line 92, in <module>
    plt.plot(time_a,(core_mass/shr_calculator.revised_theo_c_FDM_Ms(0,halo_mass,m22)),'-.', color = 'black', label = 'ms sim/theo')
  File "/usr/local/lib/python3.10/dist-packages/pandas/core/ops/common.py", line 81, in new_method
    return method(self, other)
  File "/usr/local/lib/python3.10/dist-packages/pandas/core/arraylike.py", line 210, in __truediv__
    return self._arith_method(other, operator.truediv)
  File "/usr/local/lib/python3.10/dist-packages/pandas/core/series.py", line 6113, in _arith_method
    return base.IndexOpsMixin._arith_method(self, other, op)
  File "/usr/local/lib/python3.10/dist-packages/pandas/core/base.py", line 1350, in _arith_method
    return self._construct_result(result, name=res_name)
  File "/usr/local/lib/python3.10/dist-packages/pandas/core/series.py", line 3106, in _construct_result
    out = self._constructor(result, index=self.index, dtype=dtype)
  File "/usr/local/lib/python3.10/dist-packages/pandas/core/series.py", line 509, in __init__
    data = sanitize_array(data, index, dtype, copy)
  File "/usr/local/lib/python3.10/dist-packages/pandas/core/construction.py", line 607, in sanitize_array
    subarr = _sanitize_ndim(subarr, data, dtype, index, allow_2d=allow_2d)
  File "/usr/local/lib/python3.10/dist-packages/pandas/core/construction.py", line 666, in _sanitize_ndim
    raise ValueError(
ValueError: Data must be 1-dimensional, got ndarray of shape (8, 8) instead

[hyschive]

• Add a note in README.md describing how to set-up the SHR module.
• Also check my comment in plot__all.py about an incompatibility issue.

[hyschive]

Suggested change

	ratio_shcive = []
	ratio_schive = []