[mesa-users] identity of 'v' in 'wd_ignite.mod'

Sumit Sarbadhicary s2ksumit at gmail.com
Sat Apr 4 22:31:30 EDT 2015


Robert,

Thanks a lot for clarifying the identity of 'v'. That makes sense!

Do you, or any of the MESA users, have an idea why the WD_IGNITE.MOD have
two parts? The first few lines of the first part looks like this :-

! note: initial lines of file can contain comments

!

            12 -- two consecutive models for mesa/star, with velocity
variables (v). cgs units. lnd=ln(density), lnT=ln(temperature),
lnR=ln(radius), L=lumin\

osity, dq=fraction of total mstar in cell; remaining cols are mass
fractions.


     year_month_day_when_created   20110717

                          M/Msun      1.0000000009434213E+00

                        star_age      9.4342161767044841E-01

                       initial_z      2.0000000000000000E-02

                        n_shells                            1273

                        net_name   'co_burn.net'

                         species                               9

                      tau_factor      1.5000000000000000E+02


                lnd                        lnT                        lnR
                        L                         dq                       *
v*   \

                h1                         he3                        he4
                      c12                        n14
o16   \

                     ne20                       mg24
si28

    1    -1.0783032764354498E+01     1.1640541494135320E+01
2.0135894106074399E+01     4.8571750122419371E+32
1.4075930508271764E-14     1.1361104977\

517480E-02     0.0000000000000000E+00     0.0000000000000000E+00
0.0000000000000000E+00     2.5000000000000000E-01
0.0000000000000000E+00     7.50000\

00000000000E-01     0.0000000000000000E+00     0.0000000000000000E+00
0.0000000000000000E+00

    2    -8.2810274680451421E+00     1.2428434759429717E+01
2.0135268008479237E+01     4.8571750981602734E+32
1.1250008967110072E-14    -4.6362609302\

080995E+02     0.0000000000000000E+00     0.0000000000000000E+00
0.0000000000000000E+00     2.5000000000000000E-01
0.0000000000000000E+00     7.50000\

00000000000E-01     0.0000000000000000E+00     0.0000000000000000E+00
0.0000000000000000E+00

....(more lines)

-------------------------------------------------------------------------------------------------------


and then the second part continues right after model #1273, as follows :-

-------------------------------------------------------------------------------------------------------


        previous model


               previous n_shells                            1273

           previous mass (grams)      1.9892000016223856E+33

              timestep (seconds)      4.0339007298227884E+06

               dt_next (seconds)      4.8406808757873457E+06


                lnd                        lnT                        lnR
                        L                         dq
  v   \

                h1                         he3                        he4
                      c12                        n14
o16   \

                     ne20                       mg24
si28

    1    -1.0783032764354498E+01     1.1640541494135320E+01
2.0135894106074399E+01     4.8571750122419371E+32
1.2052633511987234E-14    -3.2940374214\

366158E-02     1.0000000000000000E-99     1.0000000000000000E-99
1.0000000000000000E-99     2.5000000000000000E-01
1.0000000000000000E-99     7.50000\

00000000000E-01     1.0000000000000000E-99     1.0000000000000000E-99
1.0000000000000000E-99

    2    -8.2810274680451421E+00     1.2428434759429717E+01
2.0135268008479237E+01     4.8571750981602734E+32
9.6329144995044247E-15    -4.6808781344\

515745E+02     1.0000000000000000E-99     1.0000000000000000E-99
1.0000000000000000E-99     2.5000000000000000E-01
1.0000000000000000E-99     7.50000\

00000000000E-01     1.0000000000000000E-99     1.0000000000000000E-99
1.0000000000000000E-99

    3    -7.8621441999255000E+00     1.2633241562558258E+01
2.0135226974337208E+01     4.8522180474116445E+32
1.0000000000000000E-14    -5.5665419796\

314268E+01     1.0000000000000000E-99     1.0000000000000000E-99
1.0000000000000000E-99     2.5000000000000000E-01
1.0000000000000000E-99     7.50000\

00000000000E-01     1.0000000000000000E-99     1.0000000000000000E-99
1.0000000000000000E-99

...(more lines)

-------------------------------------------------------------------------------------------------------

I'm curious since I'm trying to use other models as input for wd_ignite
(e.g. a cooled CO WD with final Teff = 5000K) and those models are not
created by MESA in this two part format.


Sumit

On Fri, Apr 3, 2015 at 8:26 PM, Robert FARMER <rjfarmer at asu.edu> wrote:

