Estimate disk scale height, disk radius, and disk mass using plons
following the method described in Sect. 4.4. of Malfait et al. 2024a (https://ui.adsabs.harvard.edu/abs/2024arXiv240813158M/abstract)
Short overview of the method:
To estimate the disk quantities, we integrate over an array of radii (r_array), with step r_step (set to 0.01 au).
For each radius r_i in this array, we estimate the disk scale height as the median of the scale heights at certain coordinates (r_i,theta_j)
Then using the values of r_i and the median scale height, the added mass w.r.t. to total disk mass of the previous integration step is calculated.
The radius is then defined as the radius where this added mass w.r.t. the total mass is significantly small (see Eq. (9) of Malfait+ 2024a)
Because the disks are asymmetric, we can estimate these quantities in function of radius in different quadrants, to quantify the assymetry of the disk
Therefore, this script contains one function were calculations are performed over the full theta range (fullThetaCalculations()), and one over 4 different theta quadrants (thetaRegionsCalculations())
Importing packages
[1]:
%load_ext autoreload
%autoreload 2
[2]:
import plons
import os
import numpy as np
import matplotlib.pyplot as plt
# import necessary plons scripts
import plons.ConversionFactors_cgs as cgs
import plons.GeometricalFunctions as gf
import plons.AccrDisk as ad
Function to calculate scale height and disk mass ifo radius, using the median scale height in all theta directions
[3]:
'''
Input:
dumpData input data, used for density and position
run modelName, used to save in correct locations
dump dumpNumber, used to save under correct name
xH calculate quantities at first or second scale height (see Eq.7 and 8 in Malfait 2024),
'2H' used for disk radius and mass estimates
Output:
r_array array of radii starting at 'lowerR', ending at estimated disk radius, and increasing with step 'r_step'
SH_array array of scale height estimates at each radius
RhoMax_array array of maximum density value, which is probably at orbital plane, at each radius
totalMassDisk array of total disk mass assuming that disk radius is r_i of corresponding r_array
rel_mass_added array of mass added in this r_step, divided by total mass at previous r_i
Note that the surface density and optical depth values have not been in used in our paper, mainly because it is a too approximate approach
Sigma_array array of estimated surface density Sigma
SigmaTheor_array array of theoretical surface density Sigma
Tau_array array of estimated optical depth tau
TauTheor_array array of theoretical optical depth tau
! The output data is saved in a txt file, and can be used to produce plots using the accrDisks_plotsHRM_Plons.ipynb notebook
'''
def fullThetaCalculations(dumpData,run,dump,xH, printOut = False):
# Choose which and how many r and theta values you want to use to calculate scale heights and disk mass and for plot of rho vs h
n_grid = 10000
maxH = 1.5 *cgs.au
numberOfThetas = 10
thetaMin = 0
thetaMax = 2*np.pi
thetaArray = np.linspace(thetaMin,(numberOfThetas-1)/numberOfThetas*thetaMax,numberOfThetas)
crit = 0.3 #smalles constant in our case towards which the relative added mass ocnsistently keeps decreasing for all models
lowerR = 0.02
r_step = 0.01
# Calculate arrays with radii, disk mass, relative addes mass, and scale heights
r_array,totalMassDisk,rel_mass_added,SH_array,RhoMax_array,Sigma_array,SigmaTheor_array,Tau_array,TauTheor_array = ad.get_SH_diskMass_radius(lowerR,r_step,thetaArray,dumpData,maxH,n_grid,run,crit,xH,phiQuadr = False, printOut = printOut)
'''
Write data for full theta region in txt file infoAccrDisk_wind_00....txt
'''
testLimit = crit
# make directory to save results in
wdir = os.path.join(run,'plotsAnalysis')
!mkdir -p $wdir
# construct txt file with calculated arrays of data
title = os.path.join(run,'plotsAnalysis/infoAccrDisk_wind_00'+str(dump)+'_'+xH+'.txt')
with open (title,'w') as f:
f.write('Model '+str(run)+'wind_00'+str(dump)+'\n')
f.write('Data analysis accretion disks')
f.write('\n')
f.write('\n')
f.write('for rstep '+ str(r_step)+ ' au (=racc) and limit of '+ str(testLimit)+':'+'\n')
f.write('the estimated radius is: '+ str(np.round(r_array[np.where(rel_mass_added/r_step<testLimit)][0]/cgs.au,3))+' au'+'\n')
f.write('the estimated total mass is: '+ str(np.round(totalMassDisk[np.where(rel_mass_added/r_step<testLimit)][0]/cgs.Msun,11))+' Msun'+'\n')
f.write('the estimated max scale height is: '+ str(np.round(SH_array[np.where(rel_mass_added/r_step<testLimit)][0]/cgs.au,3))+' au'+'\n')
f.write('\n')
f.write('\n')
f.write('for rstep '+ str(r_step)+ ' au (=racc):'+'\n')
names = ['r [au]', 'SH [au]', 'Rho_max [g/cm^3]','Mtot [Msun]', 'Mrel/rstep []', 'Sigma []','Sigma Theor []','OptDepth []','OptDepthTheor []']
f.write("{: <30} {: <30} {: <30} {: <30} {: <30} {: <30} {: <30} {: <30} {: <30}".format(*names))
col_format = "{:<31}" * 9 + "\n" # 5 left-justfied columns with 15 character width
f.write('\n')
for i in zip(np.round(r_array/cgs.au,2),np.round(SH_array/cgs.au,4),np.round(RhoMax_array,16),np.round(totalMassDisk/cgs.Msun,13),np.round(rel_mass_added/r_step,3),np.round(Sigma_array,4),np.round(SigmaTheor_array,4),np.round(Tau_array,3),np.round(TauTheor_array,3)):
f.write(col_format.format(*i))
'''
If you want, you can immediately construct some plots here:
Make plots of relative added mass, total disk mass, midplane rho, and scale heights ifo r
'''
if xH == '1H':
fig, axs = plt.subplots(nrows = 1, ncols= 6 , figsize=(35, 5))
else:
fig, axs = plt.subplots(nrows = 1, ncols= 3 , figsize=(14, 5))
ad.plot_relAddedMass(axs[0], r_array,rel_mass_added,r_step,crit)
ad.plot_totalDiskMass(axs[1], r_array,totalMassDisk)
ad.plot_scaleHeights(axs[2], r_array,SH_array)
if xH == '1H':
ad.plot_surfDensity(axs[3],r_array,Sigma_array,SigmaTheor_array)
ad.plot_opticalDepth(axs[4],r_array,Tau_array,TauTheor_array)
ad.plot_rhoMax(axs[5],r_array,RhoMax_array)
return fig
Function to calculate scale height and disk mass ifo radius, using the median scale height in different theta quadrants
[4]:
'''
Same function as before, but in 4 different theta quadrants
Saves data in 4 seperate txt files
input
dumpData input data, used for density and position
run modelName, used to save in correct locations
dump dumpNumber, used to save under correct name
xH calculate quantities at first or second scale height (see Eq.7 and 8 in Malfait 2024)
'2H' used for disk radius and mass estimates
'''
def thetaRegionsCalculations(dumpData,run,dump,xH):
'''
Calculations of scale height, mass estimates, radius
For 4 different theta regions
'''
# Choose which and how many r and theta values you want to use to calculate scale heights and disk mass and for plot of rho vs h
n_grid = 10000
# maxH = 1.5 * cgs.au
maxH = 3 * cgs.au
numberOfThetas = 10 #Number of thetas withing one array
'''
Divide into 4 different theta regions
'''
# Away from AGB
thetaMin1 = -1./4. *np.pi
thetaMax1 = 1./4. *np.pi
thetaArray1 = np.linspace(thetaMin1,thetaMax1,numberOfThetas)
# Make sure all theta's are positive
i = 0
for theta in thetaArray1:
if theta < 0:
thetaArray1[i] = thetaArray1[i]+2*np.pi
i = i+1
# Positive y --> bow shock side (in front of companion)
thetaMin2 = 1./4. * np.pi
thetaMax2 = 3./4. * np.pi
thetaArray2 = np.linspace(thetaMin2,thetaMax2,numberOfThetas)
# Towards AGB
thetaMin3 = 3./4. * np.pi
thetaMax3 = 5./4. * np.pi
thetaArray3 = np.linspace(thetaMin3,thetaMax3,numberOfThetas)
# Negative y --> flow behind companion
thetaMin4 = 5./4. * np.pi
thetaMax4 = 7./4. * np.pi
thetaArray4 = np.linspace(thetaMin4,thetaMax4,numberOfThetas)
crit = 0.3
lowerR = 0.02
# sometimes you will see that you need to vary the step a little bit to make sure the algorithm does not get stuck in the first radii steps
r_step1 = 0.01
r_step2 = r_step1
r_step3 = r_step1
r_step4 = r_step1
'''
Calculate scale heights and mass estimates at each r; and estimate final radius
Make text files with important info to make plots
for different theta regimes
'''
testLimit = crit
# calculate scale height and mass estimates at each r for this setup
r_array1,totalMassDisk1,rel_mass_added1,SH_array1,RhoMax_array1,Sigma_array1,SigmaTheor_array1,Tau_array1,TauTheor_array1 = ad.get_SH_diskMass_radius(lowerR,r_step1,thetaArray1,dumpData,maxH,n_grid,run,crit,xH,phiQuadr=True,printOut=True)
title = os.path.join(run,'plotsAnalysis/infoAccrDisk_theta~0_wind_00'+str(dump)+'_'+xH+'.txt')
ad.writeFile(title,r_step1,testLimit,r_array1,rel_mass_added1,totalMassDisk1,SH_array1,RhoMax_array1,Sigma_array1,SigmaTheor_array1,Tau_array1,TauTheor_array1)
r_array2,totalMassDisk2,rel_mass_added2,SH_array2,RhoMax_array2,Sigma_array2,SigmaTheor_array2,Tau_array2,TauTheor_array2 = ad.get_SH_diskMass_radius(lowerR,r_step2,thetaArray2,dumpData,maxH,n_grid,run,crit,xH,phiQuadr=True,printOut=True)
title = os.path.join(run,'plotsAnalysis/infoAccrDisk_theta~pi:2_wind_00'+str(dump)+'_'+xH+'.txt')
ad.writeFile(title,r_step2,testLimit,r_array2,rel_mass_added2,totalMassDisk2,SH_array2,RhoMax_array2,Sigma_array2,SigmaTheor_array2,Tau_array2,TauTheor_array2)
r_array3,totalMassDisk3,rel_mass_added3,SH_array3,RhoMax_array3,Sigma_array3,SigmaTheor_array3,Tau_array3,TauTheor_array3 = ad.get_SH_diskMass_radius(lowerR,r_step3,thetaArray3,dumpData,maxH,n_grid,run,crit,xH,phiQuadr=True,printOut=True)
title = os.path.join(run,'plotsAnalysis/infoAccrDisk_theta~pi_wind_00'+str(dump)+'_'+xH+'.txt')
ad.writeFile(title,r_step3,testLimit,r_array3,rel_mass_added3,totalMassDisk3,SH_array3,RhoMax_array3,Sigma_array3,SigmaTheor_array3,Tau_array3,TauTheor_array3)
r_array4,totalMassDisk4,rel_mass_added4,SH_array4,RhoMax_array4,Sigma_array4,SigmaTheor_array4,Tau_array4,TauTheor_array4 = ad.get_SH_diskMass_radius(lowerR,r_step1,thetaArray4,dumpData,maxH,n_grid,run,crit,xH,phiQuadr=True,printOut=True)
title = os.path.join(run,'plotsAnalysis/infoAccrDisk_theta~pi3:2_wind_00'+str(dump)+'_'+xH+'.txt')
ad.writeFile(title,r_step4,testLimit,r_array4,rel_mass_added4,totalMassDisk4,SH_array4,RhoMax_array4,Sigma_array4,SigmaTheor_array4,Tau_array4,TauTheor_array4)
'''
If you want you can already construct some plots to see the results
'''
if xH == '1H':
fig, axs = plt.subplots(nrows = 1, ncols= 6 , figsize=(35, 5))
else:
fig, axs = plt.subplots(nrows = 1, ncols= 3 , figsize=(14, 5))
ad.plot_relAddedMass_4theta(axs[0],r_array1,r_array2,r_array3,r_array4,rel_mass_added1,rel_mass_added2,rel_mass_added3,rel_mass_added4,r_step1,r_step3,crit)
ad.plot_totalDiskMass_4theta(axs[1],r_array1,r_array2,r_array3,r_array4,totalMassDisk1,totalMassDisk2,totalMassDisk3,totalMassDisk4)
ad.plot_scaleHeights_4theta(axs[2],r_array1,r_array2,r_array3,r_array4,SH_array1,SH_array2,SH_array3,SH_array4)
if xH == '1H':
ad.plot_surfDensity_4theta(axs[3],r_array1,r_array2,r_array3,r_array4,Sigma_array1,Sigma_array2,Sigma_array3,Sigma_array4)
ad.plot_opticalDepth_4theta(axs[4],r_array1,r_array2,r_array3,r_array4,Tau_array1,Tau_array2,Tau_array3,Tau_array4)
ad.plot_rhoMax_4theta(axs[5], r_array1,r_array2,r_array3,r_array4,RhoMax_array1,RhoMax_array2,RhoMax_array3,RhoMax_array4)
return fig
Main function for the actual calculations
[5]:
'''
input:
run modelName, used to save in correct locations
dump dumpNumber, used to save under correct name
scaleHeightList list with 1H, 2H or one of both, to know for which scale height calculations need to be done
full, quadrants boolians to know if calculations should be done for full thetaregion and/or quadrants
'''
def main(run,dump,scaleHeightList,full=True,quadrants=True,printOut = False):
'''
Load in data
'''
setup = plons.LoadSetup(run, "wind")
dumpData = plons.LoadFullDump(os.path.join(run, f"wind_%05d" % dump), setup)
# Coordinate transformation: translation to rcomp + rotation such that AGB is on positive x-axis
# Coordinate transformation such that comp is at (0,0) and AGB at -x
dumpData['new_x'],dumpData['new_y'],dumpData['new_r'] = ad.calc_new_position(dumpData['position'][:,0],dumpData['position'][:,1],dumpData)
phi = gf.calcPhi(dumpData['new_x'],dumpData['new_y'])
# print(phi)
dumpData['new_Phi'] = np.where(phi < 0, phi + 2 * np.pi, phi)
'''
You can use this to test if coordinate transformation is correct
'''
# testPositionAndTheta(dumpData,run,dump)
'''
Perform calculations for full theta region and/or 4 different theta regions
make plots and write txt files with results
'''
for xH in scaleHeightList:
if full==True:
if printOut: print('Calculations full theta region:')
fig = fullThetaCalculations(dumpData,run,dump,xH, printOut)
fig.savefig(run+'plotsAnalysis/Disk_wind_%05d'%(dump)+'_'+xH+'.png', bbox_inches='tight')
if quadrants==True:
if printOut: print('Calculations for different theta regions:')
fig = thetaRegionsCalculations(dumpData,run,dump,xH)
fig.savefig(run+f'plotsAnalysis/Theta_rstep_wind_%05d'%(dump)+'_'+xH+'.png', bbox_inches='tight')
Examples
[6]:
'''
Examples, 2 models, different dumps for the eccentric, different quadrants
for the circular one
'''
model = 'model-v10e50_AD/'
!mkdir -p $model
infile = os.path.join(model, f"wind.in")
if not os.path.exists(infile):
!wget -q "https://raw.githubusercontent.com/Ensor-code/phantom-models/refs/heads/main/Malfait%2B2024a/v10e50_AD/wind.in" --output-document $infile
setupfile = os.path.join(model, f"wind.setup")
if not os.path.exists(setupfile):
!wget -q "https://raw.githubusercontent.com/Ensor-code/phantom-models/refs/heads/main/Malfait%2B2024a/v10e50_AD/wind.setup" --output-document $setupfile
dumps = [555,584]
for dumpnumber in dumps:
dump = os.path.join(model, f"wind_{dumpnumber:05d}")
if not os.path.exists(dump):
wgetlink = f"https://github.com/Ensor-code/phantom-models/raw/refs/heads/main/Malfait+2024a/v10e50_AD/wind_{dumpnumber:05d}?download="
!wget -q $wgetlink --output-document $dump
main(model,dumpnumber,['2H'],full=True,quadrants=False,printOut=True)
model = 'model-v10e00_AD/'
!mkdir -p $model
infile = os.path.join(model, f"wind.in")
if not os.path.exists(infile):
!wget -q "https://raw.githubusercontent.com/Ensor-code/phantom-models/refs/heads/main/Malfait%2B2024a/v10e00_AD/wind.in" --output-document $infile
setupfile = os.path.join(model, f"wind.setup")
if not os.path.exists(setupfile):
!wget -q "https://raw.githubusercontent.com/Ensor-code/phantom-models/refs/heads/main/Malfait%2B2024a/v10e00_AD/wind.setup" --output-document $setupfile
dumpnumber = 584
dump = os.path.join(model, f"wind_{dumpnumber:05d}")
if not os.path.exists(dump):
wgetlink = f"https://github.com/Ensor-code/phantom-models/raw/refs/heads/main/Malfait+2024a/v10e00_AD/wind_{dumpnumber:05d}?download="
!wget -q $wgetlink --output-document $dump
main(model,dumpnumber,['2H'],full=True,quadrants=True,printOut=True)
Calculations full theta region:
---------
r 0.02 au
scale height is 0.011401140114011388
mass added [Msun]: 3.4540161426105943e-11
relMass_added/rstep (rico) = 100.0
---------
r 0.03 au
/home/matse/codes/anaconda3/envs/plons/lib/python3.10/site-packages/numpy/core/fromnumeric.py:3504: RuntimeWarning: Mean of empty slice.
return _methods._mean(a, axis=axis, dtype=dtype,
/home/matse/codes/anaconda3/envs/plons/lib/python3.10/site-packages/numpy/core/_methods.py:129: RuntimeWarning: invalid value encountered in scalar divide
ret = ret.dtype.type(ret / rcount)
scale height is 0.009450945094509448
mass added [Msun]: 1.554307264174767e-10
relMass_added/rstep (rico) = 81.818
---------
r 0.04 au
scale height is 0.009000900090009034
mass added [Msun]: 5.200213192485949e-10
relMass_added/rstep (rico) = 73.243
---------
r 0.05 au
scale height is 0.008025802580258031
mass added [Msun]: 8.980441970787545e-10
relMass_added/rstep (rico) = 55.847
---------
r 0.06 au
scale height is 0.008475847584758509
mass added [Msun]: 1.3604985806171728e-09
relMass_added/rstep (rico) = 45.831
---------
r 0.07 au
scale height is 0.010201020102010244
mass added [Msun]: 1.780737211301462e-09
relMass_added/rstep (rico) = 37.495
---------
r 0.08 au
scale height is 0.0118511851185118
mass added [Msun]: 2.1837057612726974e-09
relMass_added/rstep (rico) = 31.497
---------
r 0.09 au
scale height is 0.012976297629762962
mass added [Msun]: 2.477297133394598e-09
relMass_added/rstep (rico) = 26.325
---------
r 0.1 au
scale height is 0.015076507650765061
mass added [Msun]: 2.788158586229551e-09
relMass_added/rstep (rico) = 22.857
---------
r 0.11 au
scale height is 0.017551755175517592
mass added [Msun]: 3.0030751462142103e-09
relMass_added/rstep (rico) = 19.755
---------
r 0.12 au
scale height is 0.02017701770177015
mass added [Msun]: 3.1585058726316876e-09
relMass_added/rstep (rico) = 17.203
---------
r 0.13 au
scale height is 0.023927392739273936
mass added [Msun]: 3.350395658332276e-09
relMass_added/rstep (rico) = 15.432
---------
r 0.14 au
scale height is 0.026252625262526308
mass added [Msun]: 3.3983681047574233e-09
relMass_added/rstep (rico) = 13.535
---------
r 0.15 au
scale height is 0.027977797779777977
mass added [Msun]: 3.442502755468559e-09
relMass_added/rstep (rico) = 12.057
---------
r 0.16 au
scale height is 0.03127812781278128
mass added [Msun]: 3.519258669748794e-09
relMass_added/rstep (rico) = 10.973
---------
r 0.17 au
scale height is 0.03412841284128418
mass added [Msun]: 3.565312218316935e-09
relMass_added/rstep (rico) = 10.005
---------
r 0.18 au
scale height is 0.03735373537353744
mass added [Msun]: 3.5902578904580127e-09
relMass_added/rstep (rico) = 9.153
---------
r 0.19 au
scale height is 0.04215421542154214
mass added [Msun]: 3.5825822990299886e-09
relMass_added/rstep (rico) = 8.369
---------
r 0.2 au
scale height is 0.04657965796579655
mass added [Msun]: 3.6190413583131006e-09
relMass_added/rstep (rico) = 7.795
---------
r 0.21 au
scale height is 0.05018001800180018
mass added [Msun]: 3.5633933204599293e-09
relMass_added/rstep (rico) = 7.128
---------
r 0.22 au
scale height is 0.05318031803180317
mass added [Msun]: 3.446340551182571e-09
relMass_added/rstep (rico) = 6.449
---------
r 0.23 au
scale height is 0.05610561056105611
mass added [Msun]: 3.473205121180653e-09
relMass_added/rstep (rico) = 6.103
---------
r 0.24 au
scale height is 0.059855985598559834
mass added [Msun]: 3.3139365990491646e-09
relMass_added/rstep (rico) = 5.503
---------
r 0.25 au
scale height is 0.0633063306330633
mass added [Msun]: 3.3216121904771878e-09
relMass_added/rstep (rico) = 5.227
---------
r 0.26 au
scale height is 0.06570657065706571
mass added [Msun]: 3.219910604055876e-09
relMass_added/rstep (rico) = 4.823
---------
r 0.27 au
scale height is 0.06818181818181812
mass added [Msun]: 3.1009389369215114e-09
relMass_added/rstep (rico) = 4.438
---------
r 0.28 au
scale height is 0.07043204320432038
mass added [Msun]: 3.008831839785228e-09
relMass_added/rstep (rico) = 4.129
---------
r 0.29 au
scale height is 0.07410741074107417
mass added [Msun]: 2.935913721219005e-09
relMass_added/rstep (rico) = 3.873
---------
r 0.3 au
scale height is 0.07838283828382836
mass added [Msun]: 2.8188609519416454e-09
relMass_added/rstep (rico) = 3.585
---------
r 0.31 au
scale height is 0.08175817581758171
mass added [Msun]: 2.6749436126662046e-09
relMass_added/rstep (rico) = 3.29
---------
r 0.32 au
scale height is 0.08393339333933393
mass added [Msun]: 2.5694042305308802e-09
relMass_added/rstep (rico) = 3.063
---------
r 0.33 au
scale height is 0.08700870087008697
mass added [Msun]: 2.477297133394598e-09
relMass_added/rstep (rico) = 2.869
---------
r 0.34 au
scale height is 0.09023402340234023
mass added [Msun]: 2.346812079118197e-09
relMass_added/rstep (rico) = 2.646
---------
r 0.35 au
scale height is 0.09525952595259521
mass added [Msun]: 2.2796506541229915e-09
relMass_added/rstep (rico) = 2.506
---------
r 0.36 au
scale height is 0.098934893489349
mass added [Msun]: 2.1376522127045564e-09
relMass_added/rstep (rico) = 2.296
---------
r 0.37 au
scale height is 0.10261026102610261
mass added [Msun]: 2.0225183412842035e-09
relMass_added/rstep (rico) = 2.126
---------
r 0.38 au
scale height is 0.10658565856585654
mass added [Msun]: 1.968789201288039e-09
relMass_added/rstep (rico) = 2.027
---------
r 0.39 au
scale height is 0.11131113111311126
mass added [Msun]: 1.8402230448686444e-09
relMass_added/rstep (rico) = 1.86
---------
r 0.4 au
scale height is 0.11626162616261626
mass added [Msun]: 1.7519537434463736e-09
relMass_added/rstep (rico) = 1.74
---------
r 0.41 au
scale height is 0.12256225622562256
mass added [Msun]: 1.6310631784550027e-09
relMass_added/rstep (rico) = 1.594
---------
r 0.42 au
scale height is 0.13261326132613252
mass added [Msun]: 1.5811718341728498e-09
relMass_added/rstep (rico) = 1.522
---------
r 0.43 au
scale height is 0.15729072907290725
mass added [Msun]: 1.5946041191718909e-09
relMass_added/rstep (rico) = 1.511
---------
r 0.44 au
scale height is 0.1707920792079208
mass added [Msun]: 1.469875758466508e-09
relMass_added/rstep (rico) = 1.374
---------
r 0.45 au
scale height is 0.18879387938793873
mass added [Msun]: 1.4852269413225552e-09
relMass_added/rstep (rico) = 1.369
---------
r 0.46 au
scale height is 0.21482148214821478
mass added [Msun]: 1.3528229891891491e-09
relMass_added/rstep (rico) = 1.232
---------
r 0.47 au
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
scale height is 0.22307230723072308
mass added [Msun]: 1.2779859727659196e-09
relMass_added/rstep (rico) = 1.15
---------
r 0.48 au
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
scale height is 0.22517251725172505
mass added [Msun]: 1.2434458113398136e-09
relMass_added/rstep (rico) = 1.107
---------
r 0.49 au
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
scale height is 0.22907290729072904
mass added [Msun]: 1.1225552463484427e-09
relMass_added/rstep (rico) = 0.989
---------
r 0.5 au
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
scale height is 0.2113711371137113
mass added [Msun]: 1.0304481492121606e-09
relMass_added/rstep (rico) = 0.9
---------
r 0.51 au
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
scale height is 0.1960696069606961
mass added [Msun]: 9.939890899290487e-10
relMass_added/rstep (rico) = 0.861
---------
r 0.52 au
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
scale height is 0.19966996699669967
mass added [Msun]: 9.191520735058192e-10
relMass_added/rstep (rico) = 0.79
---------
r 0.53 au
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
scale height is 0.2046204620462046
mass added [Msun]: 8.615851377956427e-10
relMass_added/rstep (rico) = 0.735
---------
r 0.54 au
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
scale height is 0.21227122712271213
mass added [Msun]: 7.752347342303777e-10
relMass_added/rstep (rico) = 0.657
---------
r 0.55 au
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
scale height is 0.21707170717071697
mass added [Msun]: 7.119111049491835e-10
relMass_added/rstep (rico) = 0.599
---------
r 0.56 au
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
scale height is 0.2256225622562256
mass added [Msun]: 7.291811856622365e-10
relMass_added/rstep (rico) = 0.61
---------
r 0.57 au
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
scale height is 0.2301230123012301
mass added [Msun]: 6.044528249568538e-10
relMass_added/rstep (rico) = 0.503
---------
r 0.58 au
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
scale height is 0.23507350735073504
mass added [Msun]: 5.83344948529789e-10
relMass_added/rstep (rico) = 0.483
---------
r 0.59 au
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
scale height is 0.23732373237323737
mass added [Msun]: 5.488047871036832e-10
relMass_added/rstep (rico) = 0.453
---------
r 0.6 au
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
scale height is 0.2427242724272427
mass added [Msun]: 5.18102421391589e-10
relMass_added/rstep (rico) = 0.426
---------
r 0.61 au
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
scale height is 0.2547254725472548
mass added [Msun]: 4.85481157822489e-10
relMass_added/rstep (rico) = 0.397
---------
r 0.62 au
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
scale height is 0.25067506750675067
mass added [Msun]: 4.202386306842889e-10
relMass_added/rstep (rico) = 0.343
---------
r 0.63 au
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
scale height is 0.2593759375937594
mass added [Msun]: 3.6650949068812407e-10
relMass_added/rstep (rico) = 0.298
-------------------------------------------------------------
estimate of r of the accretion disk is 0.63 au
Calculations full theta region:
---------
r 0.02 au
scale height is 0.016951695169516955
mass added [Msun]: 3.837795714011771e-12
relMass_added/rstep (rico) = 100.0
---------
r 0.03 au
scale height is 0.014401440144014442
mass added [Msun]: 3.070236571209417e-11
relMass_added/rstep (rico) = 88.889
---------
r 0.04 au
scale height is 0.012451245124512502
mass added [Msun]: 8.059370999424718e-11
relMass_added/rstep (rico) = 70.0
---------
r 0.05 au
scale height is 0.010651065106510722
mass added [Msun]: 2.398622321257357e-10
relMass_added/rstep (rico) = 67.568
---------
r 0.06 au
scale height is 0.010201020102010178
mass added [Msun]: 5.161835235345832e-10
relMass_added/rstep (rico) = 59.251
---------
r 0.07 au
scale height is 0.012151215121512118
mass added [Msun]: 9.345032563618662e-10
relMass_added/rstep (rico) = 51.753
---------
r 0.08 au
scale height is 0.013951395139513965
mass added [Msun]: 1.3106072363350195e-09
relMass_added/rstep (rico) = 42.057
---------
r 0.09 au
scale height is 0.01590159015901597
mass added [Msun]: 1.6349009741690143e-09
relMass_added/rstep (rico) = 34.41
---------
r 0.1 au
scale height is 0.01860186018601864
mass added [Msun]: 1.867087614866727e-09
relMass_added/rstep (rico) = 28.211
---------
r 0.11 au
scale height is 0.020102010201020103
mass added [Msun]: 2.078166379137374e-09
relMass_added/rstep (rico) = 23.897
---------
r 0.12 au
scale height is 0.022877287728772885
mass added [Msun]: 2.258542777695927e-09
relMass_added/rstep (rico) = 20.617
---------
r 0.13 au
scale height is 0.025352535253525352
mass added [Msun]: 2.3640821598312503e-09
relMass_added/rstep (rico) = 17.75
---------
r 0.14 au
scale height is 0.02835283528352834
mass added [Msun]: 2.4178112998274153e-09
relMass_added/rstep (rico) = 15.364
---------
r 0.15 au
scale height is 0.03067806780678071
mass added [Msun]: 2.540620762675792e-09
relMass_added/rstep (rico) = 13.9
---------
r 0.16 au
scale height is 0.03390339033903391
mass added [Msun]: 2.521431784105733e-09
relMass_added/rstep (rico) = 12.123
---------
r 0.17 au
scale height is 0.037128712871287106
mass added [Msun]: 2.586674311243933e-09
relMass_added/rstep (rico) = 11.061
---------
r 0.18 au
scale height is 0.039003900390038995
mass added [Msun]: 2.559809741245851e-09
relMass_added/rstep (rico) = 9.866
---------
r 0.19 au
scale height is 0.04260426042604262
mass added [Msun]: 2.515675090534715e-09
relMass_added/rstep (rico) = 8.839
---------
r 0.2 au
scale height is 0.04560456045604561
mass added [Msun]: 2.5694042305308802e-09
relMass_added/rstep (rico) = 8.28
---------
r 0.21 au
scale height is 0.04830483048304828
mass added [Msun]: 2.456189256967533e-09
relMass_added/rstep (rico) = 7.335
---------
r 0.22 au
scale height is 0.052430243024302374
mass added [Msun]: 2.540620762675792e-09
relMass_added/rstep (rico) = 7.052
---------
r 0.23 au
scale height is 0.05573057305730574
mass added [Msun]: 2.3794333426872975e-09
relMass_added/rstep (rico) = 6.195
---------
r 0.24 au
scale height is 0.058055805580557986
mass added [Msun]: 2.40629791268538e-09
relMass_added/rstep (rico) = 5.896
---------
r 0.25 au
scale height is 0.06300630063006298
mass added [Msun]: 2.3276231005481387e-09
relMass_added/rstep (rico) = 5.395
---------
r 0.26 au
scale height is 0.06518151815181512
mass added [Msun]: 2.2201648205558094e-09
relMass_added/rstep (rico) = 4.894
---------
r 0.27 au
scale height is 0.06833183318331834
mass added [Msun]: 2.2527860841249094e-09
relMass_added/rstep (rico) = 4.731
---------
r 0.28 au
scale height is 0.07275727572757275
mass added [Msun]: 2.1357333148475502e-09
relMass_added/rstep (rico) = 4.293
---------
r 0.29 au
scale height is 0.07710771077107717
mass added [Msun]: 2.1146254384204855e-09
relMass_added/rstep (rico) = 4.077
---------
r 0.3 au
scale height is 0.08003300330032999
mass added [Msun]: 2.0033293627141446e-09
relMass_added/rstep (rico) = 3.719
---------
r 0.31 au
scale height is 0.08370837083708366
mass added [Msun]: 1.968789201288039e-09
relMass_added/rstep (rico) = 3.526
---------
r 0.32 au
scale height is 0.08775877587758776
mass added [Msun]: 1.899708878435827e-09
relMass_added/rstep (rico) = 3.29
---------
r 0.33 au
scale height is 0.09225922592259228
mass added [Msun]: 1.8517364320106797e-09
relMass_added/rstep (rico) = 3.108
---------
r 0.34 au
scale height is 0.0967596759675968
mass added [Msun]: 1.7308458670193087e-09
relMass_added/rstep (rico) = 2.823
---------
r 0.35 au
scale height is 0.1011851185118512
mass added [Msun]: 1.780737211301462e-09
relMass_added/rstep (rico) = 2.822
---------
r 0.36 au
scale height is 0.10553555355535549
mass added [Msun]: 1.6540899527390734e-09
relMass_added/rstep (rico) = 2.554
---------
r 0.37 au
scale height is 0.11033603360336033
mass added [Msun]: 1.5888474256008733e-09
relMass_added/rstep (rico) = 2.395
---------
r 0.38 au
scale height is 0.11453645364536459
mass added [Msun]: 1.554307264174767e-09
relMass_added/rstep (rico) = 2.289
---------
r 0.39 au
scale height is 0.11761176117611756
mass added [Msun]: 1.489064737036567e-09
relMass_added/rstep (rico) = 2.146
---------
r 0.4 au
scale height is 0.12136213621362128
mass added [Msun]: 1.423822209898367e-09
relMass_added/rstep (rico) = 2.011
---------
r 0.41 au
scale height is 0.1256375637563756
mass added [Msun]: 1.3547418870461552e-09
relMass_added/rstep (rico) = 1.877
---------
r 0.42 au
scale height is 0.13193819381938196
mass added [Msun]: 1.4103899248993256e-09
relMass_added/rstep (rico) = 1.917
---------
r 0.43 au
scale height is 0.1359135913591359
mass added [Msun]: 1.2453647091968194e-09
relMass_added/rstep (rico) = 1.664
---------
r 0.44 au
scale height is 0.1401890189018902
mass added [Msun]: 1.20314895634269e-09
relMass_added/rstep (rico) = 1.583
---------
r 0.45 au
scale height is 0.14416441644164418
mass added [Msun]: 1.2549591984818489e-09
relMass_added/rstep (rico) = 1.624
---------
r 0.46 au
scale height is 0.15024002400240022
mass added [Msun]: 1.1283119399194605e-09
relMass_added/rstep (rico) = 1.439
---------
r 0.47 au
scale height is 0.1556405640564056
mass added [Msun]: 1.1283119399194605e-09
relMass_added/rstep (rico) = 1.419
---------
r 0.48 au
scale height is 0.15939093909390933
mass added [Msun]: 1.070745004209284e-09
relMass_added/rstep (rico) = 1.328
---------
r 0.49 au
scale height is 0.1633663366336634
mass added [Msun]: 1.0534749234962311e-09
relMass_added/rstep (rico) = 1.29
---------
r 0.5 au
scale height is 0.16921692169216915
mass added [Msun]: 9.767190092159958e-10
relMass_added/rstep (rico) = 1.182
---------
r 0.51 au
scale height is 0.1761926192619262
mass added [Msun]: 9.268276649338427e-10
relMass_added/rstep (rico) = 1.109
---------
r 0.52 au
scale height is 0.1836183618361836
mass added [Msun]: 9.575300306459369e-10
relMass_added/rstep (rico) = 1.133
---------
r 0.53 au
scale height is 0.18856885688568858
mass added [Msun]: 9.421788477898898e-10
relMass_added/rstep (rico) = 1.102
---------
r 0.54 au
scale height is 0.19306930693069305
mass added [Msun]: 8.443150570825898e-10
relMass_added/rstep (rico) = 0.978
---------
r 0.55 au
scale height is 0.20012001200120014
mass added [Msun]: 8.001804063714544e-10
relMass_added/rstep (rico) = 0.919
---------
r 0.56 au
scale height is 0.206045604560456
mass added [Msun]: 8.404772613685778e-10
relMass_added/rstep (rico) = 0.956
---------
r 0.57 au
scale height is 0.2121212121212121
mass added [Msun]: 8.232071806555249e-10
relMass_added/rstep (rico) = 0.927
---------
r 0.58 au
scale height is 0.21932193219321935
mass added [Msun]: 7.003977178071481e-10
relMass_added/rstep (rico) = 0.783
---------
r 0.59 au
scale height is 0.22742274227422737
mass added [Msun]: 7.829103256584014e-10
relMass_added/rstep (rico) = 0.867
---------
r 0.6 au
scale height is 0.23604860486048604
mass added [Msun]: 7.253433899482246e-10
relMass_added/rstep (rico) = 0.797
---------
r 0.61 au
scale height is 0.2418241824182418
mass added [Msun]: 7.234244920912188e-10
relMass_added/rstep (rico) = 0.789
---------
r 0.62 au
scale height is 0.24542454245424544
mass added [Msun]: 6.332362928119422e-10
relMass_added/rstep (rico) = 0.686
---------
r 0.63 au
scale height is 0.2518001800180017
mass added [Msun]: 6.965599220931364e-10
relMass_added/rstep (rico) = 0.749
---------
r 0.64 au
scale height is 0.2581008100810081
mass added [Msun]: 5.948583356718244e-10
relMass_added/rstep (rico) = 0.635
---------
r 0.65 au
scale height is 0.2677767776777678
mass added [Msun]: 6.428307820969715e-10
relMass_added/rstep (rico) = 0.682
---------
r 0.66 au
scale height is 0.2727272727272727
mass added [Msun]: 6.082906206708655e-10
relMass_added/rstep (rico) = 0.641
---------
r 0.67 au
scale height is 0.2775277527752775
mass added [Msun]: 5.718315613877538e-10
relMass_added/rstep (rico) = 0.599
---------
r 0.68 au
scale height is 0.2817281728172817
mass added [Msun]: 5.641559699597302e-10
relMass_added/rstep (rico) = 0.588
---------
r 0.69 au
scale height is 0.292004200420042
mass added [Msun]: 5.910205399578126e-10
relMass_added/rstep (rico) = 0.612
---------
r 0.7 au
scale height is 0.3031803180318032
mass added [Msun]: 5.296158085336244e-10
relMass_added/rstep (rico) = 0.545
---------
r 0.71 au
scale height is 0.3153315331533153
mass added [Msun]: 4.6053548568141256e-10
relMass_added/rstep (rico) = 0.472
---------
r 0.72 au
scale height is 0.32845784578457854
mass added [Msun]: 5.02751238535542e-10
relMass_added/rstep (rico) = 0.512
---------
r 0.73 au
scale height is 0.34248424842484243
mass added [Msun]: 5.718315613877538e-10
relMass_added/rstep (rico) = 0.58
---------
r 0.74 au
scale height is 0.35546054605460536
mass added [Msun]: 5.334536042476361e-10
relMass_added/rstep (rico) = 0.538
---------
r 0.75 au
scale height is 0.3641614161416142
mass added [Msun]: 4.3367091568333013e-10
relMass_added/rstep (rico) = 0.435
---------
r 0.76 au
scale height is 0.3769126912691269
mass added [Msun]: 4.4710320068237135e-10
relMass_added/rstep (rico) = 0.447
---------
r 0.77 au
scale height is 0.39056405640564057
mass added [Msun]: 4.106441413992595e-10
relMass_added/rstep (rico) = 0.409
---------
r 0.78 au
scale height is 0.39941494149414936
mass added [Msun]: 4.85481157822489e-10
relMass_added/rstep (rico) = 0.481
---------
r 0.79 au
scale height is 0.4067656765676567
mass added [Msun]: 4.6245438353841847e-10
relMass_added/rstep (rico) = 0.456
---------
r 0.8 au
scale height is 0.41704170417041697
mass added [Msun]: 3.837795714011771e-10
relMass_added/rstep (rico) = 0.377
---------
r 0.81 au
scale height is 0.4273177317731773
mass added [Msun]: 3.9145516282920064e-10
relMass_added/rstep (rico) = 0.383
---------
r 0.82 au
scale height is 0.4353435343534353
mass added [Msun]: 3.9721185640021826e-10
relMass_added/rstep (rico) = 0.387
---------
r 0.83 au
scale height is 0.4459945994599461
mass added [Msun]: 3.837795714011771e-10
relMass_added/rstep (rico) = 0.373
---------
r 0.84 au
scale height is 0.4522952295229522
mass added [Msun]: 3.8953626497219473e-10
relMass_added/rstep (rico) = 0.377
---------
r 0.85 au
scale height is 0.45829582958295834
mass added [Msun]: 3.626716949741123e-10
relMass_added/rstep (rico) = 0.35
---------
r 0.86 au
scale height is 0.4654215421542154
mass added [Msun]: 3.5883389926010053e-10
relMass_added/rstep (rico) = 0.345
---------
r 0.87 au
scale height is 0.4739723972397239
mass added [Msun]: 3.415638185470476e-10
relMass_added/rstep (rico) = 0.327
---------
r 0.88 au
scale height is 0.4795979597959795
mass added [Msun]: 3.262126356910005e-10
relMass_added/rstep (rico) = 0.311
---------
r 0.89 au
scale height is 0.4873987398739874
mass added [Msun]: 3.204559421199828e-10
relMass_added/rstep (rico) = 0.305
---------
r 0.9 au
!! At this theta direction, you are no longer in the disk or you set your maxH too small!!
scale height is 0.4578457845784578
mass added [Msun]: 3.089425549779475e-10
relMass_added/rstep (rico) = 0.293
-------------------------------------------------------------
estimate of r of the accretion disk is 0.9 au
Calculations full theta region:
---------
r 0.02 au
scale height is 0.01680168016801686
mass added [Msun]: 3.837795714011771e-12
relMass_added/rstep (rico) = 100.0
---------
r 0.03 au
scale height is 0.014776477647764873
mass added [Msun]: 2.302677428407063e-11
relMass_added/rstep (rico) = 85.714
---------
r 0.04 au
scale height is 0.013126312631263187
mass added [Msun]: 1.0362048427831777e-10
relMass_added/rstep (rico) = 79.412
---------
r 0.05 au
scale height is 0.01222622262226223
mass added [Msun]: 2.552134149817827e-10
relMass_added/rstep (rico) = 66.169
---------
r 0.06 au
scale height is 0.010576057605760608
mass added [Msun]: 5.02751238535542e-10
relMass_added/rstep (rico) = 56.587
---------
r 0.07 au
scale height is 0.011626162616261659
mass added [Msun]: 7.52207959946307e-10
relMass_added/rstep (rico) = 45.848
---------
r 0.08 au
scale height is 0.014101410141014058
mass added [Msun]: 1.0822583913513193e-09
relMass_added/rstep (rico) = 39.746
---------
r 0.09 au
scale height is 0.01590159015901597
mass added [Msun]: 1.3125261341920253e-09
relMass_added/rstep (rico) = 32.525
---------
r 0.1 au
scale height is 0.01762676267626764
mass added [Msun]: 1.5523883663177615e-09
relMass_added/rstep (rico) = 27.782
---------
r 0.11 au
scale height is 0.020327032703270374
mass added [Msun]: 1.74427815201835e-09
relMass_added/rstep (rico) = 23.79
---------
r 0.12 au
scale height is 0.022352235223522428
mass added [Msun]: 1.8901143891507975e-09
relMass_added/rstep (rico) = 20.495
---------
r 0.13 au
scale height is 0.025202520252025126
mass added [Msun]: 2.0378695241402504e-09
relMass_added/rstep (rico) = 18.098
---------
r 0.14 au
scale height is 0.0285028502850285
mass added [Msun]: 2.1702734762736564e-09
relMass_added/rstep (rico) = 16.159
---------
r 0.15 au
scale height is 0.03127812781278128
mass added [Msun]: 2.2086514334137737e-09
relMass_added/rstep (rico) = 14.123
---------
r 0.16 au
scale height is 0.033753375337533746
mass added [Msun]: 2.3141908155490972e-09
relMass_added/rstep (rico) = 12.89
---------
r 0.17 au
scale height is 0.03772877287728774
mass added [Msun]: 2.3276231005481387e-09
relMass_added/rstep (rico) = 11.477
---------
r 0.18 au
scale height is 0.040654065406540686
mass added [Msun]: 2.3487309769752034e-09
relMass_added/rstep (rico) = 10.379
---------
r 0.19 au
scale height is 0.04305430543054304
mass added [Msun]: 2.425486891255439e-09
relMass_added/rstep (rico) = 9.681
---------
r 0.2 au
scale height is 0.04755475547554749
mass added [Msun]: 2.3794333426872975e-09
relMass_added/rstep (rico) = 8.673
---------
r 0.21 au
scale height is 0.05070507050705077
mass added [Msun]: 2.369838853402268e-09
relMass_added/rstep (rico) = 7.951
---------
r 0.22 au
scale height is 0.05460546054605465
mass added [Msun]: 2.3621632619742445e-09
relMass_added/rstep (rico) = 7.344
---------
r 0.23 au
scale height is 0.056780678067806796
mass added [Msun]: 2.3640821598312503e-09
relMass_added/rstep (rico) = 6.846
---------
r 0.24 au
scale height is 0.06023102310231026
mass added [Msun]: 2.425486891255439e-09
relMass_added/rstep (rico) = 6.563
---------
r 0.25 au
scale height is 0.06405640564056403
mass added [Msun]: 2.3276231005481387e-09
relMass_added/rstep (rico) = 5.925
---------
r 0.26 au
scale height is 0.06653165316531649
mass added [Msun]: 2.3410553855471798e-09
relMass_added/rstep (rico) = 5.624
---------
r 0.27 au
scale height is 0.06945694569456944
mass added [Msun]: 2.3218664069771213e-09
relMass_added/rstep (rico) = 5.283
---------
r 0.28 au
scale height is 0.07148214821482142
mass added [Msun]: 2.2163270248417974e-09
relMass_added/rstep (rico) = 4.801
---------
r 0.29 au
scale height is 0.07485748574857484
mass added [Msun]: 2.1798679655586854e-09
relMass_added/rstep (rico) = 4.509
---------
r 0.3 au
scale height is 0.07935793579357936
mass added [Msun]: 2.0628151962813267e-09
relMass_added/rstep (rico) = 4.092
---------
r 0.31 au
scale height is 0.08153315331533151
mass added [Msun]: 2.0877608684224034e-09
relMass_added/rstep (rico) = 3.977
---------
r 0.32 au
scale height is 0.0870087008700869
mass added [Msun]: 1.9649514055740272e-09
relMass_added/rstep (rico) = 3.608
---------
r 0.33 au
scale height is 0.09173417341734175
mass added [Msun]: 1.880519899865768e-09
relMass_added/rstep (rico) = 3.338
---------
r 0.34 au
scale height is 0.09623462346234614
mass added [Msun]: 1.880519899865768e-09
relMass_added/rstep (rico) = 3.23
---------
r 0.35 au
scale height is 0.09818481848184821
mass added [Msun]: 1.7538726413033794e-09
relMass_added/rstep (rico) = 2.924
---------
r 0.36 au
scale height is 0.10141014101410134
mass added [Msun]: 1.6790356248801497e-09
relMass_added/rstep (rico) = 2.723
---------
r 0.37 au
scale height is 0.10471047104710465
mass added [Msun]: 1.6790356248801497e-09
relMass_added/rstep (rico) = 2.651
---------
r 0.38 au
scale height is 0.10876087608760876
mass added [Msun]: 1.5197671027486611e-09
relMass_added/rstep (rico) = 2.343
---------
r 0.39 au
scale height is 0.11326132613261328
mass added [Msun]: 1.561982855602791e-09
relMass_added/rstep (rico) = 2.352
---------
r 0.4 au
scale height is 0.11941194119411941
mass added [Msun]: 1.4314978013263903e-09
relMass_added/rstep (rico) = 2.11
---------
r 0.41 au
scale height is 0.12473747374737477
mass added [Msun]: 1.4161466184703432e-09
relMass_added/rstep (rico) = 2.045
---------
r 0.42 au
scale height is 0.13088808880888092
mass added [Msun]: 1.3317151127620844e-09
relMass_added/rstep (rico) = 1.886
---------
r 0.43 au
scale height is 0.13576357635763578
mass added [Msun]: 1.287580462050949e-09
relMass_added/rstep (rico) = 1.791
---------
r 0.44 au
scale height is 0.1425142514251425
mass added [Msun]: 1.2300135263407724e-09
relMass_added/rstep (rico) = 1.682
---------
r 0.45 au
scale height is 0.14926492649264927
mass added [Msun]: 1.222337934912749e-09
relMass_added/rstep (rico) = 1.644
---------
r 0.46 au
scale height is 0.15286528652865283
mass added [Msun]: 1.1417442249185018e-09
relMass_added/rstep (rico) = 1.513
---------
r 0.47 au
scale height is 0.1563156315631563
mass added [Msun]: 1.070745004209284e-09
relMass_added/rstep (rico) = 1.399
---------
r 0.48 au
scale height is 0.16201620162016203
mass added [Msun]: 1.0457993320682075e-09
relMass_added/rstep (rico) = 1.348
---------
r 0.49 au
scale height is 0.16771677167716775
mass added [Msun]: 1.0457993320682075e-09
relMass_added/rstep (rico) = 1.33
---------
r 0.5 au
scale height is 0.17259225922592258
mass added [Msun]: 9.038008906497721e-10
relMass_added/rstep (rico) = 1.136
---------
r 0.51 au
scale height is 0.17649264926492647
mass added [Msun]: 9.306654606478546e-10
relMass_added/rstep (rico) = 1.157
---------
r 0.52 au
scale height is 0.181968196819682
mass added [Msun]: 8.884497077937251e-10
relMass_added/rstep (rico) = 1.092
---------
r 0.53 au
scale height is 0.18804380438043808
mass added [Msun]: 8.308827720835484e-10
relMass_added/rstep (rico) = 1.011
---------
r 0.54 au
scale height is 0.19306930693069305
mass added [Msun]: 7.829103256584014e-10
relMass_added/rstep (rico) = 0.944
---------
r 0.55 au
scale height is 0.1970447044704471
mass added [Msun]: 7.579646535173248e-10
relMass_added/rstep (rico) = 0.905
---------
r 0.56 au
scale height is 0.20124512451245125
mass added [Msun]: 6.831276370940952e-10
relMass_added/rstep (rico) = 0.809
---------
r 0.57 au
scale height is 0.2069456945694569
mass added [Msun]: 6.217229056699068e-10
relMass_added/rstep (rico) = 0.731
---------
r 0.58 au
scale height is 0.2113711371137113
mass added [Msun]: 6.485874756679892e-10
relMass_added/rstep (rico) = 0.757
---------
r 0.59 au
scale height is 0.2166966696669666
mass added [Msun]: 5.468858892466773e-10
relMass_added/rstep (rico) = 0.634
---------
r 0.6 au
scale height is 0.2275727572757276
mass added [Msun]: 5.35372502104642e-10
relMass_added/rstep (rico) = 0.617
---------
r 0.61 au
scale height is 0.2357485748574857
mass added [Msun]: 5.00832340678536e-10
relMass_added/rstep (rico) = 0.574
---------
r 0.62 au
scale height is 0.241074107410741
mass added [Msun]: 4.6437328139542433e-10
relMass_added/rstep (rico) = 0.529
---------
r 0.63 au
scale height is 0.2487998799879988
mass added [Msun]: 4.240764263983007e-10
relMass_added/rstep (rico) = 0.481
---------
r 0.64 au
scale height is 0.2587758775877588
mass added [Msun]: 4.3942760925434776e-10
relMass_added/rstep (rico) = 0.496
---------
r 0.65 au
scale height is 0.2688268826882688
mass added [Msun]: 3.185370442629769e-10
relMass_added/rstep (rico) = 0.358
---------
r 0.66 au
scale height is 0.27955295529552954
mass added [Msun]: 3.6075279711710644e-10
relMass_added/rstep (rico) = 0.404
---------
r 0.67 au
scale height is 0.2903540354035403
mass added [Msun]: 3.185370442629769e-10
relMass_added/rstep (rico) = 0.356
---------
r 0.68 au
scale height is 0.3020552055205521
mass added [Msun]: 3.089425549779475e-10
relMass_added/rstep (rico) = 0.344
---------
r 0.69 au
scale height is 0.31113111311131114
mass added [Msun]: 2.955102699789063e-10
relMass_added/rstep (rico) = 0.328
---------
r 0.7 au
scale height is 0.3185568556855686
mass added [Msun]: 2.3410553855471805e-10
relMass_added/rstep (rico) = 0.259
-------------------------------------------------------------
estimate of r of the accretion disk is 0.7000000000000002 au
Calculations for different theta regions:
---------
r 0.02 au
scale height is 0.013651365136513776
mass added [Msun]: 0.0
relMass_added/rstep (rico) = 1.0
---------
r 0.03 au
scale height is 0.010801080108010685
mass added [Msun]: 5.756693571017657e-12
relMass_added/rstep (rico) = 100.0
---------
r 0.04 au
scale height is 0.009300930093009484
mass added [Msun]: 2.302677428407063e-11
relMass_added/rstep (rico) = 80.0
---------
r 0.05 au
scale height is 0.00885088508850881
mass added [Msun]: 4.6053548568141254e-11
relMass_added/rstep (rico) = 61.538
---------
r 0.06 au
scale height is 0.008700870087008846
mass added [Msun]: 8.635040356526482e-11
relMass_added/rstep (rico) = 53.571
---------
r 0.07 au
scale height is 0.00990099009900986
mass added [Msun]: 1.285661564193943e-10
relMass_added/rstep (rico) = 44.371
---------
r 0.08 au
scale height is 0.011101110111011003
mass added [Msun]: 1.49674032846459e-10
relMass_added/rstep (rico) = 34.061
---------
r 0.09 au
scale height is 0.012601260126012596
mass added [Msun]: 1.6310631784550021e-10
relMass_added/rstep (rico) = 27.07
---------
r 0.1 au
scale height is 0.014551455145514602
mass added [Msun]: 1.6694411355951196e-10
relMass_added/rstep (rico) = 21.696
---------
r 0.11 au
scale height is 0.016801680168016794
mass added [Msun]: 1.9764647927160615e-10
relMass_added/rstep (rico) = 20.437
---------
r 0.12 au
scale height is 0.018751875187518932
mass added [Msun]: 1.938086835575944e-10
relMass_added/rstep (rico) = 16.694
---------
r 0.13 au
scale height is 0.02160216021602163
mass added [Msun]: 2.2067325355567684e-10
relMass_added/rstep (rico) = 15.972
---------
r 0.14 au
scale height is 0.02445244524452446
mass added [Msun]: 2.2259215141268272e-10
relMass_added/rstep (rico) = 13.876
---------
r 0.15 au
scale height is 0.02685268526852688
mass added [Msun]: 2.2259215141268272e-10
relMass_added/rstep (rico) = 12.185
---------
r 0.16 au
scale height is 0.029852985298529804
mass added [Msun]: 2.475378235537592e-10
relMass_added/rstep (rico) = 11.933
---------
r 0.17 au
scale height is 0.03270327032703276
mass added [Msun]: 2.417811299827416e-10
relMass_added/rstep (rico) = 10.439
---------
r 0.18 au
scale height is 0.036003600360036005
mass added [Msun]: 2.3218664069771216e-10
relMass_added/rstep (rico) = 9.111
---------
r 0.19 au
scale height is 0.039903990399039885
mass added [Msun]: 2.456189256967534e-10
relMass_added/rstep (rico) = 8.791
---------
r 0.2 au
scale height is 0.04455445544554463
mass added [Msun]: 2.475378235537592e-10
relMass_added/rstep (rico) = 8.139
---------
r 0.21 au
scale height is 0.04950495049504943
mass added [Msun]: 2.417811299827416e-10
relMass_added/rstep (rico) = 7.364
---------
r 0.22 au
scale height is 0.053255325532553216
mass added [Msun]: 2.5905121069579446e-10
relMass_added/rstep (rico) = 7.313
---------
r 0.23 au
scale height is 0.05520552055205509
mass added [Msun]: 2.5905121069579446e-10
relMass_added/rstep (rico) = 6.815
---------
r 0.24 au
scale height is 0.058505850585058465
mass added [Msun]: 2.609701085528004e-10
relMass_added/rstep (rico) = 6.424
---------
r 0.25 au
scale height is 0.06180618061806171
mass added [Msun]: 2.648079042668122e-10
relMass_added/rstep (rico) = 6.12
---------
r 0.26 au
scale height is 0.06600660066006604
mass added [Msun]: 2.628890064098063e-10
relMass_added/rstep (rico) = 5.727
---------
r 0.27 au
scale height is 0.0694569456945695
mass added [Msun]: 2.5329451712477684e-10
relMass_added/rstep (rico) = 5.23
---------
r 0.28 au
scale height is 0.07455745574557453
mass added [Msun]: 2.456189256967534e-10
relMass_added/rstep (rico) = 4.827
---------
r 0.29 au
scale height is 0.07785778577857803
mass added [Msun]: 2.5329451712477684e-10
relMass_added/rstep (rico) = 4.741
---------
r 0.3 au
scale height is 0.08055805580558051
mass added [Msun]: 2.2642994712669446e-10
relMass_added/rstep (rico) = 4.066
---------
r 0.31 au
scale height is 0.08370837083708366
mass added [Msun]: 2.3218664069771216e-10
relMass_added/rstep (rico) = 4.003
---------
r 0.32 au
scale height is 0.0862586258625863
mass added [Msun]: 2.2067325355567684e-10
relMass_added/rstep (rico) = 3.665
---------
r 0.33 au
scale height is 0.09015901590159005
mass added [Msun]: 1.9188978570058852e-10
relMass_added/rstep (rico) = 3.088
---------
r 0.34 au
scale height is 0.09420942094209415
mass added [Msun]: 2.1683545784166507e-10
relMass_added/rstep (rico) = 3.372
---------
r 0.35 au
scale height is 0.09780978097809784
mass added [Msun]: 1.938086835575944e-10
relMass_added/rstep (rico) = 2.926
---------
r 0.36 au
scale height is 0.10201020102010191
mass added [Msun]: 1.6886301141651792e-10
relMass_added/rstep (rico) = 2.486
---------
r 0.37 au
scale height is 0.10696069606960684
mass added [Msun]: 1.8997088784358264e-10
relMass_added/rstep (rico) = 2.721
---------
r 0.38 au
scale height is 0.11191119111911191
mass added [Msun]: 1.6694411355951196e-10
relMass_added/rstep (rico) = 2.335
---------
r 0.39 au
scale height is 0.11671167116711674
mass added [Msun]: 1.5926852213148845e-10
relMass_added/rstep (rico) = 2.179
---------
r 0.4 au
scale height is 0.1203120312031203
mass added [Msun]: 1.707819092735238e-10
relMass_added/rstep (rico) = 2.283
---------
r 0.41 au
scale height is 0.12301230123012304
mass added [Msun]: 1.3432284999041193e-10
relMass_added/rstep (rico) = 1.764
---------
r 0.42 au
scale height is 0.12631263126312628
mass added [Msun]: 1.362417478474178e-10
relMass_added/rstep (rico) = 1.758
---------
r 0.43 au
scale height is 0.13381338133813372
mass added [Msun]: 1.4583623713244726e-10
relMass_added/rstep (rico) = 1.847
---------
r 0.44 au
scale height is 0.14161416141614175
mass added [Msun]: 1.2280946284837665e-10
relMass_added/rstep (rico) = 1.531
---------
r 0.45 au
scale height is 0.14806480648064788
mass added [Msun]: 1.419984414184355e-10
relMass_added/rstep (rico) = 1.74
---------
r 0.46 au
scale height is 0.1558655865586559
mass added [Msun]: 1.0170158642131189e-10
relMass_added/rstep (rico) = 1.231
---------
r 0.47 au
scale height is 0.16291629162916293
mass added [Msun]: 1.1897166713436486e-10
relMass_added/rstep (rico) = 1.419
---------
r 0.48 au
scale height is 0.16696669666966704
mass added [Msun]: 1.2280946284837665e-10
relMass_added/rstep (rico) = 1.444
---------
r 0.49 au
scale height is 0.16876687668766868
mass added [Msun]: 1.0745827999232958e-10
relMass_added/rstep (rico) = 1.248
---------
r 0.5 au
scale height is 0.17056705670567046
mass added [Msun]: 1.0170158642131189e-10
relMass_added/rstep (rico) = 1.167
---------
r 0.51 au
scale height is 0.17221722172217227
mass added [Msun]: 1.0745827999232958e-10
relMass_added/rstep (rico) = 1.218
---------
r 0.52 au
scale height is 0.17551755175517553
mass added [Msun]: 1.1705276927735897e-10
relMass_added/rstep (rico) = 1.31
---------
r 0.53 au
scale height is 0.17581758175817572
mass added [Msun]: 1.1321497356334722e-10
relMass_added/rstep (rico) = 1.251
---------
r 0.54 au
scale height is 0.1779177917791778
mass added [Msun]: 1.0553938213532366e-10
relMass_added/rstep (rico) = 1.153
---------
r 0.55 au
scale height is 0.18166816681668146
mass added [Msun]: 1.285661564193943e-10
relMass_added/rstep (rico) = 1.385
---------
r 0.56 au
scale height is 0.18706870687068694
mass added [Msun]: 1.1897166713436486e-10
relMass_added/rstep (rico) = 1.265
---------
r 0.57 au
scale height is 0.19396939693969387
mass added [Msun]: 9.210709713628251e-11
relMass_added/rstep (rico) = 0.97
---------
r 0.58 au
scale height is 0.20177017701770175
mass added [Msun]: 1.1321497356334722e-10
relMass_added/rstep (rico) = 1.178
---------
r 0.59 au
scale height is 0.21002100210020994
mass added [Msun]: 9.402599499328838e-11
relMass_added/rstep (rico) = 0.969
---------
r 0.6 au
scale height is 0.21677167716771678
mass added [Msun]: 1.0745827999232958e-10
relMass_added/rstep (rico) = 1.095
---------
r 0.61 au
scale height is 0.22262226222622267
mass added [Msun]: 7.867481213724131e-11
relMass_added/rstep (rico) = 0.795
---------
r 0.62 au
scale height is 0.2296729672967297
mass added [Msun]: 6.908032285221189e-11
relMass_added/rstep (rico) = 0.694
---------
r 0.63 au
scale height is 0.23552355235523545
mass added [Msun]: 7.483701642322955e-11
relMass_added/rstep (rico) = 0.746
---------
r 0.64 au
scale height is 0.24602460246024593
mass added [Msun]: 6.7161424995206e-11
relMass_added/rstep (rico) = 0.665
---------
r 0.65 au
scale height is 0.2554755475547554
mass added [Msun]: 4.413465071113537e-11
relMass_added/rstep (rico) = 0.435
---------
r 0.66 au
scale height is 0.26432643264326433
mass added [Msun]: 3.837795714011771e-11
relMass_added/rstep (rico) = 0.377
---------
r 0.67 au
scale height is 0.26957695769576956
mass added [Msun]: 4.2215752854129486e-11
relMass_added/rstep (rico) = 0.413
---------
r 0.68 au
scale height is 0.2764776477647764
mass added [Msun]: 3.837795714011771e-11
relMass_added/rstep (rico) = 0.374
---------
r 0.69 au
scale height is 0.2838283828382837
mass added [Msun]: 2.6864569998082404e-11
relMass_added/rstep (rico) = 0.261
-------------------------------------------------------------
estimate of r of the accretion disk is 0.69 au
---------
r 0.02 au
scale height is 0.019201920192019214
mass added [Msun]: 1.9188978570058856e-12
relMass_added/rstep (rico) = 100.0
---------
r 0.03 au
scale height is 0.01770177017701775
mass added [Msun]: 3.837795714011771e-12
relMass_added/rstep (rico) = 66.667
---------
r 0.04 au
scale height is 0.015301530153015202
mass added [Msun]: 1.5351182856047085e-11
relMass_added/rstep (rico) = 72.727
---------
r 0.05 au
scale height is 0.010801080108010815
mass added [Msun]: 3.645905928311183e-11
relMass_added/rstep (rico) = 63.333
---------
r 0.06 au
scale height is 0.009600960096009671
mass added [Msun]: 7.867481213724131e-11
relMass_added/rstep (rico) = 57.746
---------
r 0.07 au
scale height is 0.011551155115511415
mass added [Msun]: 1.2472836070538253e-10
relMass_added/rstep (rico) = 47.794
---------
r 0.08 au
scale height is 0.013651365136513646
mass added [Msun]: 1.938086835575944e-10
relMass_added/rstep (rico) = 42.616
---------
r 0.09 au
scale height is 0.016051605160516062
mass added [Msun]: 2.686456999808239e-10
relMass_added/rstep (rico) = 37.135
---------
r 0.1 au
scale height is 0.019201920192019214
mass added [Msun]: 3.6075279711710644e-10
relMass_added/rstep (rico) = 33.274
---------
r 0.11 au
scale height is 0.022052205220522043
mass added [Msun]: 4.106441413992595e-10
relMass_added/rstep (rico) = 27.471
---------
r 0.12 au
scale height is 0.0250525052505251
mass added [Msun]: 4.70129974966442e-10
relMass_added/rstep (rico) = 23.926
---------
r 0.13 au
scale height is 0.028502850285028566
mass added [Msun]: 5.334536042476361e-10
relMass_added/rstep (rico) = 21.352
---------
r 0.14 au
scale height is 0.03285328532853286
mass added [Msun]: 5.622370721027243e-10
relMass_added/rstep (rico) = 18.37
---------
r 0.15 au
scale height is 0.03570357035703582
mass added [Msun]: 6.044528249568538e-10
relMass_added/rstep (rico) = 16.492
---------
r 0.16 au
scale height is 0.040354035403540434
mass added [Msun]: 6.581819649530187e-10
relMass_added/rstep (rico) = 15.224
---------
r 0.17 au
scale height is 0.04275427542754272
mass added [Msun]: 6.581819649530187e-10
relMass_added/rstep (rico) = 13.213
---------
r 0.18 au
scale height is 0.04620462046204619
mass added [Msun]: 7.119111049491835e-10
relMass_added/rstep (rico) = 12.504
---------
r 0.19 au
scale height is 0.051155115511551115
mass added [Msun]: 7.38775674947266e-10
relMass_added/rstep (rico) = 11.486
---------
r 0.2 au
scale height is 0.053855385538553985
mass added [Msun]: 7.598835513743306e-10
relMass_added/rstep (rico) = 10.566
---------
r 0.21 au
scale height is 0.057155715571557096
mass added [Msun]: 7.52207959946307e-10
relMass_added/rstep (rico) = 9.469
---------
r 0.22 au
scale height is 0.05970597059705974
mass added [Msun]: 7.426134706612776e-10
relMass_added/rstep (rico) = 8.549
---------
r 0.23 au
scale height is 0.06390639063906393
mass added [Msun]: 7.598835513743306e-10
relMass_added/rstep (rico) = 8.044
---------
r 0.24 au
scale height is 0.06675667566756664
mass added [Msun]: 8.251260785125309e-10
relMass_added/rstep (rico) = 8.033
---------
r 0.25 au
scale height is 0.07110711071107119
mass added [Msun]: 7.771536320873835e-10
relMass_added/rstep (rico) = 7.034
---------
r 0.26 au
scale height is 0.07530753075307539
mass added [Msun]: 7.963426106574424e-10
relMass_added/rstep (rico) = 6.723
---------
r 0.27 au
scale height is 0.08010801080107996
mass added [Msun]: 7.963426106574424e-10
relMass_added/rstep (rico) = 6.299
---------
r 0.28 au
scale height is 0.08445844584458452
mass added [Msun]: 7.944237128004365e-10
relMass_added/rstep (rico) = 5.913
---------
r 0.29 au
scale height is 0.08895889588958904
mass added [Msun]: 7.886670192294189e-10
relMass_added/rstep (rico) = 5.544
---------
r 0.3 au
scale height is 0.09375937593759374
mass added [Msun]: 7.426134706612776e-10
relMass_added/rstep (rico) = 4.962
---------
r 0.31 au
scale height is 0.09900990099009886
mass added [Msun]: 7.809914278013955e-10
relMass_added/rstep (rico) = 4.959
---------
r 0.32 au
scale height is 0.10456045604560468
mass added [Msun]: 7.598835513743306e-10
relMass_added/rstep (rico) = 4.603
---------
r 0.33 au
scale height is 0.10861086108610866
mass added [Msun]: 7.138300028061893e-10
relMass_added/rstep (rico) = 4.145
---------
r 0.34 au
scale height is 0.11371137113711381
mass added [Msun]: 7.675591428023542e-10
relMass_added/rstep (rico) = 4.267
---------
r 0.35 au
scale height is 0.11926192619261924
mass added [Msun]: 6.773709435230776e-10
relMass_added/rstep (rico) = 3.629
---------
r 0.36 au
scale height is 0.12436243624362428
mass added [Msun]: 6.927221263791246e-10
relMass_added/rstep (rico) = 3.578
---------
r 0.37 au
scale height is 0.13141314131413143
mass added [Msun]: 6.735331478090657e-10
relMass_added/rstep (rico) = 3.362
---------
r 0.38 au
scale height is 0.13636363636363624
mass added [Msun]: 6.485874756679892e-10
relMass_added/rstep (rico) = 3.136
---------
r 0.39 au
scale height is 0.14206420642064202
mass added [Msun]: 6.735331478090657e-10
relMass_added/rstep (rico) = 3.154
---------
r 0.4 au
scale height is 0.1474647464746475
mass added [Msun]: 5.756693571017656e-10
relMass_added/rstep (rico) = 2.625
---------
r 0.41 au
scale height is 0.1534653465346535
mass added [Msun]: 6.082906206708655e-10
relMass_added/rstep (rico) = 2.699
---------
r 0.42 au
scale height is 0.15976597659765965
mass added [Msun]: 6.082906206708655e-10
relMass_added/rstep (rico) = 2.628
---------
r 0.43 au
scale height is 0.16591659165916584
mass added [Msun]: 5.392102978186538e-10
relMass_added/rstep (rico) = 2.276
---------
r 0.44 au
scale height is 0.17266726672667268
mass added [Msun]: 5.468858892466773e-10
relMass_added/rstep (rico) = 2.257
---------
r 0.45 au
scale height is 0.18076807680768064
mass added [Msun]: 5.123457278205714e-10
relMass_added/rstep (rico) = 2.07
---------
r 0.46 au
scale height is 0.18736873687368738
mass added [Msun]: 5.430480935326655e-10
relMass_added/rstep (rico) = 2.147
---------
r 0.47 au
scale height is 0.19561956195619556
mass added [Msun]: 4.3558981354033604e-10
relMass_added/rstep (rico) = 1.693
---------
r 0.48 au
scale height is 0.2049204920492049
mass added [Msun]: 4.3942760925434776e-10
relMass_added/rstep (rico) = 1.68
---------
r 0.49 au
scale height is 0.2133213321332133
mass added [Msun]: 4.662921792524302e-10
relMass_added/rstep (rico) = 1.751
---------
r 0.5 au
scale height is 0.22202220222022204
mass added [Msun]: 4.106441413992595e-10
relMass_added/rstep (rico) = 1.519
---------
r 0.51 au
scale height is 0.22907290729072904
mass added [Msun]: 3.6842838854513003e-10
relMass_added/rstep (rico) = 1.344
---------
r 0.52 au
scale height is 0.23507350735073515
mass added [Msun]: 4.0296854997123594e-10
relMass_added/rstep (rico) = 1.449
---------
r 0.53 au
scale height is 0.2403240324032403
mass added [Msun]: 3.262126356910005e-10
relMass_added/rstep (rico) = 1.159
---------
r 0.54 au
scale height is 0.2475247524752474
mass added [Msun]: 3.204559421199828e-10
relMass_added/rstep (rico) = 1.126
---------
r 0.55 au
scale height is 0.2568256825682569
mass added [Msun]: 3.3005043140501225e-10
relMass_added/rstep (rico) = 1.146
---------
r 0.56 au
scale height is 0.26537653765376523
mass added [Msun]: 2.571323128387886e-10
relMass_added/rstep (rico) = 0.885
---------
r 0.57 au
scale height is 0.2736273627362735
mass added [Msun]: 2.897535764078886e-10
relMass_added/rstep (rico) = 0.988
---------
r 0.58 au
scale height is 0.28172817281728174
mass added [Msun]: 2.686456999808239e-10
relMass_added/rstep (rico) = 0.907
---------
r 0.59 au
scale height is 0.289078907890789
mass added [Msun]: 2.2834884498370035e-10
relMass_added/rstep (rico) = 0.765
---------
r 0.6 au
scale height is 0.2938793879387938
mass added [Msun]: 2.2642994712669446e-10
relMass_added/rstep (rico) = 0.753
---------
r 0.61 au
scale height is 0.29987998799880006
mass added [Msun]: 2.2259215141268272e-10
relMass_added/rstep (rico) = 0.735
---------
r 0.62 au
scale height is 0.30333033303330337
mass added [Msun]: 2.0148427498561797e-10
relMass_added/rstep (rico) = 0.661
---------
r 0.63 au
scale height is 0.3102310231023102
mass added [Msun]: 2.1107876427064736e-10
relMass_added/rstep (rico) = 0.688
---------
r 0.64 au
scale height is 0.3156315631563157
mass added [Msun]: 1.8997088784358264e-10
relMass_added/rstep (rico) = 0.615
---------
r 0.65 au
scale height is 0.32073207320732056
mass added [Msun]: 1.4583623713244726e-10
relMass_added/rstep (rico) = 0.47
---------
r 0.66 au
scale height is 0.3249324932493249
mass added [Msun]: 1.6118741998849433e-10
relMass_added/rstep (rico) = 0.517
---------
r 0.67 au
scale height is 0.3286828682868287
mass added [Msun]: 1.4775513498945312e-10
relMass_added/rstep (rico) = 0.471
---------
r 0.68 au
scale height is 0.3313831383138314
mass added [Msun]: 1.419984414184355e-10
relMass_added/rstep (rico) = 0.451
---------
r 0.69 au
scale height is 0.3358835883588359
mass added [Msun]: 1.49674032846459e-10
relMass_added/rstep (rico) = 0.473
---------
r 0.7 au
scale height is 0.3393339333933393
mass added [Msun]: 1.1513387142035309e-10
relMass_added/rstep (rico) = 0.363
---------
r 0.71 au
scale height is 0.34233423342334235
mass added [Msun]: 1.1513387142035309e-10
relMass_added/rstep (rico) = 0.361
---------
r 0.72 au
scale height is 0.34683468346834684
mass added [Msun]: 1.0170158642131189e-10
relMass_added/rstep (rico) = 0.318
---------
r 0.73 au
scale height is 0.3510351035103509
mass added [Msun]: 9.402599499328838e-11
relMass_added/rstep (rico) = 0.293
-------------------------------------------------------------
estimate of r of the accretion disk is 0.73 au
---------
r 0.02 au
scale height is 0.02010201020102017
mass added [Msun]: 0.0
relMass_added/rstep (rico) = 1.0
---------
r 0.03 au
scale height is 0.018451845184518482
mass added [Msun]: 3.837795714011771e-12
relMass_added/rstep (rico) = 100.0
---------
r 0.04 au
scale height is 0.01665166516651657
mass added [Msun]: 1.5351182856047085e-11
relMass_added/rstep (rico) = 80.0
---------
r 0.05 au
scale height is 0.01560156015601565
mass added [Msun]: 2.1107876427064743e-11
relMass_added/rstep (rico) = 52.381
---------
r 0.06 au
scale height is 0.01665166516651657
mass added [Msun]: 7.291811856622365e-11
relMass_added/rstep (rico) = 64.407
---------
r 0.07 au
scale height is 0.018151815181518163
mass added [Msun]: 1.49674032846459e-10
relMass_added/rstep (rico) = 56.934
---------
r 0.08 au
scale height is 0.020102010201020037
mass added [Msun]: 2.648079042668122e-10
relMass_added/rstep (rico) = 50.182
---------
r 0.09 au
scale height is 0.022202220222022267
mass added [Msun]: 3.8186067354417114e-10
relMass_added/rstep (rico) = 41.983
---------
r 0.1 au
scale height is 0.02475247524752491
mass added [Msun]: 4.835622599654832e-10
relMass_added/rstep (rico) = 34.711
---------
r 0.11 au
scale height is 0.02760276027602761
mass added [Msun]: 5.603181742457184e-10
relMass_added/rstep (rico) = 28.684
---------
r 0.12 au
scale height is 0.03045304530453044
mass added [Msun]: 6.082906206708655e-10
relMass_added/rstep (rico) = 23.745
---------
r 0.13 au
scale height is 0.03255325532553254
mass added [Msun]: 6.639386585240362e-10
relMass_added/rstep (rico) = 20.583
---------
r 0.14 au
scale height is 0.034203420342034226
mass added [Msun]: 6.965599220931364e-10
relMass_added/rstep (rico) = 17.759
---------
r 0.15 au
scale height is 0.038103810381038106
mass added [Msun]: 6.888843306651127e-10
relMass_added/rstep (rico) = 14.94
---------
r 0.16 au
scale height is 0.0397539753975398
mass added [Msun]: 7.195866963772071e-10
relMass_added/rstep (rico) = 13.499
---------
r 0.17 au
scale height is 0.04230423042304231
mass added [Msun]: 6.869654328081069e-10
relMass_added/rstep (rico) = 11.416
---------
r 0.18 au
scale height is 0.04575457545754577
mass added [Msun]: 7.023166156641541e-10
relMass_added/rstep (rico) = 10.451
---------
r 0.19 au
scale height is 0.04830483048304828
mass added [Msun]: 7.138300028061893e-10
relMass_added/rstep (rico) = 9.602
---------
r 0.2 au
scale height is 0.05145514551455144
mass added [Msun]: 6.773709435230776e-10
relMass_added/rstep (rico) = 8.351
---------
r 0.21 au
scale height is 0.05400540054005408
mass added [Msun]: 6.831276370940952e-10
relMass_added/rstep (rico) = 7.768
---------
r 0.22 au
scale height is 0.05760576057605764
mass added [Msun]: 6.389929863829598e-10
relMass_added/rstep (rico) = 6.774
---------
r 0.23 au
scale height is 0.06090609060906088
mass added [Msun]: 6.274795992409243e-10
relMass_added/rstep (rico) = 6.237
---------
r 0.24 au
scale height is 0.06255625562556257
mass added [Msun]: 5.967772335288302e-10
relMass_added/rstep (rico) = 5.6
---------
r 0.25 au
scale height is 0.06525652565256518
mass added [Msun]: 5.603181742457184e-10
relMass_added/rstep (rico) = 4.995
---------
r 0.26 au
scale height is 0.06810681068106814
mass added [Msun]: 5.564803785317067e-10
relMass_added/rstep (rico) = 4.726
---------
r 0.27 au
scale height is 0.0715571557155716
mass added [Msun]: 5.679937656737421e-10
relMass_added/rstep (rico) = 4.602
---------
r 0.28 au
scale height is 0.07395739573957402
mass added [Msun]: 4.969945449645243e-10
relMass_added/rstep (rico) = 3.871
---------
r 0.29 au
scale height is 0.07680768076807672
mass added [Msun]: 4.893189535365008e-10
relMass_added/rstep (rico) = 3.671
---------
r 0.3 au
scale height is 0.0793579357935795
mass added [Msun]: 4.70129974966442e-10
relMass_added/rstep (rico) = 3.407
---------
r 0.31 au
scale height is 0.08145814581458133
mass added [Msun]: 4.797244642514713e-10
relMass_added/rstep (rico) = 3.36
---------
r 0.32 au
scale height is 0.0841584158415842
mass added [Msun]: 4.0872524354225357e-10
relMass_added/rstep (rico) = 2.783
---------
r 0.33 au
scale height is 0.08880888088808894
mass added [Msun]: 4.1256303925626533e-10
relMass_added/rstep (rico) = 2.732
---------
r 0.34 au
scale height is 0.09390939093909384
mass added [Msun]: 3.626716949741123e-10
relMass_added/rstep (rico) = 2.345
---------
r 0.35 au
scale height is 0.09645964596459634
mass added [Msun]: 4.0872524354225357e-10
relMass_added/rstep (rico) = 2.575
---------
r 0.36 au
scale height is 0.09960996099609963
mass added [Msun]: 3.262126356910005e-10
relMass_added/rstep (rico) = 2.014
---------
r 0.37 au
scale height is 0.10486048604860487
mass added [Msun]: 3.4732051211806523e-10
relMass_added/rstep (rico) = 2.099
---------
r 0.38 au
scale height is 0.1098109810981098
mass added [Msun]: 3.185370442629769e-10
relMass_added/rstep (rico) = 1.889
---------
r 0.39 au
scale height is 0.11386138613861378
mass added [Msun]: 2.974291678359122e-10
relMass_added/rstep (rico) = 1.733
---------
r 0.4 au
scale height is 0.11881188118811883
mass added [Msun]: 2.83996882836871e-10
relMass_added/rstep (rico) = 1.628
---------
r 0.41 au
scale height is 0.12736273627362732
mass added [Msun]: 2.724834956948357e-10
relMass_added/rstep (rico) = 1.538
---------
r 0.42 au
scale height is 0.1308130813081308
mass added [Msun]: 2.552134149817827e-10
relMass_added/rstep (rico) = 1.42
---------
r 0.43 au
scale height is 0.1344134413441345
mass added [Msun]: 2.4945672141076507e-10
relMass_added/rstep (rico) = 1.369
---------
r 0.44 au
scale height is 0.13801380138013805
mass added [Msun]: 2.1107876427064736e-10
relMass_added/rstep (rico) = 1.145
---------
r 0.45 au
scale height is 0.14296429642964312
mass added [Msun]: 2.3218664069771216e-10
relMass_added/rstep (rico) = 1.244
---------
r 0.46 au
scale height is 0.15046504650465042
mass added [Msun]: 2.0724096855663562e-10
relMass_added/rstep (rico) = 1.098
---------
r 0.47 au
scale height is 0.15871587158715886
mass added [Msun]: 1.6310631784550021e-10
relMass_added/rstep (rico) = 0.857
---------
r 0.48 au
scale height is 0.16381638163816375
mass added [Msun]: 1.9572758141460027e-10
relMass_added/rstep (rico) = 1.018
---------
r 0.49 au
scale height is 0.17011701170117005
mass added [Msun]: 1.5734962427448256e-10
relMass_added/rstep (rico) = 0.812
---------
r 0.5 au
scale height is 0.17581758175817583
mass added [Msun]: 1.3240395213340604e-10
relMass_added/rstep (rico) = 0.678
---------
r 0.51 au
scale height is 0.18136813681368127
mass added [Msun]: 1.5159293070346488e-10
relMass_added/rstep (rico) = 0.771
---------
r 0.52 au
scale height is 0.1858685868586858
mass added [Msun]: 1.2472836070538253e-10
relMass_added/rstep (rico) = 0.63
---------
r 0.53 au
scale height is 0.1915691569156916
mass added [Msun]: 1.1897166713436486e-10
relMass_added/rstep (rico) = 0.597
---------
r 0.54 au
scale height is 0.19846984698469852
mass added [Msun]: 1.0362048427831777e-10
relMass_added/rstep (rico) = 0.518
---------
r 0.55 au
scale height is 0.20747074707470742
mass added [Msun]: 1.1129607570634133e-10
relMass_added/rstep (rico) = 0.553
---------
r 0.56 au
scale height is 0.2152715271527153
mass added [Msun]: 8.443150570825893e-11
relMass_added/rstep (rico) = 0.418
---------
r 0.57 au
scale height is 0.2223222322232222
mass added [Msun]: 7.675591428023542e-11
relMass_added/rstep (rico) = 0.378
---------
r 0.58 au
scale height is 0.22907290729072904
mass added [Msun]: 9.594489285029426e-11
relMass_added/rstep (rico) = 0.471
---------
r 0.59 au
scale height is 0.23432343234323444
mass added [Msun]: 6.7161424995206e-11
relMass_added/rstep (rico) = 0.328
---------
r 0.6 au
scale height is 0.24122412241224123
mass added [Msun]: 5.3729139996164796e-11
relMass_added/rstep (rico) = 0.262
-------------------------------------------------------------
estimate of r of the accretion disk is 0.6 au
---------
r 0.02 au
scale height is 0.01560156015601552
mass added [Msun]: 1.9188978570058856e-12
relMass_added/rstep (rico) = 100.0
---------
r 0.03 au
scale height is 0.012601260126012596
mass added [Msun]: 7.675591428023542e-12
relMass_added/rstep (rico) = 80.0
---------
r 0.04 au
scale height is 0.012601260126012596
mass added [Msun]: 3.262126356910005e-11
relMass_added/rstep (rico) = 77.273
---------
r 0.05 au
scale height is 0.01170117011701164
mass added [Msun]: 1.1321497356334722e-10
relMass_added/rstep (rico) = 72.84
---------
r 0.06 au
scale height is 0.012451245124512502
mass added [Msun]: 1.9764647927160615e-10
relMass_added/rstep (rico) = 55.978
---------
r 0.07 au
scale height is 0.012901290129012914
mass added [Msun]: 2.648079042668122e-10
relMass_added/rstep (rico) = 42.857
---------
r 0.08 au
scale height is 0.013651365136513646
mass added [Msun]: 3.4732051211806523e-10
relMass_added/rstep (rico) = 35.984
---------
r 0.09 au
scale height is 0.013951395139513965
mass added [Msun]: 3.8761736711518887e-10
relMass_added/rstep (rico) = 28.652
---------
r 0.1 au
scale height is 0.01590159015901584
mass added [Msun]: 4.0872524354225357e-10
relMass_added/rstep (rico) = 23.203
---------
r 0.11 au
scale height is 0.017851785178517845
mass added [Msun]: 4.3558981354033604e-10
relMass_added/rstep (rico) = 19.825
---------
r 0.12 au
scale height is 0.0187518751875188
mass added [Msun]: 4.509409963963831e-10
relMass_added/rstep (rico) = 17.029
---------
r 0.13 au
scale height is 0.021002100210020993
mass added [Msun]: 4.4902209853937726e-10
relMass_added/rstep (rico) = 14.498
---------
r 0.14 au
scale height is 0.02340234023402354
mass added [Msun]: 4.85481157822489e-10
relMass_added/rstep (rico) = 13.551
---------
r 0.15 au
scale height is 0.024902490249024873
mass added [Msun]: 4.70129974966442e-10
relMass_added/rstep (rico) = 11.6
---------
r 0.16 au
scale height is 0.026702670267026916
mass added [Msun]: 4.874000556794949e-10
relMass_added/rstep (rico) = 10.735
---------
r 0.17 au
scale height is 0.029552955295529617
mass added [Msun]: 4.950756471075184e-10
relMass_added/rstep (rico) = 9.832
---------
r 0.18 au
scale height is 0.0313531353135314
mass added [Msun]: 4.547787921103949e-10
relMass_added/rstep (rico) = 8.284
---------
r 0.19 au
scale height is 0.034053405340534
mass added [Msun]: 4.720488728234478e-10
relMass_added/rstep (rico) = 7.918
---------
r 0.2 au
scale height is 0.03735373537353724
mass added [Msun]: 4.4902209853937726e-10
relMass_added/rstep (rico) = 7.004
---------
r 0.21 au
scale height is 0.03945394539453934
mass added [Msun]: 4.3558981354033604e-10
relMass_added/rstep (rico) = 6.362
---------
r 0.22 au
scale height is 0.043054305430543165
mass added [Msun]: 4.4902209853937726e-10
relMass_added/rstep (rico) = 6.155
---------
r 0.23 au
scale height is 0.04620462046204619
mass added [Msun]: 4.4134650711135367e-10
relMass_added/rstep (rico) = 5.704
---------
r 0.24 au
scale height is 0.048604860486048604
mass added [Msun]: 4.547787921103949e-10
relMass_added/rstep (rico) = 5.552
---------
r 0.25 au
scale height is 0.05205520552055194
mass added [Msun]: 4.2791422211231246e-10
relMass_added/rstep (rico) = 4.964
---------
r 0.26 au
scale height is 0.054755475547554805
mass added [Msun]: 4.144819371132713e-10
relMass_added/rstep (rico) = 4.588
---------
r 0.27 au
scale height is 0.05880588058805878
mass added [Msun]: 4.3942760925434776e-10
relMass_added/rstep (rico) = 4.638
---------
r 0.28 au
scale height is 0.06105610561056111
mass added [Msun]: 3.780228778301594e-10
relMass_added/rstep (rico) = 3.837
---------
r 0.29 au
scale height is 0.06345634563456352
mass added [Msun]: 3.799417756871653e-10
relMass_added/rstep (rico) = 3.713
---------
r 0.3 au
scale height is 0.06675667566756677
mass added [Msun]: 3.5883389926010053e-10
relMass_added/rstep (rico) = 3.388
---------
r 0.31 au
scale height is 0.06930693069306929
mass added [Msun]: 3.262126356910005e-10
relMass_added/rstep (rico) = 2.988
---------
r 0.32 au
scale height is 0.07260726072607265
mass added [Msun]: 3.089425549779475e-10
relMass_added/rstep (rico) = 2.752
---------
r 0.33 au
scale height is 0.07455745574557453
mass added [Msun]: 3.051047592639357e-10
relMass_added/rstep (rico) = 2.646
---------
r 0.34 au
scale height is 0.0777077707770778
mass added [Msun]: 2.705645978378298e-10
relMass_added/rstep (rico) = 2.293
---------
r 0.35 au
scale height is 0.08100810081008092
mass added [Msun]: 2.417811299827416e-10
relMass_added/rstep (rico) = 2.008
---------
r 0.36 au
scale height is 0.08430843084308429
mass added [Msun]: 2.4945672141076507e-10
relMass_added/rstep (rico) = 2.029
---------
r 0.37 au
scale height is 0.0889588958895889
mass added [Msun]: 2.3026774284070628e-10
relMass_added/rstep (rico) = 1.839
---------
r 0.38 au
scale height is 0.09270927092709269
mass added [Msun]: 2.1107876427064736e-10
relMass_added/rstep (rico) = 1.658
---------
r 0.39 au
scale height is 0.09555955595559565
mass added [Msun]: 2.0532207069962974e-10
relMass_added/rstep (rico) = 1.587
---------
r 0.4 au
scale height is 0.09810981098109803
mass added [Msun]: 1.784575007015473e-10
relMass_added/rstep (rico) = 1.36
---------
r 0.41 au
scale height is 0.1015601560156015
mass added [Msun]: 1.9572758141460027e-10
relMass_added/rstep (rico) = 1.47
---------
r 0.42 au
scale height is 0.10471047104710465
mass added [Msun]: 1.4391733927544137e-10
relMass_added/rstep (rico) = 1.069
---------
r 0.43 au
scale height is 0.11011101110111012
mass added [Msun]: 1.6886301141651792e-10
relMass_added/rstep (rico) = 1.239
---------
r 0.44 au
scale height is 0.11491149114911482
mass added [Msun]: 1.4583623713244726e-10
relMass_added/rstep (rico) = 1.059
---------
r 0.45 au
scale height is 0.1182118211821182
mass added [Msun]: 1.4775513498945312e-10
relMass_added/rstep (rico) = 1.061
---------
r 0.46 au
scale height is 0.12256225622562263
mass added [Msun]: 1.2472836070538253e-10
relMass_added/rstep (rico) = 0.888
---------
r 0.47 au
scale height is 0.12736273627362746
mass added [Msun]: 1.3816064570442375e-10
relMass_added/rstep (rico) = 0.974
---------
r 0.48 au
scale height is 0.1329132913291329
mass added [Msun]: 1.1129607570634133e-10
relMass_added/rstep (rico) = 0.779
---------
r 0.49 au
scale height is 0.13771377137713786
mass added [Msun]: 1.0937717784933546e-10
relMass_added/rstep (rico) = 0.759
---------
r 0.5 au
scale height is 0.1431143114311431
mass added [Msun]: 1.0170158642131189e-10
relMass_added/rstep (rico) = 0.701
---------
r 0.51 au
scale height is 0.14866486648664878
mass added [Msun]: 9.210709713628251e-11
relMass_added/rstep (rico) = 0.631
---------
r 0.52 au
scale height is 0.1539153915391539
mass added [Msun]: 8.251260785125306e-11
relMass_added/rstep (rico) = 0.562
---------
r 0.53 au
scale height is 0.1578157815781579
mass added [Msun]: 9.402599499328838e-11
relMass_added/rstep (rico) = 0.636
---------
r 0.54 au
scale height is 0.1627662766276627
mass added [Msun]: 8.443150570825893e-11
relMass_added/rstep (rico) = 0.568
---------
r 0.55 au
scale height is 0.16756675667566753
mass added [Msun]: 6.908032285221189e-11
relMass_added/rstep (rico) = 0.463
---------
r 0.56 au
scale height is 0.17476747674767465
mass added [Msun]: 7.675591428023542e-11
relMass_added/rstep (rico) = 0.512
---------
r 0.57 au
scale height is 0.18271827182718264
mass added [Msun]: 6.908032285221189e-11
relMass_added/rstep (rico) = 0.458
---------
r 0.58 au
scale height is 0.18901890189018894
mass added [Msun]: 6.7161424995206e-11
relMass_added/rstep (rico) = 0.444
---------
r 0.59 au
scale height is 0.1968196819681967
mass added [Msun]: 6.332362928119422e-11
relMass_added/rstep (rico) = 0.417
---------
r 0.6 au
scale height is 0.2046204620462046
mass added [Msun]: 6.7161424995206e-11
relMass_added/rstep (rico) = 0.44
---------
r 0.61 au
scale height is 0.20987098709870972
mass added [Msun]: 5.756693571017657e-11
relMass_added/rstep (rico) = 0.376
---------
r 0.62 au
scale height is 0.21782178217821782
mass added [Msun]: 6.908032285221189e-11
relMass_added/rstep (rico) = 0.449
---------
r 0.63 au
scale height is 0.22682268226822686
mass added [Msun]: 4.6053548568141254e-11
relMass_added/rstep (rico) = 0.298
-------------------------------------------------------------
estimate of r of the accretion disk is 0.63 au
Other possibilities:
Do calculations for different dumps throughout simulation to see how disk mass evolves
Do calculations for different dumps within one eccentric orbit to see how disk properties change depending on orbital phase
If you want something else, check if there is not yet a function for it in AccrDisk.py in the source code