[mesa-users] irradiation

[ztt <ztt...@12...>]

  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
            atm_grey_irradiated_max_tries = 50 ! limit on iterations
            trace_atm_grey_irradiated = .false.

Thank you very much!

Best wishes, Tingtao