spectro Variables de type beta Cephei: sigma Sco

[Olivier Thizy]

Bonjour,

 

 

Après BW Vul, j'ai voulu faire une autre étoile pulsante. J'ai choisi sigma Sco et je l'ai suivi pendant 4h (soit les deux tiers de la période de pulsation). Comme BW Vul, j'ai fait des spectrogrammes 2D dont voici un exemple autour de H-alpha:

 

 

 

Ces spectrogrammes sont réalisés en deux temps:

 

1/ je sélectionne une partie du spectre échelle, je mets le continuum à 1 et je converti en format texte (DAT). J'utilise pour cela les fonctions de SpcAudace dans AudeLA. Voic le script utilisé pour ceux que cela intéresse:

 

#
# Script pour les fichiers FITS de sigma Sco dans le répertoire courant en DAT
# SpcAudace doit être lancé...
#
# v1.10 (Olivier Thizy): version "batch" simple mais à éditer à chaque fois...
#

global audace

::console::affiche_resultat "Traitement du répertoire: "
::console::affiche_resultat $audace(rep_images)
::console::affiche_resultat "\n"

set liste_fichier [glob -dir $audace(rep_images) -tails *.fit]

set nbre_fichier [llength $liste_fichier]

if {$nbre_fichier > 0} {
 for {set i 0} {$i<$nbre_fichier} {incr i} {
   set nom_fichier [file rootname [string tolower [lindex $liste_fichier $i]]]

### define continuum value to set all spectra to continuum=1 in this area
# Sodium doublet:
#    set continuum_value [spc_icontinuum $nom_fichier 5840 5860]
# Halpha:
#    set continuum_value [spc_icontinuum $nom_fichier 6690 6710]
# Hbeta:
   set continuum_value [spc_icontinuum $nom_fichier 4970 4990]

### Crop le spectre dans la zone qui nous intéresse
# Sodium doublet:
#    spc_select $nom_fichier 5600 6000
# Halpha:
#    spc_select $nom_fichier 6450 6700
# Hbeta:
   spc_select $nom_fichier 4775 5090

# mets le continuum à 1
   loadima "${nom_fichier}_sel"
   mult [expr 1/$continuum_value]
   saveima "${nom_fichier}_sel"

# Conversion du spectre FITS en DAT
   spc_fits2dat "${nom_fichier}_sel"

#    ::console::affiche_resultat $nom_fichier
#    ::console::affiche_resultat "\n"
   }
 } else  { ::console::affiche_resultat "pas de fichier / no file to process.\n"}

 

 

2/ j'utilise IRIS pour créer le spectrogramme:

 

>data_anim IRIS_filelist output 5800 5950 0.05
>data2image r 20000 70
>scale 4 1 5

 

Le fichier IRIS_filelist.lst contient la liste des fichiers DAT avec la phase. Le fichier output contient une liste de fichiers (R1, R2, R3...) avec la phase correspondante.

 

 

Plus de spectrogrammes de cette étoile sont sur: http://www.shelyak.com/en/sigsco.html

 

 

Cordialement,

Olivier Thizy

Vous ne verrez plus des étoiles comme avant !

http://www.shelyak.com/fr/

[gabal]

Pas de limite à l'observation des atmosphère stellaire maintenant !

Elle pulse aussi de façon plus régulière, sans chocs et sans dédoublement de raies (?).

Pourquoi mets tu le continuum à 1 ?

[Olivier Thizy]

Pas de limite à l'observation des atmosphère stellaire maintenant !

Elle pulse aussi de façon plus régulière, sans chocs et sans dédoublement de raies (?).

Pourquoi mets tu le continuum à 1 ?

 

Bonjour,

 

en fait, je fais en sorte que tous les spectres soit au mêem niveau pour ^tre comparable dans le spectrogramme 2D. J'ai pris 1 comme cela.

 

Cordialement,

Olivier Thizy

Vous ne verrez plus des étoiles comme avant !

http://www.shelyak.com/fr/