NAME:
transf
PURPOSE: (one line)
Determine transformation coefficients from instrumental to standard mags.
DESCRIPTION:
This will take a set of all-sky standard star measurements and derive
the best fitting set of transformation coefficients from instrumental
to standard magnitudes.
Normally, you give this routine ALL the data taken with a single filter.
Only those objects that are standards in this filter and the alternate
color will be used in the fit. Only those that should be standards
can end up marked as bad from this program.
My conventions for the photometric fit are begin with Hardie's
nomenclature. However, there are three additional terms I've added:
the first order expansion term of extinction with respect to time
and the full second-order expansion in airmass and color.
The basic formula for the transformation is:
m0 = m - kX - n(t-t0)X - k"CX + eC + fX^2 + gC^2 + Z
where
m = instrumental magnitude
k = extinction coefficient, mag/airmass
X = airmass
n = coefficient of the 1st order expansion of extinction as a
function of time
t = Time of observation (in hours)
t0 = Reference time for n, time dependent correction is zero at
this time, usually is the middle of the observation set.
k" = second order extinction coefficient
C = Standard system color of the object
e = color term
f = second-order term in airmass squared
g = second-order term in color squared
Z = zero point
m0 = Standard magnitude
The fitting process takes advantage of the uncertainties on the instrumental
magnitude and the uncertainties on the forced terms. The uncertainties
on the catalog standard magnitudes are ignored. The coefficients, k",
e, and n are considered to be auxillary coefficients and can be either
fit for or forced. k", e, f, and g can be forced to explicit values.
n is either fit for or set to 0. The extinction and zero-point are always
fit for.
CATEGORY:
Photometry
CALLING SEQUENCE:
transf,names,jd,filter,color1,color2,airm,mag,err,bad,
[ RESID=resid,FILE=landfile,NOPLOT=noplot,OTHER=other,
OLAB=olab,TITLE=title ]
INPUTS:
names - String array of standard names.
jd - Julian date of observations.
filter - Filter number of data. (0-U, 1-B, 2-V, 3-R, 4-I)
color1 - First filter for color.
color2 - Second filter for color. (eg., B-V is 1,2)
airm - Vector of airmasses for observations.
mag - Vector of instrumental magnitudes.
err - Vector of uncertainties.
bad - Vector of flags, 0=good, 1=bad, default is good
OPTIONAL INPUT PARAMETERS:
KEYWORD INPUT PARAMETERS:
These keywords affect how the fit is done. All except KTIME must be
two-elements vectors, the first element is the value of the coefficient
and the second element is the uncertainty on the coefficient. The default
is to fit for these terms.
CTERM - Optional, fixed color term. If provided, the color term is
assumed, not fitted.
K2 - Second order extinction coefficient. If provided, this term
is assumed, not fitted.
COLSQ - Term proportional to the standard color squared
AMSQ - Term proportional to the airmass squared
KTIME - Flag, enable time dependent extinction.
REFCOLOR - Reference color to use during the fit. The old default value for
this was implicitly zero. Now it is the mean of all the good
values used.
REFAM - Reference airmass to use during the fit. The old default value
for this was implicitly zero. Now it is the middle of the range
of airmass for all good values used.
-------------------
These keywords affect the plotting output
NOPLOT - Flag, if true inhibits the summary plots.
OLAB - Label for the plot of residuals vs. OTHER.
OTHER - If provided, allows plotting the residuals against some other
variable, eg., exposure time. Must match the length of
the data input vectors.
-------------------
Other control options
FILE - Optional name for the Landolt "system" Standards catalog.
NOEDIT - Flag, if set inhibits final interactive editing of fitted points.
This keyword has no effect and is not necessary if the current
plotting device is 'PS'.
RESID - Flag, if true enables printing a full residual table.
TITLE - Title for plots.
FILES - Flag, if true requests saving information to summary files.
This feature is deprecated in favor of saving to the database.
NOSAVE - Flag, if true supresses saving information to databases or
summary files. NOSAVE implies NOFILES.
CHISCALE -Flag, if true rescales final solution to force reduced Chi=1.0
TAGDATE -String identifying date of data.
TAGINST -String identifying instrument.
HISTFILE -File name of history file to update extintion values in.
Only one line for TAGDATE/TAGINST is allowed. If a duplicate
is found, the appropriate line is replaced.
MAGRESID - Edit bad values by mag residual, not sigma residuals
NODB - flag, if set will suppress saving data to the database.
DATABASE - Name of MYSQL database for transformation fit. Passed
to puttran. The default is 'phot'.
DVERBOSE - flag passed to gettran and puttran to set db screen verbosity.
TABLENAME- Name of table in MYSQL database for transformation search.
Passed to puttran. The default is 'transf'.
OUTPUTS:
No explicit vector outputs. All output is graphical, printed to the
screen, or saved to external files.
KEYWORD OUTPUT PARAMETERS:
COMMON BLOCKS:
None
SIDE EFFECTS:
RESTRICTIONS:
PROCEDURE:
MODifICATION HISTORY:
93/08/05, Initial version completed, Marc W. Buie, Lowell Observatory.
93/10/05, MWB, added support for 1993 catalog.
95/02/22, MWB, modified plots.
96/10/16, MWB, added summary save to file, NOFILES, CHISCALE, TAG*, and
HISTFILE keywords actions added.
96/11/22, MWB, changed to new master photometry catalog
97/01/24, MWB, added MAGRESID flag
97/2/7, MWB, massive rewrite
2004/02/09, MWB, changed path to transf files.
2006/08/02, MWB, added DB, DATABASE, TABLENAME keywords
2006/10/16, PLC, added NOSAVE keyword
2006/12/28, PLC, added DVERBOSE flag.
2014/03/06, MWB, added REFCOLOR and REFAM keywords, added new transformation
terms. Saving to files is deprecated, it still happens
but the new terms are not included in the save. If you
want all of the information you must use the database.
The NOFILES keyword has been changed to FILES.
The DB keyword has been changed to NODB.