Max-Planck-Institut für extraterrestrische Physik

Recent changes in EXSAS spatial analysis software

Deutsche Version

MPE

HEG

X-Ray Astronomy

Wave

Exsas

X-Ray Astronomy

Wave

ROSAT

XMM-Newton

Chandra

Other Projects

Data Center

Exsas

Search

Impressum

Contact

Spatial Analysis

Recent changes in EXSAS spatial analysis software

Recent changes in EXSAS (* means: not copied to LOCAL yet).

1995
Date	Code	Execution	Remark
950206	ul3.fom	ul3.exe	Moderated error message NOT ANY SOURCES FOUND
950207	mkpsfg.fom	mkpsfg.exe	Internal write error and message output fixed
950207	corcts.for	sp2.exe	Average of vignetting correction factor
		ul3.exe	instead of vign. value itself
950217	cspio.inc	sp2.exe	New parameter VIGNET_TABLE
	cspio.cmn	ul3.exe	passed to spextr.for via common cspio
	spini.for	spipar.exe	If not present in parfile, then set to
	ulini.for	ulipar.exe	EXSAS_CAL:vignet.tbl
	sppari.for
	ulpari.for
	spipar.fom
	ulipar.fom
	spparc.for
	ulparc.for
	spextr.for
950220		crepar.prg	introduced keyword [extended]
		wripar.prg	(but disabled as parameter yet)
950220	spextr.for	sp2.exe	vign. table via vigcor only for PSPC
		ul3.exe	HRI and WFC are still done via offaxa
950220	spipar.fom	spipar.exe	new style of parameter file
	spparc.for	sp2.exe
		sppari.for
950221	cspio.inc	sp2.exe	support of experiment ROSAT,SURVEY
	cspio.cmn	spipar.exe
	spini.for
	sp2.fom
	sppari.for
	spipar.fom
	spparc.fom
		detsou.prg
		crepar.prg
950221	ul3.fom	ul3.exe	support of experiment ROSAT,SURVEY
	ulipar.fom	ulipar.exe	and new style of parameter file
	ulini.for
	ulpari.for
	ulparc.for
950227	exsmil.for	exsm.exe	bug in closing/opening tables
	exsmim.for		(if lslst1 or mlslt1 empty)
950229	csmio.cmn	exsm.exe	support of experiment ROSAT,SURVEY
	csmio.inc	smipar.exe	and new style of parameter file
	exsm.fom	detsou.prg
	smipar.fom
	smini.for
	smparc.for
	smpari.for
950229	cdmio.cmn	exdm.exe	support of experiment ROSAT,SURVEY
	cdmio.inc	dmipar.exe	and new style of parameter file
	exdm.fom	detsou.prg
	dmipar.fom

	dmini.for
	dmparc.for
	dmpari.for
950302	getfnc.for	sp2.exe	avoided real overflow
950302	cdlio.cmn	exdl.exe	support of experiment ROSAT,SURVEY
	cdlio.inc	dlipar.exe	and new style of parameter file
	exdl.fom	detsou.prg
	dlipar.fom
	dlini.for
	dlparc.for
	dlpari.for
950305	fo_wri.for	foldo.exe	bug: divison by zero
950306	csmdat.cmn	exsm.exe	checks for maximum number of input
csmio.cmn	smipar.exe	or output source list entries
	csmio.inc
	exsmil.for
	exsmim.for
	tsmlis.for
950307	cdbio.cmn	exdb.exe	support of experiment ROSAT,SURVEY
	cdbio.inc	dbipar.exe	and new style of parameter file
	exdb.fom	detsou.prg
	dbipar.fom
	dbini.for
	dbparc.for
	dbpari.for
950307		detsou.prg	mask defaults for PSPC,HRI,WFC
		detectsour.hlq
950309	imini.for	imipar.exe	defaults for ROSAT,SURVEY
	dlimag.fom	dlimag.exe	bug in calling imini from dlimag
950309	blind.for	exdb.exe	removed old common csool, replaced by common csoio.
	edbout.for
	exdb.fom		sooli.for, soolri.for,csool.cmn are now obsolete
950309		detsou.prg	simplified
		creima.prg
950309	csoio.inc	soipar.exe	removed obsolete parameters
	csoio.cmn	crepar.prg	flagout, fwtbl, sbpar, sdpar.
	soini.for	wripar.prg	sopar.par is now detector-independent
	sopari.for		fillfw.prg now obsolete
	soparc.for
	soipar.fom
950309		detsou.prg	simplified
		creima.prg
