Tags:
create new tag
, view all tags

GSC4SIM

depuis une position (historiquement le catalogue positionnel GSC) ou un nom : génère des commandes pour l'entrée dans Simbad

voir les réunions : Réunions Equipe Biblio

Options d'entrée sur la ligne de commande:

  • -source= : liste de catalogues séparés par des virgules (pour limiter à ces catalogues)
  • -c= : coordonnées (toujours en dernier) les espaces sont autorisés
  • -c.rs = : rayon en seconde (par défaut 1.5)
  • -v : ajoute quelques affichages pour mieux comprendre le déroulement
  • -v2 : affiche encore plus de détails (asu)
  • -beta : utilise la version en cours de test
  • -h : affiche l'aide
  • -a : équivalent à -source=.......tous les catalogues de gsc4sim..... (utile uniquement pour la recherche par nom) sinon, c'est déjà le comportement par défaut
  • -list : montre les catalogues par défaut "gaia2,tycho2,ucac2,ucac3,ucac4,2mass,wise_ok,sdss7,sdss9,usnob,denis,gsc1.2,gsc2.2,gsc2.3,ppmx,ukidss,tdsc"
BEGIN :
  • Prépare les commandes "asu" pour chaque catalogue
EXECUTION :
  • si la ligne commence par "-list" : affiche les sources par défaut

  • si la ligne commence par: (recherche par nom)

    • général pour tous les catalogues
      1. cat_init
      2. appelle exec_asu avec :
        • -oc=S (affichage des coordonnées en sexagésimal)
        • -source=...
        • -out.add=RAJ,DEJ
        • -out.all(sauf Gaia, où les colonnes sont précisées)
        • mode=1 pour SDSS>3
      3. définit le nom principal ($gsc) selon le catalogue et les chiffres dans le nom donné (après le 1er espace)
      4. si une réponse est trouvée, et option '-a' (search_all), alors lance la recherche par position (from_pos avec la position trouvée) sauf ??? pourquoi erreur ?? 2QZ, NVSS, SDSS>3 (faute de frappe:defalut_sources), ULAS, UGCS, UDXS, UKIDSS(déjà from_pos par défaut), Gaia2 (oubli?)
    • USNOB /[Uu][SNOsno-]*[Bb][1.0]* *[0-9]/
      1. définit le bibcode à "2003AJ....125..984M" et la qualité de magnitude à "E"
      2. exec_asu USNO-B1.0 avec USNO-B1.0=chiffres
    • USNOA2 /^[Uu][SNOAsnoa-]*2[.0]* *[0-9]/ (idem USNOB1)
      1. définit le bibcode à "1998yCat.1252....0M" et la qualité de magnitude à "E"
      2. exec_asu USNO-A2.0 avec USNO-A2.0=chiffres
    • 2QZ /^2[Qq][Zz] *[Jj][0-2][0-9]*/
      1. définit le bibcode à "2004MNRAS.349.1397C" et la qualité de magnitude à "E"
      2. exec_asu VII/241/2qz avec 2QZ=J+chiffres
    • UCAC2 /^[Uu][Cc][Aa][Cc]2 *[0-9]/
      1. définit le bibcode à "2004AJ....127.3043Z" et la qualité de magnitude à "D" (non utilisé ???)
      2. exec_asu I/289 avec 2UCAC=chiffres
    • UCAC3 /^[Uu][Cc][Aa][Cc]3 *[0-9]/
      1. définit le bibcode à "2009yCat.1315....0Z" et la qualité de magnitude à "D"
      2. exec_asu I/315 avec 3UC=chiffres
    • UCAC4 /^[Uu][Cc][Aa][Cc]4 *[0-9]/
      1. définit le bibcode à "2012yCat.1322....0Z" et la qualité de magnitude à "D"
      2. exec_asu I/322 avec UCAC4=chiffres
    • CMC14 /^[Cc][Mm][Cc]14 *J[0-9]/
      1. qualité de magnitude à "C
      2. exec_asu 1304 CMC14=chiffres
    • PPMX /^[Pp][Pp][Mm][Xx] *J[0-9]/
      • exec_asu 1322 PPMX=chiffres
    • MSX6C /^[Mm][Ss][Xx][6C]* *G[0-3]/
      • exec_asu 5114 MSX6C==chiffres
    • NSV /^[Nn][Ss][Vv] *[0-9]/
      • exec_asu B/gcvs/nsv_cat MSV===chiffres
    • DENIS /^DENIS *J[0-2][0-9.]*[+-][0-8][0-9]*/
      • exec_asu B/denis DENIS=chiffres
    • NVSS /^[Nn][Vv][Ss][Ss] *J[012]/ sub(/^[NVSnvs ]*/, ""); Jname = $0;
      • exec_asu VIII/65/nvss NVSS=chiffres
    • 2MASS /^2[Mm][asAS]* *J[0-9][0-9][0-9][0-9][0-9][0-9][0-9]/
      • exec_asu 2246 2MASS=chiffres
    • 2MASSI /^2[Mm][asAS-]*[Ii] *J[0-9][0-9][0-9][0-9][0-9][0-9][0-9]/
      • exec_asu 2241 Name=chiffres
    • 2MASX /^2[Mm][Aa][Ss][Ss]*[Xx] *J[0-9][0-9][0-9][0-9][0-9][0-9][0-9]/
      • exec_asu 7233 2MASX===chiffres
    • SDSS[3-9] JHHMMSS /^ [Ss][Dd][Ss][Ss][3456789] [Jj][0-2][0-9]/ (et les dr 10/11/12...?)
      • extrait le numéro de la DR utilisée dans le nom, et la position
      • définit la qualité de magnitude à "C"
      • définit le bibcode selon le numéro de la DR:
        • 3 : "2005AJ....129.1755A"
        • 4 : "2006ApJS..162...38A"
        • 5 : "2007yCat.2276....0A"
        • 6 : "2008ApJS..175..297A"
        • 7 : "2009yCat.2294....0A"
        • 8 : "2011yCat.2306....0A"
        • 9 : "2012ApJS..203...21A"
      • exec_asu :
        • rayon : -c.rs=0.2
        • -sort=mode,_r
        • mode=1
        • catalogue source :
          • sdds3 : 2259
          • sdss4: 2267
          • sdss5 : 2276
          • sdss6 : 2282
          • sdss7 : 2294
          • sdss8 : 2306
          • sdss9 : 5139
        • -c=position
      • si rien trouvé : recherche par position dans SDSS7
        • si toujours rien, "#===Looking for possible nearby SDSS stars" : recherche par position dans le catalogue de la DR donnée en entrée + rayon à 1sec (-c.rs=1)
      • (si une réponse est trouvée, et option '-a' (search_all) et (??estimated_vmag??), faute de frappe:defalut_sources)
    • SDSS NNNNN /^ [Ss][Dd][Ss][Ss][1-9] [0-9][0-9][0-9][0-9][0-9][0-9][0-9][0-9][0-9][0-9]*/
      • exec_asu source et bibcode selon la DR comme SDSS précédent et objID=chiffres
    • SDSS NNNN.NNN.NN /^ [Ss][Dd][Ss][Ss][1-9] [0-9][0-9][0-9]*[.][0-9]*[.][0-9]*[.][0-9]*/
      • exec_asu source et bibcode selon la DR comme SDSS précédent et SDSS-ID=chiffres
    • ULAS/UGCS/UDXS /^ [Uu][LlGgDd][AaCcXx][Ss][0-9]* *[Jj][0-2][0-9]*/ (valide que pour les DR 6/8/9 ???)
      • exec_asu :
        • ULAS8 : 2310/las8 et ULAS=
        • UGCS8: 2310/gcs8 et UGCS=
        • UDXS8: 2310/dxs8 et UDXS=
        • ULAS9: 2319/las9 et ULAS=
        • UGCS9: 2319/gcs9 et UGCS=
        • UDXS9: 2319/dxs9 et UDXS=
      • si rien trouvé recherche par position dans UKIDSS : from_pos
    • UGPS /^ [Uu][Gg][Pp][Ss][0-9] [Jj][0-2][0-9]*/
      • extrait le nombre après "[Uu][LlGgDd][AaCcXxPp][Ss][0-9]" (???erreur de copié/collé : (ne fonctionne jamais par nom ???)
      • exec_asu UGPS6 : 2316/gps6 et UGPS= (valide que pour la DR6 explicite ???)
      • si rien trouvé recherche par position dans UKIDSS : from_pos
    • UKIDSS /^ [Uu][Kk][Ii][Dd][Ss] *[Jj][0-2][0-9]*/
      • recherche par position dans UKIDSS : from_pos
    • GSC 2.3 /^[Gg][Ss][Cc]2[.]3 *[NS]/ et finit par /[NS][0-9A-Z][0-9A-Z][0-9A-Z][0-9][0-9][0-9][0-9][0-9][0-9]/
      • exec_asu de GSC2.3 avec GSC2.3=NS+chiffres
      • (définit le bibcode + qmag, mais écrasé plus tard)
    • GSC2.2 /^[Gg][Ss][Cc]2[.2]* [NS]/ et finit par /[NS][0-9]*/
      • exec_asu de GSC2.2 avec GSC2.2=NS+chiffres
      • (définit le bibcode + qmag, mais écrasé plus tard)
    • WISE JHHMM.. /^[Ww][Ii][Ss][Ee] *J[012]/
      • exec_asu 2311 et WISE= (avec l'option "ccf=0000" ???pour wise_ok)
      • (définit le bibcode + qmag, mais écrasé plus tard)
    • Gaia DR2 /^[Gg][Aa][Ii][Aa] [Dd][Rr]2 *[0-9][0-9]*/
      • enlève ",_r" de la commande asu (pas nécessaire ???)
    • GSC /^[GSC ]*[0-9]/
      • compte le nombre de chiffres entre les '-':
        • si <2 ou >3 -> erreur "****Not a GSC number:"
        • si ==2 : ajoute une dernier "-1"
      • définit le nom gsc = "GSC " +2 premiers chiffres
      • définit le nom tyc = "TYC " +3 groupes de chiffres
      • définit le code GSC=avec les 2 premiers chiffres alignés à 5 chiffres en complétant avec des 0 # as written in GSC1.2
      • appelle exec_asu Tycho-2 avec TYC1=.. TYC2=... TYC3=...
        • si plusieurs résultats :
          • message ***Blend:"
          • gsc=tyc
        • si une position est trouvée: recherche par position sur les catalogues par défaut sauf tycho2
      • appelle exec_asu gsc1.2 et GSC=code GSC
      • définit bibcode et qmag mais écrasé plus tard
      • si une position trouvée, lance la recherche autour de la position sur les catalogues par défaut sauf gsc1.2 + ajoute denis (si la position a une valeur négative)

  • recherche par position : (-c=....) pour chacun des catalogues définis (dans l'option -source= ou par défaut)

    • Affiche "!# Look around"
    • Appelle exec_asu en récupérant les options définies pour la recherche par nom et enlève le dernier critère (supposé xxxACROxxx=)
Fonction cat_init : appelée au début à chaque recherche par nom
  1. initialise le tableau cat_done à 0
  2. initialise le tableau cat_qual à vide
  3. initialise DR à vide
  4. initialise writo à 0 (pour savoir si on met "o " en début du paragraphe)
Fonction exec_asu (execute asu) : appelle de Vizier sur un catalogue précis + critère de colonne/position
  • tt : fichier contenant les arguments de commande
  • position du résultat
  • flag si on utilise la position(utilisé uniquement pour SDSS)
  1. lance "asu -b" avec les options données
  2. extrait les données: 
        phase=0: debut
        phase=1: (nouvelle)table ("^----Table")
                 lecture de l'entete: nom de la table, catid, nom et indice des colonnes
        phase=2: lecture des valeurs ("^---*\t---")
                 un bloc par type de catalogue en se basant sur catid ou sur un nom de colonne (ex UCAC4, Vtmag (pour tycho)...)
                    -> utilisation de from_SDSS, from_2MASS ou from_UKIDSS 
                Note: toute les lignes de chaque table est lue - dans le cas d'une recherche positionnelle, seule la 1ere ligne est prise en compte 
                      pour la generation de la commande Simbad (voir fonction from_pos)
  1. Cas particuliers:
    1. Tycho : (vérifie l'existance du nom de colonne "VTmag")
      1. init_Tycho2
      2. switch(catid){
        case 1259: catref = "2000A&A...355L..27H" qmag = "D"
        case 1239: catref = "2000A&A...355L..27H" qmag = "D"
        case 1276: catref = "2002A&A...384..180F" qmag = "D"
        }
      3. recalcule l'erreur de position : e_RAdeg = e_RAdeg + 8.75*e_pmRA; // idem DE
      4. écrit "o TYC" n1-n2-n3
      5. Si c'est le 1er trouvé alors prépare aussi le nom "GSC" n1-n2
      6. Si l'erreur de position:
        • existe : qualité = "B"
        • est vide : qualité = "C"
      7. stocke les positions pour les renvoyer
      8. error_ellipse avec e_RAdeg et e_DEdeg
      9. affiche la commande des positions (en fonction de la qualité) : prcom avec longueur d'onde "(O)" et l'erreur de position calculée avec err2d(e_DEdeg, e_RAdeg)
      10. Si l'erreur de position existe :
        • lance error_ellipse avec les erreurs du mouvement propre
        • affiche les commentaire de l'ellipse d'erreur
        • la commande de mouvement propre avec qualité "B" (sans passer par prcom !!???)
      11. Si VTmag n'est pas vide ET si BTmag n'est pas vide:
        • mag = VTmag + interpol(BTmag-VTmag, aTc, aVT, nTmag)
        • affiche la commande de la magnitude "V" : amag avec mag et e_VTmag
        • mag = BTmag + interpol(BTmag-VTmag, aTc, aBT, nTmag);
        • affiche la commande de la magnitude "B" : amag avec mag et e_BTmag
      13. affiche l'identificateur supplémentaire GSC s'il y en a un
      14. note que le catalogue tycho2 a été fait (cat_done["tycho2"] = 1)
    2. 2MASS/2MASSI : catid==2246 ou catid==2241
      1. appelle from_2MASS
        1. définit le bibcode, selon le catalogue :
          • 2246 : "2003yCat.2246....0C"
          • sinon : "2000yCat.2241....0S"
        2. définit le nom principal ($gsc) avec "2MASS"+J+num de la colonne "2MASS"
        3. écrit "o "+nom 2MASS
        4. si le nom principal ($gsc) avait été définit avant et différent de celui-ci alors affiche "!Name changed from '"
        5. si la colonne "prox" existe, alors l'affiche en commentaire : "!+++Source has prox=" prox2M "arcsec";
        6. affiche la commande d'ajout des positions avec prcom avec:
          • qualité="C",
          • (NIR),
          • err2d avec l'ellipse d'erreur avec error_ellipse et les colonnes errMaj,errMin,errPA ou MajAxis,MajAxis,Theta (2 fois MajAxis ????) *1000 et flag -1
        7. affiche les commandes d'ajout de magnitudes Jmag,Hmag,Kmag avec amag2mass et la qualité venant de la colonne "Qflg" (1/2/3e caractère) (### Warning quality "U" ==> undetected !)
        8. récupère les magnitudes Bmag et Rmag et si elles sont bonnes (l'une ou l'autre =0 ou vides), affiche la commande d'ajout de magnitude avec amag, et la qualité "!E"
      2. note que le catalogue 2MASS a été fait (cat_done["2mass"]=1)
    3. 2MASX : catid == 7233
      1. catref="2006AJ....131.1163S"
      2. définit le nom principal ($gsc) avec le nom "2MASX J" + num de la colonne "2MASX"
      3. écrit "o "+ nom principal
      4. stocke les positions pour les renvoyer
      5. qualité ="C"
      6. affiche les commandes de positions (en fonction de la qualité) : prcom avec longueur d'onde "(NIR)" sans l'erreur
      7. si rS3sig et Sb/a et Spa ne sont pas vides : affiche la commande des dimensions, avec :
        • majaxis = rS3sig / 30.0 (multiply by 2 (diam), divide by 60 (arcmin))
        • minaxis = majaxis * Sb/a
        • si Spa<0 : pos_angle = 180-Spa
        • (NIR)
        • qualité=C
      8. affiche la commande de magnitude des colonnes J.ext, H.ext, K.ext ( amagt)
      9. note que le catalogue 2masx a été fait (cat_done["2masx"] = 1)
    4. SDSS : catalogue parmi : 2259,2267,2276,2282,2294,2306,5139,2294
      1. appelle from_SDSS :
        1. récupère le numéro dans la colonne SDSS sinon SDSS<N>
        2. définit le nom principal ($gsc) avec le nom SDSS
        3. si mode=2 alors sort avec found2=1
        4. si on trouve q_mode alors l'affiche
        5. écrit "o "+ nom SDSS trouvé si 1er
        6. affiche les mag en commentaire (!)
        7. définit le bibcode selon :
          • ref_sdss[2259] = "2005AJ....129.1755A";
          • ref_sdss[2267] = "2006ApJS..162...38A";
          • ref_sdss[2276] = "2007yCat.2276....0A";
          • ref_sdss[2282] = "2008ApJS..175..297A";
          • ref_sdss[2294] = "2009yCat.2294....0A";
          • ref_sdss[2306] = "2011yCat.2306....0A";
          • ref_sdss[5139] = "2012ApJS..203...21A";
        8. affiche la commande de position avec prcom et :
          • si e_RAJ2000 existe (SDSS-3 has no error...) : calcule l'ellipse avec erro_ellipse et err2d et valeurs*1000 sinon sans erreur
          • (O)
          • qualité ="C"
        9. Si on trouve q_mode: (# SDSS>=7 has a flag for 'clean photometry')
          • et ="+" alors qmag="B", (# in case of clean photometry)
          • sinon qmag="C"
        10. Affiche en commentaire (!) la commande d'ajout de type d'objet avec une qualité "E" (juste pour éviter le "c o", et forcer le "a o" ???) avec, si la colonne cl :
          • =3 : Galaxie,
          • =6 : et z>0.2 : QSO?,
          • sinon Star
        11. Affiche les commandes d'ajout de magnitudes umag, gmag, rmag imag, zmag avec amagt et l'option"AB".
        12. # Convert magnitudes from Fukugita 1996AJ....111.1748F
          # g' = V + 0.56(B-V) - 0.12
          # r' = V - 0.49(B-V) + 0.11
          #==> V = 0.47g' + 0.53r'
          #==> B = 1.42g' - 0.42r' -0.22
        13. et affiche B et V en commentaire
        14. S'il y a un redshift : colonne zsp présente (# SDSS-6 has redshift), sinon z=-9, affiche la commande (sans prcom ???) d'ajout de redshift avec:
          • z:s
          • (O)
          • Les critères pour déterminer la qualité du redshift ne sont
            applicables que si z >= 0.01. Pour les redshifts plus faibles, et en
            particulier pour les étoiles, la qualité du redshift sdss doit être
            automatiquement "E". Sinon, "C"
          • et en commentaire (!) la commande des "data velocities"
        15. Si le mouvement propre existe, ajoute la commande d'ajout de pm en commentaire (!) avec:
          • un angle d'ellipse d'erreur à 90,
          • qualité="D"
        16. Affiche les autres magnitudes en commentaire
      2. note que le catalogue sdss7 a été fait (cat_done["sdss7"] = 1)
      3. précise si au moins un mode=2 a été trouvé (si on a trouvé found2=1)
    5. Gaia2 : catid=1345
      1. définit le bibcode ="2018yCat.1345....0G"
      2. récupère le numéro dans la colonne Source
      3. définit le nom principal = "Gaia DR2 "+num
      4. écrit la commande "o "+nom Gaia
      5. affiche la commande d'ajout des coordonnées avec prcom avec :
        • si le mouvement propre :
          • existe : qualité ="A",
          • est vide : qualité ="C",
        • err2d sur les erreur position ICRS,
        • (O),
        • position_angle de l'ellipse d'erreur = 90
      6. Si la paralaxe>0 :
        1. affiche la commande de correction de paralaxe (sans prcom ???),
        2. affiche la commande d'ajout des data plx,
        3. affiche la commande d'ajout des data distance avec:
          • 1000.0/Plx
          • ( 1000.0*e_Plx ) / (e_Plx ^2)
      7. Si le mouvement propre existe, affiche la commande d'ajout via prcom avec err2d et qualité="A"
      8. Si la magnitude G existe, affiche la commande de correction de magnitude G avec qualité ="C" (sans passer par prcom ???)
      9. Si la colonne RV existe, affiche la commande d'ajout de de vitesse avec :
        • selon l'erreur:
          • <0.9 qualité = "A",
          • <2 qualité ="B",
          • <5 qualité ="C",
          • <10 quaité ="D",
          • sinon qualité="E"
        • (O)
      10. Affiche la commande d'ajout de data velocities avec en plus o_RV
    6. UCAC2 : si 2UCAC présent dans les colonnes
      1. récupère le numéro dans la colonne 2UCAC
      2. définit le nom principal =UCAC2 + num
      3. écrit "o " +nom UCAC2
      4. définit le bibcode="2003yCat.1289....0Z"
      5. calcule l'ellipse d'erreur avec error_ellipse et un angle à 90 et flag à 0
      6. affiche la commande d'ajout de positions avec prcom avec :
        • si pmRA=0 et pmDE=0:
          • qualité="C"
        • si pm ok :
          • qualité ="B"
        • err2d avec sqrt(e_RAdeg^2 + (EpRA-2000)^2 * e_pmRA^2) et idem pour DE
        • (O)
      7. affiche le commentaire pmok="!" (=> ne devrait pas utiliser prcom si en commentaire ????)
      8. calcule l'ellipse d'erreur du mouvement propre avec error_ellipse, et un agle à 90 et le flag à 1
      9. affiche la commande d'ajout du mouvement propre avec :
        • si pmRA=0 et pmDE=0:
          • qualité="D"
        • si pm ok :
          • qualité ="C"
          • err2d
      10. note ucac2 comme fait
    7. UCAC4 : Si UCAC4 présent dans les colonnes
      1. récupère le numéro dans la colonne UCAC4
      2. définit le nom principal = UCAC4 + num
      3. écrit "o "+nom UCAC4
      4. définit le bibcode="2012yCat.1322....0Z"
      5. met en commentaire : pmok="!", dès que :
        • pmRA ou pmDE ou e_pmRA ou e_pmDE sont vides
        • la somme du carré des erreurs <0.05
        • |pmRA|<0.0001 et |pmDE|<0.0001
      6. affiche "!# pm=" et "!# dt="
      7. calcule l'ellipse d'erreur avec error_ellipse et un angle à 90 et flag à 0
      8. affiche la commande d'ajout de positions avec prcom avec :
        • si pmok="!"
          • qualité="C"
        • si pm ok :
          • qualité ="B"
        • err2d avec sqrt(e_RAdeg^2 + (EpRA-2000)^2 * e_pmRA^2) et idem pour DE
        • (O)
      9. affiche le commentaire pmok="!" (=> ne devrait pas utiliser prcom si en commentaire ????)
      10. calcule l'ellipse d'erreur du mouvement propre avec error_ellipse, et un agle à 90 et le flag à 1
      11. affiche la commande d'ajout du mouvement propre avec :
        • si pmok="!"
          • qualité="D"
        • si pm ok :
          • qualité ="C"
          • err2d
      12. affiche la commande d'ajout de magnitude R avec amag avec:
        • valeur=f.mag
        • erreur=e_a.mag
        • qualité=E ou qualité=D si erreur<0.05
      13. ajoute d'autres magnitudes B,V,g,r,i avec amag avec :
        • valeur=colonne name+"mag"
        • erreur=e_XXX
          • si l'erreur n'est pas en décimales :
            • et =99 : affiche en commentaire "!"
            • erreur = erreur/100 + 2 chiffres après la virgule
        • qualité = D, sauf si f_XXX contient '-' alors qualité=E
      14. note ucac4 comme fait
    8. UCAC3
    9. DENIS
    10. GSC2.2 : si "GSC2.2" dans les colonnes
      1. récupère le numéro dans la colonne GSC2.2
      2. définit le nom principal = GSC2 + num
      3. écrit "o "+nom GSC2
      4. définit le bibcode="2001yCat.1271....0S" et qualité de magnitude="E"
      5. affiche la commande d'ajout de positions avec prcom avec :
        • qualité = "D"
        • bande (O)
        • ellipse d'erreur avec error_ellipse et 1000*e_RAdeg 1000*e_DEdeg, angle à 90
      6. affiche la commande d'ajout de type d'objet en commentaire: "!" devant (avec prcom, avec une qualité E)
        • en fonction de la colonne "Class" 0 =>Star ; 1 => Galaxy
      7. affiche les magnitude avec amagt sur les colonnes : Bjmag, Rmag, Vmag, Imag
      8. note gsc2.2 comme fait
    11. GSC2.3 : si "GSC2.3" dans les colonnes
      1. récupère le numéro dans la colonne GSC2.3
      2. définit le nom principal = GSC2.3 + num
      3. écrit "o "+nom GSC2.3
      4. définit le bibcode="2008AJ....136..735L" et qualité de magnitude="E"
      5. affiche la commande d'ajout de positions avec prcom avec :
        • qualité = "D"
        • bande (O)
        • ellipse d'erreur avec error_ellipse et 1000*e_RAdeg 1000*e_DEdeg, angle à 90 et flag à -2 #at least 1mas error
      6. affiche la commande d'ajout de type d'objet en commentaire: "!" devant (avec prcom, avec une qualité E)
      7. affiche les magnitudes avec la commande amagt des colonnes : Umag, Bmag, Vmag
      8. affiche les magnitudes avec la commande amagt2 des colones "jmag" -> B, "Nmag" -> I, "Fmag" -> R
      9. note gsc2.3 comme fait
    12. USNO-A2 : si "USNO-A2.0" dans les colonnes
      1. récupère le numéro dans la colonne USNO-A2.0
      2. définit le nom principal = USNO-A2.0 + num
      3. écrit "o "+nom USNOA
      4. définit le bibcode="1998yCat.1252....0M" et qualité de magnitude="E"
      5. affiche la commande d'ajout de positions avec prcom avec :
        • qualité "D"
        • bande (O)
        • sans erreur
      6. affiche les commandes d'ajout de magnitudes Bmag et Rmag avec amag (si >0), sans erreur, et avec 1 précision .N
      7. note usnoa comme fait
    13. USNO-B1 : si "USNO-B1.0" dans les colonnes
      1. récupère le numéro dans la colonne USNO-B1.0
      2. définit le nom principal = USNO-B1.0 + num
      3. écrit "o "+nom USNOB
      4. définit le bibcode="2003AJ....125..984M" et qualité de magnitude="E"
      5. si muPr<=5
        1. affiche le mouvement propre ignoré pmRA,pmDE
        2. affiche la commande d'ajout de positions avec prcom avec :
          • qualité "D"
          • bande (O)
          • les positions calculées en 2000 : RAJ2000+( (Epoch-2000)*pmRA/3.6e6 )/cos(DEJ2000*P1/180) et DEJ2000+( (Epoch-2000)*pmDE/3.6e6 )
          • error_ellipse avec e_RA, e_DE (sauf si =0 : transformé en 300) angle à 90, flag 0 #1mas
      6. si muPr>5 (#True error)
        1. affiche la commande d'ajout de positions avec prcom avec :
          • qualité "C"
          • bande (O)
          • error_ellipse avec sqrt(e_RAJ2000^2 + (Epoch-2000)*(Epoch-2000)*e_pmRA*e_pmRA) (idem DE) avec e_XX transformé en 300 si =0 et e_pmXX transformé en 1 si <1
          • error_ellipse sur
        2. affiche la commande d'ajout de mouvement propre pmRA,pmDE avec une qualité =D et avec e_pmXX transformé en 1 si <1
      7. affiche les commandes d'ajout de magnitudes Bmag et Rmag avec amag, sans erreur, et avec 1 précision .N
        • si Bmag=0 : Bmag=B1mag (idem pour R)
      8. note usnob comme fait
    14. ...
Fonctions from_2MASS, from_SDSS ... voir directement dans exec_asu.

Fonction prcom (print command) : affiche la commande d'ajout ou correction (commenté ou non) en fonction de la qualité

  • type de donnée
  • qualité
  • erreur
  1. S'il s'agit de coordonnées ou mouvement propres : prépare la commande "C" (change), sinon "A" (add)
  2. Si la qualité commence par un "!', alors affichage en commentaire avec "[bad qual]"
  3. Si aucune donnée de ce type a été entré avant, on garde en mémoire la qualité pour ce catalogue pour cette donnée.
  4. Si on avait déjà une donnée :
    • et que sa lettre de qualité est meilleure alors la commande est changé en "C", sinon commenté "![qual<]"
    • et à qualité égale, regarde l'erreur, si elle est meilleur : la commande est changée en "C", sinon commenté "![!err>]"
Fonction prmag (print magnitude) : affiche la commande d'ajout de magnitude pour SIMBAD (générique pour toutes les magnitudes de tous les catalogues, sauf cas exceptionnel généré en amont dans les catalogues)
  • command (a/c) # not used ?
  • lettre de la couleur
  • valeur
  • erreur
  1. Garde la lettre de qualité ($qmag) si définie avant, sauf si à A/B/C, et qu'on a :
    • une erreur<0.1 => qualité="C"
    • une erreur<0.5 => qualité="D"
    • sinon qualité="E"
  2. si la valeur et l'erreur sont vides => return
  3. Si il y a aucune erreur ou erreur=0 ou valeur=0 => [bad qual], qualité = "!E"
  4. Prépare la commande simbad "m " + "AB" ou non définie avant +" couleur"
  5. Affiche la commande (prcom) avec :
    • si l'erreur est vide
      • et si la valeur a plus de 4 charactères, affiche la valeur tronquée à 2 décimales
      • sinon, affiche la valeur telque
    • si l'erreur <= 4 charactères
      • et au moins 3 décimales tronque la valeur à 2 décimales
      • en dessous de 3 décimales, affiche la valeur tel que
    • sinon (l'erreur > 5 charactères), affiche la valeur et l'erreur tronquées à 3 décimales
Fonction amag (add magnitude) : appelle la commande d'ajout de magnitude pour SIMBAD (appelle la méthode générique, peu utile)
  • lettre de la couleur
  • valeur
  • erreur
  1. appelle prmag avec le paramètre "commande" = "a m" + variable AB si définie avant (SDSS+CMC14)
Fonction amagt2 (add magnitude) : prépare la commande d'ajout de magnitude depuis un nom de colonne (utilisée uniquement dans GSC2.3)
  • nom de colonne de la table Vizier
  • lettre de la couleur
  1. récupère la valeur depuis la table en cours, et prend la valeur correspondant à la colonne donnée
  2. Si la valeur n'existe pas, ou est vide : ressort
  3. Si cette valeur existe récupère l'erreur dans la colonne "e_"+nom_de_colonne,
    • si la colonne n'existe pas, mets une erreur vide
  4. appelle amag avec ces données
Fonction amagt (add magnitude) : prépare la commande d'ajout de magnitude depuis un nom de colonne
  • nom de colonne de la table Vizier
  1. récupère la valeur depuis la table en cours, et prend la valeur correspondant à la colonne donnée
  2. Si la valeur n'existe pas, ou est vide : ressort
  3. Si cette valeur existe récupère l'erreur dans la colonne "e_"+nom_de_colonne,
    • si la colonne n'existe pas, mets une erreur vide
  4. appelle amag avec ces données et la bande=1er caratère du nom
Fonction amag2mass (add magnitude for 2mass only) : prépare le commande d'ajout de magnitude avec un changement de qualité spécial pour 2MASS

  • nom de colonne de la table Vizier
  • lettre de qualité venant du catalogue
  1. si la lettre de qualité du catalogue est :
    • A ou B : qualité="C"
    • C ou D : qualité = "D"
    • E ou F : qualité = "E"
    • sinon : affiche "!+++magnitude has quality , ignored" et sort
  2. appelle amagt avec le nom de colonne
Fonction error_ellipse :

Fonction err2d :

Fonction from_pos (recherche dans asu selon une position) :

gsc4sim

NOTE: commandes simbad: a: addition, c: correction

gsc4sim-c.awk: filtre pour la sortie Simbad avec des listes de priorites

  1. la position: Tycho-2,UCAC2,2MASS,2Mass,SDSS,GSC2.2,GSC1.2,DENIS
  2. la magnitude optic: SDSS,Tycho-2,UCAC2,GSC2.2,GSC1.2
  3. la magnitude IR: 2MASS,2Mass,SDSS,DENIS,GSC2.2,GSC1.2
  4. le type d'objet: SDSS,GSC2.2,USNO-B1.0
Le script commente les ligne (voir en dessous) puis ajoute une unique sortie qui est definie par la liste des priorites

Note: le script n'est utilise qu'avec l'option -c= (et pas -c ...!!) => plus utilisé depuis la version du 7/11/2018

transformation sortie:
  • "o (.*)" => "a i $1"
  • suppression "v .*"
  • suppression ligne vide
  • c c (.*) => utilise la priorite pour memoriser valeur position $1 (1)
  • a p (.*) => utilise la priorite pour memoriser valeur position $1 (1)
  • c o (.*) => utilise la priorite pour memoriser valeur IR magnitude $1 (3) !!! TRES bizare!!!! ca devrait etre (4) non?
  • a m[AB]* [UBVRugr] => utilise la priorite pour memoriser valeur optical magnitude (2)
  • a m[AB]* [IJHKiz] => utilise la priorite pour memoriser valeur IR magnitude (3)

gsc4sim.awk: script principal

  • ReadMe: gsc4sim notes (Gilles)
-- AnaisOBERTO - 2018-12-17
Attachments Attachments
Topic attachments
I Attachment Action Size Date Who Comment
Unknown file formatEXT ReadMe manage 8.6 K 2019-01-14 - 15:55 GillesLandais gsc4sim notes (Gilles)
Unknown file formatawk gsc4sim.beta.awk manage 86.7 K 2020-09-15 - 08:44 AnaisOBERTO Code AWK de gsc4sim
Edit | Attach | Print version | History: r20 < r19 < r18 < r17 < r16 | Backlinks | Raw View | Raw edit | More topic actions
Topic revision: r20 - 2021-04-06 - AnaisOBERTO
 
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2024 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback