Réunion VizieR du 04/02/2014

Présents : Marianne B., Sylvain G., Gilles L., François O., Emmanuelle P., Patricia V.

Questions Patricia :

Problèmes techniques liés à la nouvelle machine :

- xpdf ne fonctionne plus : renvoit le message "segmentation fault (core dumped)"
=> le test avec evince (programme plus récent que xpdf pour lire les PDF) ne fonctionne pas non plus : "connection rejected (wrong authentification)"

- gv n'est pas installé (pour lire postscripts (?) et PDF).

- show_images pour afficher une liste d'images. En fait, n'est pas utilisable tel quel, il s'agit d'un programme local pour le cat. J/A+A/501/539

- dans les groupes d'appartenance pour cats, par rapport à l'ancienne version il manque "man". => Gilles dit que tant que ça ne pose pas de problème, il n'est pas nécessaire de le rajouter. Pour info : "nogroup" est le groupe qui permet l'accès au FTP, "www-data" = accès au submit.

Problème de Fabienne avec un viz4sim pour la ref. J/MNRAS/419/3505 alors que la col. Simbad est bien là et fonctionne...

Pour rappel, le programme viz4sim (suivi d'une référence) permet de sortir les commandes de MAJ SIMBAD pour ajouter une réf. à la liste des objets de la col. Simbad (ou SimbadName).
Ici, la col. est retrouvée mais ne récupère pas les noms d'objets... En testant avec une autre ref. dont le viz4sim a été effectué hier (le même jour que pour Fabienne) (J/AJ/144/42), le viz4sim fonctionne bien...

Rep. François : En fait le problème du viz4sim est un problème un peu général, un flag (VO_OSELECT) qui n'est pas mis proprement dans vizin
pour les colonnes qui sont utilisées dans \vizSimbad (ou vizNED).
J'ai corrigé manuellement, mais je regarde de plus près dans vizin...

Cat. Jan/29 : des tables de courbes de lumières

Les tables 2, 3 et 4 contiennent, respectivement, 16, 9 et 34 bandes différentes. Cela fait beaucoup pour un seul menu... Comment faire le graphe ?
=> En fait, la table 1 avec le graphe en lumière blanche suffira...

Cat. Jan/27 (Gerner et al. The chemical evolution in the early phases of massive star formation. I.) : spectres FITS avec des Fréquences à transformer en vitesses et des valeurs à récupérer dans les headers de chaque fichier...

Le .graph :

set v0 = `fits2a -TOC fits/$file | fgrep "VELO_LSR" | acut -c10-20`
set nu0 = `fits2a -TOC fits/$file | fgrep "RESTFREQ" | acut -c10-30`

fits2a fits/$file | gawk -v v0=$v0  -v nu0=$nu0 '\
  BEGIN {c=299792458; f=c/nu0/1000; v0=37 }\
  { v = v0-$1*f;  printf "%11.3f %10s\n", v, $2 }}' \
| tee /tmp/Jan.27      \
| graph -T $Vgraph -C $plotarg -X 'v(LSR) [km/s]' -Y 'Tmb [K]' -L "$Name" \
  --pen-colors 1=red4:2=red $argv:q

devrait fonctionner...
set valeur = cmde permet de récupérer les valeurs à partir des fichiers FITS.

Cat. SDSS quasars DR10 en approche...

Il y a 4 codes d'erreurs pour le redshift. Il n'y en avait qu'un dans les version précédentes, est-ce normal ?
=> Laisser ces 4 codes et voir à quoi ils correspondent si possible...

Questions Emmanuelle :

Ref. 2012ApJ.752.53 : table 1 liste d'id de KOI qui n'ont pas d'excès IR (donc pas utilisés dans l'analyse)... A prendre ??

