SKIRT / SKIRT9

SKIRT version 9 -- advanced radiative transfer in dusty systems
http://www.skirt.ugent.be
GNU Affero General Public License v3.0
37 stars 31 forks source link

Add observer info to output of distant instruments #214

Closed petercamps closed 5 months ago

petercamps commented 5 months ago

Description With this update, parallel-projected instruments such as Frame and SED instruments add information on the observer sight line and model-observer distance to the output files, as a convenience to the user.

Specifically, FITS file headers now include additional keywords in the file header as in the following example:

    CROTA1  =      4.500000000E+01 / Inclination angle, in deg                      
    CROTA2  =      4.500000000E+01 / Azimuth angle, in deg                          
    CROTA3  =      3.000000000E+01 / Roll angle, in deg                             
    REDSHIFT=      2.000000000E+00 / Redshift (if zero, distances are equal)        
    DISTLUMI=      1.607026089E+04 / Luminosity distance                            
    DISTANGD=      1.785584544E+03 / Angular diameter distance                      
    DISTUNIT= 'Mpc     '           / Units of distances                             

Similarly, SED column text files now include a header line as in the following example:

    # SED at inclination 45 deg, azimuth 45 deg, redshift 2, luminosity distance 16070.26089 Mpc

Motivation This information may facilitate automated processing of the SKIRT output without the need for referring to the ski file. The feature was suggested by user @djsavic.

Compatibility The extra header line in SED column text files will break code that skips a predefined number of header lines rather than relying on the # prefix to skip all comment lines.