indexamajig: Use static wavelength from DataTemplate as estimate

If the geometry file contains a literal wavelength, then let's use it.
That's better than nagging the user for a --wavelength-estimate.

1 files changed, 11 insertions, 0 deletions

diff --git a/src/indexamajig.c b/src/indexamajig.c
index e4da883a..6980c4b7 100644
--- a/src/indexamajig.c
+++ b/src/indexamajig.c
@@ -601,6 +601,7 @@ int main(int argc, char *argv[])
 	FelixOptions *felix_opts = NULL;
 	XGandalfOptions *xgandalf_opts = NULL;
 	PinkIndexerOptions *pinkindexer_opts = NULL;
+	double wl_from_dt;
 
 	/* Defaults for "top level" arguments */
 	args.filename = NULL;
@@ -930,6 +931,16 @@ int main(int argc, char *argv[])
 
 	}
 
+	wl_from_dt = data_template_get_wavelength_if_possible(args.iargs.dtempl);
+	if ( !isnan(wl_from_dt) ) {
+		if ( !isnan(args.iargs.wavelength_estimate) ) {
+			ERROR("WARNING: Ignoring your value for "
+			      "--wavelength-estimate because the geometry file "
+			      "already contains a static value.\n");
+		}
+		args.iargs.wavelength_estimate = wl_from_dt;
+	}
+
 	/* Prepare the indexing system */
 	if ( args.indm_str == NULL ) {