950309		crepar.prg	output of experiment suppressed for task SO
950309	cdbio.cmn	exdb.exe	new parameters, removed old ones
	cdbio.inc	dbipar.exe	db_fof, dbfwhm
	dbini.for
	dbpari.for
	dbparc.for
	dbipar.fom
	edbout.for
	blind.for
	csmio.cmn	exsm.exe
	csmio.inc	smipar.exe	sm_fof, smfwhm
	smini.for
	smpari.for
	smparc.for
	smipar.fom
	esmout.for
	identf.for
	cspio.cmn	sp2.exe
	cspio.inc	spipar.exe	sp_fof, spfwhm
	spini.for
	sppari.for
	spparc.for
	spipar.fom
	spextr.for
	sp2.fom
	ulini.for	ul3.exe
	ulpari.for	ulipar.exe
	ulparc.for
	ulipar.fom
	ul3.fom
950314		creima.prg	automatic recognition of survey mode fill parameter file correctly
950314		createparf.hlq	help cre/par: [extended] option
		create.hlc	update examples section
950315	exdbin.for	exdb.exe	Parameters SKYPIX_PER_IMAPIX and
	dbipar.fom	dbipar.exe	SKY_PIX_SIZE now read from image
	dbini.for		descriptors BIN_SIZE and SKY_PIX_SIZE
	dbparc.for
	dbpari.for
950315	exsmin.for	exsm.exe	Parameters SKYPIX_PER_IMAPIX and
	smipar.fom	smipar.exe	SKY_PIX_SIZE now read from bacmp1
	smparc.for		descriptors STEP(1) and
	smpari.for		SKY_PIX_SIZE
950315	sp2.fom	sp2.exe	CSPIO Parameters PIXIM, PIXEL, DETPIX,
	spipar.fom	spipar.exe	DETCEN, and SKYCEN now read from
	spparc.for		descriptors of background image bacmp1
	sppari.for		instead of reading from parameter file
950315	ul3.fom	ul3.exe	CSPIO Parameters PIXIM, PIXEL, DETPIX,
	ulipar.fom	ulipar.exe	DETCEN, and SKYCEN now read from
	ulparc.for		descriptors of background image bacmp1
	ulpari.for		instead of reading from parameter file
950315		detsou.prg	Further simplified. Because of the above described changes it is not necessary to pass these instrument-specific parameters to the parameter files
950315		detloc.prg	Simplified call of CREATE/PARFIL SO
		crebgi.prg
		detmap.prg
		mersou.prg
		detmax.prg
		comupp.prg
950317	spextr.for	sp2.exe	Introduced vignet_hri.tbl for HRI vignetting
	spini.for	spipar.exe	Changed vignet.tbl to vignet_pspc.tbl
	ulini.for	ul3.exe
		ulipar.exe
950320	sp2.fom	sp2.exe	Removed extraction of mis_id and det_id
	ul3.fom	ul3.exe	from FWHM table (because name of FWHM is not defined at this time).
950320	sppari.for	sp2.exe	Warning message if old-style parameter
	ulpari.for	ul3.exe	file is encountered
	smpari.for	exsm.exe
	dmpari.for	exdm.exe
	dbpari.for	exdb.exe
	dlpari.for	exdl.exe
		spipar.exe
		ulipar.exe
		smipar.exe
		dmipar.exe
		dbipar.exe
		dlipar.exe
950321	sp2.fom	sp2.exe	Removed forgotten write statement
950321	getvig.for	sp2.exe	HRI compatibility: vignetting table need not
		ul3.exe	have at least 256 rows. Close vigtab correctly.
950323		detsou.prg	Defaults for OFF_AX_MAX now depending on detector.
		detectsour.hlq
950323	spini.for	sp2.exe	Amplitude defaults for HRI 0-15 instead of 1-16
	ulini.for	spipar.exe
		ul3.exe
		ulipar.exe
950323	dlcmd.for	dlimag.exe	Lowest allowed value for HRI lower amplitude is zero instead of 1 (parfile reading)
950323		crepar.prg	Removed warning for task EXUL / ROSAT,WFC. WFC defaults are available.
950328	ul3.fom	ul3.exe	Text output for debugging, actvated by flag within program
950329	spcs01.for	prisma.exe	(Spectral package) Integrated photon flux is now written to keyword OUTPUTR(13)
950405	exdbin.for	exdb.exe	Fixed bug: read double prec descriptors SKY_PIX_SIZE and BIN_SIZE using single prec variables. MIDAS doesn't like that.
	edbout.for
950405	exsmin.for	exsm.exe	changed some statements
	sp2.fom	sp2.exe	"single prec var" = "double prec var" to "spv" = SNGL("dpv")
	ul3.fom	ul3.exe	No problems for our compiler, but maybe for others.
950405	boxpro.fom	boxpro.exe	Removed unneccessary input parameters
	linout.for	drabox.prg	[bins] and [limit]. Removed nonsense
		drawbox.hlq	column :STEP_SIZE in output table
		draw.hlc	linprof.tbl.
