smmills / Phodymm

Photodynamics Code v1
2 stars 2 forks source link

running with out mass prio #28

Closed petigura closed 5 years ago

petigura commented 5 years ago

Hey @smmills,

Are you sure the instructions here are right?

// Do you wanto enforce a Gaussian mass ratio prior on the planets (1=yes, 2=no)
// If yes, the following line specifies first the center of the prior for each planets, and the next the uncertainty (sigma) of the Gaussian
// Thus each each line should have Nplanet entries, and both lines must be present even if a mass prior is not used (but will be ignored) 
int massprior= 0
0.     112e-6   1e-7
11e-6  17e-6    34e-6

I set massprior= 2 and the code was still echoing out the mass prior. I set it to zero (also seen in some other examples) and the code did not echo out any mass prior.

smmills commented 5 years ago

Yes, that's a typo, should be 0=no. I will fix that in the .in file descriptions.

smmills commented 5 years ago

Fixed in 9951f3fd78f90177f141ee05f1a9d60b4328e4fc. It will also not let you proceed if massprior is not 0 or 1, so setting it equal to 2 will fail.