Commit 145762e2 authored by Roberto Susino's avatar Roberto Susino
Browse files

Fix ephemerides calculation

parent 57bad12a

metis_l2_prep_uv.pro

+1 −3
Original line number Diff line number Diff line
@@ -121,13 +121,11 @@ pro metis_l2_prep_uv 
	; append wcs keywords



	wcs = metis_wcs(header, cal_pack, ref_detector = ref_detector)



	foreach element, wcs do fxaddpar, primary_header, element.name, element.value, element.comment, before = 'DATATYPE'



	; append solar ephemeris keywords



	ephemerides = solo_get_ephemerides(header.date_avg, constants = cal_pack.constants)



	ephemerides = solo_get_ephemerides(header, constants = cal_pack.constants)

	foreach element, ephemerides do fxaddpar, primary_header, element.name, element.value, element.comment, before = 'DATATYPE'



	history = [history, 'Solar ephemerides and WCS:', '  SKD version = ' + kernel_version]

metis_l2_prep_vl_generic.pro

+1 −2
Original line number Diff line number Diff line
@@ -123,8 +123,7 @@ pro metis_l2_prep_vl_generic 


	; append solar ephemeris keywords



	ephemerides = solo_get_ephemerides(header.date_avg, constants = cal_pack.constants)



	ephemerides = solo_get_ephemerides(header, constants = cal_pack.constants)

	foreach element, ephemerides do fxaddpar, primary_header, element.name, element.value, element.comment, before = 'DATATYPE'



	history = [history, 'Solar ephemerides and WCS:', '  SKD version = ' + kernel_version]

metis_l2_prep_vl_polariz.pro

+1 −2
Original line number Diff line number Diff line
@@ -235,8 +235,7 @@ pro metis_l2_prep_vl_polariz 


	; append solar ephemeris keywords



	ephemerides = solo_get_ephemerides(header.date_avg, constants = cal_pack.constants)



	ephemerides = solo_get_ephemerides(header, constants = cal_pack.constants)

	foreach element, ephemerides do fxaddpar, primary_header, element.name, element.value, element.comment, before = 'DATATYPE'



	history = ['Solar ephemerides and WCS:', '  SKD version = ' + kernel_version]

solo_get_ephemerides.pro

+231 −232
Original line number Diff line number Diff line 
function solo_get_ephemerides, utc, constants = constants

function solo_get_ephemerides, header, constants = constants



	utc = header.date_avg



    if ~ keyword_set(constants) then begin

        rsun = 695508000.0 ; metres

        au = 149597870691.0 ; metres

    endif else begin

	rsun = constants.rsun.value

	au = constants.au.value

    endelse



	; compute and add the wcs keyword
@@ -216,6 +213,8 @@ function solo_get_ephemerides, utc, constants = constants 


	; corrections of the acquisition date



	utc = header.date_beg



	jul_utc = date_conv(utc, 'julian')

	date_earth = date_conv(jul_utc + tdel / 86400.d0, 'fits')

	date_sun = date_conv(jul_utc - sun_time / 86400.d0, 'fits')