> Hi
> So if you look in the star/private folder you'll find a file called
> read_model.f90 that shows what mesa does with model files.
>
> If you look through the function read1_model, you'll see we load the v
> data into something to do with i_v. On a hunch if you grep inside the star
> folder
> grep -rin i_v star/p*/*
>
> You'll find this:
> public/star_data.inc:102:         integer :: i_v ! Lagrangian velocity at
> outer face of cell
>
> So it looks like its the Lagrangian velocity at outer face of a cell
>
> I think,
> Rob
>
>
>
>
>
> On Thu, Apr 2, 2015 at 9:14 AM, Bill Paxton <paxton at kitp.ucsb.edu> wrote:
>
>> Hi,
>>
>> Those are good questions, but they need to go to the mesa-users email
>> list where there are lots of folks willing and able to help you.
>>
>> If you haven't already done so, please sign-up for mesa-users and repost
>> your message there.
>>
>> Cheers,
>> Bill
>>
>>
>>
>> On Apr 1, 2015, at 5:31 PM, Sumit Sarbadhicary wrote:
>>
>> Hello Dr. Paxton,
>>
>> I am a graduate student and a new MESA user, working with version r7184
>> on a class project. I am a bit confused about the identity of a variable
>> 'v', defined in wd_ignite.mod contained in the MESA test_suite folder,
>> wd_ignite/. Is this the velocity of each Lagrangian mass shell? Or the
>> convective velocity? Or something else?
>>
>> I'm also a bit confused why there are *two* consecutive models in
>> 'wd_ignite.mod', instead of one (as it says in the very first line below).
>> I understand the purpose is to follow the evolution of an initial WD
>> profile through the carbon ignition phase, but the requirement of two
>> initial profiles eludes me.
>>
>> The first few lines of 'wd_ignite.mod' are as follows. 'v' appears as one
>> of the column headers (right of 'dq') :-
>>
>> ------------------------------------------------------------------------------------------------------------
>>
>> ! note: initial lines of file can contain comments
>>
>> !
>>
>>             12 -- two consecutive models for mesa/star, with velocity
>> variables (v). cgs units. lnd=ln(density), lnT=ln(temperature),
>> lnR=ln(radius), L=lumin\
>>
>> osity, dq=fraction of total mstar in cell; remaining cols are mass
>> fractions.
>>
>>
>>      year_month_day_when_created   20110717
>>
>>                           M/Msun      1.0000000009434213E+00
>>
>>                         star_age      9.4342161767044841E-01
>>
>>                        initial_z      2.0000000000000000E-02
>>
>>                         n_shells                            1273
>>
>>                         net_name   'co_burn.net'
>>
>>                          species                               9
>>
>>                       tau_factor      1.5000000000000000E+02
>>
>>
>>                 lnd                        lnT
>> lnR                          L                         dq
>>      *   v*   \
>>
>>                 h1                         he3
>> he4                        c12                        n14
>>       o16   \
>>
>>                      ne20                       mg24
>>   si28
>>
>>     1    -1.0783032764354498E+01     1.1640541494135320E+01
>> 2.0135894106074399E+01     4.8571750122419371E+32
>> 1.4075930508271764E-14     1.1361104977\
>>
>> 517480E-02     0.0000000000000000E+00     0.0000000000000000E+00
>> 0.0000000000000000E+00     2.5000000000000000E-01
>> 0.0000000000000000E+00     7.50000\
>>
>> 00000000000E-01     0.0000000000000000E+00     0.0000000000000000E+00
>> 0.0000000000000000E+00
>>
>>     2    -8.2810274680451421E+00     1.2428434759429717E+01
>> 2.0135268008479237E+01     4.8571750981602734E+32
>> 1.1250008967110072E-14    -4.6362609302\
>>
>> 080995E+02     0.0000000000000000E+00     0.0000000000000000E+00
>> 0.0000000000000000E+00     2.5000000000000000E-01
>> 0.0000000000000000E+00     7.50000\
>>
>> 00000000000E-01     0.0000000000000000E+00     0.0000000000000000E+00
>> 0.0000000000000000E+00
>>
>>     3    -7.8621441999255000E+00     1.2633241562558258E+01
>> 2.0135226974337208E+01     4.8522180474116445E+32
>> 1.1678717762613631E-14    -5.3005735647\
>>
>> 722219E+01     0.0000000000000000E+00     0.0000000000000000E+00
>> 0.0000000000000000E+00     2.5000000000000000E-01
>> 0.0000000000000000E+00     7.50000\
>>
>> 00000000000E-01     0.0000000000000000E+00     0.0000000000000000E+00
>> 0.0000000000000000E+00
>>
>>     4    -7.6227340198235742E+00     1.2703597628288353E+01
>> 2.0135198951381092E+01     4.8509545544014497E+32
>> 6.1656608267007517E-15    -3.9122855811\
>>
>> 425765E+01     0.0000000000000000E+00     0.0000000000000000E+00
>> 0.0000000000000000E+00     2.5000000000000000E-01
>> 0.0000000000000000E+00     7.50000\
>>
>> 00000000000E-01     0.0000000000000000E+00     0.0000000000000000E+00
>> 0.0000000000000000E+00
>>
>>
>> ------------------------------------------------------------------------------------------------------------------------------------------------------------
>>
>> Apologies if the questions are too trivial or ill-defined. Looking
>> forward to your reply!
>>
>> Yours sincerely,
>>
>> Sumit Sarbadhicary
>>
>> University of Pittsburgh
>>
>>
>>
>>
>> ------------------------------------------------------------------------------
>> Dive into the World of Parallel Programming The Go Parallel Website,
>> sponsored
>> by Intel and developed in partnership with Slashdot Media, is your hub
>> for all
>> things parallel software development, from weekly thought leadership
>> blogs to
>> news, videos, case studies, tutorials and more. Take a look and join the
>> conversation now. http://goparallel.sourceforge.net/
>> _______________________________________________
>> mesa-users mailing list
>> mesa-users at lists.sourceforge.net
>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150404/39e23293/attachment.html>


More information about the Mesa-users mailing list