[mesa-users] Fwd: I/O in MESA

[Bill Paxton]

Anyone want to help Manos on this?

Thanks,
Bill



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.
> 
> Manos
> 
> -- 
> Dr. Manos Chatzopoulos
> Department of Astronomy
> University of Texas at Austin
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20130712/844cb888/attachment.html>