Hi there,
I found something I'm confused about in the function set_T_rho_args(T, rho, T_rho_type_id, T_rho_args, mat_id) in woma/eos/T_rho.py. In the T_rho.py file, line 119 is used to assign T_rho_args with a fixed entropy, while it recalculates the surface entropy here. For SESAME and ANEOS, this would be adiabatic instead of fix entropy. The entropy of placed particles would be different from what you set initially in the planet profile if you use random surface T and P.
Hi there,
I found something I'm confused about in the function
set_T_rho_args(T, rho, T_rho_type_id, T_rho_args, mat_id)inwoma/eos/T_rho.py. In the T_rho.py file, line 119 is used to assignT_rho_argswith a fixed entropy, while it recalculates the surface entropy here. For SESAME and ANEOS, this would be adiabatic instead of fix entropy. The entropy of placed particles would be different from what you set initially in theplanetprofile if you use random surface T and P.