def get_accretion_rate(nout=0): data = visu_ramses.load_snapshot(nout) yr = 86400.0 * 365.25 time = data["info"]["time"] * data["info"]["unit_t"] / yr / 1000.0 u_m = data["info"]["unit_d"] * data["info"]["unit_l"]**3 msun = 1.98847e33 rate = data["sinks"]["acc_rate"] * u_m * yr / msun return time, rate
import matplotlib.pyplot as plt import visu_ramses fig = plt.figure() ratio = 0.8 sizex = 12.0 fig.set_size_inches(sizex, ratio * sizex) ax1 = plt.subplot(221) ax2 = plt.subplot(222) ax3 = plt.subplot(223) ax4 = plt.subplot(224) ax5 = ax4.twinx() # Load RAMSES output data = visu_ramses.load_snapshot(2) time = data["time"] ngr = data["ngrp"] ar = 7.56591469318689378e-015 xshift = 1013.67 order = data["x"].argsort() x = data["x"][order] - xshift amrlev = data["level"][order] rho = data["density"][order] u = data["velocity_x"][order] p = data["thermal_pressure"][order] T = data["temperature"][order] Tr = (data["radiative_energy_1"][order] / ar)**0.25 # Read analytical solution