include convective driving in gyre
Posted: Fri Jun 20, 2025 4:21 pm
Hi,
I am working with the driving of g-modes in DAVs.. now, the way that nad works is by 'frozen convection' where FROZEN_PESNELL_4 gives the driving but it is ignoring the lagrangian perturbation of the convective flux.
while scheme1
conv_term = lambda*c_rad*(3._WP + dc_rad)/(c_1*alpha_omg*omega_c**2)
completely ignores the convective flux,
scheme4
conv_term = lambda*(c_lum*(3._WP + dc_lum) - (c_lum - c_rad))/(c_1*alpha_omg*omega_c**2)
just ignores the lagrangian variation of the convective flux.
I was wondering, if there was any way to include the lagrangian variation of the convective flux to the equations. I tried by modifying the second equation in gyre_nad_eqns.fpp, to
conv_term = lambda*c_lum*(3._WP + dc_lum)/(c_1*alpha_omg*omega_c**2)
but that is not producing any different result than the original one (scheme4)..
can I get some help on how to go about this one so that I can include that term in calculation.
Thank you!
I am working with the driving of g-modes in DAVs.. now, the way that nad works is by 'frozen convection' where FROZEN_PESNELL_4 gives the driving but it is ignoring the lagrangian perturbation of the convective flux.
while scheme1
conv_term = lambda*c_rad*(3._WP + dc_rad)/(c_1*alpha_omg*omega_c**2)
completely ignores the convective flux,
scheme4
conv_term = lambda*(c_lum*(3._WP + dc_lum) - (c_lum - c_rad))/(c_1*alpha_omg*omega_c**2)
just ignores the lagrangian variation of the convective flux.
I was wondering, if there was any way to include the lagrangian variation of the convective flux to the equations. I tried by modifying the second equation in gyre_nad_eqns.fpp, to
conv_term = lambda*c_lum*(3._WP + dc_lum)/(c_1*alpha_omg*omega_c**2)
but that is not producing any different result than the original one (scheme4)..
can I get some help on how to go about this one so that I can include that term in calculation.
Thank you!