From: Yevgeni K. <ki...@as...> - 2017-08-02 16:36:20
|
Hi Frank, Josiah, and Bill, Frank, thank you for the suggestion. The inlist ran as expected after I made that change. Josiah, thank you for clarifying what the fundamental issue was. It was unfortunate not to receive an error message, which is why I was confused. But it all makes sense now. Bill, I don't imagine anybody is blaming you. Thank you for all your work creating and maintaining MESA. Best, Yevgeni Kissin PhD Canadian Institute for Theoretical Astrophysics University of Toronto On 31 July 2017 at 14:08, Frank Timmes <fx...@ma...> wrote: > most excellent catch on "/" being a reserved character! > > fxt > > > > > On Jul 31, 2017, at 9:47 AM, Josiah Schwab <jws...@uc...> wrote: > > > > Hi All, > > > > I wanted to take a minute and pedagogically describe what is going on > > here for future reference. > > > > The first, and most important point is that you can't do math in > > inlists. So as Frank says, the solution is to write out your numbers. > > > > Now if you were to try to do some other math operation, for example, > > > > am_D_mix_factor = 1d0 + 1d0 > > > > you'd receive an error. > > > > However, the character '/' ends a fortran namelist. That is, one > > normally writes something like > > > > &controls > > > > am_D_mix_factor = 0.03333 > > > > / ! end of controls namelist > > > > where the comment is informative, but syntactically irrelevant. > > > > So when you write > > > > &controls > > > > am_D_mix_factor = 1d0/30d0 > > some_other_option = 'something' > > > > / ! end of controls namelist > > > > the namelist ends immediately at the first '/'. So everything beyond > > (i.e. the 30d0 and any subsequent options) will be ignored, because it > > is located outside of the controls namelist. > > > > As you found, this mistake can be insidious because it will not produce > > an error message. MESA will simply run with the truncated inlist. > > > > Josiah > > > ------------------------------------------------------------ > ------------------ > Check out the vibrant tech community on one of the world's most > engaging tech sites, Slashdot.org! http://sdm.link/slashdot > _______________________________________________ > mesa-users mailing list > mes...@li... > https://lists.sourceforge.net/lists/listinfo/mesa-users > |