SPectro-Interferometric Data Analysis Software Tool (SPIDAST)

Overview

SPIDAST is a spectro-interferometric calibration and interpretation software tool suite. It can be used to :

• estimate the angular diameters of the calibrator targets from indirect methods;
• calibrate the spectro-interferometric measurements, i.e. visibilities, fluxes, coherent fluxes, closure phases, and/or complex bispectral elements;
• interprete the measurements with wavelength-dependent spatial 1D and/or 2D (still to come) parametric models.

SPIDAST is a product of the ESOM team of the ​OCA-Fizeau Laboratory (UMR-CNRS 6525). SPIDAST has been developed using ​PV-WAVE since May 2006. Some procedures has been translated from Fortran or IDL.

Installing SPIDAST

1. Create the main directory SPIDAST/ and put the compilation files in it.
2. Create the subdirectories :
   • CALIB/ for spectrophotometric data calibration procedures,
   • DIAM/ for angular diameter estimation procedures,
   • FIT/ for model fitting procedures,
   • MANIP/ for spectrophotometric data manipulation procedures,
   • MODEL/ for modelling procedures,
   • PHOT/ for photometric procedures,
   • STAT/ for statistical procedures,
   • UTIL/ for utility procedures,

and put the program files in them as following :

SPIDAST/	SPIDAST/CALIB/	SPIDAST/DIAM/	SPIDAST/FIT/	SPIDAST/MANIP/	SPIDAST/MODEL/	SPIDAST/PHOT/	SPIDAST/STAT/	SPIDAST/UTIL/
compile	calib_coh.pro	average_diameter.pro	average_parameters.pro	average_bisp.pro	abs_eff.pro	average_magnitudes.pro	cov2cor.pro	angular_dist.pro
compile_calib	calib_flux.pro	get_broadband_diameter.pro	boot_fit.pro	average_data.pro	cent_path.pro	complete_flux_data.pro	get_outlier_fences.pro	bisp_proj.pro
compile_common	calibrate_bisp_source.pro	get_color_diameter.pro	create_fitting_batch.pro	average_observations.pro	central_depth.pro	convert_lim_data.pro	loca_scale_chi2.pro	complex2real.pro
compile_diam	calibrate_source.pro	get_MK_diameter.pro	fit_model.pro	concatenate_observations.pro	core_depth.pro	create_spec.pro	loca_scale_norm.pro	ct2lst.pro
compile_enviro	create_calibration_batch.pro	get_VK_diameter.pro	fit_sed.pro	get_broadband_flux.pro	core_path.pro	extinction.pro	mem_stdev.pro	date_conv.pro
compile_fit	create_data.pro	ld_poly_fit.pro		get_flux_scale.pro	delta_ring.pro	interstellar_correct.pro	plot_distri.pro	daycnv.pro
compile_instr	create_data_spec.pro	ld_poly_funct.pro		merge_bands.pro	engelke.pro	Jy2Wm2mu.pro	rect_stdev.pro	deg2dms.pro
compile_manip	cross_calibrate_bisp_source.pro	lm_fit.pro		merge_bisp_bands.pro	gaussian_disk.pro	magnitude2flux.pro	tri_stdev.pro	deg2hms.pro
compile_model	cross_calibrate_source.pro	poly_boot.pro		replace_bisp_mod.pro	inner_temp.pro	make_sed.pro	wavg.pro	euler.pro
compile_phot	get_bisp_response.pro	polybootw.pro		scale_flux.pro	marcs_disk.pro	scale_spectrum.pro	wstdev.pro	extract_proheader.pro
compile_stat	get_observations.pro	select_bisp.pro			model_data_funct.pro	wsterr.pro		functproj.pro
compile_util	get_response.pro	select_data.pro			model_data_funct_proj.pro	gaussian.pro
interpol_bisp_response.pro		shift_vis.pro			model_spectrum.pro	get_syn_ref.pro
interpol_observations.pro		nb_dens.pro				get_wave_list.pro
interpol_response.pro		planck.pro				helio_jd.pro
print_night_stat.pro		shell_depth.pro				heliocor.pro
shell_path.pro						juldate.pro
uniform_disk.pro						lorentzian.pro
make_tab3.pro
min_pdr_radius.pro
plot_spectra.pro
precess.pro
premat.pro
regutrap.pro
regutrap2.pro
reset_plot.pro
sunpos.pro
sym.pro
trap.pro
trap2.pro
trapez_fourier.pro
trapez_hankel.pro
trigo_deg_lib.pro
value_locate.pro
xyz.pro
ydn2md.pro
zparcheck.pro

SPIDAST runtime parameters

The SPIDAST suite of programs use the following parameters, configurable via the system variables defined in the 3 following compilation files :

compile_enviro

Parameter	Description	Default
PLOT_DEV	device used for plot	Win32 (WinXP) / X (Unix)
SYNTHE_PATH	absolute path to synthetic flux/intensity files (Planck, Engelke, Marcs...)	tbd by user
ALL_DATA_PATH	absolute path to general purpose data files	tbd by user
SPIDAST_PATH	absolute path to SPIDAST routines	tbd by user
OBS_PATH	absolute path to main directory containing observation data files	tbd by user

