[mesa-users] Brunt-Vaisala frequency in convective regions.

[Daniel Kagan]

I am trying to calculate the value of N^2 in convective regions of a stellar 
model; while this will be negative, both of the formulas for the Brunt-Vaisala 
frequency ought to work there. However, Mesa doesn't seem to use these formulae 
in convective regions, and N^2 (as opposed to log(N^2))
is not a standard output. Is there a way to force Mesa to do this calculation? 
I have tried to do it manually (without success), but is there an easier way?

Thanks,
Daniel Kagan