For the final assignment of this series you will assemble all of your EB code to form eb_reduce.pro. This procedure will accept as input a string variable specifying the directory holding the raw data. For instance, for the first night of observations, the input will be “2003June18”. The procedure will then call the necessary subroutines (e.g. flats.pro, etc) needed to reduce each observation of OO Aql. Copy the working versions of all your subroutines into your hw10/ directory.
eb_reduce will then write to disk a postscript file showing the light curve of OO Aql for the specified night. Make a plot of FA/FB + FA/FC as a function of JD – min(JD), where FA/FB is the flux measured from OO Aql with respect to star B. Use psopen to create the postscript figures (use doc_library, ‘psopen’ for the program documentation).
Repeat the reduction process for three additional nights: 2002June19, 2002July02, and 2002July03.
The following material will be due in class Wednesday, Dec 7:
1) The 4 postscript plots
2) A master light curve that combines the data from all 4 nights.
3) A printout of each routine called by eb_reduce. Use ctrl-c ctrl-h to insert a standard comments header at the top of each routine. Use escape ctrl-q to ensure your code has the proper, standard indentation.
4) Attach a cover sheet to each subroutine printout with your name and the subroutine’s name. Center this information on the page, both vertically and horizontally. Please staple the pages for each subroutine separately and then paperclip the entire stack together.
Here is an example plot (note, this is NOT the light curve of OO Aql!)
Here is an example of a standard program comments header and standard indentation.
HINTS:
1) Use textoidl to generate your axis labels. For example, textoidl(‘F_{A}’) produces FA.
2) When creating the plot, make sure to make nice, thick axes using the xthick and ythick keywords.
3) Want to learn more about the plot procedure? Check out the…wait for it…online help file!
