----- Documentation for Diam\average_diameter.pro ----- NAME: AVERAGE_DIAMETER AUTHOR: pierre.cruzalebes@oca.eu PURPOSE: average diameter estimations using the weighted mean CATEGORY: angular diameter estimation INPUTS: need text input file (in OBS_PATH/INP/ directory) with: - input file number - input files name(s) REQUIREMENTS: data read in input file(s) (in OBS_PATH/OUT/ directory) must be in ASCII format and contain (in this order) : - model flag, flux flag, synthetic exitance+radiance file names, data format flag, starting parameters, free parameter flags, goodness-of-fit parameter, input parameters, inferior errors, superior errors, confidence level OUTPUTS: data written in output file (in OBS_PATH/OUT/ directory) must be in ASCII format and contain (in this order) : - model flag, flux flag, synthetic exitance+radiance file name, data format flag, starting parameters, free parameter flags, goodness-of-fit parameter, output parameters, inferior errors, superior errors, nber of degrees of freedom LOCAL PROCEDURE CALLED: RESET_PLOT LOCAL FUNCTION USED: EXTRACT_PROHEADER WSTDEV LOCAL SYSTEM VARIABLE USED: FORMAT_OUT NB_MAX_PARAM OBS_PATH PLOT_DEV SPIDAST_PATH REVISION HISTORY: Written by pcr 2008/04/09 last modification by pcr 2010/02/11 ----- Documentation for Diam\get_broadband_diameter.pro ----- NAME: GET_BROADBAND_DIAMETER AUTHOR: pierre.cruzalebes@oca.eu PURPOSE: measure broadband photometric diameter (IRFM) by scaling synthetic spectral radiant exitance with broadband-photometry flux measurement CATEGORY: angular diameter estimation INPUTS: need text input file (in OBS_PATH/INP/ directory) with : - input broadband photometric flux measurement file name - number of filter file(s) - input broadband filter transmission file name(s) (in same row) - synthetic spectral radiant exitance file name ('none' if not available) - if needed : effective temperature (for Planck and Engelke) REQUIREMENTS: data read in broadband-photometry flux measurement file (in subdirectory of OBS_PATH/) must be in ASCII format and contain (in this order) : - mean wavelength and equivalent width [mum], flux and absolute error on flux [W/m^2/mum] data read in filter file (in ALL_DATA_PATH/ directory) must be in ASCII format and contain (in this order) : - wavelength [mum] and transmission data read in input synthetic exitance file (in SYNTHE_PATH/ directory) must be in ASCII format and contain : - wavelength [AA = 10^-4 mum] and exitance [erg/cm^2/s/AA = 10 W/m^2/mum] OUTPUTS: data written in output file (of name input_file_stem.out in OBS_PATH/OUT/ directory) must be in ASCII format and contain (in this order, in same row) : - model flag , flux flag, initial synthetic exitance+radiance file names, data format flag, starting parameters, free flags, goodness-of-fit parameter, output parameters, inferior errors, superior errors, nber of degrees of freedom data written in covariance file (of name input_file_stem.cova in OBS_PATH/OUT/ directory) must be in ASCII format and contain the covariance matrix (of size NB_MAX_PARAM*NB_MAX_PARAM) LOCAL PROCEDURE CALLED: ENGELKE PLANCK RESET_PLOT LOCAL FUNCTION USED: EXTRACT_PROHEADER TRAP TRI_STDEV LOCAL SYSTEM VARIABLE USED: ALL_DATA_PATH FORMAT_OUT MAS2RADIAN MAX_REL_ERR NB_LAMBDA_MAX NB_MAX_PARAM OBS_PATH SPIDAST_PATH SYNTHE_PATH PLOT_DEV REVISION HISTORY: Written by pcr 2006/11/27 last modification by pcr 2010/06/27 ----- Documentation for Diam\get_color_diameter.pro ----- NAME: GET_COLOR_DIAMETER AUTHOR: pierre.cruzalebes@oca.eu PURPOSE: calculate angular diameters from interstellar corrected Johnson color indexes (surface brightness method) according to : - Bonneau et al., A&A 456, 789 (2006), with (V,B-V), (V,V-R), and (V,V-K) - Groenewegen, MNRAS 353, 903 (2004), with (V,V-K), (V,V-R), and (K,J-K) - van Belle et al., AJ 117, 521 (1999), with (PLX,V-K), where PLX is the parallax - di Benedetto, A&A 339, 858 (1998), with (V,V-K) - van Belle, PASP 111, 1515 (1999), with (V,V-K) and (B,V-K), for so-called "evolved" (giants and supergiants) and "variable" (Miras, carbon...) stars CATEGORY: angular diameter estimation INPUTS: need text input file with: - magnitude file name - basic data file name REQUIREMENTS: data read in magnitude file (in subdirectory of OBS_PATH/) must be in ASCII format and contain (in this order) : - effective wavelengths and widths [mum], magnitudes and errors basic data file (in subdirectory of OBS_PATH/) must be in standard FITS format and contain the following ASCII lines : PARALLAX = ppp.pp / Parallax (deg) PARA_ERR = e.ee / Error in parallax (deg) TEFF = tttt. / Effective temperature (K) TEFF_ERR = ttt. / Error in effective temperature (K) OUTPUTS: data written in output file(s) (of name(s) input_file_stem_BV.out, input_file_stem_VR.out, input_file_stem_VK.out, input_file_stem_BK.out in OBS_PATH/OUT/ directory) must be in ASCII format and contain (in this order, in same row) : - model flag, flux flag, initial synthetic exitance+radiance file names, data format flag, starting parameters, free flags, chi-square, output parameters, inferior errors, superior errors, number of degrees of freedom data written in covariance file(s) (of name(s) input_file_stem_BV.cova, input_file_stem_VR.cova, and/or input_file_stem_VK.cova, in OBS_PATH/OUT/ directory) must be in ASCII format and contain the covariance matrix(ces) (of size NB_MAX_PARAM*NB_MAX_PARAM) LOCAL PROCEDURE CALLED: RESET_PLOT LOCAL FUNCTION USED: EXTRACT_PROHEADER SYM TRI_STDEV LOCAL SYSTEM VARIABLE USED: FORMAT_OUT MAX_HEADER MAX_REL_ERR NB_MAX_PARAM OBS_PATH PHISUN_1PC PLOT_DEV SPIDAST_PATH REVISION HISTORY: Written by pcr 2009/03/25 modified by pcr 2010/07/16 ----- Documentation for Diam\get_MK_diameter.pro ----- NAME: GET_MK_DIAMETER AUTHOR: pierre.cruzalebes@oca.eu PURPOSE: calculate linear diameter from absolute luminosity and effective temperative derived from Morgan-Keenan spectral type according to - de Jager & Nieuwenhuijzen A&A 177, 217 (1987) - Straizys & Kuriliene ApSS 80, 353 (1981) - Allen's Astrophysical Quantities (4th ed.) with revised solar constants in Table 2.1 of - Amsler et al. Phys. Lett B, 667 (2008) deduce angular diameter from parallax if diameter file given, recalculate stellar parameters CATEGORY: angular diameter estimation. INPUTS: need text input file with: - basic data file name - log(Teff), log(L/L_sun), and log(g) file names (in same row) - diameter file name ('none' if not needed) REQUIREMENTS: basic data file (in subdirectory of OBS_PATH/) must be in standard FITS format and contain the following ASCII lines with: PARALLAX = ppp.ppp / Parallax (deg) PARA_ERR = e.eee / Error in parallax (deg) SPECTYP = 'sssss' / Morgan-Keenan spectral type data read in temperature/luminosity/gravity files (in ALL_DATA_PATH/ directory) must be in ASCII format and contain (tempe from de Jager, lumin from de Jager or Allen, gravity from Straizys or Allen) : - first row : values of luminosity continuous analytical b-variable - first column : values of spectral class continuous analytical s-variable - tabulated values of log(Teff)(s,b) or log(L/L_sun)(s,b) or log(g)(s,b) data read in diameter file (in subdirectory of OBS_PATH/ ) must be in ASCII format and contain (in this order) : - model flag, flux flag, synthetic exitance+radiance file names, data format flag, starting parameters, free parameter flags, goodness-of-fit parameter, input parameters, inferior errors, superior errors, confidence level OUTPUTS: data written in temperature file (of name input_file_stem.tempe in the same subdirectory of OBS_PATH/) must be in ASCII format and contain (in this order) : - effective temperature and error [K] data written in output file (of name input_file_stem.out in OBS_PATH/OUT/ directory) must be in ASCII format and contain (in this order, in same row) : - model flag, flux flag, initial synthetic exitance+radiance file names, data format flag, starting parameters, free flags, chi-square, output parameters, inferior errors, superior errors, number of degrees of freedom data written in covariance file (of name input_file_stem.cova in OBS_PATH/OUT/ directory) must be in ASCII format and contain the covariance matrix (of size NB_MAX_PARAM*NB_MAX_PARAM) LOCAL PROCEDURE CALLED: RESET_PLOT LOCAL FUNCTION USED: EXTRACT_PROHEADER TRI_STDEV LOCAL SYSTEM VARIABLE USED: ALL_DATA_PATH FORMAT_OUT MAS2RADIAN MAX_HEADER MAX_REL_ERR NB_MAX_PARAM OBS_PATH PHISUN_1PC PLOT_DEV SPIDAST_PATH REVISION HISTORY: Written by pcr 2009/04/17 last modification by pcr 2010/12/07 ----- Documentation for Diam\get_VK_diameter.pro ----- NAME: GET_VK_DIAMETER AUTHOR: pierre.cruzalebes@oca.eu PURPOSE: calculate angular diameter and effective temperature thanks to V-K color index according to Van Belle et al., AJ 117, 521 (1999) CATEGORY: angular diameter estimation INPUTS: need text input file with: - magnitude file name - basic data file name REQUIREMENTS: data read in magnitude file (in subdirectory of OBS_PATH/) must be in ASCII format and contain (in this order) : - effective wavelengths and widths [mum], magnitudes and errors basic data file (in subdirectory of OBS_PATH/) must be in standard FITS format and contain the following ASCII lines with: PARALLAX = ppp.pp / Parallax (deg) PARA_ERR = e.ee / Error in parallax (deg) OUTPUTS: data written in output file (of name input_file_stem.out in OBS_PATH/OUT/ directory) must be in ASCII format and contain (in this order, in same row) : - model flag, flux flag, initial synthetic exitance+radiance file names, data format flag, starting parameters, free flags, chi-square, output parameters, inferior errors, superior errors, number of degrees of freedom data written in covariance file (of name input_file_stem.cova in OBS_PATH/OUT/ directory) must be in ASCII format and contain the covariance matrix (of size NB_MAX_PARAM*NB_MAX_PARAM) LOCAL PROCEDURE CALLED: RESET_PLOT LOCAL FUNCTION USED: EXTRACT_PROHEADER TRI_STDEV LOCAL SYSTEM VARIABLE USED: FORMAT_OUT MAS2RADIAN MAX_REL_ERR NB_MAX_PARAM OBS_PATH PHISUN_1PC PLOT_DEV SPIDAST_PATH REVISION HISTORY: Written by pcr 2008/03/09 last modification by pcr 2010/06/20