[Mesa-users] More resolution at core boundary ?
Conny Aerts
conny.aerts at kuleuven.be
Mon Feb 1 13:37:24 UTC 2021
Dear Giovanni,
regarding your question on how to treat the convective core boundary in mesa for massive stars and your reference to a thread by Ehsan Moravveji, I disadvise you to go back so far in time. Mesa has improved a lot in the treatment of the boundaries and major tests have been done to set the resolutions optimally for asteroseismology of stars with a convective core. I suggest you read the papers by
Pedersen et al. (2018): https://www.aanda.org/articles/aa/full_html/2018/06/aa32317-17/aa32317-17.html
Michielsen et al. (2019): https://www.aanda.org/articles/aa/full_html/2019/08/aa35754-19/aa35754-19.html
You can download the inlist of the Michielsen et al. paper from the mesa market place: http://cococubed.asu.edu/mesa_market/inlists.html
Although things have moved on again meanwhile compared to v10398 used by Mathias, he will surely be prepared to help you with the current versions of the software to set things optimally for your case (p modes or g modes have different requirements).
Best regards,
Conny
On 1 Feb 2021, at 13:00, mesa-users-request at lists.mesastar.org<mailto:mesa-users-request at lists.mesastar.org> wrote:
Send Mesa-users mailing list submissions to
mesa-users at lists.mesastar.org<mailto:mesa-users at lists.mesastar.org>
To subscribe or unsubscribe via the World Wide Web, visit
https://lists.mesastar.org/mailman/listinfo/mesa-users
or, via email, send a message with subject or body 'help' to
mesa-users-request at lists.mesastar.org
You can reach the person managing the list at
mesa-users-owner at lists.mesastar.org
When replying, please edit your Subject line so it is more specific
than "Re: Contents of Mesa-users digest..."
Today's Topics:
1. About the warning issue during installing MESA15140 (Hongwei Ge)
2. Re: About the warning issue during installing MESA15140
(Rob Farmer)
3. More resolution at core boundary ? (GM Mirouh)
----------------------------------------------------------------------
Message: 1
Date: Mon, 1 Feb 2021 17:12:01 +0800
From: Hongwei Ge <hongwei.ge at gmail.com>
To: mesa-users at lists.mesastar.org
Subject: [Mesa-users] About the warning issue during installing
MESA15140
Message-ID:
<CAFL9L2sRnddJHiJA5B1Kd0AmeRb6Z2jVrm8_PUr3Gx4YY-AxDw at mail.gmail.com>
Content-Type: text/plain; charset="utf-8"
Hi MESA developers and users:
As a beginner, I am recently reading through MESA15140.
I just found some warning message in my basic.log file, such as:
../private/pre_conditioners.f90:81:81:
71 | do j=1,nblocks
| 2
......
81 | pre_conditioner(k,j) = pre_conditioner(k,j) +
abs(ublk(k2,k,j-1) / right_pre_conditioner(k2,j-1))
|
1
Warning: Array reference at (1) out of bounds (0 < 1) in loop beginning at
(2) [-Wdo-subscript]--
There are about 6 places for such kind of warning about the Array bounds.
The code can run smoothly actually, but should we just neglect these
warnings or should we modify the array?
With best regards,
Hongwei
Yunnan Observatories, Chinese Academy of Sciences,
NO. 396, YangFangWang, Guandu District,
Kunming, 650216, China
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210201/2a68a486/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: build.log
Type: application/octet-stream
Size: 44349 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210201/2a68a486/attachment-0001.obj>
------------------------------
Message: 2
Date: Mon, 1 Feb 2021 10:23:54 +0100
From: Rob Farmer <r.j.farmer at uva.nl>
To: Hongwei Ge <hongwei.ge at gmail.com>
Cc: "mesa-users at lists.mesastar.org" <mesa-users at lists.mesastar.org>
Subject: Re: [Mesa-users] About the warning issue during installing
MESA15140
Message-ID:
<CAPvTSybODEGZFuo7efqZOCfwabx2DGbxgg3Ea7qMTT2p8cKVBA at mail.gmail.com>
Content-Type: text/plain; charset="utf-8"
Hi
Yes they can be ignored:
find . -name pre_conditioners.f90
./mtx/private/pre_conditioners.f90
Lines 80-82
if (j > 1) then
pre_conditioner(k,j) = pre_conditioner(k,j) +
abs(ublk(k2,k,j-1) / right_pre_conditioner(k2,j-1))
end if
Note the if block that means j must be at least 2 and thus j-1 >=1 keeping
the index in bounds.
Rob
On Mon, 1 Feb 2021 at 10:13, Hongwei Ge via Mesa-users <
mesa-users at lists.mesastar.org> wrote:
Hi MESA developers and users:
As a beginner, I am recently reading through MESA15140.
I just found some warning message in my basic.log file, such as:
../private/pre_conditioners.f90:81:81:
71 | do j=1,nblocks
| 2
......
81 | pre_conditioner(k,j) = pre_conditioner(k,j) +
abs(ublk(k2,k,j-1) / right_pre_conditioner(k2,j-1))
|
1
Warning: Array reference at (1) out of bounds (0 < 1) in loop beginning at
(2) [-Wdo-subscript]--
There are about 6 places for such kind of warning about the Array bounds.
The code can run smoothly actually, but should we just neglect these
warnings or should we modify the array?
With best regards,
Hongwei
Yunnan Observatories, Chinese Academy of Sciences,
NO. 396, YangFangWang, Guandu District,
Kunming, 650216, China
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210201/443665bd/attachment-0001.htm>
------------------------------
Message: 3
Date: Mon, 1 Feb 2021 11:07:19 +0000
From: GM Mirouh <g.mirouh at surrey.ac.uk>
To: mesa-users at lists.mesastar.org
Subject: [Mesa-users] More resolution at core boundary ?
Message-ID: <033c28db-1fe4-789a-f34c-e196a13e0b16 at surrey.ac.uk>
Content-Type: text/plain; charset="utf-8"; Format="flowed"
Dear all,
I am computing MESA models of upper main-sequence stars and I want to
locate and describe the interface between the convective core and the
radiative envelope precisely. What is the best way to increase the
number of mesh points locally at the interface?
I have MESA 11701 and 12115 installed. I have read through the default
files, but this thread here has got me confused:
https://sourceforge.net/p/mesa/mailman/message/34628504/ .
Here is an example of my work: with MESA 11701, I have computed a 17Msun
model with default mesh settings: my model has 1064 shells and the
core-envelope interface is between the mesh points at m=6.46 and m=6.54.
I would like to resolve this 0.08 Msun shell much better than this, with
about 100 mesh points there. I don?t need more resolution elsewhere
though, and ideally I would like to simply restart from my model file,
simply rederive my variables on a new mesh and save that in a profile file.
Thank you for any help you can provide,
Best.
--
Gio.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210201/d1b342e0/attachment-0001.htm>
End of Mesa-users Digest, Vol 43, Issue 1
*****************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210201/d63a06dc/attachment.htm>
More information about the Mesa-users
mailing list