[mesa-users] .osc output format

Mon Feb 6 04:11:47 EST 2017

Warrick,

Thanks a lot for your so-well documented mail of last Friday. It gives me a good starting point for my investigations. I must tell you that in the meantime I managed to obtain .osc output format files from MESA. Now my current problem is how can I make MESA to produce smaller files, since one single run of MESA produce more than 1.5 GB of file size.

Thanks a lot, Warrick!
Jesús
________________________________________
De: Warrick Ball [wball at bison.ph.bham.ac.uk]
Enviado: viernes, 03 de febrero de 2017 10:05
Para: URIBE MARCHANTE Jesús
CC: RICHARD H D TOWNSEND; mesa-users at lists.sourceforge.net
Asunto: Re: [mesa-users] .osc output format

Jesús,

I have some experience in writing/re-writing data between different file
formats so I though I'd offer a few tips.

First, several formats (including .OSC) are described in this document:

https://www.astro.up.pt/corot/ntools/docs/CoRoT_ESTA_Files.pdf

which is actually mentioned in the default files around where the
pulsation format is set.  That will at least tell you which variables are
going into the .OSC files (although many aren't actually used in the
oscillation calculation).

I also thought it's worth mentioning that writing the oscillation data is
slightly delicate. First, the oscillation codes want some values that
aren't obviously available in the profile data (e.g. the second
derivatives of pressure and density at the centre).  Second, MESA has some
cell-averaged variables which have to be interpolated to get the values on
the cell boundaries, because that's what the oscillation codes expect.
(Certainly the density, not sure about any others.)

I've personally been messing around with the oscillation model files and
found MESA's subroutines very helpful, so this serves as another example
of how to find out what MESA does.

For example, try searching for the control 'pulse_info_format' in the
code.