compile_instr

Parameter	Description	Default
GEOELEV	observing site elevation (in m)	2681. (VLTI)
GEOLAT	observing site latitude (in decimal deg)	-24.62742941 (VLTI)
GEOLON	observing site longitude (in decimal deg)	-70.40498688 (VLTI)
MAX_FLUX	maximal number of photoelectrons measured in raw flux	1000000
MAX_HEADER	maximal length of each header file	10000
MIN_FLUX	minimum number of photoelectrons measured in raw flux	0.001
NB_BAND_MAX	maximal number of spectral bands simultaneously measured	5
NB_LAMBDA_MAX	maximal number of wavelengths per spectral band	10000
NB_NIGHT_MAX	maximal number of nights allocated to each osberving run	10
NB_OB_MAX	maximal number of observing blocks per target	100
NB_OUV_MAX	maximal number of apertures simultaneously available	4
NB_PAIR_MAX	maximal number of pairs of apertures simultaneously available	NB_OUV_MAX*(NB_OUV_MAX-1)/2
NB_PART_MAX	maximal number of parts per night (each part defined by a given instrumental setting)	10
NB_TARG_MAX	maximal number of targets per part	20
NB_TRIPLE_MAX	maximal number of triplets of apertures simultaneously available	NB_OUV_MAX*(NB_OUV_MAX-1)*(NB_OUV_MAX-2)/6
SPEC_RES	instrumental spectral resolution	1500. (AMBER MR-K) / 35. (AMBER LR-JHK)

compile_common

Parameter	Description	Default
CHI2_LM_TOLE	chi-square tolerance for convergence of Levenberg-Marquardt procedure	0.001
DEL_FLAMBDA_LM	percentage of lambda to increase/decrease used in Levenberg-Marquardt procedure	10.
EPSILON	smallest positive quantity	1.e-15
FLAMBDA_LM	starting value of constant used in Levenberg-Marquardt procedure	0.001
FORMAT_OUT	format of data contained in the output parameter files	(1I8, 2A40, 1I8, 8F15.7, 8I8, 1E15.7, 8F15.7, 8F15.7, 8F15.7, 1I8)
LUMIN_STEP	spectral step of luminance data (in microns)	0.01
MARCS_RES	MARCS-model spectral resolution	20000.
MAS2RADIAN	conversion factor milliseconds of arc to radians	PI/180./3600.*0.001
MAX_LONG	highest unsigned long integer	30000
MAX_REL_ERR	maximal relative error for validity of the 2nd-order error propagation with algebraic approach	0.3
NB_BASE_STEP	number of steps per projected baseline element	0
NB_BOOT_LOOPS	number of bootstrap loops	1000
NB_BOOT_MIN	minimum number of bootstrap loops	10
NB_CORE	number of optical depth values in core direction	31
NB_DATA_BOOT	number of minimum data required for bootstrap calculation	5
NB_DATA_MAX	maximal number of measured data, used in fit	1000
NB_LM_ITER	maximal number of Levenberg-Marquardt iterations	100
NB_MAX_BOOTS	maximal number of bootstrap procedures	3
NB_MAX_PARAM	maximal number of model parameters	8
NB_RAD_MAX	maximal number of radii per wavelength	500
NB_RES_STEP	number of steps per spectral resolution element	2
NB_SHELL	number of optical depth values in shell direction	51
NB_STAR	number of optical depth values in stellar direction	11
PARAM_REL_DIST	defining parameter range [param-rel_dist*param,param+rel_dist*param] used for fit or test	0.2
REF_ANG_DIST	reference angular distance between science and calibrator (in deg), used for calibration weightening	15.
REF_FWHM	reference seeing fwhm (in arcsecond), used for calibration weightening	1.
REF_TAU0	reference coherence time (in second), used for calibration weightening	0.005
REF_TIME_DELAY	reference time delay between science and calibrator observation (in hour), used for calibration weightening	1.
RSUN2PARSEC	solar radius in parsec unit	2.254e-8
PHISUN_1PC	solar angular diameter (in mas) seen at 1 parsec distance	2*RSUN2PARSEC/MAS2RADIAN

Compiling SPIDAST

still to come

Running SPIDAST

still to come

Copyright

This Software "SPIDAST" is copyrighted OCA-Fizeau.

OCA-Fizeau holds all the ownership rights on the Software. The scientific community is asked to use this Software in order to test and evaluate it.

OCA-Fizeau freely grants the right to use this Software for research purpose only. Any use or reproduction of the Software to obtain profit or for commercial ends being subject to obtaining the prior express authorization of OCA-Fizeau.

Any commercial use made without obtaining the prior express agreement of OCA-Fizeau would therefore constitute a fraudulent imitation.

The Software being currently developed, OCA-Fizeau is assuming no liability, and should not be responsible, in any manner or any case, for any direct or indirect dammages sustained by the user.

Any user of the Software shall notify OCA-Fizeau of any comments concerning the use of the Software (e-mail: pierre.cruzalebes@…).