.. _1.4M_ms_op_mono:

***************
1.4M_ms_op_mono
***************

.. tags:: star, kap, opacity, op-mono, low-mass, main-sequence

The test checks the functionality of OP mono opacities. The test vehicle is a
1.4 |Msun| solar metallicity model.


This test case has three parts. Click to see a larger view of a plot.

* Part 1 (``inlist_1.4M_ms_initial_model``) creates the pre-main-sequence model and evolves the model until the central hydrogen mass fraction falls below 0.5. During this evolution the OP mono opacities are off.

.. image:: ../../../star/test_suite/1.4M_ms_op_mono/docs/initial_model.svg
   :width: 100%


The upper plot shows the fraction of the opacity that is from OP mono (yellow curve) and the opacity (blue curve)
as a function of the temperature profile - center on right, surface on left.
The lower plot shows the opacity derivatives.


* Part 2 (``inlist_1.4M_ms_op_mono``) continues the evolution of ``inlist_1.4M_ms_initial_model`` for 10 timesteps using the OP mono opacities, as shown by the fraction of the opacity that is from OP mono (yellow curve, upper plot)

.. image:: ../../../star/test_suite/1.4M_ms_op_mono/docs/op_mono.svg
   :width: 100%


* Part 3 (``inlist_1.4M_ms_op_mono``) continues the evolution of ``inlist_1.4M_ms_initial_model`` for 10 timesteps using an alternative formulation of the OP mono opacities, as shown by the smoother opacity derivatives in the lower plot. See $MESA_DIR/kap/public/kap_lib.f.90 for the differences between the regular and alternative OP mono formulations.

.. image:: ../../../star/test_suite/1.4M_ms_op_mono/docs/op_mono_alt.svg
   :width: 100%


pgstar commands used for the plots above:

.. code-block:: console

 &pgstar

   pgstar_interval = 10

   file_white_on_black_flag = .true.
   file_device = 'vcps'

  ! two profile panels
    Profile_Panels1_win_flag = .true.

    Profile_Panels1_win_width = 15
    Profile_Panels1_win_aspect_ratio = 0.75

    Profile_Panels1_xleft = 0.15
    Profile_Panels1_xright = 0.85
    Profile_Panels1_ybot = 0.15
    Profile_Panels1_ytop = 0.85
    Profile_Panels1_txt_scale = 1.0
    Profile_Panels1_title = '1.4_ms_op_mono'

    Profile_Panels1_num_panels = 2

    Profile_Panels1_xaxis_name = 'logT'
    Profile_Panels1_xaxis_reversed = .false.
    Profile_Panels1_xmin = -101d0
    Profile_Panels1_xmax = -101d0

    Profile_Panels1_yaxis_name(1) = 'kap_frac_op_mono'
    Profile_Panels1_other_yaxis_name(1) = 'log_opacity'

    Profile_Panels1_yaxis_name(2) = 'dlnkap_dlnd'
    Profile_Panels1_other_yaxis_name(2) = 'dlnkap_dlnT'

    Profile_Panels1_file_flag = .true.
    Profile_Panels1_file_dir = 'pgstar_out'
    Profile_Panels1_file_prefix = 'profile_panels1_'
    Profile_Panels1_file_interval = 10
    Profile_Panels1_file_width = 15
    Profile_Panels1_file_aspect_ratio = -1

 / ! end of pgstar namelist


Last-Updated: 27May2021 (MESA ebecc10) by fxt