=> Pas pour VizieR. Voir avec Cécile si on veut ajouter la ref. sur les objets dans Simbad et récupérer les positions le cas échéant... (Marianne rappelle qu'on ne va pas systématiquement ajouter les KOI dans Simbad lorsque les coo ne sont pas données...)

Ref. J/ApJ/745/173 : problème dont la solution a été oubliée... 2 tables avec un fichier .sim par table ; les 2 tables sont mergées dans VizieR. La 2è liste de nom est plus large que la première donc certains noms sont tronqués dans VizieR... Comment faire ?

=>Utiliser le vizSet pour indiquer la largeur de colonne à la Col. vizSimbadName...
\vizSet{ table[12] }{ SimbadName }{ type=a22 } (par exemple).

Ref. J/AJ/144/62 : table 2, une absolute J magnitude entre 4 et 15, normal ?

=> Cela dépend des types d'objets. Pour les étoiles, la magnitude absolue est la magnitude apparente de l'étoile à 10pc.
Une magnitude absolue de 15 correspond à des objets de luminosité très faible comme les brown dwarfs par exemple.
Le soleil a une magnitude apparente de -27 mais une magnitude absolue de 5, tandis que Sirius a une magnitude absolue de 0 pour une magnitude apparente de -1.5. Les supergéantes, par contre, ont une magnitude absolue autour de -5, etc.

Complément de François : Pour les magnitudes absolues: la magnitude absolue ne dépend pas de la distance (c'est la luminosité intrinsèque de l'étoile); on peut rajouter la formule :

m - M = 5 log(r/10) + absorption
si m = magnitude apparente,
M = magnitude absolue
r = distance en parsecs de la source

c'est-à-dire que la magnitude apparente augmente de 5 magnitudes si la distance augment d'un facteur 10, par ex. Sirius serait à la magnitude 0 si elle était à 10pc, à la magnitude 10 si elle était à 1000pc (1kpc), à la magnitude 15 si elle était à 10kpc, etc.

Ref. J/ApJS/209/39 : FITS image cubes - pour exemple : quelles informations à donner ?

Le fits.list donne les principales informations (fits.list fits/*) :

277.33335+00.50000 10       305x 721x 133 fits/Serpens_12CO.fits
277.33335+00.50000 10       305x 721x 133 fits/Serpens_13CO.fits

Où on a les RAdeg+DEdeg, l'échelle en arcsec/pix, le Nx x Ny (voir description qui suit la table fournit par le fits.list).
Il manque dans la description le "133" qui correspond au Nz (Number of pixels along Z-axis)
En faisant le fits2a -toc fits/*, la dernière partie de chaque ligne nous fournit le CTYPE3 (ce Nz), ici "[VELOCITY=-10925/28675(300)]" que l'on retrouve avec l'option -TOC également. Nous n'avons pas les unités mais comme il s'agit de nuages moléculaires galactiques, elles doivent être en m/s !
Pour les images, le RESTFREQ serait un bon paramètre à ajouter, ici on ne l'a pas dans le header mais on peut l'ajouter au titre de chaque fichier.

Peut-être y aurait-il moyen d'améliorer le fits.list pour avoir directement le 3è axe et le RESTFREQ quand ils sont disponibles ?

Ref. J/ApJ/751/50 : Les auteurs (Gerke et al.) disent "the catalogs developed here replace those presented in G05" (Cat. J/ApJ/625/6) : rendre ce dernier cat. obsolète ?

=> Oui !
La procédure :
1. Rendre le vieux cat. (J/ApJ/625/6) obsolète dans le FTP :
Faire un cdcat J/ApJ/625/6 pour se retrouver dans le répertoire concerné.
Utiliser la commande modcat -old . et dérouler les différents fichiers édités :
- Keep it public? => y
- .status édité. Se rendre tout en bas pour remplir le \cRemoved avec le nom du Cat. qui le remplace ( \cRemoved{ 04-Feb-2014 }{ J/ApJ/751/50 } ).
- C'est bon? > y (ou "o" c'est pareil...)
- Dans le fichier ~/lib/Obsolete qui s'ouvre : tout en bas ajouter le commentaire : "# See J/ApJ/751/50" à la ligne de mise à jour qui s'est créée automatiquement (comme pour le fichier ~/lib/Public).
- Un fichier "=obsolete= " est créé dans le répertoire. Ajouter également le commantaire "See J/ApJ/751/50" (ce sera un commentaire publique).
- Dans le fichier ".history", on peut également ajouter un commentaire (celui-ci est pour nous).
- Les autres fichiers ne sont pas à modifier.

Pour vérifier, dans le FTP, on voit le fichier "=obsolete=" et dans le Browse, on voit le commentaire :

*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
J/ApJ/625/6   ==OBSOLETE version of Catalogue==
04-Feb-2014: (by Marianne Brouty) See J/ApJ/751/50
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*

2. Rendre le catalogue obsolète dans VizieR
Pour cela, il faut refaire un 2v depuis le répertoire de J/ApJ/625/6 mais uniquement lorsque le cat. J/ApJ/751/50 est public (pour l'instant, il n'est que dans le VizieR local donc on ne le fait pas).
François doit aussi vérifier si le 2v sur un catalogue rendu obsolète fonctionne correctement...

La semaine suivante, le catalogue se retrouvera dans la liste des catalogues obsolètes. Il faudra donc faire un findcat -old si on veut le retrouver...

Remarque François: En fait, il existe la commande make_obsolete, plus rapide !

Elle est rapidement décrite dans la section "Pipe-line d'introduction" http://cdsarc.u-strasbg.fr/doc/viz/index.htx#pipeline

Remarque Gilles :

De la documentation pour la création de graphes en langage pseudoSQL est disponible dans la doc VizieR "Génération d'un graphe VizieR à partir d'un fichier SQL"

L'URL : http://cdsarc.u-strasbg.fr/doc/viz/vizgraph.htx
Il y a notamment un beau schéma synthétique et des exemples.
Dans la section 4. Description des outils, on a la commande vizgraph qui permet de voir ce que récupère le graphe comme données :
Exemple Cat. III/185 :
vizgraph : tout seul => donne la liste des graphes SQL
vizgraph graphdiaghr : liste les valeurs du graphe générées par le graph_sql "graphdiaghr". (Sauf que ça ne marche pas chez moi).
Le fichier temporaire donnée par l'URL /tmp/... à la fin permet de voir la syntaxe du graph.
vizgraph graphdiaghr -d : est l'option "debug" qui permet de voir les erreurs...

Problème de synchronisation du B/wds sur les miroirs à cause d'un caractère bizarre... ? A voir.

-- EmmanuellePerret - 2014-02-04