Dear colleagues,

Here you can find the code to calculate corrections 'Dm' 
according to Eq. (6) in my paper: "On the effect of emission 
lines on the UBVR photometry", 2007, New Astronomy, Vol 12, 597-604.

Available also at: http://arxiv.org/abs/0704.3856
                  (here as "neaspa164.pdf")
--------------------------------------------------------------

The files are: 
     - the code: "dml",
     - examples of the input files: "data.1" and "data.2". 

"dml" - the code calculates Dm corrections DU, DB, DV, DRc, DRj 
        to the true continuum due to emission lines. The input 
        file includes line fluxes. The magnitude of the true 
        continuum (i.e. emission line free continuum) is then
                     m(cont.) = m(obs.) - Dm.

"dml.f" - a fortran source file of the code.

"data.1" - example of emission line spectrum in the form: 
           [lambda, line flux], where the 'line flux' is in units 
           of the local continuum. Example concerns to RR Tel, 
           the HST/STIS (18/10/2000) spectrum. 

"data.2" - the same as for the "data.1" file, but in the form: 
           [lambda, I/Ic, sigma], where 'I/Ic' is the height of 
           the line in units of the local continuum and 'sigma'
           is the half-width of its Gaussian fit. Thus the line 
           flux is expressed as its approximation by the Gauss curve. 
           This form of the input data may be convenient if one 
           reconstructs the emission line spectrum from a figure. 

Output is displayed as: 

  Type of the input file:
         1 = lambda,F
         2 = lambda,I,sigma
  Your case: 1
 Input file [lambda,F]: data.1
    DU =    -1.67              ! in mag
    DB =    -1.74              ! in mag
    DV =    -0.96              ! in mag
    DRc=    -1.04              ! in mag
    DRj=    -0.89              ! in mag
-------------------------------------------------------------
 
   With best wishes, 

              Augustin Skopal                 18th June 2007