950407		ul3.exe	Compiled with -O0 option. Compiler optimisation does something strange with sub laguer.for.
950410		conong.prg	Replaced call of home/ajv/exe/sendplot.prg
		conspl.prg	with EXSAS_PRG:conong.prg
950410	conrad.fom	conrad.exe	New command CONTOUR/RADEC_GRID by ajv.
	makefile	conrad.prg	Contour plot onto RA/DEC-grid.
		conong2.prg	Output in new_plot.ps.
		conracg.prg	Not exportable.
		contourade.hlq	Sourcecodes in /spatial/contour.
		contou.hlc
		spat.prg
950410		crenor.prg	New command CREATE/NORTH_ALIGNED_IMAGE
		createnort.hlq	by ajv. Create north-aligned image from
		create.hlc	digitzed sky image.
		spat.prg	Not exportable.
950410		conracg.prg	CONTOUR/RADEC now (probably) compatible with all non-512x512 images.
950411	*commandsum.tex		Updated: new commands CREATE/NORTH, CONTOUR/RADEC_GRID
950412	*chapter5.tex		Updated: new-style parameter files, vignetting correction, some other items
950419		distrib.prg	Corrected decision SCC (switch from PSPC C to B) to 20530000
950419		distrib.prg	Put under SCCS
		dist_rib.prg
950420		pripro.prg	Added -h option to lpr: no header pages
950427	traequ.fom	traequ.exe	Simplified. No limit of row-# anymore.
950509		allrfp.prg	GET/PROT now tests for additional files, and converts them to add1.ps ... add3.ps if they exist. This is only a temporary workaround for RODs with nonstandard file contents.
950510	boxpro.fom	boxpro.exe	Internal filename for use in error-output was wrong.
950601	*addmap.for	expma4.exe	Complete reworking of the shift, scaling
	*ex_att.for		rotating, and wobbling of the instrument
	*ex_acc.for		maps. Now all transformations use reals
	*ex_cre.for		and are combined in improt.for, except
	*ex_che.for		the (whole-pixel) adding of the wobbled
	*ex_cmd.for		instruments.
	*imprea.for		This was only possible for a fixed attitude
	*improt.for		rebin in x and y, set now to the size of
		create_exposure_map.par	the reference image pixels.
		create_exposure_map.par_e	Not HRI compatible anymore.
950602	* crewob.fom	crewob.exe	New command CREATE/WOBBLED_MASK
		create.hlc
		crewob.prg
		createwobb.hlq
		*makefile
		spat.prg
950602	*sodraw.fom	sodraw.exe	Hidden feature: If INPUTI/I/6/1 = 12345
		drasou.prg	then the screen output is suppressed. Implemented on request of Thomas Boller. (Changed already in March, but erroneously not implemented in /LOCAL then).
950731		detloc.prg	Problems with WRITE/PARFIL 'SOPAR':
		crebgi.prg	the short form by addressing the parameter
		detmap.prg	by its internal variable name is no
		mersou.prg	longer allowed by Carlo's new WRITE/PARF.
		detmax.prg	Inserted the true parameter names
		comupp.prg	(This only happened when writing to sopar).
		creima.prg
		seasou.prg
