[mesa-users] optical_depth -- and a message for MESA BLACK BELT USERS
Aaron Dotter
aaron.dotter at gmail.com
Fri Mar 23 22:11:05 EDT 2012
Hi Jing,
I took a look at the model you sent and the numbers that you're getting for
tau are reasonable; there's nothing wrong with the calculation. The
physical explanation for why tau grows so quickly is that you're looking at
a white dwarf with log(g)~8.
The reason your "tau_phot" is different from what MESA reports in s% tau(:)
is that you're not taking into account that some variables are "face"
variables and some are "cell" mass-average, but MESA is. Also, if you look
at the way s% photosphere_r is defined, you'll see that it's s% r(1)/Rsun
unless you intentionally modify tau_factor in the star_job inlist. That is
to say that cell #1 is the "photosphere." In your model, if you output s%
tau(1) you should get 2/3 (or whatever is appropriate for your atmosphere
BC).
Aaron
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120323/f037dd5f/attachment.html>
More information about the Mesa-users
mailing list