41 #include "gravi_data.h" 42 #include "gravi_pfits.h" 43 #include "gravi_dfs.h" 45 #include "gravi_utils.h" 47 #include "gravi_metrology.h" 48 #include "gravi_calib.h" 49 #include "gravi_preproc.h" 50 #include "gravi_wave.h" 51 #include "gravi_p2vm.h" 53 #include "gravi_p2vmred.h" 60 static int gravity_p2vm_create(cpl_plugin *);
61 static int gravity_p2vm_exec(cpl_plugin *);
62 static int gravity_p2vm_destroy(cpl_plugin *);
63 static int gravity_p2vm(cpl_frameset *,
const cpl_parameterlist *);
69 static char gravity_p2vm_short[] =
"Calibrate the instrument bad pixels, wavelength table, interferometric contrast and phase.";
70 static char gravity_p2vm_description[] =
71 "This recipe reduces the internal calibrations. As a special sequence of shutter opening is required, it is advised to always build the SOF with a complete sequence of files obtained within a single execution of the p2vm calibration template. However it is still possible to input a SOF with DARK_RAW only, or DARK_RAW and FLAT_RAW only. It is also possible to input a SOF with some already processed calibration (e.g WAVE).\n" 72 GRAVI_RECIPE_INPUT
"\n" 73 GRAVI_DARK_RAW
" : raw dark, all shutters closed (DPR.TYPE=DARK)\n" 74 GRAVI_FLAT_RAW
" x4 : raw flats, one sutter open (DPR.TYPE=FLAT)\n" 75 GRAVI_P2VM_RAW
" x6 : raw p2vms, two sutters open (DPR.TYPE=P2VM)\n" 76 GRAVI_WAVE_RAW
" : raw wavelength calibration for FT (DPR.TYPE=WAVE)\n" 77 GRAVI_WAVESC_RAW
" : raw wavelength calibration for SC (DPR.TYPE=WAVESC)\n" 78 GRAVI_RECIPE_OUTPUT
"\n" 79 GRAVI_DARK_MAP
" : dark calibration\n" 80 GRAVI_FLAT_MAP
" : flat calibration\n" 81 GRAVI_BAD_MAP
" : badpixel calibration\n" 82 GRAVI_WAVE_MAP
" : wave calibration\n" 83 GRAVI_P2VM_MAP
" : p2vm calibration\n" 101 int cpl_plugin_get_info(cpl_pluginlist * list)
103 cpl_recipe * recipe = cpl_calloc(1,
sizeof *recipe );
104 cpl_plugin * plugin = &recipe->interface;
106 if (cpl_plugin_init(plugin,
108 GRAVI_BINARY_VERSION,
109 CPL_PLUGIN_TYPE_RECIPE,
112 gravity_p2vm_description,
113 "Nabih Azouaoui, Vincent Lapeyrere, JB. Le Bouquin",
118 gravity_p2vm_destroy)) {
119 cpl_msg_error(cpl_func,
"Plugin initialization failed");
120 (void)cpl_error_set_where(cpl_func);
124 if (cpl_pluginlist_append(list, plugin)) {
125 cpl_msg_error(cpl_func,
"Error adding plugin to list");
126 (void)cpl_error_set_where(cpl_func);
142 static int gravity_p2vm_create(cpl_plugin * plugin)
148 if (cpl_error_get_code() != CPL_ERROR_NONE) {
149 cpl_msg_error(cpl_func,
"%s():%d: An error is already set: %s",
150 cpl_func, __LINE__, cpl_error_get_where());
151 return (
int)cpl_error_get_code();
154 if (plugin == NULL) {
155 cpl_msg_error(cpl_func,
"Null plugin");
156 cpl_ensure_code(0, (
int)CPL_ERROR_NULL_INPUT);
160 if (cpl_plugin_get_type(plugin) != CPL_PLUGIN_TYPE_RECIPE) {
161 cpl_msg_error(cpl_func,
"Plugin is not a recipe");
162 cpl_ensure_code(0, (
int)CPL_ERROR_TYPE_MISMATCH);
166 recipe = (cpl_recipe *)plugin;
169 recipe->parameters = cpl_parameterlist_new();
170 if (recipe->parameters == NULL) {
171 cpl_msg_error(cpl_func,
"Parameter list allocation failed");
172 cpl_ensure_code(0, (
int)CPL_ERROR_ILLEGAL_OUTPUT);
178 gravi_parameter_add_static_name (recipe->parameters);
181 gravi_parameter_add_debug_file (recipe->parameters);
182 gravi_parameter_add_preproc_file (recipe->parameters);
185 gravi_parameter_add_badpix (recipe->parameters);
186 gravi_parameter_add_profile (recipe->parameters);
189 gravi_parameter_add_wave (recipe->parameters);
192 p = cpl_parameter_new_enum (
"gravity.calib.phase-calibration", CPL_TYPE_STRING,
193 "The relative phase of the P2VM are defined by \n " 194 "NONE defines phiA(lbd) at zero for all baselines;\n " 195 "CLOSURE defines phiA(lbd) at zero for baselines (01,02,03);\n " 196 "DISP defines phiA(lbd) to have zero mean and minimum GD for baselines (01,02,03);\n " 197 "FULL defines phiA(lbd) to have zero-GD for baselines (01,02,03)",
198 "gravi.p2vm",
"CLOSURE", 4,
"NONE",
"CLOSURE",
"DISP",
"FULL");
199 cpl_parameter_set_alias (p, CPL_PARAMETER_MODE_CLI,
"phase-calibration");
200 cpl_parameter_disable (p, CPL_PARAMETER_MODE_ENV);
201 cpl_parameterlist_append (recipe->parameters, p);
213 static int gravity_p2vm_exec(cpl_plugin * plugin)
218 cpl_errorstate initial_errorstate = cpl_errorstate_get();
221 if (cpl_error_get_code() != CPL_ERROR_NONE) {
222 cpl_msg_error(cpl_func,
"%s():%d: An error is already set: %s",
223 cpl_func, __LINE__, cpl_error_get_where());
224 return (
int)cpl_error_get_code();
227 if (plugin == NULL) {
228 cpl_msg_error(cpl_func,
"Null plugin");
229 cpl_ensure_code(0, (
int)CPL_ERROR_NULL_INPUT);
233 if (cpl_plugin_get_type(plugin) != CPL_PLUGIN_TYPE_RECIPE) {
234 cpl_msg_error(cpl_func,
"Plugin is not a recipe");
235 cpl_ensure_code(0, (
int)CPL_ERROR_TYPE_MISMATCH);
239 recipe = (cpl_recipe *)plugin;
242 if (recipe->parameters == NULL) {
243 cpl_msg_error(cpl_func,
"Recipe invoked with NULL parameter list");
244 cpl_ensure_code(0, (
int)CPL_ERROR_NULL_INPUT);
246 if (recipe->frames == NULL) {
247 cpl_msg_error(cpl_func,
"Recipe invoked with NULL frame set");
248 cpl_ensure_code(0, (
int)CPL_ERROR_NULL_INPUT);
252 recipe_status = gravity_p2vm(recipe->frames, recipe->parameters);
255 if (cpl_dfs_update_product_header(recipe->frames)) {
256 if (!recipe_status) recipe_status = (int)cpl_error_get_code();
259 if (!cpl_errorstate_is_equal(initial_errorstate)) {
262 cpl_errorstate_dump(initial_errorstate, CPL_FALSE, NULL);
265 return recipe_status;
275 static int gravity_p2vm_destroy(cpl_plugin * plugin)
279 if (plugin == NULL) {
280 cpl_msg_error(cpl_func,
"Null plugin");
281 cpl_ensure_code(0, (
int)CPL_ERROR_NULL_INPUT);
285 if (cpl_plugin_get_type(plugin) != CPL_PLUGIN_TYPE_RECIPE) {
286 cpl_msg_error(cpl_func,
"Plugin is not a recipe");
287 cpl_ensure_code(0, (
int)CPL_ERROR_TYPE_MISMATCH);
291 recipe = (cpl_recipe *)plugin;
293 cpl_parameterlist_delete(recipe->parameters);
307 static int gravity_p2vm(cpl_frameset * frameset,
308 const cpl_parameterlist * parlist)
310 cpl_frameset * p2vm_frameset=NULL, * wavecalib_frameset=NULL, * darkcalib_frameset=NULL,
311 * flatcalib_frameset=NULL, * dark_frameset=NULL, * wave_frameset=NULL, * wavesc_frameset=NULL,
312 * badcalib_frameset=NULL, * flat_frameset=NULL, * used_frameset=NULL, * current_frameset=NULL;
314 cpl_frame * frame=NULL, * frame_p2vm=NULL;
316 gravi_data * p2vm_map=NULL, * data=NULL, * dark_map=NULL, * wave_map=NULL,
317 * profile_map=NULL, * badpix_map=NULL, * wave_data=NULL;
318 gravi_data * spectrum_data=NULL;
319 gravi_data * preproc_data=NULL;
320 gravi_data ** raw_data=NULL;
321 gravi_data * p2vmred_data = NULL;
323 int nb_frame, i, nb_frame_gain = 0;
325 cpl_propertylist * met_plist=NULL;
326 int ** valid_trans = cpl_malloc (2 *
sizeof (
int*));
327 int ** valid_CP = cpl_malloc (2 *
sizeof (
int*));
330 for (i = 0 ; i < 2; i++){
331 valid_trans[i] = cpl_calloc (4,
sizeof (
int));
332 valid_CP[i] = cpl_calloc (6,
sizeof (
int));
336 gravity_print_banner ();
337 cpl_msg_set_time_on();
338 cpl_msg_set_component_on();
339 gravi_msg_function_start(1);
343 cpl_ensure_code(gravi_dfs_set_groups(frameset) == CPL_ERROR_NONE, cpl_error_get_code()) ;
346 used_frameset = cpl_frameset_new ();
349 dark_frameset = gravi_frameset_extract_dark_data (frameset);
350 darkcalib_frameset = gravi_frameset_extract_dark_map (frameset);
353 flat_frameset = gravi_frameset_extract_flat_data (frameset);
354 flatcalib_frameset = gravi_frameset_extract_flat_map (frameset);
357 badcalib_frameset = gravi_frameset_extract_bad_map (frameset);
360 wave_frameset = gravi_frameset_extract_wave_data (frameset);
361 wavesc_frameset = gravi_frameset_extract_wavesc_data (frameset);
362 wavecalib_frameset = gravi_frameset_extract_wave_map (frameset);
365 p2vm_frameset = gravi_frameset_extract_p2vm_data (frameset);
373 if (!cpl_frameset_is_empty (dark_frameset)) {
374 cpl_msg_info (cpl_func,
" ***** Compute DARK map ***** ");
377 frame = cpl_frameset_get_position (dark_frameset, 0);
378 data = gravi_data_load_rawframe (frame, used_frameset);
379 gravi_data_detector_cleanup (data, parlist);
382 dark_map = gravi_compute_dark (data);
383 FREE (gravi_data_delete, data);
385 CPLCHECK_CLEAN (
"Cannot compute the DARK map");
388 gravi_data_save_new (dark_map, frameset, NULL, parlist,
389 NULL, frame,
"gravity_p2vm",
390 NULL, GRAVI_DARK_MAP);
392 CPLCHECK_CLEAN (
"Could not save the DARK map");
394 else if (!cpl_frameset_is_empty (darkcalib_frameset)) {
395 cpl_msg_info (cpl_func,
" ***** Get DARK map ***** ");
398 frame = cpl_frameset_get_position (darkcalib_frameset, 0);
399 dark_map = gravi_data_load_frame (frame, used_frameset);
401 CPLCHECK_CLEAN (
"Could not load the DARK map");
404 cpl_frameset_insert (dark_frameset, cpl_frame_duplicate (frame));
412 if (!cpl_frameset_is_empty (badcalib_frameset)) {
413 cpl_msg_info (cpl_func,
" ***** Get BAD pixel map ***** ");
416 frame = cpl_frameset_get_position (badcalib_frameset, 0);
417 badpix_map = gravi_data_load_frame (frame, used_frameset);
419 else if (!cpl_frameset_is_empty (dark_frameset) &&
420 !cpl_frameset_is_empty (flat_frameset)) {
421 cpl_msg_info (cpl_func,
" ***** Compute BAD pixel map from DARK and FLAT_RAW ***** ");
424 nb_frame_gain = cpl_frameset_get_size (flat_frameset);
425 raw_data = cpl_calloc (nb_frame_gain,
sizeof(gravi_data *));
428 for (
int i = 0; i < nb_frame_gain; i++) {
429 frame = cpl_frameset_get_position (flat_frameset, i);
430 raw_data[i] = gravi_data_load_rawframe (frame, NULL);
431 gravi_data_detector_cleanup (raw_data[i], parlist);
435 badpix_map = gravi_compute_badpix (dark_map, raw_data,
436 nb_frame_gain, parlist);
438 CPLCHECK_CLEAN(
"Cannot compute the BAD pixel from DARK and FLAT");
441 frame = cpl_frameset_get_position (dark_frameset, 0);
442 gravi_data_save_new (badpix_map, frameset, NULL, parlist,
443 NULL, frame,
"gravity_p2vm",
444 NULL, GRAVI_BAD_MAP);
446 CPLCHECK_CLEAN (
"Could not save the BAD pixel map");
448 FREELOOP (gravi_data_delete, raw_data, nb_frame_gain);
455 if (!cpl_frameset_is_empty(flatcalib_frameset)) {
456 cpl_msg_info (cpl_func,
" ***** Get FLAT map ***** ");
459 frame = cpl_frameset_get_position (flatcalib_frameset, 0);
460 profile_map = gravi_data_load_frame (frame, used_frameset);
462 else if (!cpl_frameset_is_empty (flat_frameset)) {
463 cpl_msg_info (cpl_func,
" ***** Compute FLAT map ***** ");
467 if (!badpix_map || !dark_map) {
468 ERROR_CLEAN (CPL_ERROR_ILLEGAL_INPUT,
"Missing DARK or BAD in frameset");
472 nb_frame_gain = cpl_frameset_get_size (flat_frameset);
473 raw_data = cpl_malloc (nb_frame_gain *
sizeof(gravi_data *));
476 for (i = 0; i < nb_frame_gain; i++) {
477 frame = cpl_frameset_get_position (flat_frameset, i);
478 raw_data[i] = gravi_data_load_rawframe (frame, used_frameset);
479 gravi_data_detector_cleanup (raw_data[i], parlist);
484 profile_map = gravi_compute_profile (raw_data, dark_map, badpix_map, nb_frame_gain, parlist);
485 CPLCHECK_CLEAN (
"Cannot compute the FLAT profile");
488 cpl_propertylist * gain_header;
489 gain_header = gravi_compute_gain (raw_data, nb_frame_gain, dark_map);
490 CPLCHECK_CLEAN (
"Cannot compute the GAIN");
493 cpl_propertylist * profile_header = gravi_data_get_header (profile_map);
494 cpl_propertylist_append (profile_header, gain_header);
495 FREE (cpl_propertylist_delete, gain_header);
498 frame = cpl_frameset_get_position (flat_frameset, 0);
500 gravi_data_save_new (profile_map, frameset, NULL, parlist,
501 used_frameset, frame,
"gravity_p2vm",
502 NULL, GRAVI_FLAT_MAP);
504 CPLCHECK_CLEAN (
"Could not save the FLAT profile_map");
507 FREELOOP (gravi_data_delete, raw_data, nb_frame_gain);
514 if (!cpl_frameset_is_empty (wavecalib_frameset)) {
515 cpl_msg_info (cpl_func,
" ***** Get wave map ***** ");
518 frame = cpl_frameset_get_position (wavecalib_frameset, 0);
519 wave_map = gravi_data_load_frame (frame, used_frameset);
521 else if (!cpl_frameset_is_empty (wave_frameset)) {
522 cpl_msg_info (cpl_func,
" ***** Process the WAVE_RAW ***** ");
525 if ( !badpix_map || !profile_map || !dark_map) {
526 ERROR_CLEAN (CPL_ERROR_ILLEGAL_INPUT,
"Missing DARK, FLAT, or BAD in frameset");
530 wave_map = gravi_data_new (0);
533 frame = cpl_frameset_get_position (wave_frameset, 0);
534 data = gravi_data_load_rawframe (frame, used_frameset);
535 gravi_data_detector_cleanup (data, parlist);
538 cpl_msg_info (cpl_func,
"Compute the P2VM_MET from WAVE_RAW");
539 cpl_table * met_table = gravi_data_get_table (data, GRAVI_METROLOGY_EXT);
540 cpl_table * p2vm_met = gravi_metrology_compute_p2vm (met_table);
543 gravi_data_add_table (wave_map, NULL, GRAVI_P2VM_MET_EXT, p2vm_met);
545 CPLCHECK_CLEAN (
"Cannot compute P2VM_MET");
548 cpl_msg_info (cpl_func,
"Extract SPECTRUM for WAVE_RAW");
549 spectrum_data = gravi_extract_spectrum (data, profile_map, dark_map,
550 badpix_map, NULL, parlist);
552 cpl_msg_info (cpl_func,
"Compute OPDs for WAVE_RAW");
553 gravi_wave_compute_opds (spectrum_data, gravi_data_get_table (data, GRAVI_METROLOGY_EXT));
554 FREE (gravi_data_delete, data);
556 CPLCHECK_CLEAN (
"Cannot process the WAVE_RAW");
559 if (gravi_param_get_bool (parlist,
"gravity.dfs.debug-file")) {
560 gravi_data_save_new (spectrum_data, frameset, NULL, parlist,
561 used_frameset, frame,
"gravity_p2vm",
566 gravi_compute_wave (wave_map, spectrum_data, GRAVI_FT, parlist);
568 CPLCHECK_CLEAN (
"Cannot compute wave for FT");
571 if (!cpl_frameset_is_empty (wavesc_frameset)) {
573 cpl_msg_info (cpl_func,
"Delete WAVE_RAW data");
574 FREE (gravi_data_delete, spectrum_data);
576 cpl_msg_info (cpl_func,
" ***** Process the WAVESC_RAW ***** ");
579 frame = cpl_frameset_get_position(wavesc_frameset, 0);
580 data = gravi_data_load_rawframe (frame, used_frameset);
581 gravi_data_detector_cleanup (data, parlist);
583 cpl_msg_info (cpl_func,
"Extract SPECTRUM for WAVESC_RAW");
584 spectrum_data = gravi_extract_spectrum (data, profile_map, dark_map,
585 badpix_map, NULL, parlist);
587 cpl_msg_info (cpl_func,
"Compute OPDs for WAVESC_RAW");
588 gravi_wave_compute_opds (spectrum_data, gravi_data_get_table (data, GRAVI_METROLOGY_EXT));
589 FREE (gravi_data_delete, data);
591 CPLCHECK_CLEAN (
"Cannot process the WAVESC_RAW");
594 if (gravi_param_get_bool (parlist,
"gravity.dfs.debug-file")) {
595 gravi_data_save_new (spectrum_data, frameset, NULL, parlist,
596 used_frameset, frame,
"gravity_p2vm",
601 cpl_msg_warning (cpl_func,
"No WAVESC_RAW in the SOF:" 602 " SC wavelength will be inaccurate");
606 gravi_compute_wave (wave_map, spectrum_data, GRAVI_SC, parlist);
608 CPLCHECK_CLEAN (
"Cannot compute wave for SC");
611 FREE (gravi_data_delete, spectrum_data);
614 gravi_wave_qc (wave_map, profile_map);
617 gravi_data_save_new (wave_map, frameset, NULL, parlist,
618 used_frameset, frame,
"gravity_p2vm",
619 NULL, GRAVI_WAVE_MAP);
621 CPLCHECK_CLEAN (
"Could not save the WAVE map");
629 if ( cpl_frameset_is_empty (p2vm_frameset) ) {
630 cpl_msg_info (cpl_func,
"All RAW data reduced... stop recipe");
635 if ( !badpix_map || !profile_map || !wave_map || !dark_map) {
636 ERROR_CLEAN (CPL_ERROR_ILLEGAL_INPUT,
"Missing DARK, FLAT, BAD, or WAVE in frameset");
640 if ( !cpl_frameset_is_empty (flat_frameset) )
641 cpl_frameset_join (p2vm_frameset, flat_frameset);
643 if ( !cpl_frameset_is_empty (wave_frameset) )
644 cpl_frameset_join (p2vm_frameset, wave_frameset);
647 nb_frame = cpl_frameset_get_size (p2vm_frameset);
650 if (nb_frame != 11) {
651 ERROR_CLEAN (CPL_ERROR_ILLEGAL_INPUT,
"Missing P2VM in frameset");
659 cpl_msg_info (cpl_func,
" ***** Create the P2VM ***** ");
660 p2vm_map = gravi_create_p2vm (wave_map);
661 CPLCHECK_CLEAN (
"Cannot create the P2VM data");
664 for (i = 0; i < nb_frame; i++) {
666 cpl_msg_info (cpl_func,
" ***** file %d over %d ***** ", i+1, nb_frame );
667 current_frameset = cpl_frameset_duplicate (used_frameset);
670 frame = cpl_frameset_get_position (p2vm_frameset, i);
671 cpl_frameset_insert (used_frameset, cpl_frame_duplicate (frame));
672 data = gravi_data_load_rawframe (frame, current_frameset);
673 gravi_data_detector_cleanup (data, parlist);
674 CPLCHECK_CLEAN (
"Cannot load data");
677 cpl_msg_info (cpl_func,
"Shutters: %d-%d-%d-%d",
678 gravi_data_get_shutter (data, 0), gravi_data_get_shutter (data, 1),
679 gravi_data_get_shutter (data, 2), gravi_data_get_shutter (data, 3));
684 if ( frame_p2vm==NULL && gravi_data_check_shutter (data, 1,1,0,0) )
686 cpl_msg_info (cpl_func,
"Use this frame for the P2VM product");
693 if ( gravi_data_check_shutter (data, 1,1,1,1) ) {
696 cpl_msg_info (cpl_func,
"Nothing to be done with the file... yet.");
697 wave_data = data; data = NULL;
699 FREE (cpl_frameset_delete, current_frameset);
705 preproc_data = gravi_extract_spectrum (data, profile_map, dark_map,
706 badpix_map, NULL, parlist);
707 FREE (gravi_data_delete, data);
708 CPLCHECK_CLEAN (
"Cannot extract spectrum");
711 gravi_align_spectrum (preproc_data, wave_map, p2vm_map);
712 CPLCHECK_CLEAN (
"Cannot re-interpolate spectrum");
715 gravi_compute_p2vm (p2vm_map, preproc_data, valid_trans, valid_CP);
716 CPLCHECK_CLEAN(
"Cannot compute the P2VM");
719 if (gravi_param_get_bool (parlist,
"gravity.dfs.preproc-file")) {
721 gravi_data_save_new (preproc_data, frameset, NULL, parlist,
722 current_frameset, frame,
"gravity_p2vm",
723 NULL, GRAVI_PREPROC);
728 FREE (gravi_data_delete,preproc_data);
729 cpl_msg_info(cpl_func,
"Execution time to delete preproc_data : %f s",
730 (clock() - start) / (
double)CLOCKS_PER_SEC);
733 FREE (cpl_frameset_delete, current_frameset);
741 gravi_p2vm_normalisation (p2vm_map, valid_trans, valid_CP);
742 CPLCHECK_CLEAN(
"Cannot normalise the p2vm_map");
750 cpl_msg_info (cpl_func,
" ***** Analyse the WAVE file to calibrate the internal closure and transmission ***** ");
753 preproc_data = gravi_extract_spectrum (wave_data, profile_map, dark_map,
754 badpix_map, NULL, parlist);
755 CPLCHECK_CLEAN (
"Cannot extract spectrum");
758 gravi_align_spectrum (preproc_data, wave_map, p2vm_map);
759 CPLCHECK_CLEAN (
"Cannot re-interpolate spectrum");
762 gravi_data_move_ext (preproc_data, wave_data, GRAVI_ARRAY_GEOMETRY_EXT);
763 gravi_data_move_ext (preproc_data, wave_data, GRAVI_OPTICAL_TRAIN_EXT);
764 FREE (gravi_data_delete, wave_data);
765 CPLCHECK_CLEAN (
"Cannot move ext");
768 p2vmred_data = gravi_compute_p2vmred (preproc_data, p2vm_map,
"gravi_single", parlist);
769 FREE (gravi_data_delete, preproc_data);
770 CPLCHECK_CLEAN (
"Cannot apply p2vm");
773 if (!strcmp (gravi_param_get_string (parlist,
"gravity.calib.phase-calibration"),
"CLOSURE")) {
774 gravi_p2vm_phase_correction (p2vm_map, p2vmred_data, 0);
776 else if (!strcmp (gravi_param_get_string (parlist,
"gravity.calib.phase-calibration"),
"DISP")) {
777 gravi_p2vm_phase_correction (p2vm_map, p2vmred_data, 1);
779 else if (!strcmp (gravi_param_get_string (parlist,
"gravity.calib.phase-calibration"),
"FULL")) {
780 gravi_p2vm_phase_correction (p2vm_map, p2vmred_data, 2);
783 cpl_msg_info (cpl_func,
"P2VM phases are kept to zero (option phase-calibration=NONE)");
785 CPLCHECK_CLEAN (
"Cannot recalibrate the P2VM phases");
788 gravi_p2vm_transmission (p2vm_map, p2vmred_data);
789 FREE (gravi_data_delete, p2vmred_data);
791 CPLCHECK_CLEAN (
"Cannot compute the transmission");
799 gravi_data_save_new (p2vm_map, frameset, NULL, parlist,
800 used_frameset, frame_p2vm,
"gravity_p2vm",
801 NULL, GRAVI_P2VM_MAP);
803 CPLCHECK_CLEAN(
"Could not save the P2VM on the output file");
807 cpl_msg_info (cpl_func,
"Cleanup memory");
809 FREE (cpl_frameset_delete, dark_frameset);
810 FREE (cpl_frameset_delete, darkcalib_frameset);
811 FREE (cpl_frameset_delete, badcalib_frameset);
812 FREE (gravi_data_delete, spectrum_data);
813 FREE (gravi_data_delete, p2vmred_data);
814 FREE (gravi_data_delete, wave_data);
815 FREE (gravi_data_delete, dark_map);
816 FREELOOP (gravi_data_delete, raw_data, nb_frame_gain);
817 FREE (gravi_data_delete, badpix_map);
818 FREE (cpl_propertylist_delete, met_plist);
819 FREE (cpl_frameset_delete, wavecalib_frameset);
820 FREE (cpl_frameset_delete, flatcalib_frameset);
821 FREE (cpl_frameset_delete, flat_frameset);
822 FREE (cpl_frameset_delete, wave_frameset);
823 FREE (cpl_frameset_delete, wavesc_frameset);
824 FREE (cpl_frameset_delete, used_frameset);
825 FREE (gravi_data_delete, data);
826 FREE (gravi_data_delete, preproc_data);
827 FREE (gravi_data_delete, dark_map);
828 FREE (gravi_data_delete, wave_map);
829 FREE (gravi_data_delete, profile_map);
830 FREELOOP (cpl_free, valid_CP, 2);
831 FREELOOP (cpl_free, valid_trans, 2);
832 FREE (gravi_data_delete, p2vm_map);
833 FREE (cpl_frameset_delete, p2vm_frameset);
834 FREE (cpl_frameset_delete, current_frameset);
837 CPLCHECK_INT (
"Could not cleanup memory");
839 gravi_msg_function_exit(1);
840 return (
int)cpl_error_get_code();