950824	*sodraw.fom	sodraw.exe	Replaced write(6,*) with call sttput(... Makes change of 950602 obsolete. Screen output can now be disabled with LOG(4)=2.
950905	*restor.for	exdb.exe	Converted bgr image output of CREATE/BG
	*edbout.for		from integer to real (only array iback).
	*cdbimg.cmn
950913		trans3.prg	TRANSF/COORD: New default for parameter
		transfcoor.hlq	[ampl] (channel 30) to make it compatible with DETECT/MAXLIK and COMPUTE/UPPER.
951006		crenor.prg	CREATE/NORTH_ALIGNED_IMAGE fixed division-by-zero error
951017	*bdfr4w.for		Error occurred when image to write to didn't exist before. Fixed now.
951017	*cdbimg.cmn	exdb.exe	CREATE/BG_IMAGE
	*exdb.fom		Full conversion of integer to real arrays.
	*bkgimg.for		Variable binsizes for better handling of
	*blind.for		low-intensity images.
	*cheese.for
	*edbout.for
	*exdbin.for
	*restor.for
951018	*edmout.for	exdm.exe	DETECT/MAP Variable 'flagout' in edmout was no longer defined by the sopar common.
951018	*tdmlis.for	exdm.exe	DETECT/MAP Changed format of column BACKG in mslst from 9.2 to 9.4. Necessary to be able to work with very low backgrounds.
951120	*uslstr.for	sp2.exe	Serious bug: BLUR was initialized always
		ul3.exe	with fwrops.tbl (PSPC), if uslstr was called (always if START_POSITION=user_list). This caused every subsequent call of blur to deliver PSPC blur radii.
951121	*getexp.fom	getexp.exe	READ/EXPOSURE
		read_expo.hlq	Streamlined, removed some arrays which are not necessarily needed, and set array dimensions (mxrow) to 262200.
951205	*sp2.fom	sp2.exe	DETECT/MAXLIK
	*ul3.fom	ul3.exe	diagnosis output in tables usphot.tbl
	*spexcl.for		and diagnose.tbl, if parameter FULL_OUT_FLAG
	*maxlik.for		in parameter file is set to TRUE.
	*sourex.for		usphot.tbl contains source number, xpos, and ypos of photons for each source diagnose.tbl contains various important values for each source, for test issues (only works for DET/MAX, not for COMP/UPPER)
951206	*sp2.fom	sp2.exe	DETECT/MAXLIK, COMPUTE/UPPER
	*ul3.fom	ul3.exe	New parameter EXCLUDE_SIGMA ( = variable
	*spini.for		SP_EXC). Determines the threshold above
	*ulini.for		which sectors with eventual neighbouring
	*sppari.for		sources within the cut radius are excluded.
	*ulpari.for		Up to this time this was hardcoded as
	*spexcl.for		5 sigma in spexcl.for.
		*cspio.cmn
		*cspio.inc
		*sppar_h.par
		*sppar_he.par
		*sppar_p.par
		*sppar_pe.par
		*sppar_s.par
		*sppar_se.par
		*sppar_w.par
		*sppar_we.par
		*ulpar_h.par
		*ulpar_he.par
		*ulpar_p.par
		*ulpar_pe.par
		*ulpar_s.par
		*ulpar_se.par
		*ulpar_w.par
		*ulpar_we.par
951211	*bkgimg.for	exdb.exe	CREATE/BG_IMAGE
	*blind.for		calculations/maps now in REAL instead
	*cheese.for		of INTEGER
	*edbout.for		variable rebinning for bad photon
	*exdb.fom		statistics
	*exdbin.for		set to constant for very bad photon
	*restor.for		statistics
	*cdbimg.cmn
951211	*exdmin.for	exdm.exe	DETECT/MAP
	*mdetec.for	dmipar.exe	all background calculations done now
	*msrbkg.for		in REAL instead of INTEGER
	*msnrmp.for
	*tdmlis.for
	*cdmimg.cmn
951212	*esmout.for	exsm.exe	MERGE/SOURCES
	*exsmin.for		background map now in REAL instead of INTEGER
	*tsmlis.for		Max. no. of rows of output table now 10000 instead of 1000. Format of BACKG column now f11.5
951214		sppar_h.par	DETECT/MAXLIK, COMPUTE/UPPER_LIMIT
		sppar_he.par	Parameter EXCLUDE_SIGMA moved to
		sppar_p.par	the extended parameter set.
		sppar_pe.par
		sppar_s.par
		sppar_se.par
		sppar_w.par
		sppar_we.par
		ulpar_h.par
		ulpar_he.par
		ulpar_p.par
		ulpar_pe.par
		ulpar_s.par
		ulpar_se.par
		ulpar_w.par
		ulpar_we.par
951214		creima.prg	CREATE/SOURCE amplitude default set to '1-256' instead of 'all'. Reason: projection interpretes 'all' as 1-512.
951215	*mlcnt1.fom	mlcnt1.exe	FIT/MULTI_SOURCE
	*cgausn.for		map arrays now REAL instead of INTEGER
	*mlcout.for
951215	*tdllis.for	exdl.exe	background columns in source list output
	*tdmlis.for	exdm.exe	now set to unified format : f10.5
	*mlcout.for	mlcnt1.exe	(to be able to deal with very low backgrounds)
	*tsmlis.for	exsm.exe
	*sp2.fom	sp2.exe
	*tmlli1.for	ul3.exe
	*tsoli1.for
951215		dbpar_he.par	CREATE/BG, FIT/MULTI, DETECT/MAP,
		dbpar_pe.par	MERGE/SOURCES, DETECT/MAXLIK,
		dbpar_se.par	COMPUTE/UPPER
		dbpar_we.par	defaults for BG_SCALE in parameter
		mcpar.par	files set to 1.0
		dmpar_he.par
		dmpar_pe.par
		dmpar_se.par
		dmpar_we.par
		smpar_he.par
		smpar_pe.par
		smpar_se.par
		smpar_we.par
		sppar_he.par
		sppar_pe.par
		sppar_se.par
		sppar_we.par
		ulpar_he.par
		ulpar_pe.par
		ulpar_se.par
		ulpar_we.par
951215	*dbini.for	exdb.exe	default values for background scale factor set to 1
	*mcini.for	mlcnt1.exe
	*dmini.for	dbipar.exe
	*smini.for	mcipar.exe
	*spini.for	exdm.exe
	*ulini.for	dmipar.exe
		exsm.exe
		smipar.exe
		sp2.exe
		spipar.exe
		ul3.exe
		ulipar.exe
		detectsour.hlq
951215		backma.prg	CREATE/SMOOTHED_BG
		createsmoo.hlq	Background scale factor = 1.0
			Add descriptor BG_SCALE to output image
951215		creima.prg	CREATE/SOURCE_DETECT_IMAGE
		createsour.hlq	Parameter pet now defaulted to events.tbl,
		create.hlc	parameter ampl new default 1-256
951215		detsou.prg DETECT/SOURCES
		detectsour.hlq Parameter ima now defaulted to image1.bdf
		detect.hlc parameter pem now defaulted to events.tbl

1996
960110	*spipar.fom	spipar.exe	Fixed syntax error which made it
	*ulipar.fom	ulipar.exe	uncompilable on certain machines. 'READ(VALUE,*) (DUMMY(1))'
960115		detsou.prg	DETECT/SOURCES Fixed error (wrong binsize) when called as DET/SOU create events for HRI data.
960116		sppar_h.par	DETECT/MAXLIK, COMPUTE/UPPER
		sppar_he.par	HRI lower amplitude default set to 1
		ulpar_h.par	instead of 0
		ulpar_he.par
960129		crenor.prg	CREATE/NORTH_ALIGNED_IMAGE Bugfixes by Andreas Vogler: Works now also for small rotation angles (bypasses REBIN/ROTATE). Format problem (introduced with the new MIDAS SET/FORMAT convention) during filling of descr. POINT_LONG and POINT_LAT fixed.
960129	*ex_att.for	expma4.exe	CREATE/EXPOSURE
		expma3.prg	Add X-offset and Y-offset to attitude X and Y-values, to compensate for off-center reference images. This had worked in the old version, but was not possible to implement in the same way as before due to the new integer binning of the attitude histogram (see changes of 950601).
960201		detsou.prg	DETECT/SOURCES Change amplitudes to INTEGER before writing them to parameter file. (REAL amplitudes in parfile caused trouble when run on a HP)
960207	*getvig.for	sp2.exe	DETECT/MAXLIK, COMPUTE/UPPER
		ul3.exe	Inserted missing SAVE statement. Had no consequences on OSF, but was running wild on HP machines.
960221	*spextr.for	sp2.exe	DETEXT/MAXLIK, COMPUTE/UPPER
		ul3.exe	rewrote algorithm which selects the photons to accept in case the photon count exceeds the array sizes IDM_SP. Fixed a bug in calculating the correction factor FA_PH (was too small for sources with 4000-8000 counts within the extraction circle).
960228	*meratt.fom	meratt.exe	MERGE/FILES rewrote roll angle transformation in merge-atttitude part - the old version using hpdet.for was in error.
960313	*spextr.fom	sp2.exe	DETECT/MAXLIK, COMPUTE/UPPER
		ul3.exe	Fixed bug: since modification on 960221 the background was corrected with FA_PH before FA_PH was even calculated. The consequence was that the background of each source was corrected with the FA_PH of the one preceding in the source list. Influence on results: the first source in each source list was thrown out, the bgr of each very bright source and the following was wrong. Further modification now: the correction is only applied if photons have been filtered out.
960412	*uslstr.for	sp2.exe	DETECT/MAXLIK, COMPUTE/UPPER If a userlist contained the column :EXT, the extraction circles were supposed to be enlarged by its values. However, the default for EXTENT_THR of 0 threw out those sources without warning. If one e.g. used a output source list as user-list, it contained this column. Reading of EXT_MPX from user-list disabled now. Initialization of SRC_MPX fixed.
960425	*ex_att.for	expma4.exe	CREATE/EXPOSURE Deadtime correction factor was applied wrong (multiplied instead of divided). Fixed.
960425	*sp2.fom	sp2.exe	DETECT/MAX, DETECT/SOURCES, COMPUTE/UPPER
	*ul3.fom	ul3.exe	Check if descriptor SC_MODE in events table and parameter OBS_MODE in parameter file match. If not, error message and stop. Necessary because SC_MODE overrides OBS_MODE.
960430	*sp2.fom	sp2.exe	DETECT/SOURCES, COMPUTE/UPPER
	*ul3.fom	ul3.exe	Up to now the number of continuation lines to skip in the merged_list (not in the user_list) was determined by the value of :COFL of the last valid entry. This forbid sorting or selecting of the merged list. Now simply each line with :COFL=0 is skipped.
960531	*simsou.fom	simsou.exe	New command SIMULATE/SOURCE
	*simatt.for		Creates simulated observation dataset
	*simcat.for		containing one source.
	*simcer.for		Command is for test issues and will be
	*simcev.for		not exported.
	*simevr.for		Up to now only available in SCCS environment.
	*simexh.for
	*simrbg.for
	*simsrc.for
		simsou.prg
		spat.prg
		simula.hlc
		simulasour.hlq
960613	*exsm.fom	exsm.exe	MERGE/SOURCES
	*getnew.for		The output for SIGNAL of the merged
	*polish.for		source is now the maximum of the sources
	*compar.for		with the highest priority which had been
	*esmout.for		merged.
960617	*getvig.for	sp2.exe	DET/MAX, COM/UPP, SIM/SOU
	*simsrc.for	ul3.exe	getvig.for did only initialization at first
		simsou.exe	call, no actual calculation. That means, the first photon got a wrong vignetting. Fixed - now calculation already at first call.
960617	*spextr.for	sp2.exe	DET/MAX, COM/UPP
	*makefile	ul3.exe	replaced PSF-sigma subroutine psfex.for with new fast sigsq_g.for. psfex.for had also too small telescope resolutions compared to the formulae documented in blur.for or in sigsq_g.for.
960621	*compar.for	exsm.exe	MERGE/SOURCES
	*exsmil.for		Sorting order for list priorities and SNR
	*exsmim.for		was wrong. Fixed. Small bug in sorini.for
	*sorini.for		for array dimension = 1. Bug was in since 1991, SASS software is ok. Should now be faster, some positions more reliable.
960625		crenor.prg	CREATE/NORTH_ALIGNED_IMAGE Inserted abs function of rotation angle in checking for very small angles
960702	*sp2.fom	sp2.exe	DET/MAX, COMP/UPPER, DET/SOURCE
	*ul3.fom	ul3.exe	Flux error is now corrected for FA_PH (photon selection factor). Up to now the y-error was erroneously corrected with FA_PH! Only sources with cts>4000 affected.
960718		expo.par	CREATE/EXPOSURE
		createexpo.hlq	Help text changed to warn that image size of reference image must be 512 x 512 pixels.
960718	*clslst.cmn	exdb.exe	CREATE/BG_IMAGE, FIT/MULTI,
	*db/blind.for	mlcnt1.exe	MERGE/SOURCES
	*db/exdb.fom	exsm.exe	Replaced integer ISN_LSL with
	*db/exdbin.for		real SNR_LSL (value of column
	*general/lslstr.for		SIGNA in local source lists).
*	mlcnts/mlcnt1.fom		No changes in DL yet.
*	sm/exsmil.for
960718	*cmslst.cmn exsm.exe	MERGE/SOURCES
	*general/mslstr.for	Replaced integer ISN_MSL with real
	*sm/exsmim.for	SNR_MSL (value of column SIGNA in map detect source lists). No changes in DM yet.
960718	*cdlimg.cmn	exdl.exe	DETECT/LOCAL
	*rclnr.for		All variables and calculations dealing
	*exdl.fom		with SNR now done in real instead of
	*ldetec.for		integer. Fixes overflow problem with
	*snrmap.for		bright sources.
	*match.for		Table output of SIGNA in lslst1.tbl
	*cmass.for		now in f9.2 format.
	*edlout.for
	*tdllis.for
960718	*cdmimg.cmd	exdm.exe	DETECT/MAP
	*edmout.for		All variables and calculations dealing
	*tdmlis.for		with SNR now done in real instead of
	*exdm.fom		integer. Fixes overflow problem with
	*mdetec.for		bright sources.
	*msrbkg.for		Table output of SIGNA in mslst1.tbl
	*msnrmp.for		now in f9.2 format.
960722		crenor.prg	CREATE/NORTH_ALIGNED_IMAGE Put in some checks for the proper range of right ascension, and difference in right ascensions. Near 0h there were cases where the difference in RA between object and rotation reference point was >360. deg, leading to wrong rotation angles.
960724	*diskcnt.fom	diskcnt.prg	PSF/CUT
	*dprob.for	spat.prg	preliminary installation of PSF-overlap
		diskcnt.exe	calculation PSF/CUT by G.Boese
960902	*intpsf.fom	intpsf.exe	CALCULATE/PSF corrected misleading text output
		calculpsf.hlq	added warning about non-existent selection for accepted time intervals for the attitude table
960905		intrin.prg	INTEGRATE/RING fixed cursor circle size for shrunken (scale negative) images.
960906	*srcbkg.for	exdl.exe	DETECT/LOCAL corrected spurious detections with background < 0 (has no effect on :SIGNAL, but on :COUNT in the lslst1 table)
961024	*mdetec.for	exdm.exe	DETECT/MAP
	*msnrmp.for		Faulty likelihood calculation rewritten.
	*msrbkg.for		Replaced subroutines POISLK, SNRFIT and
	-poislk.for		ERRFC with POISL2, GAMMP and GAMMLN.
	-snrfit.for		Output positions are now much more
	-errfc.for		accurate.
	+poisl2.for
	+gammp.for
	+gammln.for
961111	*gammp.for	exdm.exe	DETECT/MAP ITMAX increased to 1000, to be able to deal with extremely high backgrounds. High numbers of iterations for calculating the Poissonian likelihood are necessary for high backgrounds with the measured counts approaching the background itself.
961112	*trans5.fom	*trans5.exe	TRANSFORM/COORDINATES
		*transfcoor.hlq	Fixed calculation of distance to next source (faulty unit conversions). Max. number of table rows now limited to 100,000 instead of 300,000. Rewrote some less-than-clear explanations in the on-line help.
961115	*trarad.fom	*trarad.prg	TRANS/COORD,TRANS/RADEC
	*tracoo.fom	+tracoo.prg	names change to clear up the confusion a bit/TD>
	-trans5.fom	-trans3.prg	transy.prg --> trarad.prg
		-transy.prg	tracoo.fom --> trarad.fom
			trans3.prg --> tracoo.prg
			trans5.fom --> tracoo.fom
			Older versions of trarad.prg and trarad.fom were still existing. transy.prg, trans3.prg, and trans5.fom are now obsolete.
961121	*tracoo.fom	tracoo.prg	TRANSFORM/COORD
		tracoo.exe	use now subroutine FWHM_G instead of BLUR for calculation of FWHM. Fully compatible with HRI, WFC and SURVEY. Interpreting the mission- and detector-related descriptors is now done in the FORTRAN code instead of the MIDAS prg. Some further minor changes.
961210		creima.prg	CREATE/SOURCE_DETECT_IMAGE
		createsour.hlq	PSPC pointing amplitude defaults changed from 1-256 to 8-256.
961210	*sp2.fom	sp2.exe	DETECT/MAXLIK (,COMPUTE/UPPER)
	*spini.for	ul3.exe	PSPC pointing amplitude defaults
	*ul3.fom	sppar_p.par	changed from 1-256 to 8-256.
		sppar_pe.par	For COMP/UPP those defaults were already active, the correction was necessary only for fixing a small inconsistency.
961210	*sp2.fom	sp2.exe	DETECT/MAXLIK, COMPUTE/UPPER
	*ul3.fom	ul3.exe	Check energy ranges (selected vs background energy ranges): no check for PSPC (pointing) amplitudes below channel 8

1997
970113	*simsou.fom	simsou.exe	SIMULATE/SOURCE
		simulasour.hlq	Added some error messages
		simsou.prg	(ERROR creating table xxx...),
		spat.prg	HRI completely defaulted now. Released.
970113	*aveoff.fom	aveoff.prg	AVERAGE/OFF_AXIS_ANGLE
	*offavm.for	aveoff.exe	New command. Released.
		averagoff_.hlq
		averag.hlc
		spagen.a
		spat.prg
970115	*gammln.for	exdm.exe	Fixed error in application of reflection
		dm.a	formula for calculation of ln(gamma). No consequences for source detection since all treated counts are >=1 and reflection formula would only be used for 0<a<1.
970115	*conpro.fom	conpro.exe	CONVERT/PROBABILITY
	*gammq.for	conpro.prg	New command. Released.
	*makefile	prob.a
		conver.hlc
		converprob.hlq
		spat.prg
970116	*conpro.fom	conpro.exe	CONVERT/PROBABILITY Added warning about applicability reliability of conversions for EXSAS source detection) to screen printout
970117		averagoff_.hlq	AVERAGE/OFF_AXIS_ANGLE Added warning note to on-line help.
970121		crenor.prg	CREATE/NORTH_ALIGNED_IMAGE Bug in converting declinations of image center for declinations between 0 and -1 deg. Both alignment and descriptor POINT_LAT of the output image were in error. Fixed now.
970203		creima.prg	CREATE/SOURCE in command-line mode:
		createsour.hlq	Up to now there was no special HRI amplitude default in commandline mode. This had caused no problems until 961210 when the PSPC default had been changed from 1-256 to 8-256. From then on using CRE/SOU for HRI data resulted in a amplitude default 8-15, if no amplitude was explicitely specified. Fixed: CRE/SOU now uses detector-specific amplitude defaults.
