Star formation sub grid update

I changed the temperature criteria to an entropy criteria to be faster, furthermore I added the birth temperature to the snapshots and added a histogram plot to the solutions.

So I run tests using the lowres8 isolated galaxy without feedback, and I don't seem to see any difference in the results: Old: New:

Same for the SFH: Old: New:

I will run the one with feedback to check if it gives a similar result.

One thing, I assumed that in general we want to keep the functionality to call the cooling function in star_formation_is_star_forming function for other subgrid models, so that function still has the required arguments.

examples/IsolatedGalaxy/IsolatedGalaxy_starformation/plotSolution.py

@@ -123,6 +123,7 @@ stars_pos = f["/PartType4/Coordinates"][:, :]
 stars_BirthDensity = f["/PartType4/BirthDensity"][:]
 stars_BirthTime = f["/PartType4/BirthTime"][:]
 stars_XH = f["/PartType4/ElementAbundance"][:,0]
+stars_BirthTemperature = f["/PartType4/BirthTemperature"][:]
 
 # Centre the box
 gas_pos[:, 0] -= centre[0]
@@ -207,9 +208,17 @@ figure()
 subplot(111, xscale="linear", yscale="linear")
 hist(np.log10(stars_BirthDensity),density=True,bins=20,range=[-2,5])
 xlabel("${\\rm Stellar~birth~density}~n_{\\rm H}~[{\\rm cm^{-3}}]$", labelpad=0)
-ylabel("${\\rm Probability}$", labelpad=-7)
+ylabel("${\\rm Probability}$", labelpad=3)
 savefig("BirthDensity.png", dpi=200)
 
+# Histogram of the birth temperature 
+figure()
+subplot(111, xscale="linear", yscale="linear")
+hist(np.log10(stars_BirthTemperature),density=True,bins=20,range=[3.5,5.0])
+xlabel("${\\rm Stellar~birth~temperature}~[{\\rm K}]$", labelpad=0)
+ylabel("${\\rm Probability}$", labelpad=3)
+savefig("BirthTemperature.png", dpi=200)
+
 # Plot of the specific star formation rate in the galaxy
 rhos = 10**np.linspace(-1,np.log10(KS_high_den_thresh),100)
 rhoshigh = 10**np.linspace(np.log10(KS_high_den_thresh),5,100)