[mesa-users] Fwd: I/O in MESA
paxton at kitp.ucsb.edu
Fri Jul 12 17:58:50 EDT 2013
Anyone want to help Manos on this?
Begin forwarded message:
> From: Manos Chatzopoulos <manolis at astro.as.utexas.edu>
> Date: July 12, 2013 1:27:43 PM PDT
> To: Bill Paxton <paxton at kitp.ucsb.edu>
> Subject: I/O in MESA
> Hi Bill,
> I have a quick question on MESA basic coding. I want to create a new variable
> that is a simple function of other variables calculated in the rotation_mix_info.f
> subroutine (specifically, eta, kappa and N2_T). I want, upon re-compilation, this variable
> to be calculated for every model, everywhere within the star and then eventually saved
> as an output in the profileXX.data files.
> How do I do that?
> Here is a relevant piece of code:
> kap = 16*boltz_sigma*T(k)**3/(3*opacity(k)*rho(k)**2*Cp(k)) ! thermal diffusivity
> q = max(1d-30,min(1d30,abs(domega_dlnr(k)/omega(k)))) ! shear
> s% omega_shear(k) = q
> eta = 7d11*lnLambda/T(k)**1.5d0 ! magnetic diffusivity -- eqn 5
> eta_div_r2_omega = eta/(r(k)**2*omega(k))
> How do I add my own variable that is then recognized by the whole MESA tree (where do I make
> the definitions etc).
> Thanks a lot! This is important help.
> Dr. Manos Chatzopoulos
> Department of Astronomy
> University of Texas at Austin
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Mesa-users