$ grep pulse_info_format $MESA_DIR/star/p*/*.f90
/home/wball/mesa/7624/star/private/ctrls_io.f90:         pulse_info_format, &
/home/wball/mesa/7624/star/private/ctrls_io.f90:         s% pulse_info_format = pulse_info_format
/home/wball/mesa/7624/star/private/ctrls_io.f90:         pulse_info_format = s% pulse_info_format
/home/wball/mesa/7624/star/private/profile.f90:               fname_out = trim(fname) // '.' // trim(s% pulse_info_format)
/home/wball/mesa/7624/star/private/profile.f90:               call save_pulsation_info(s% id, s% pulse_info_format, &
/home/wball/mesa/7624/star/public/star_lib.f90:            id, pulse_info_format, &
/home/wball/mesa/7624/star/public/star_lib.f90:         character (len=*), intent(in) :: pulse_info_format
/home/wball/mesa/7624/star/public/star_lib.f90:         which_format = get_pulsation_format(pulse_info_format, ierr)

Looks like the call to save_pulsation_info star/private/profile.f90 is a
good place to start.  Let's find that subroutine directly.

$ grep subroutine\ save_pulsation_info $MESA_DIR/star/p*/*.f90
/home/wball/mesa/7624/star/private/pulsation_info.f90:      subroutine save_pulsation_info( &
/home/wball/mesa/7624/star/private/pulsation_info.f90:      end subroutine save_pulsation_info

The (well-named) file star/private/pulsation_info.f90 sounds right. I've
just looked through it myself and its various function calls will show you
how MESA writes the pulsation info.  You have to follow a few function
calls back and forth but ultimately nearly everything is there.

Hope that helps!

Cheers,
Warrick

------------
Warrick Ball
Postdoc, School of Physics and Astronomy
University of Birmingham, Edgbaston, Birmingham B15 2TT
wball at bison.ph.bham.ac.uk
+44 (0)121 414 4552

On Wed, 1 Feb 2017, URIBE MARCHANTE Jesús wrote:

> Thank you, Rich!
> One additional question, do you know if it is possible to migrate the .data files into .osc format?
>
> Jesús
> ________________________________________
> De: RICHARD H D TOWNSEND [townsend at astro.wisc.edu]
> Enviado: miércoles, 01 de febrero de 2017 21:51
> Para: URIBE MARCHANTE Jesús
> CC: mesa-users at lists.sourceforge.net
> Asunto: Re: [mesa-users] .osc output format
>
> Hi Jesús —
>
> You need to set the following flags to the &controls namelist:
>
> pulse_info_format = ‘OSC’
> write_pulse_info_with_profile = .TRUE.
>
> cheers,
>
> Rich
>
>> On Feb 1, 2017, at 2:08 PM, URIBE MARCHANTE Jesús <Jesus.URIBE at euipo.europa.eu> wrote:
>>
>> Hello!
>>
>> I'm using MESA to produce typical .data output. Does anyone know how to get .osc format output? transform .osc into .data files?
>>
>> thank you!
>>
>> Jesús Uribe
>> ------------------------------------------------------------------------------
>> AMENDING EU TRADE MARK REGULATION: The amending Regulation on the EU trade mark entered into force on March 23 2016.
>>
>> As a result, the Office is now called the European Union Intellectual Property Office (EUIPO).
>>
>> Please update your contacts and bookmarks. Our web address changed to www.euipo.europa.eu on March 23 2016.
>>
>> Our email addresses changed to name.surname at euipo.europa.eu on March 23 2016, although the old name.surname at oami.europa.eu email addresses will run in parallel to the new addresses for several months.
>>
>> Learn about changes to fees, technical changes and institutional changes here: https://euipo.europa.eu/ohimportal/en/eu-trade-mark-regulation
>> ------------------------------------------------------------------------------
>> IMPORTANT: This message is intended exclusively for information purposes. It cannot be considered as an official EUIPO communication concerning procedures laid down in the European Union Trade Mark Regulations and Community Designs Regulations. It is therefore not legally binding on the EUIPO for the purpose of those procedures.
>> The information contained in this message and attachments is intended solely for the attention and use of the named addressee and may be confidential. If you are not the intended recipient, you are reminded that the information remains the property of the sender. You must not use, disclose, distribute, copy, print or rely on this e-mail. If you have received this message in error, please contact the sender immediately and irrevocably delete or destroy this message and any copies.
>> ------------------------------------------------------------------------------
>>
>>
>> ------------------------------------------------------------------------------
>> Check out the vibrant tech community on one of the world's most
>> engaging tech sites, SlashDot.org! http://sdm.link/slashdot
>> _______________________________________________
>> mesa-users mailing list
>> mesa-users at lists.sourceforge.net
>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>
> ------------------------------------------------------------------------------
> AMENDING EU TRADE MARK REGULATION: The amending Regulation on the EU trade mark entered into force on March 23 2016.
>
> As a result, the Office is now called the European Union Intellectual Property Office (EUIPO).
>
> Please update your contacts and bookmarks. Our web address changed to www.euipo.europa.eu on March 23 2016.
>
> Our email addresses changed to name.surname at euipo.europa.eu on March 23 2016, although the old name.surname at oami.europa.eu email addresses will run in parallel to the new addresses for several months.
>
> Learn about changes to fees, technical changes and institutional changes here: https://euipo.europa.eu/ohimportal/en/eu-trade-mark-regulation
> ------------------------------------------------------------------------------
> IMPORTANT: This message is intended exclusively for information purposes. It cannot be considered as an official EUIPO communication concerning procedures laid down in the European Union Trade Mark Regulations and Community Designs Regulations. It is therefore not legally binding on the EUIPO for the purpose of those procedures.
> The information contained in this message and attachments is intended solely for the attention and use of the named addressee and may be confidential. If you are not the intended recipient, you are reminded that the information remains the property of the sender. You must not use, disclose, distribute, copy, print or rely on this e-mail. If you have received this message in error, please contact the sender immediately and irrevocably delete or destroy this message and any copies.
> ------------------------------------------------------------------------------
>
>
>
>
------------------------------------------------------------------------------
AMENDING EU TRADE MARK REGULATION: The amending Regulation on the EU trade mark entered into force on March 23 2016.

As a result, the Office is now called the European Union Intellectual Property Office (EUIPO).

Please update your contacts and bookmarks. Our web address changed to www.euipo.europa.eu on March 23 2016.

Our email addresses changed to name.surname at euipo.europa.eu on March 23 2016, although the old name.surname at oami.europa.eu email addresses will run in parallel to the new addresses for several months.

Learn about changes to fees, technical changes and institutional changes here: https://euipo.europa.eu/ohimportal/en/eu-trade-mark-regulation 
------------------------------------------------------------------------------
IMPORTANT: This message is intended exclusively for information purposes. It cannot be considered as an official EUIPO communication concerning procedures laid down in the European Union Trade Mark Regulations and Community Designs Regulations. It is therefore not legally binding on the EUIPO for the purpose of those procedures.
The information contained in this message and attachments is intended solely for the attention and use of the named addressee and may be confidential. If you are not the intended recipient, you are reminded that the information remains the property of the sender. You must not use, disclose, distribute, copy, print or rely on this e-mail. If you have received this message in error, please contact the sender immediately and irrevocably delete or destroy this message and any copies.
------------------------------------------------------------------------------