[Mesa-users] Post-processing trajectories in MESA networks

[Rob Farmer]

tPlease keep mesa-users in any replies

large_mtx_decsol = ‘klu’ sets the choice of the matrix solver. So lets find
out what options mesa has available:

First find out where large_mtx_decsol is used
cd net/test/
grep -rinI large_mtx_decsol src/*
src/mod_one_zone_burn.f:78:      character (len=32) :: small_mtx_decsol,
large_mtx_decsol
src/mod_one_zone_burn.f:506:            decsol_choice =
decsol_option(large_mtx_decsol, ierr)
src/mod_one_zone_burn.f:508:               write(*,*) 'ERROR: unknown
large_mtx_decsol ' // trim(large_mtx_decsol)
src/mod_one_zone_burn.f:1506:         decsol_switch, small_mtx_decsol,
large_mtx_decsol, &
src/mod_one_zone_burn.f:1561:            ! else use large_mtx_decsol.
src/mod_one_zone_burn.f:1563:      large_mtx_decsol = ''

If you look over each instance that uses large_mtx_decsol then you find:

         if (species >= decsol_switch) then
            decsol_choice = decsol_option(large_mtx_decsol, ierr)
            if (ierr /= 0) then
               write(*,*) 'ERROR: unknown large_mtx_decsol ' //
trim(large_mtx_decsol)
               return
            end if
         else

Is where your error message is coming from. Let's find the decsol_option
function and see what options it allows. Searching in the
src/mod_one_zone_burn.f file shows this line:

use mtx_lib, only: decsol_option

Which means the function is defined in the mtx_lib file so we need to
switch files and find the mtx_lib file. That will be in
mtx/public/mtx_lib.f90. In that file you'll find:

      integer function decsol_option(which_decsol_option, ierr)
         use mtx_def
         character (len=*), intent(in) :: which_decsol_option
         integer, intent(out) :: ierr
         character (len=64) :: option
         ierr = 0
         option = which_decsol_option

         if (option == 'lapack') then
            decsol_option = lapack

         else if (option == 'bcyclic_dble') then
            decsol_option = bcyclic_dble

         else
            ierr = -1
            decsol_option = -1
         end if
      end function decsol_option

So it seems that large_mtx_decsol can only be 'lapack' or 'bcyclic_dble'
and not 'klu'. So klu was an old option that got must of been deleted and
not cleaned up everywhere. So you have a choice of setting
large_mtx_decsol to either 'lapack' or 'bcyclic_dble'

What you are finding here is that the one zone burn code in mesa is not
well tested. So you may find other tools better suited, and better tested,
for instance frank timmes's torch code
http://cococubed.asu.edu/code_pages/net_torch.shtml

Rob


On Thu, 1 Apr 2021 at 17:24, Khanak Bhargava <kbhargava at umassd.edu> wrote:

> Hello Rob,
>
>
>
> Sorry for the delayed response. I tried to use the default
> inlist_one_zone_burn with read_T_Rho_History = .true. and mesa_495.net. I
> was finally able to get it to read my temperature-density history file but
> I am getting an error with a parameter called large_mtx_decsol = ‘klu’.
> It says unknown large_mtx_decsol for some reason.
>
>
>
> I am not sure what this parameter does since I am new to MESA. Would you
> happen to know anything about this?
>
>
>
> Thanks so much,
>
> Khanak
>
>
>
> *From: *Rob Farmer <r.j.farmer at uva.nl>
> *Date: *Wednesday, March 31, 2021 at 4:09 AM
> *To: *Khanak Bhargava <kbhargava at umassd.edu>,
> mesa-users at lists.mesastar.org <mesa-users at lists.mesastar.org>
> *Subject: *Re: Post-processing trajectories in MESA networks
>
> [EXTERNAL SENDER]
>
> Hi,
>
> There is not enough information here to help you. How are you using
> inlist_one_zone_burn? Did you change anything? How did you set up the
> inlist to work with your data? What file does it say is unavailable?
>
>
>
> Rob
>
>
> ------------------------------
>
> *From:* Mesa-users <mesa-users-bounces at lists.mesastar.org> on behalf of
> Khanak Bhargava via Mesa-users <mesa-users at lists.mesastar.org>
> *Sent:* 28 March 2021 20:41
> *To:* mesa-users at lists.mesastar.org <mesa-users at lists.mesastar.org>
> *Subject:* [Mesa-users] Post-processing trajectories in MESA networks
>
>
>
> Hello MESA users,
>
>
>
> Could someone guide me on how to run a temperature-density history file in
> MESA’s networks?
>
>
>
> I am new to MESA and I am trying to post-process temperature-density
> trajectories using the 495 isotope network in MESA. I tried to use the
> inlist_one_zone_burn with read_T_Rho_history = .true. but it gives an error
> that says unable to open file.
>
>
>
> Best,
>
> Khanak Bhargava
>
> Department of Physics
>
> University of Massachusetts Dartmouth
>
>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210407/273c0f8c/attachment.htm>