maubin/sbmy_control
0
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

parameter card

Browse source
This commit is contained in:
Mayeul Aubin 2025-06-24 11:47:24 +02:00
parent 5eefe1061c
commit 3e93c3fec3
1 changed files with 29 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

29
sbmy_control/parameters_card.py
View file

@ -452,6 +452,35 @@ def create_parameter_card_dict(
MeshDensity=MeshDensity,
LoadPhiBCs=int(LoadPhiBCs),
InputPhiBCsBase=InputPhiBCsBase,
# P3M default parameters
n_Tiles = 60,
# Number of tiles in each dimension for the P3M code
nPairsForceDiagnostic =0,
# Number of particle-hole pairs to be used for the force diagnostic
nBinsForceDiagnostic =0,
maxTrialsForceDiagnostic =0,
OutputForceDiagnostic ='',
PrintOutputTimestepsLog = 0,
OutputTimestepsLog = "",
fac_dyn_custom = 0.03,
da_max_early_custom = 0.0013,
fac_H_custom = 0.05,
fac_bend = 0.3125,
sub_bend1 = 0.012,
sub_bend2 = 0.014,
fac_p3m_fit = 0.15,
da_max_late_custom = 0.02,
# Polynomial coefficients for the P3M time step limiter. Only used if use_p3m_fit=1
use_p3m_fit = 1,
p3m_fit_coeff0 = -3.19544528,
p3m_fit_coeff1 = -0.01948223,
p3m_fit_coeff2 = 0.5180241,
p3m_fit_coeff3 = -1.08372045,
p3m_fit_coeff4 = -0.71808255,
p3m_fit_coeff5 = -0.14370312,
RunForceDiagnostic = 0,
h=h,
Omega_m=Omega_m,