970204	*bkgimg.for	exdb.exe	CREATE/BG_IMAGE
	*splin2.for		Systematic errors eliminated by setting statistical weights to 1.0 and 0.5 (replaced bin). No minimum limit for number of photons per bin anymore. Log-mode disabled. This is only a temporary fix until we have tested a more reliable spline fit routine.
970214	*offavm.for		General routine for averaging off-axis angles. Fixed bug: attitude table was not closed. The logics of error handling is still in need of a major rework. Only preliminarily fixed now.
970225		createbg_i.hlq	CREATE/BG_IMAGE Updated on-line help file. Added notes about the use of masks and exposure maps.
970225		sppar_he.par	DETECT/MAXLIK parameter files Changed values of extended parameters BG_SPECTRUM and BG_RATE from 'bacsp.tbl' and 'bacrt.tbl' to 'none'.
970307	*edbout.for	exdb.exe	CREATE/BG Copy descriptor TIM_SEL from input image to background map
970320	*sp2.fom	sp2.exe	DETECT/MAX, COMPUTE/UPPER
	*ul3.fom	ul3.exe	Copy descriptor TIM_SEL from events table to source lists
970320	*cresim.fom	cresim.exe	CREATE/SIMULATED_IMAGE
	*errf.for	cresim.prg	New command: create a simulated image
	*poidev.for	reco.a	from a source list (ML-detection output).
	*recadd.for	spat.prg	Useful for testing the results of a
	*recbac.for	create.hlc	completed source detection task.
	*recgau.for	createsimu.hlq
	*recima.for
	*recmsk.for
	*recout.for
	*recpoi.for
	*recsol.for
	*makefile
