Re: Phantom Binary Accretion Rates

[Daniel Price]

Hi Jeremy,

I recently completely rewrote all the ev2mdot (ev2mdot, ev2mdotsink, ev2mdotbinary) utilities and combined them into a single working version (ev2mdot). The problem is that they needed to know which column is which and we frequently change this :)

Now the utility scans the column labels to identify the correct columns, and the same program works for all .ev files.

Best wishes,

Daniel

> On 7 Sep 2017, at 12:19 pm, Jeremy Smallwood <xxxu\> wrote:
>
> I have a question regarding the calculation of accretion rates in a system with a circumbinary disc (I have attached my disc.setup). Following your instructions on the wiki (using ev2mdot), I only obtain the first three columns (a sample shown before) instead of seven columns (time mdot macc mdot1 macc1 mdot2 macc2). Plus the macc is a negative value at times. I used discSink0001N01.ev and discSink0002N01.ev.
>
> # t mdot macc
> 1.0001554740E+03 2.2993177999E-05 2.2996752840E-02
> 2.0002189090E+03 0.0000000000E+00 -6.8004137020E-01
> 3.0002363250E+03 6.9895698097E-05 -6.1014445480E-01
> 4.0002537410E+03 6.1703203544E-04 6.8983268720E-03
> 5.0002864970E+03 0.0000000000E+00 -2.2882316670E-01
> 6.0003345920E+03 1.2854242545E-04 -1.0027455900E-01
>
> We have used an old version of ev2mdotbinary.f90, but this also has negative values for macc, macc1, and macc2, plus the command for obtaining smoother accretion rates does not work (example: ./ev2mdotbinary {value} discSink000*.ev). This removes all points less than {value} and plots all remaining points greater than {value}.
>
>
>
> <disc.setup>
> <ev2mdotbinary.f90>
>
> Let me know if you need any further information.
>
> Best,
>
> Jeremy Smallwood
> UNLV