[mesa-users] .osc output format

RICHARD H D TOWNSEND townsend at astro.wisc.edu
Fri Feb 3 10:51:18 EST 2017


Hi Warrick —

Thanks for sharing your wisdom!

One think I’ll add is that as of more-recent releases of MESA, the pulsation output code now lives in the mesa/star/private/pulse.f90 file and other files (see especially pulse_osc.f90) files. It’s also a little tidier than before!

cheers,

Rich

> On Feb 3, 2017, at 3:05 AM, Warrick Ball <wball at bison.ph.bham.ac.uk> wrote:
> 
> 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.
>> ------------------------------------------------------------------------------
>> 
>> 
>> 
>> 



More information about the Mesa-users mailing list