DOI: 10.5281/zenodo.8051249
title: Three-component modelling of O-rich AGB star winds I. Effects of drift using forsterite – dataset

description: The data provided here include all parameter files, log files, and a set of the<br> binary output files that are the basis for the publication in A&A. The file 'file_listing.txt' contains a complete list of files and directories<br> in all gzipped tar files. Each individual gzipped tar file is formatted as<br> follows: Mm.m_Ll.ll_Ttttt.tar.gz where<br> m.m :: the assumed mass of the model, in solar masses<br> l.ll :: The assumed luminosity, in log10(solar luminosities)<br> tttt :: The effective temperature of the star, in Kelvin. <br> The contents of the tar files vary according to the model, but here is the<br> general directory structure: nodr/ :: non-drift / PC models<br> drift/ :: drift models nodr/init<br> drift/init :: Initial model files created using John Connor. <br> File suffixes are the following: .par :: Plain-text parameter file that contains all parameters that are<br> different from the respective default value in the model.<br> Consequently, to see what parameters were actually used, it is<br> necessary to look in the log file (see below). .bin :: Binary file that contains output of converged models. Each model is<br> stored in two versions, first the previous time step and then the<br> current time step (having access to the model code T-800, data of both<br> time steps are needed to restart model calculations at that time<br> step). The initial model file only contains one model; where the previous<br> time step data are the same as the current time step data. We provide a tool to read this file, see below. Note! These files can get pretty large and are therefore only<br> available for a smaller number of the models in the Zenodo dataset.<br> Please ask the corresponding author for the missing files is the<br> need should appear. .log :: Plain-text log file that shows the used model parameters and a number<br> of key properties for each converged model.<br> The encoding of this file is UTF-8. .inf :: Plain-text secondary log file that contains the header of the<br> [primary] log file as well as timing information.<br> The encoding of this file is UTF-8. .tpb :: Secondary binary file that contains a number of properties specified<br> at the outer boundary, typically for each consecutive time step. We provide a tool to read this file, see below. .lis :: Plain-text file with the iteration history. Unavailable here. .liv :: Plain-text file with values specified for a number of properties at<br> each gridpoint. Unavailable here. .inp :: Plain-text file that is used to launch a model; some are present.<br> This file is automatically generated by the tool that launches T-800<br> and is typically removed when T-800 launches. Unavailable here. .eps :: Encapsulated PostScript file created by John Connor when calculating<br> the initial model. <br> Model evolution structure - file endings before the suffix: _rlx :: Files related to relaxing the T-800 calculations on the initial model<br> created by John Connor. _exp :: Files related to expanding the initially compact model to using the<br> full radial domain. _fix :: Files related to the intermediate stage where calculations are changed<br> from expansion to outflow.<br> <br> _out :: Files related to the outflow stage of the calculations; this is what<br> you want to look at to see the wind evolution. Results in the paper<br> are calculated using these data. Note! Some outflow stage calculations continue the evolution of the previous<br> set of files. The underlying reason for continued calculations is typically<br> that the calculated time interval is too short. Such files are typically<br> given the extension '_cont.lin_out', '_cont2.lin_out', etc. <br> Stored data in the binary files: The binary files (suffix '.bin') contain the full radial structure in the<br> following 10 (PC models) or 11 (drift models) primary variables: mr: radius<br> mm: integrated [gas] mass<br> md: gas density<br> mu: gas velocity<br> me: internal energy<br> mj: radiative energy<br> mh: radiative flux<br> n0: dust moment, forsterite (Fo)<br> nm: number density of magnesium atoms<br> ns: number density of silicon atoms<br> v0: dust velocity, forsterite (only drift models) Other properties are derived from these primary variables using auxiliary code<br> that isn't part of this dataset. <br> Load files: Two tools are provided here that can load the binary data files using the<br> Interactive Data Language (IDL): sc_load_bin (for files with the suffix '.bin'): Loads the full content of a T-800 binary file and returns a structure<br> with the data. <br> sc_load_tpb (for files with the suffix '.tpb'): Loads the full content of a T-800 'tpb' binary file and returns a<br> structure with the data. Note! Due to the way models run on clusters, this file is sometimes<br> incomplete; this happens when the model code T-800 is stopped as the<br> cluster-specific walltime is reached. If this is the case, it is<br> necessary to use the binary file instead, where data are saved<br> typically every 20:th time step. Alternative tools for use with Python and Julia could be considered for<br> writing, but where not yet available when this dataset was made public.<br> Please contact the corresponding author for a current status on this issue.