feat: Update make_particles.py to include disk mass calculation

This commit modifies the make_particles.py script to include the calculation of the disk mass. The disk mass is now calculated based on the sum of the masses of the particles in the output data. This ensures that the disk mass is accurately represented in the simulation.

Author: Gianuzzi

tools/make_particles.py

@@ -81,12 +81,15 @@ def norm_dist(x0=0.3, scale=0.01, xmin=0, xmax=1.0):
 data_in = {}
 names = ["mass", "a", "e", "M", "w", "R"]
 # PARÁMETROS A VARIAR (Poner unidades de ser necesario)
-data_in["mass"] = [rndm(12, 15, lambda x: 10 ** x)]
+data_in["mass"] = [rndm(0.0, 1.)]
 data_in["a"] = [0.0]
 data_in["e"] = [0.0]  # Podría ser: rayleigh_dist(0, 0.1, 1)
 data_in["M"] = [rndm(0.0, 360.0)]
 data_in["w"] = [0.0]
-data_in["R"] = [n_steps(0.5, 3.5, 4)]
+data_in["R"] = [n_steps(0.5, 3.5, 4000)]
+
+# Disk mass (in kg)
+disk_mass = 6.3e15
 
 # -----------------
 # OUTPUT
@@ -228,6 +231,9 @@ def shuffle_matrix(matrix):
     # Remove mass if not needed
     if not use_mass:
         data_out = data_out[:, 1:]
+    elif disk_mass > 0:
+        temp = np.sum(data_out[:, 0])
+        data_out[:, 0] = disk_mass * data_out[:, 0] / temp
 
     # Guardamos
     if not check_continue(filename):