Joined: 10 Feb 2006 Posts: 135 Location: Flagstaff
Posted: Wed Dec 31, 2008 9:39 pm Post subject: Adding in sidereal time to image headers
I'm just back from a Magellan run, and discovered to my chagrin that my headers all have a bogus value for the sidereal time. Everything else is fine, including the airmass, UT time, date, etc. I actually discovered this after running setjd and doslit, and somewhere along the line some IRAF task recomputed the airmasses based on my bogus sidereal times, with the result that my airmasses are now wacky as well.
Given that I know the longitude, time, and data, it should be possible to put valid sidereal time values into my headers, and then run setairmass. I've gotten as far as determining that either astedit or astcalc might be able to do this, but I confess that after multiple readings of the help pages, I'm still scratching my head.
Any additional guidance would be greatly appreciated.
cheers,
phil
ASTHEDIT is probably the task you need, specifically the 'mst' function will help: Create a command file called 'cmd' that contains the following
Code:
st = mst (@'date-obs', ut, obsdb ("lco", "longitude"))
Then running as
Code:
cl> asthedit *.fits cmd
should update the 'ST' keyword (assuming the default parameters and that the obsdb value for Las Campanas is good enough for you, if not you can enter a more precise longitude). Let me know if you still have questions.
You cannot post new topics in this forum You cannot reply to topics in this forum You cannot edit your posts in this forum You cannot delete your posts in this forum You cannot vote in polls in this forum