From 7a4420be5a080207b0c3897bab7a1224b2b6adab Mon Sep 17 00:00:00 2001
From: Roberto Susino <roberto.susino@inaf.it>
Date: Mon, 5 Sep 2022 09:42:21 +0200
Subject: [PATCH] Add patch to correct for missing SEQ_NUM keyword

---
 metis_l2_prep_vl_generic.pro | 15 ++++++++++++++-
 metis_l2_prep_vl_polariz.pro | 16 ++++++++++++++++
 2 files changed, 30 insertions(+), 1 deletion(-)

diff --git a/metis_l2_prep_vl_generic.pro b/metis_l2_prep_vl_generic.pro
index cde56d3..fca9edb 100644
--- a/metis_l2_prep_vl_generic.pro
+++ b/metis_l2_prep_vl_generic.pro
@@ -34,7 +34,20 @@ pro metis_l2_prep_vl_generic
 	; read the primary hdu
 
 	data = mrdfits(input.file_name, 0, primary_header, /silent)
-
+	
+	; patch to fix the lack of the keyword SEQ_NUM
+	
+	pol_id = fxpar(primary_header, 'POL_ID', missing = 0)
+	if pol_id ge 1 and pol_id le 4 then begin
+		seq_num = fxpar(primary_header, 'SEQ_NUM', missing = 0)
+		if seq_num eq 0 then begin
+			obj_cnt = fxpar(primary_header, 'OBJ_CNT')
+			n_pol = fxpar(primary_header, 'N_POL', missing = 4)
+			seq_num = ((obj_cnt - 1)/n_pol + 1)
+			fxaddpar, primary_header, 'SEQ_NUM', seq_num, before = 'POL_ID'
+		endif
+	endif
+	
 	; read the quality matrix
 
 	quality_matrix = mrdfits(input.file_name, 'quality matrix', /silent)
diff --git a/metis_l2_prep_vl_polariz.pro b/metis_l2_prep_vl_polariz.pro
index 03d203b..9b6c791 100755
--- a/metis_l2_prep_vl_polariz.pro
+++ b/metis_l2_prep_vl_polariz.pro
@@ -55,6 +55,22 @@ pro metis_l2_prep_vl_polariz
 
 		image = mrdfits(input.file_name[k], 0, primary_header, /silent)
 
+		; patch to fix the lack of the keyword SEQ_NUM
+		
+		seq_num = fxpar(primary_header, 'SEQ_NUM', missing = 0)
+		if seq_num eq 0 then begin
+			obj_cnt = fxpar(primary_header, 'OBJ_CNT')
+			n_pol = fxpar(primary_header, 'N_POL', missing = 4)
+			seq_num = ((obj_cnt - 1)/n_pol + 1)
+			fxaddpar, primary_header, 'SEQ_NUM', seq_num, before = 'POL_ID'
+		endif
+
+		if k eq 0 then comment = fxpar(primary_header, 'COMMENT') else begin
+			sxdelpar, primary_header, 'COMMENT'
+			foreach line, comment do $
+				fxaddpar, primary_header, 'COMMENT', line
+		endelse
+
 		header = fits_hdr2struct(primary_header)
 
 		; ====================================
-- 
GitLab