From 9fbb4f7f3ce82fcb7743fb4e821d0a9692cac1f1 Mon Sep 17 00:00:00 2001 From: Andrew Aquila Date: Fri, 2 Dec 2011 18:56:23 +0100 Subject: add strem_grep and old powder_plot --- src/powder_plot.c | 18 ++++++++++++++---- 1 file changed, 14 insertions(+), 4 deletions(-) (limited to 'src') diff --git a/src/powder_plot.c b/src/powder_plot.c index 1769c949..4cc2d1f8 100644 --- a/src/powder_plot.c +++ b/src/powder_plot.c @@ -887,7 +887,7 @@ int main(int argc, char *argv[]) if ( is_stream(filename) == 1 ) { file_type = FILE_STREAM; - + } else if ( H5Fis_hdf5(filename) > 0 ) { file_type = FILE_H5; @@ -949,6 +949,15 @@ int main(int argc, char *argv[]) if ( file_type == FILE_HKL ) { need_geometry = 0; need_beam = 0; + + /* Logic checks */ + if ( need_geometry && (image.lambda < 0.0) ) { + need_beam = 1; + } + if ( hist_info.histsize <= 0 ) { + ERROR("You need to specify a histogram with more then 0 " + "bins\n"); + return 1; } /* Get geometry, beam and pdb files and parameters as needed */ @@ -970,7 +979,7 @@ int main(int argc, char *argv[]) } free(geometry); - /* Open files to get wavelength if it exists & camera length + /* Open files to get wavelength if it exists & camera length if they are not found in the geometry file */ if (file_type == FILE_STREAM) { fh = fopen(filename, "r"); @@ -1067,11 +1076,11 @@ int main(int argc, char *argv[]) resolution_limits(image.reflections, cell, &hist_info.q_min, &hist_info.q_max); } else if (hist_info.q_min < 0.0) { - double dummy; + double dummy; resolution_limits(image.reflections, cell, &hist_info.q_min, &dummy); } else if (hist_info.q_max < 0.0) { - double dummy; + double dummy; resolution_limits(image.reflections, cell, &dummy, &hist_info.q_max); } @@ -1090,6 +1099,7 @@ int main(int argc, char *argv[]) hist_info.q_min, hist_info.q_max); return 1; } + if ( hist_info.spacing == LINEAR) { hist_info.q_delta = (hist_info.q_max - hist_info.q_min)/ hist_info.histsize; -- cgit v1.2.3