970407		integrring.hlq	INTEGRATE/RING, INTEGRATE/BOX
		integrbox.hlq	Added note: 'image pixels' refer to MIDAS image pixels, and not to the EXSAS image pixels coordinate system (X_IMA, Y_IMA) which is flipped in y-direction.
970408	*dlimag.fom	dlimag.exe	CREATE/SOURCE_DETECTION_IMAGE Added warning message, if HRI amplitude band is not the whole band (measured fluxes may be in error due to spatial gain variations over the HRI detector).
970416	expo.par		CREATE/EXPOSURE_MAP
	expo_e.par		Default for instrument map changed from the rather untypical moimp_8_19_b to the more typical moimp_42_51_b. Ideal would be a broadband IM, but such an animal is not available.
970421	*/sim/makefile		SIM/SOURCE, CREATE/SIM
	*/reco/makefile		Removed link to library spectral.a in makefiles since random number generator ran1.o is now in general.a
970507	*iditxc.for		GET/COORD, DRAW/COORD, DRAW/SOURCE
	*gridc.fom	gridc.exe	With new MIDAS version 96NOV greek
	*ploli1.fom	ploli1.exe	letters appeared in display labelling
	*sodraw.fom	sodraw.exe	(font size 3 is now special, greek font) Fixed (set to 0). Changed also definition of T_OR (text angle) from double prec to Integer, as required by MIDAS IIGTXT.C)
