Hallo zusammen,
der Befehl compute/barycorr hat mir heute bei der gleichen Beobachtung zwei leicht verschiedene Korrekturwerte für Total barycentric RV correction geliefert. Lasse ich Datum und Uhrzeit aus der bdf-Datei lesen, so ist der Wert um 0,4 km/s verschieden, als wenn ich Datum und Uhrzeit von Hand in den Befehl einsetze.

Anzeige bei Datum und Uhrzeit aus der Datei (Stationskoordinaten sind in rvcorr.prg enthalten):
Midas 005> (==>)(==>) rvcorr 18,13,45.8 -21,03,32
current dir. now: /home/jottschi/spektro/analyse
no CDELTi nor CDi_j matrix found... STEPi set to 1.0
FITS file: tmp00001.fts converted to: barycorr.bdf
Barycentric correction time: 5.817411185035E-03 day
Heliocentric correction time: 5.816522465156E-03 day

Total barycentric RV correction: -3.89611 km/s
Total heliocentric RV correction: -3.906156 km/s
(incl. diurnal RV correction of -0.0776 km/s)

Descriptor O_TIME of image barycorr.bdf used for date and UT.
current dir. now: /home/jottschi/spektro

Anzeige bei händischer Eingabe:
Midas 003> compute/barycorr 2011,07,03 23,16,53.1 18,13,45.8 -21,03,32 +07,59,33 +47,07,23
Barycentric correction time: 5.804280576029E-03 day
Heliocentric correction time: 5.803359576328E-03 day

Total barycentric RV correction: -4.324439 km/s
Total heliocentric RV correction: -4.33448 km/s
(incl. diurnal RV correction of -0.027816 km/s)

Ich habe mir dann den O_TIME descriptor angesehen. Bei der Wandlung von fts nach bdf wird aus 2011-07-03T23-16-53.096 der Wert 2011.5013698630

Gebe ich diesen Wert bei der händischen Eingabe ein, so erhalte ich wieder den Korrekturwert, als würde compute/barycorr den Wert selbst aus dem header der bdf-Datei lesen.

Läuft da bei der Wandlung was schief?

Ratlos,
Jörg!

Sali Otmar,

und vielen Dank. Tja, mein altes Leiden ist immer noch aktiv. Wieder einmal
nicht genau genug hingeschaut, um die Ursache selbst finden zu können.
Gut, bleibt aber noch die Tatsache, dass die Beobachtungszeit nicht
vollständig übernommen wird. Deshalb anhängend der Fits-Header als
Textausgabe.

Schöne Grüsse,
Jörg!

SIMPLE = T/CCDSOFT-SOFTWARE BISQUE 3
BITPIX = 16
NAXIS = 2
NAXIS1 = 1530
NAXIS2 = 1020
BSCALE = +1.000000000000e+000
BZERO = +3.276800000000e+004
BIAS = 100
FOCALLEN= +0.000000000000e+000
APTAREA = +0.000000000000e+000
APTDIA = +0.000000000000e+000
DATE-OBS= '2011-07-03T23:16:53.096'
TIME-OBS= '23:16:53.096 '
SWCREATE= 'CCDSoft Version 5.00.186'
SET-TEMP= -1.000000000000e+001
COLORCCD= 0
DISPCOLR= 1
IMAGETYP= 'Light Frame '
CCDSFPT = 1
XORGSUBF= 0
YORGSUBF= 0
CCDSUBFL= 0
CCDSUBFT= 0
XBINNING= 1
CCDXBIN = 1
YBINNING= 1
CCDYBIN = 1
EXPSTATE= 293
CCD-TEMP= -9.787817267812e+000
TEMPERAT= -9.787817267812e+000
INSTRUME= 'SBIG ST-8 3 CCD Camera'
EGAIN = +2.470000000000e+000
E-GAIN = +2.470000000000e+000
XPIXSZ = +9.000000000000e+000
YPIXSZ = +9.000000000000e+000
SBIGIMG = 6
USER_2 = 'SBIG ST-8 3 CCD Camera'
DATAMAX = 65535
SBSTDVER= 'SBFITSEXT Version 1.0'
EXPTIME = +6.000000000000e+002
EXPOSURE= +6.000000000000e+002
CBLACK = 47
CWHITE = 698
HISTORY Darkframe subtracted.
END


________________________________

Wer hat's gemacht? .... Die Schweizer ...

Hallo Jörg,

Midas hätte gerne einen Header TM-START mit der Startzeit in Sekunden.
Ich hänge mal ein kleines Perl-Skript an, dass diesen Header einfügt,
modifiziert aus meinem Baukasten

wie kann ich denn dann das Ergebnis aus compute/bary auf das Spektrum anwenden? Ist es ok, wenn man dies nach dem rebinning macht rebin/line macht?

!
! rebin file to heliocentric velocity, keep number of pixels
!
set/format i5 g15.8
define/local cvel/d/1/1 2.99792458d5
!
define/local fact/d/1/1 0.0
define/local n/i/1/1 0

defi/para p1 ? ? "input file"
defi/para p2 ? ? "output file"
defi/para p3 ? ? "helvel"
!
!
fact = 1.0/(1.0+{p3}/cvel)           ! factor (1+v/c)
!
write/out rebininng file {p1} with velocity {p3} to new file {p2}
! 
rebi/ii {p2} {p1} {p1} lin 0,{fact} spg             ! rebin