From 008a336f98e975451472a040e0688bdef8839532 Mon Sep 17 00:00:00 2001
From: Roberto Susino <roberto.susino@inaf.it>
Date: Wed, 26 Jan 2022 14:07:17 +0100
Subject: [PATCH] Add the keyword POLANGLE to VL-image FITS headers

---
 metis_l2_prep_vl_generic.pro | 26 ++++++++++++++++++++++++++
 1 file changed, 26 insertions(+)

diff --git a/metis_l2_prep_vl_generic.pro b/metis_l2_prep_vl_generic.pro
index 932a131..ceb69c9 100644
--- a/metis_l2_prep_vl_generic.pro
+++ b/metis_l2_prep_vl_generic.pro
@@ -158,6 +158,32 @@ pro metis_l2_prep_vl_generic
 
 	sxdelpar, primary_header, 'BLANK'
 	
+	; read the calibration curve to convert pmp raw voltages (dacpol) into effective polarization angles
+
+	dacpol_cal = cal_pack.vl_channel.dacpol_cal
+
+	if header.pol_id ne 5 then begin
+		if fix(header.hdr_vers) le 4 then begin
+			case header.pol_id of
+				0: dacpol = header.dac1pol1
+				1: dacpol = header.dac1pol1
+				2: dacpol = header.dac1pol2
+				3: dacpol = header.dac1pol3
+				4: dacpol = header.dac1pol4
+			endcase
+		endif
+
+		if fix(header.hdr_vers) ge 5 then begin
+			dacpol = header.dac1pol1
+		endif
+
+		k = where(dacpol_cal.dacpol eq dacpol)
+		angle = dacpol_cal.angle[k]
+		angle = float(angle[0])
+
+		fxaddpar, primary_header, 'POLANGLE', angle, '[deg] Polarization angle', after = 'POL_ID'
+	endif
+
 	; append wcs keywords
 
 	wcs = metis_wcs(header, cal_pack, ref_detector = ref_detector)
-- 
GitLab