970514	*tracoo.fom	tracoo.exe	TRANSFORM/COORDINATES Now program doesn't shut down anymore if descriptor FILTER_ID is missing (exception is WFC of course)
970514	*getnew.for	exsm.exe	MERGE/SOURCES Replaced IEOR function for merging of priority codes with IOR. E.g. 16+16 now is correctly 16 instead of 0. This improves merging within the same list or merging of lists of the same priority. IOR was originally in the software, but was replaced in 1991 with IEOR for unknown reasons.
970522	*equsky.for	trarad.exe	TRANS/RADEC, DET/MAXLIK, COMPUTE/UPPER
	*deqsky.for	sp2.exe	replaced call of BLUR with FWHM_G in
	*uslstr.for	ul3.exe	deqsky and equsky. Therefore parameters
	*trarad.fom	trarad.prg	DET_ID and FIL_ID are now required
	*trans/makefile		instead of FWTBL (FWHM table). No consequences for the user.
970606	*sp2.fom	sp2.exe	DETECT/MAX, COMPUTE/UPPER
	*ul3.fom	ul3.exe	New scaling of flux error for bright (>4000 cts) sources. Scaled now with sqrt(FA_PH). This is justified by the fact that the flux fit error is dominated by source photon statistics for bright sources. This modification supersedes the modification of 960702.
970606	*recima.for	cresim.exe	CREATE/SIMULATED_IMAGE Mask array was defined in recmsk.for, with sizes passed as parameters to this routine. Some compilers weren't too happy with this. Fixed now.
	*recmsk.for
970618	*getamp.for	dlimag.exe	CREATE/SOURCE_DETECT_IMAGE Fixed bug in counting the number of amplitude ranges in parameter input. There also was a problem for amplitude ranges like 0-12. They were set to 0-15. This bug was related with the one noted above.

Spatial Analysis

Recent changes in EXSAS spatial analysis software

1995

1996

1997