Skip to content
Snippets Groups Projects
Commit 484b9f2a authored by Roberto Susino's avatar Roberto Susino
Browse files

Prettify comments in FITS header

parent d3acf981
No related branches found
No related tags found
No related merge requests found
...@@ -169,6 +169,11 @@ pro metis_l2_prep_uv ...@@ -169,6 +169,11 @@ pro metis_l2_prep_uv
fxaddpar, primary_header, 'COMMENT', 'WARNING: UV radiometric calibration is still preliminary.' fxaddpar, primary_header, 'COMMENT', 'WARNING: UV radiometric calibration is still preliminary.'
fxaddpar, primary_header, 'COMMENT', 'Uncertainty matrix in the FITS extension is preliminary.' fxaddpar, primary_header, 'COMMENT', 'Uncertainty matrix in the FITS extension is preliminary.'
if ref_detector then begin
fxaddpar, primary_header, 'COMMENT', 'Flip vertically and rotate CROTA degrees counter-clockwise'
fxaddpar, primary_header, 'COMMENT', ' to have Solar North up.'
endif else $
fxaddpar, primary_header, 'COMMENT', 'Rotate CROTA degrees counter-clockwise to have Solar North up.' fxaddpar, primary_header, 'COMMENT', 'Rotate CROTA degrees counter-clockwise to have Solar North up.'
for k = 0, n_elements(history) - 1 do $ for k = 0, n_elements(history) - 1 do $
...@@ -185,9 +190,9 @@ pro metis_l2_prep_uv ...@@ -185,9 +190,9 @@ pro metis_l2_prep_uv
; add the extension with the quality matrix ; add the extension with the quality matrix
extension_header = !null extension_header = !null
fxaddpar, extension_header, 'PCOUNT', 0, 'Parameter count' fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'Group count' fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'Extension name' fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'extension name'
if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'UV') if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'UV')
fits_add_checksum, extension_header, quality_matrix fits_add_checksum, extension_header, quality_matrix
mwrfits, float(quality_matrix), out_file_name, extension_header, /no_comment, /silent mwrfits, float(quality_matrix), out_file_name, extension_header, /no_comment, /silent
...@@ -197,9 +202,9 @@ pro metis_l2_prep_uv ...@@ -197,9 +202,9 @@ pro metis_l2_prep_uv
; add the extension with the error matrix ; add the extension with the error matrix
extension_header = !null extension_header = !null
fxaddpar, extension_header, 'PCOUNT', 0, 'Parameter count' fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'Group count' fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'Extension name' fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'extension name'
error_matrix = data * sqrt(error) error_matrix = data * sqrt(error)
if not ref_detector then error_matrix = metis_rectify(error_matrix, 'UV') if not ref_detector then error_matrix = metis_rectify(error_matrix, 'UV')
fits_add_checksum, extension_header, float(error_matrix) fits_add_checksum, extension_header, float(error_matrix)
......
...@@ -182,7 +182,7 @@ pro metis_l2_prep_vl_generic ...@@ -182,7 +182,7 @@ pro metis_l2_prep_vl_generic
angle = dacpol_cal.angle[k] angle = dacpol_cal.angle[k]
angle = float(angle[0]) angle = float(angle[0])
fxaddpar, primary_header, 'POLANGLE', angle, '[deg] Polarization angle', after = 'POL_ID' fxaddpar, primary_header, 'POLANGLE', angle, '[deg] polarization angle', after = 'POL_ID'
endif endif
; append wcs keywords ; append wcs keywords
...@@ -218,9 +218,9 @@ pro metis_l2_prep_vl_generic ...@@ -218,9 +218,9 @@ pro metis_l2_prep_vl_generic
; add the extension with the quality matrix ; add the extension with the quality matrix
extension_header = !null extension_header = !null
fxaddpar, extension_header, 'PCOUNT', 0, 'Parameter count' fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'Group count' fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'Extension name' fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'extension name'
if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL') if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL')
fits_add_checksum, extension_header, quality_matrix fits_add_checksum, extension_header, quality_matrix
mwrfits, float(quality_matrix), out_file_name, extension_header, /no_comment, /silent mwrfits, float(quality_matrix), out_file_name, extension_header, /no_comment, /silent
...@@ -230,9 +230,9 @@ pro metis_l2_prep_vl_generic ...@@ -230,9 +230,9 @@ pro metis_l2_prep_vl_generic
; add the extension with the error matrix ; add the extension with the error matrix
extension_header = !null extension_header = !null
fxaddpar, extension_header, 'PCOUNT', 0, 'Parameter count' fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'Group count' fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'Extension name' fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'extension name'
error_matrix = data * sqrt(error) error_matrix = data * sqrt(error)
if not ref_detector then data = metis_rectify(error_matrix, 'VL') if not ref_detector then data = metis_rectify(error_matrix, 'VL')
fits_add_checksum, extension_header, float(error_matrix) fits_add_checksum, extension_header, float(error_matrix)
......
...@@ -294,14 +294,14 @@ pro metis_l2_prep_vl_polariz ...@@ -294,14 +294,14 @@ pro metis_l2_prep_vl_polariz
fxaddpar, primary_header, 'VERS_SW', input.sw_version fxaddpar, primary_header, 'VERS_SW', input.sw_version
fxaddpar, primary_header, 'VERS_CAL', cal_pack.version fxaddpar, primary_header, 'VERS_CAL', cal_pack.version
fxaddpar, primary_header, 'VERSION', version fxaddpar, primary_header, 'VERSION', version
fxaddpar, primary_header, 'DATE', date, 'Date and time of FITS file creation' fxaddpar, primary_header, 'DATE', date, 'date and time of FITS file creation'
fxaddpar, primary_header, 'DATE-BEG', header.date_beg, 'Start time of observation' fxaddpar, primary_header, 'DATE-BEG', header.date_beg, 'start time of observation'
fxaddpar, primary_header, 'DATE-OBS', header.date_obs, 'Same as DATE-BEG' fxaddpar, primary_header, 'DATE-OBS', header.date_obs, 'same as DATE-BEG'
fxaddpar, primary_header, 'DATE-AVG', header.date_avg, 'Average time of observation' fxaddpar, primary_header, 'DATE-AVG', header.date_avg, 'average time of observation'
fxaddpar, primary_header, 'DATE-END', header.date_end, 'End time of observation' fxaddpar, primary_header, 'DATE-END', header.date_end, 'end time of observation'
fxaddpar, primary_header, 'OBT_BEG', header.obt_beg, 'Start acquisition time in on-board time', format = 'F0.5' fxaddpar, primary_header, 'OBT_BEG', header.obt_beg, 'start acquisition time in on-board time', format = 'F0.5'
fxaddpar, primary_header, 'OBT_END', header.obt_end, 'End acquisition time in on-board time', format = 'F0.5' fxaddpar, primary_header, 'OBT_END', header.obt_end, 'end acquisition time in on-board time', format = 'F0.5'
fxaddpar, primary_header, 'TELAPSE', telapse, '[s] Elapsed time between beginning and end of observation' fxaddpar, primary_header, 'TELAPSE', telapse, '[s] elapsed time between beginning and end of observation'
fxaddpar, primary_header, 'WAVEBAND', cal_pack.vl_channel.name fxaddpar, primary_header, 'WAVEBAND', cal_pack.vl_channel.name
fxaddpar, primary_header, 'XPOSURE', header.xposure fxaddpar, primary_header, 'XPOSURE', header.xposure
fxaddpar, primary_header, 'NSUMEXP', header.nsumexp fxaddpar, primary_header, 'NSUMEXP', header.nsumexp
...@@ -388,9 +388,9 @@ pro metis_l2_prep_vl_polariz ...@@ -388,9 +388,9 @@ pro metis_l2_prep_vl_polariz
; add the extension with the quality matrix ; add the extension with the quality matrix
extension_header = !null extension_header = !null
fxaddpar, extension_header, 'PCOUNT', 0, 'Parameter count' fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'Group count' fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'Extension name' fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'extension name'
if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL') if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL')
fits_add_checksum, extension_header, quality_matrix fits_add_checksum, extension_header, quality_matrix
mwrfits, float(quality_matrix), out_file_name[0], extension_header, /no_comment, /silent mwrfits, float(quality_matrix), out_file_name[0], extension_header, /no_comment, /silent
...@@ -400,9 +400,9 @@ pro metis_l2_prep_vl_polariz ...@@ -400,9 +400,9 @@ pro metis_l2_prep_vl_polariz
; add the extension with the error matrix ; add the extension with the error matrix
extension_header = !null extension_header = !null
fxaddpar, extension_header, 'PCOUNT', 0, 'Parameter count' fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'Group count' fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'Extension name' fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'extension name'
error_matrix = intarr(header.naxis1, header.naxis2) error_matrix = intarr(header.naxis1, header.naxis2)
if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL') if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL')
fits_add_checksum, extension_header, error_matrix fits_add_checksum, extension_header, error_matrix
...@@ -442,9 +442,9 @@ pro metis_l2_prep_vl_polariz ...@@ -442,9 +442,9 @@ pro metis_l2_prep_vl_polariz
; add the extension with the quality matrix ; add the extension with the quality matrix
extension_header = !null extension_header = !null
fxaddpar, extension_header, 'PCOUNT', 0, 'Parameter count' fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'Group count' fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'Extension name' fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'extension name'
if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL') if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL')
fits_add_checksum, extension_header, quality_matrix fits_add_checksum, extension_header, quality_matrix
mwrfits, float(quality_matrix), out_file_name[1], extension_header, /no_comment, /silent mwrfits, float(quality_matrix), out_file_name[1], extension_header, /no_comment, /silent
...@@ -454,9 +454,9 @@ pro metis_l2_prep_vl_polariz ...@@ -454,9 +454,9 @@ pro metis_l2_prep_vl_polariz
; add the extension with the error matrix ; add the extension with the error matrix
extension_header = !null extension_header = !null
fxaddpar, extension_header, 'PCOUNT', 0, 'Parameter count' fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'Group count' fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'Extension name' fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'extension name'
error_matrix = intarr(header.naxis1, header.naxis2) error_matrix = intarr(header.naxis1, header.naxis2)
if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL') if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL')
fits_add_checksum, extension_header, error_matrix fits_add_checksum, extension_header, error_matrix
...@@ -496,9 +496,9 @@ pro metis_l2_prep_vl_polariz ...@@ -496,9 +496,9 @@ pro metis_l2_prep_vl_polariz
; add the extension with the quality matrix ; add the extension with the quality matrix
extension_header = !null extension_header = !null
fxaddpar, extension_header, 'PCOUNT', 0, 'Parameter count' fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'Group count' fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'Extension name' fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'extension name'
if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL') if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL')
fits_add_checksum, extension_header, quality_matrix fits_add_checksum, extension_header, quality_matrix
mwrfits, float(quality_matrix), out_file_name[2], extension_header, /no_comment, /silent mwrfits, float(quality_matrix), out_file_name[2], extension_header, /no_comment, /silent
...@@ -508,9 +508,9 @@ pro metis_l2_prep_vl_polariz ...@@ -508,9 +508,9 @@ pro metis_l2_prep_vl_polariz
; add the extension with the error matrix ; add the extension with the error matrix
extension_header = !null extension_header = !null
fxaddpar, extension_header, 'PCOUNT', 0, 'Parameter count' fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'Group count' fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'Extension name' fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'extension name'
error_matrix = intarr(header.naxis1, header.naxis2) error_matrix = intarr(header.naxis1, header.naxis2)
if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL') if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL')
fits_add_checksum, extension_header, error_matrix fits_add_checksum, extension_header, error_matrix
...@@ -572,9 +572,9 @@ pro metis_l2_prep_vl_polariz ...@@ -572,9 +572,9 @@ pro metis_l2_prep_vl_polariz
; add the extension with the stokes q image ; add the extension with the stokes q image
extension_header = !null extension_header = !null
fxaddpar, extension_header, 'PCOUNT', 0, 'Parameter count' fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'Group count' fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Stokes Q', 'Extension name' fxaddpar, extension_header, 'EXTNAME', 'Stokes Q', 'extension name'
fxaddpar, extension_header, 'BTYPE', 'Stokes Q' fxaddpar, extension_header, 'BTYPE', 'Stokes Q'
fxaddpar, extension_header, 'BUNIT', cal_pack.vl_channel.cal_units fxaddpar, extension_header, 'BUNIT', cal_pack.vl_channel.cal_units
fxaddpar, extension_header, 'DATAMIN', min(q, /nan) fxaddpar, extension_header, 'DATAMIN', min(q, /nan)
...@@ -589,9 +589,9 @@ pro metis_l2_prep_vl_polariz ...@@ -589,9 +589,9 @@ pro metis_l2_prep_vl_polariz
; add the extension with the stokes u image ; add the extension with the stokes u image
extension_header = !null extension_header = !null
fxaddpar, extension_header, 'PCOUNT', 0, 'Parameter count' fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'Group count' fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Stokes U', 'Extension name' fxaddpar, extension_header, 'EXTNAME', 'Stokes U', 'extension name'
fxaddpar, extension_header, 'BTYPE', 'Stokes U' fxaddpar, extension_header, 'BTYPE', 'Stokes U'
fxaddpar, extension_header, 'BUNIT', cal_pack.vl_channel.cal_units fxaddpar, extension_header, 'BUNIT', cal_pack.vl_channel.cal_units
fxaddpar, extension_header, 'DATAMIN', min(u, /nan) fxaddpar, extension_header, 'DATAMIN', min(u, /nan)
...@@ -606,9 +606,9 @@ pro metis_l2_prep_vl_polariz ...@@ -606,9 +606,9 @@ pro metis_l2_prep_vl_polariz
; add the extension with the quality matrix ; add the extension with the quality matrix
extension_header = !null extension_header = !null
fxaddpar, extension_header, 'PCOUNT', 0, 'Parameter count' fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'Group count' fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'Extension name' fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'extension name'
if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL') if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL')
fits_add_checksum, extension_header, quality_matrix fits_add_checksum, extension_header, quality_matrix
...@@ -619,9 +619,9 @@ pro metis_l2_prep_vl_polariz ...@@ -619,9 +619,9 @@ pro metis_l2_prep_vl_polariz
; add the extension with the error matrix ; add the extension with the error matrix
extension_header = !null extension_header = !null
fxaddpar, extension_header, 'PCOUNT', 0, 'Parameter count' fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'Group count' fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'Extension name' fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'extension name'
error_matrix = intarr(header.naxis1, header.naxis2) error_matrix = intarr(header.naxis1, header.naxis2)
if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL') if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL')
......
...@@ -108,23 +108,23 @@ function metis_wcs, header, cal_pack, ref_detector = ref_detector ...@@ -108,23 +108,23 @@ function metis_wcs, header, cal_pack, ref_detector = ref_detector
wcs.add, { $ wcs.add, { $
name: 'WCSNAME', $ name: 'WCSNAME', $
value: 'Helioprojective-Cartesian', $ value: 'Helioprojective-Cartesian', $
comment: 'Name of coordinate system'} comment: 'name of coordinate system'}
wcs.add, { $ wcs.add, { $
name: 'CTYPE1', $ name: 'CTYPE1', $
value: ctype1, $ value: ctype1, $
comment: ctype1 eq 'HPLT-TAN' ? 'Helioprojective latitude (Solar Y)' : 'Helioprojective longitude (Solar X)'} comment: ctype1 eq 'HPLT-TAN' ? 'helioprojective latitude (solar y)' : 'helioprojective longitude (solar x)'}
wcs.add, { $ wcs.add, { $
name: 'CTYPE2', $ name: 'CTYPE2', $
value: ctype2, $ value: ctype2, $
comment: ctype2 eq 'HPLT-TAN' ? 'Helioprojective latitude (Solar Y)' : 'Helioprojective longitude (Solar X)'} comment: ctype2 eq 'HPLT-TAN' ? 'helioprojective latitude (solar y)' : 'helioprojective longitude (solar x)'}
wcs.add, { $ wcs.add, { $
name: 'CUNIT1', $ name: 'CUNIT1', $
value: 'arcsec', $ value: 'arcsec', $
comment: 'Units along axis 1'} comment: 'units along axis 1'}
wcs.add, { $ wcs.add, { $
name: 'CUNIT2', $ name: 'CUNIT2', $
value: 'arcsec', $ value: 'arcsec', $
comment: 'Units along axis 2'} comment: 'units along axis 2'}
wcs.add, { $ wcs.add, { $
name: 'PC1_1', $ name: 'PC1_1', $
value: pc[0, 0], $ value: pc[0, 0], $
...@@ -144,47 +144,47 @@ function metis_wcs, header, cal_pack, ref_detector = ref_detector ...@@ -144,47 +144,47 @@ function metis_wcs, header, cal_pack, ref_detector = ref_detector
wcs.add, { $ wcs.add, { $
name: 'CDELT1', $ name: 'CDELT1', $
value: cdelt[0], $ value: cdelt[0], $
comment: '[arcsec] Pixel scale along axis 1'} comment: '[arcsec] pixel scale along axis 1'}
wcs.add, { $ wcs.add, { $
name: 'CDELT2', $ name: 'CDELT2', $
value: cdelt[1], $ value: cdelt[1], $
comment: '[arcsec] Pixel scale along axis 2'} comment: '[arcsec] pixel scale along axis 2'}
wcs.add, { $ wcs.add, { $
name: 'CROTA', $ name: 'CROTA', $
value: atan(pc[0, 1], pc[0, 0]) * !radeg, $ value: atan(pc[0, 1], pc[0, 0]) * !radeg, $
comment: '[deg] Rotation angle'} comment: '[deg] rotation angle'}
wcs.add, { $ wcs.add, { $
name: 'CRVAL1', $ name: 'CRVAL1', $
value: crval[0], $ value: crval[0], $
comment: '[arcsec] Value of reference pixel along axis 1'} comment: '[arcsec] value of reference pixel along axis 1'}
wcs.add, { $ wcs.add, { $
name: 'CRVAL2', $ name: 'CRVAL2', $
value: crval[1], $ value: crval[1], $
comment: '[arcsec] Value of reference pixel along axis 2'} comment: '[arcsec] value of reference pixel along axis 2'}
wcs.add, { $ wcs.add, { $
name: 'CRPIX1', $ name: 'CRPIX1', $
value: crpix[0], $ value: crpix[0], $
comment: '[pixel] Reference pixel location along axis 1'} comment: '[pixel] reference pixel location along axis 1'}
wcs.add, { $ wcs.add, { $
name: 'CRPIX2', $ name: 'CRPIX2', $
value: crpix[1], $ value: crpix[1], $
comment: '[pixel] Reference pixel location along axis 2'} comment: '[pixel] reference pixel location along axis 2'}
wcs.add, { $ wcs.add, { $
name: 'SUN_XCEN', $ name: 'SUN_XCEN', $
value: sunpix[0], $ value: sunpix[0], $
comment: '[pixel] Sun center location along axis 1'} comment: '[pixel] sun center location along axis 1'}
wcs.add, { $ wcs.add, { $
name: 'SUN_YCEN', $ name: 'SUN_YCEN', $
value: sunpix[1], $ value: sunpix[1], $
comment: '[pixel] Sun center location along axis 2'} comment: '[pixel] sun center location along axis 2'}
wcs.add, { $ wcs.add, { $
name: 'SUNPIX1', $ name: 'SUNPIX1', $
value: sunpix[0], $ value: sunpix[0], $
comment: '[pixel] Sun center location along axis 1'} comment: '[pixel] sun center location along axis 1'}
wcs.add, { $ wcs.add, { $
name: 'SUNPIX2', $ name: 'SUNPIX2', $
value: sunpix[1], $ value: sunpix[1], $
comment: '[pixel] Sun center location along axis 2'} comment: '[pixel] sun center location along axis 2'}
wcs.add, { $ wcs.add, { $
name: 'IO_XCEN', $ name: 'IO_XCEN', $
value: borpix[0], $ value: borpix[0], $
...@@ -248,11 +248,11 @@ function metis_wcs, header, cal_pack, ref_detector = ref_detector ...@@ -248,11 +248,11 @@ function metis_wcs, header, cal_pack, ref_detector = ref_detector
wcs.add, { $ wcs.add, { $
name: 'INN_FOV', $ name: 'INN_FOV', $
value: 1.6, $ value: 1.6, $
comment: '[deg] Inner Metis FOV'} comment: '[deg] inner Metis FOV'}
wcs.add, { $ wcs.add, { $
name: 'OUT_FOV', $ name: 'OUT_FOV', $
value: 3.4, $ value: 3.4, $
comment: '[deg] Outer Metis FOV'} comment: '[deg] outer Metis FOV'}
return, wcs return, wcs
end end
...@@ -16,7 +16,7 @@ function solo_get_ephemeris, header, cal_pack ...@@ -16,7 +16,7 @@ function solo_get_ephemeris, header, cal_pack
ephemerides.add, { $ ephemerides.add, { $
name: 'LONPOLE', $ name: 'LONPOLE', $
value: 180., $ value: 180., $
comment: '[deg] Native longitude of the celestial pole'} comment: '[deg] native longitude of the celestial pole'}
; spherical coordinates of the S/C in the heeq (i.e., stonyhurst) frame ; spherical coordinates of the S/C in the heeq (i.e., stonyhurst) frame
...@@ -30,33 +30,33 @@ function solo_get_ephemeris, header, cal_pack ...@@ -30,33 +30,33 @@ function solo_get_ephemeris, header, cal_pack
ephemerides.add, { $ ephemerides.add, { $
name: 'RSUN_ARC', $ name: 'RSUN_ARC', $
value: rsun_arc, $ value: rsun_arc, $
comment: '[arcsec] Apparent photospheric solar radius'} comment: '[arcsec] apparent photospheric solar radius'}
ephemerides.add, { $ ephemerides.add, { $
name: 'RSUN_REF', $ name: 'RSUN_REF', $
value: rsun, $ value: rsun, $
comment: '[m] Assumed physical solar radius'} comment: '[m] assumed physical solar radius'}
solar_angles = solo_get_solar_angles(utc) solar_angles = solo_get_solar_angles(utc)
ephemerides.add, { $ ephemerides.add, { $
name: 'SOLAR_B0', $ name: 'SOLAR_B0', $
value: solar_angles[0], $ value: solar_angles[0], $
comment: '[deg] S/C tilt of solar North pole'} comment: '[deg] S/C tilt of solar north pole'}
ephemerides.add, { $ ephemerides.add, { $
name: 'SOLAR_P0 ', $ name: 'SOLAR_P0 ', $
value: solar_angles[1], $ value: solar_angles[1], $
comment: '[deg] S/C celestial North to solar North angle'} comment: '[deg] S/C celestial north to solar north angle'}
ephemerides.add, { $ ephemerides.add, { $
name: 'SOLAR_EP', $ name: 'SOLAR_EP', $
value: solar_angles[2], $ value: solar_angles[2], $
comment: '[deg] S/C ecliptic North to solar North angle'} comment: '[deg] S/C ecliptic north to solar north angle'}
carrot = solo_get_carringrot(utc) carrot = solo_get_carringrot(utc)
ephemerides.add, { $ ephemerides.add, { $
name: 'CAR_ROT', $ name: 'CAR_ROT', $
value: carrot, $ value: carrot, $
comment: 'Carrington rotation number'} comment: 'carrington rotation number'}
ephemerides.add, { $ ephemerides.add, { $
name: 'HGLT_OBS', $ name: 'HGLT_OBS', $
...@@ -74,24 +74,24 @@ function solo_get_ephemeris, header, cal_pack ...@@ -74,24 +74,24 @@ function solo_get_ephemeris, header, cal_pack
ephemerides.add, { $ ephemerides.add, { $
name: 'CRLT_OBS', $ name: 'CRLT_OBS', $
value: coord[2], $ value: coord[2], $
comment: '[deg] S/C Carrington latitude (B0 angle)'} comment: '[deg] S/C carrington latitude (B0 angle)'}
ephemerides.add, { $ ephemerides.add, { $
name: 'CRLN_OBS', $ name: 'CRLN_OBS', $
value: coord[1], $ value: coord[1], $
comment: '[deg] S/C Carrington longitude (L0 angle)'} comment: '[deg] S/C carrington longitude (L0 angle)'}
ephemerides.add, { $ ephemerides.add, { $
name: 'DSUN_OBS', $ name: 'DSUN_OBS', $
value: solo_dist, $ value: solo_dist, $
comment: '[m] S/C distance from Sun'} comment: '[m] S/C distance from sun'}
ephemerides.add, { $ ephemerides.add, { $
name: 'DSUN_AU', $ name: 'DSUN_AU', $
value: solo_dist/au, $ value: solo_dist/au, $
comment: '[AU] S/C distance from Sun'} comment: '[AU] S/C distance from sun'}
ephemerides.add, { $ ephemerides.add, { $
name: 'AU_REF', $ name: 'AU_REF', $
value: au, $ value: au, $
comment: '[m] Assumed physical Astronomical Unit'} comment: '[m] assumed physical astronomical unit'}
; coordinates of the S/C in the hee frame ; coordinates of the S/C in the hee frame
...@@ -101,15 +101,15 @@ function solo_get_ephemeris, header, cal_pack ...@@ -101,15 +101,15 @@ function solo_get_ephemeris, header, cal_pack
ephemerides.add, { $ ephemerides.add, { $
name: 'HEEX_OBS', $ name: 'HEEX_OBS', $
value: coord[0], $ value: coord[0], $
comment: '[m] S/C Heliocentric Earth Ecliptic X'} comment: '[m] S/C heliocentric earth ecliptic x'}
ephemerides.add, { $ ephemerides.add, { $
name: 'HEEY_OBS', $ name: 'HEEY_OBS', $
value: coord[1], $ value: coord[1], $
comment: '[m] S/C Heliocentric Earth Ecliptic Y'} comment: '[m] S/C heliocentric earth ecliptic y'}
ephemerides.add, { $ ephemerides.add, { $
name: 'HEEZ_OBS', $ name: 'HEEZ_OBS', $
value: coord[2], $ value: coord[2], $
comment: '[m] S/C Heliocentric Earth Ecliptic Z'} comment: '[m] S/C heliocentric earth ecliptic z'}
; coordinates of the S/C in the hci frame ; coordinates of the S/C in the hci frame
...@@ -120,28 +120,28 @@ function solo_get_ephemeris, header, cal_pack ...@@ -120,28 +120,28 @@ function solo_get_ephemeris, header, cal_pack
ephemerides.add, { $ ephemerides.add, { $
name: 'HCIX_OBS', $ name: 'HCIX_OBS', $
value: coord[0], $ value: coord[0], $
comment: '[m] S/C Heliocentric Inertial X'} comment: '[m] S/C heliocentric inertial x'}
ephemerides.add, { $ ephemerides.add, { $
name: 'HCIY_OBS', $ name: 'HCIY_OBS', $
value: coord[1], $ value: coord[1], $
comment: '[m] S/C Heliocentric Inertial Y'} comment: '[m] S/C heliocentric inertial y'}
ephemerides.add, { $ ephemerides.add, { $
name: 'HCIZ_OBS', $ name: 'HCIZ_OBS', $
value: coord[2], $ value: coord[2], $
comment: '[m] S/C Heliocentric Inertial Z'} comment: '[m] S/C heliocentric inertial z'}
ephemerides.add, { $ ephemerides.add, { $
name: 'HCIX_VOB', $ name: 'HCIX_VOB', $
value: vel[0], $ value: vel[0], $
comment: '[m/s] S/C Heliocentric Inertial X velocity'} comment: '[m/s] S/C heliocentric inertial x velocity'}
ephemerides.add, { $ ephemerides.add, { $
name: 'HCIY_VOB', $ name: 'HCIY_VOB', $
value: vel[1], $ value: vel[1], $
comment: '[m/s] S/C Heliocentric Inertial Y velocity'} comment: '[m/s] S/C heliocentric inertial y velocity'}
ephemerides.add, { $ ephemerides.add, { $
name: 'HCIZ_VOB', $ name: 'HCIZ_VOB', $
value: vel[2], $ value: vel[2], $
comment: '[m/s] S/C Heliocentric Inertial Z velocity'} comment: '[m/s] S/C heliocentric inertial z velocity'}
; coordinates of the S/C in the hae frame ; coordinates of the S/C in the hae frame
...@@ -151,15 +151,15 @@ function solo_get_ephemeris, header, cal_pack ...@@ -151,15 +151,15 @@ function solo_get_ephemeris, header, cal_pack
ephemerides.add, { $ ephemerides.add, { $
name: 'HAEX_OBS', $ name: 'HAEX_OBS', $
value: coord[0], $ value: coord[0], $
comment: '[m] S/C Heliocentric Aries Ecliptic X'} comment: '[m] S/C heliocentric aries ecliptic x'}
ephemerides.add, { $ ephemerides.add, { $
name: 'HAEY_OBS', $ name: 'HAEY_OBS', $
value: coord[1], $ value: coord[1], $
comment: '[m] S/C Heliocentric Aries Ecliptic Y'} comment: '[m] S/C heliocentric aries ecliptic y'}
ephemerides.add, { $ ephemerides.add, { $
name: 'HAEZ_OBS', $ name: 'HAEZ_OBS', $
value: coord[2], $ value: coord[2], $
comment: '[m] S/C Heliocentric Aries Ecliptic Z'} comment: '[m] S/C heliocentric aries ecliptic z'}
; coordinates of the S/C in the heeq frame ; coordinates of the S/C in the heeq frame
...@@ -169,15 +169,15 @@ function solo_get_ephemeris, header, cal_pack ...@@ -169,15 +169,15 @@ function solo_get_ephemeris, header, cal_pack
ephemerides.add, { $ ephemerides.add, { $
name: 'HEQX_OBS', $ name: 'HEQX_OBS', $
value: coord[0], $ value: coord[0], $
comment: '[m] S/C Heliocentric Earth Equatorial X'} comment: '[m] S/C heliocentric earth equatorial x'}
ephemerides.add, { $ ephemerides.add, { $
name: 'HEQY_OBS', $ name: 'HEQY_OBS', $
value: coord[1], $ value: coord[1], $
comment: '[m] S/C Heliocentric Earth Equatorial Y'} comment: '[m] S/C heliocentric earth equatorial y'}
ephemerides.add, { $ ephemerides.add, { $
name: 'HEQZ_OBS', $ name: 'HEQZ_OBS', $
value: coord[2], $ value: coord[2], $
comment: '[m] S/C Heliocentric Earth Equatorial Z'} comment: '[m] S/C heliocentric earth equatorial z'}
; coordinates of the S/C in the gse frame ; coordinates of the S/C in the gse frame
...@@ -187,15 +187,15 @@ function solo_get_ephemeris, header, cal_pack ...@@ -187,15 +187,15 @@ function solo_get_ephemeris, header, cal_pack
ephemerides.add, { $ ephemerides.add, { $
name: 'GSEX_OBS', $ name: 'GSEX_OBS', $
value: coord[0], $ value: coord[0], $
comment: '[m] S/C Geocentric Solar Ecliptic X'} comment: '[m] S/C geocentric solar ecliptic x'}
ephemerides.add, { $ ephemerides.add, { $
name: 'GSEY_OBS', $ name: 'GSEY_OBS', $
value: coord[1], $ value: coord[1], $
comment: '[m] S/C Geocentric Solar Ecliptic Y'} comment: '[m] S/C geocentric solar ecliptic y'}
ephemerides.add, { $ ephemerides.add, { $
name: 'GSEZ_OBS', $ name: 'GSEZ_OBS', $
value: coord[2], $ value: coord[2], $
comment: '[m] S/C Geocentric Solar Ecliptic Z'} comment: '[m] S/C geocentric solar ecliptic z'}
; light travel times and radial velocity of the S/C ; light travel times and radial velocity of the S/C
...@@ -206,15 +206,15 @@ function solo_get_ephemeris, header, cal_pack ...@@ -206,15 +206,15 @@ function solo_get_ephemeris, header, cal_pack
ephemerides.add, { $ ephemerides.add, { $
name: 'OBS_VR', $ name: 'OBS_VR', $
value: rad_vel, $ value: rad_vel, $
comment: '[m/s] Radial velocity of S/C relative to Sun'} comment: '[m/s] radial velocity of S/C relative to sun'}
ephemerides.add, { $ ephemerides.add, { $
name: 'EAR_TDEL', $ name: 'EAR_TDEL', $
value: t_del, $ value: t_del, $
comment: '[s] Time(Sun to Earth) - Time(Sun to S/C)'} comment: '[s] time(sun to earth) - time(sun to S/C)'}
ephemerides.add, { $ ephemerides.add, { $
name: 'SUN_TIME', $ name: 'SUN_TIME', $
value: sun_time, $ value: sun_time, $
comment: '[s] Time(Sun to S/C)'} comment: '[s] time(sun to S/C)'}
; corrections of the acquisition date ; corrections of the acquisition date
...@@ -227,11 +227,11 @@ function solo_get_ephemeris, header, cal_pack ...@@ -227,11 +227,11 @@ function solo_get_ephemeris, header, cal_pack
ephemerides.add, { $ ephemerides.add, { $
name: 'DATE_EAR', $ name: 'DATE_EAR', $
value: date_ear, $ value: date_ear, $
comment: '[UTC] Obs. start time corrected to Earth'} comment: '[UTC] obs. start time corrected to earth'}
ephemerides.add, { $ ephemerides.add, { $
name: 'DATE_SUN', $ name: 'DATE_SUN', $
value: date_sun, $ value: date_sun, $
comment: '[UTC] Obs. start time corrected to Sun'} comment: '[UTC] obs. start time corrected to sun'}
return, ephemerides return, ephemerides
end end
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment