Re: [mesa-users] irradiation

[Jean-Claude P. <jc...@gm...>]

Hi,

no you need both. The 'relax_irradiation' helps you getting started by
relaxing the star as you turn on irradiation progressively. The section
in your inlist should look like:

----------------------------------------------------
set_irradiation = .true.

relax_irradiation = .true.
relax_irradiation_min_steps = 200
relax_irradiation_max_yrs_dt = 1d-1

new_irrad_flux = 1d7 ! erg s^-1 cm^-2
new_irrad_col_depth = 1 ! g cm^-2
----------------------------------------------------

Hope that helps.
Good luck,

JC


On 22/07/11 21:59, ztt wrote:
> Hi Jean,
>
> Thanks for your reply. My mail box didn't get your email but I saw it
> on the mailing list.
> So is this the correct procedure example for turning on irradiation??
> 1. do relax with:
> relax_irradiation = .true.
> relax_irradiation_min_steps = 200
> relax_irradiation_max_yrs_dt = 1d-1
>
> new_irrad_flux = 1d7 ! erg s^-1 cm^-2
> new_irrad_col_depth = 1 ! g cm^-2
> with 'set_irradiation = .true.' commented;
> 2. turn on irradiation with:
> set_irradiation = .true.
> new_irrad_flux = 1d7 ! erg s^-1 cm^-2
> new_irrad_col_depth = 1 ! g cm^-2
> but with all the relax_* commented now;
>
> while keep the setup for tau & atmosphere all the time:
> set_tau_factor = .true. ! change tau_factor without reconverging.
> set_to_this_tau_factor = 100
>
> which_atm_option = 'grey_and_kap'
> !which_atm_option = 'grey_irradiated'
>
> atm_grey_irradiated_atol = 1d-4
> atm_grey_irradiated_rtol = 1d-4
> atm_grey_irradiated_T_eq = 1000 ! equilibrium temperature based on
> irradiation
> ! e.g., might use T_eq = Teff_star*(R_star/(2*Distance_to_star))^(1/2)
> atm_grey_irradiated_kap_nu = 4d-3 ! opacity for irradiation
> atm_grey_irradiated_P_surf = 1d6 ! 1 bar in cgs units
> &nb! sp; ; atm_grey_irradiated_max_tries = 50 ! limit on iterations
> trace_atm_grey_irradiated = .false.
>
> Thank you very much!
>
> Best wishes, Tingtao
>
>