Commit 92de877f authored by Roberto Susino's avatar Roberto Susino
Browse files

Fix bug in output filenames definition

parent 7bb12d79

metis_l2_prep_vl_polariz.pro

+26 −22
Original line number Diff line number Diff line
@@ -324,6 +324,10 @@ pro metis_l2_prep_vl_polariz 
	date_beg_string = date_beg_string.replace('-', '')

	date_beg_string = date_beg_string.replace(':', '')



	; array of output file names



	out_file_name = strarr(4)



	; keywords specific for polarized brightness images



	primary_pb_header = primary_header
@@ -331,7 +335,7 @@ pro metis_l2_prep_vl_polariz 
	; name of the fits file



	file_name = 'solo_L2_metis-vl-pb_' + date_beg_string + '_V' + version + '.fits'

	out_file_name = 'output/' + file_name

	out_file_name[0] = 'output/' + file_name



	fxaddpar, primary_pb_header, 'FILENAME', file_name

	fxaddpar, primary_pb_header, 'BTYPE', 'VL polarized brightness'
@@ -350,10 +354,10 @@ pro metis_l2_prep_vl_polariz 


	if not ref_detector then pb_image = metis_rectify(pb_image, 'VL')

	fits_add_checksum, primary_pb_header, pb_image

	mwrfits, pb_image, out_file_name, primary_pb_header, /no_comment, /create, /silent

	mwrfits, pb_image, out_file_name[0], primary_pb_header, /no_comment, /create, /silent



	journal, 'Polarized-brightness FITS file created:'

	journal, '  file name = ' + file_basename(out_file_name)

	journal, '  file name = ' + file_basename(out_file_name[0])



	; add the extension with the quality matrix
@@ -363,7 +367,7 @@ pro metis_l2_prep_vl_polariz 
	fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'Extension name'

	if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL')

	fits_add_checksum, extension_header, quality_matrix

	mwrfits, quality_matrix, out_file_name, extension_header, /no_comment, /silent

	mwrfits, quality_matrix, out_file_name[0], extension_header, /no_comment, /silent



	journal, 'Quality-matrix extension correctly added.'
@@ -376,7 +380,7 @@ pro metis_l2_prep_vl_polariz 
	error_matrix = intarr(header.naxis1, header.naxis2)

	if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL')

	fits_add_checksum, extension_header, error_matrix

	mwrfits, error_matrix, out_file_name, extension_header, /no_comment, /silent

	mwrfits, error_matrix, out_file_name[0], extension_header, /no_comment, /silent



	journal, 'Error-matrix extension correctly added.'
@@ -387,7 +391,7 @@ pro metis_l2_prep_vl_polariz 
	; name of the fits file



	file_name = 'solo_L2_metis-vl-tb_' + date_beg_string + '_V' + version + '.fits'

	out_file_name = 'output/' + file_name

	out_file_name[1] = 'output/' + file_name



	fxaddpar, primary_tb_header, 'FILENAME', file_name

	fxaddpar, primary_tb_header, 'BTYPE', 'VL total brightness'
@@ -406,10 +410,10 @@ pro metis_l2_prep_vl_polariz 


	if not ref_detector then tb_image = metis_rectify(tb_image, 'VL')

	fits_add_checksum, primary_tb_header, tb_image

	mwrfits, tb_image, out_file_name, primary_tb_header, /no_comment, /create, /silent

	mwrfits, tb_image, out_file_name[1], primary_tb_header, /no_comment, /create, /silent



	journal, 'Total-brightness FITS file created:'

	journal, '  file name = ' + file_basename(out_file_name)

	journal, '  file name = ' + file_basename(out_file_name[1])



	; add the extension with the quality matrix
@@ -419,7 +423,7 @@ pro metis_l2_prep_vl_polariz 
	fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'Extension name'

	if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL')

	fits_add_checksum, extension_header, quality_matrix

	mwrfits, quality_matrix, out_file_name, extension_header, /no_comment, /silent

	mwrfits, quality_matrix, out_file_name[1], extension_header, /no_comment, /silent



	journal, 'Quality-matrix extension correctly added.'
@@ -432,7 +436,7 @@ pro metis_l2_prep_vl_polariz 
	error_matrix = intarr(header.naxis1, header.naxis2)

	if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL')

	fits_add_checksum, extension_header, error_matrix

	mwrfits, error_matrix, out_file_name, extension_header, /no_comment, /silent

	mwrfits, error_matrix, out_file_name[1], extension_header, /no_comment, /silent



	journal, 'Error-matrix extension correctly added.'
@@ -443,7 +447,7 @@ pro metis_l2_prep_vl_polariz 
	; name of the fits file



	file_name = 'solo_L2_metis-vl-pol-angle_' + date_beg_string + '_V' + version + '.fits'

	out_file_name = 'output/' + file_name

	out_file_name[2] = 'output/' + file_name



	fxaddpar, primary_polangle_header, 'FILENAME', file_name

	fxaddpar, primary_polangle_header, 'BTYPE', 'VL polarization angle'
@@ -462,10 +466,10 @@ pro metis_l2_prep_vl_polariz 


	if not ref_detector then pol_angle = metis_rectify(pol_angle, 'VL')

	fits_add_checksum, primary_polangle_header, pol_angle

	mwrfits, pol_angle, out_file_name, primary_polangle_header, /no_comment, /create, /silent

	mwrfits, pol_angle, out_file_name[2], primary_polangle_header, /no_comment, /create, /silent



	journal, 'Polarization-angle FITS file created:'

	journal, '  file name = ' + file_basename(out_file_name)

	journal, '  file name = ' + file_basename(out_file_name[2])



	; add the extension with the quality matrix
@@ -475,7 +479,7 @@ pro metis_l2_prep_vl_polariz 
	fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'Extension name'

	if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL')

	fits_add_checksum, extension_header, quality_matrix

	mwrfits, quality_matrix, out_file_name, extension_header, /no_comment, /silent

	mwrfits, quality_matrix, out_file_name[2], extension_header, /no_comment, /silent



	journal, 'Quality-matrix extension correctly added.'
@@ -488,7 +492,7 @@ pro metis_l2_prep_vl_polariz 
	error_matrix = intarr(header.naxis1, header.naxis2)

	if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL')

	fits_add_checksum, extension_header, error_matrix

	mwrfits, error_matrix, out_file_name, extension_header, /no_comment, /silent

	mwrfits, error_matrix, out_file_name[2], extension_header, /no_comment, /silent



	journal, 'Error-matrix extension correctly added.'
@@ -522,7 +526,7 @@ pro metis_l2_prep_vl_polariz 
	; name of the fits file



	file_name = 'solo_L2_metis-vl-stokes_' + date_beg_string + '_V' + version + '.fits'

	out_file_name = 'output/' + file_name

	out_file_name[3] = 'output/' + file_name



	fxaddpar, primary_stokes_header, 'FILENAME', file_name

	fxaddpar, primary_stokes_header, 'BTYPE', 'Stokes I'
@@ -540,10 +544,10 @@ pro metis_l2_prep_vl_polariz 
	; add checksum and datasum to the fits header

	if not ref_detector then i = metis_rectify(i, 'VL')

	fits_add_checksum, primary_stokes_header, i

	mwrfits, i, out_file_name, primary_stokes_header, /no_comment, /create, /silent

	mwrfits, i, out_file_name[3], primary_stokes_header, /no_comment, /create, /silent



	journal, 'Stokes parameters FITS file created:'

	journal, '  file name = ' + file_basename(out_file_name)

	journal, '  file name = ' + file_basename(out_file_name[3])



	; add the extension with the stokes q image
@@ -557,7 +561,7 @@ pro metis_l2_prep_vl_polariz 
	fxaddpar, extension_header, 'DATAMAX', max(q, /nan)

	if not ref_detector then q = metis_rectify(q, 'VL')

	fits_add_checksum, extension_header, q

	mwrfits, q, out_file_name, extension_header, /no_comment, /silent

	mwrfits, q, out_file_name[3], extension_header, /no_comment, /silent



	journal, 'Q parameter extension correctly added.'
@@ -573,7 +577,7 @@ pro metis_l2_prep_vl_polariz 
	fxaddpar, extension_header, 'DATAMAX', max(u, /nan)

	if not ref_detector then u = metis_rectify(u, 'VL')

	fits_add_checksum, extension_header, u

	mwrfits, u, out_file_name, extension_header, /no_comment, /silent

	mwrfits, u, out_file_name[3], extension_header, /no_comment, /silent



	journal, 'U parameter extension correctly added.'
@@ -585,7 +589,7 @@ pro metis_l2_prep_vl_polariz 
	fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'Extension name'

	if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL')

	fits_add_checksum, extension_header, quality_matrix

	mwrfits, quality_matrix, out_file_name, extension_header, /no_comment, /silent

	mwrfits, quality_matrix, out_file_name[3], extension_header, /no_comment, /silent



	journal, 'Quality-matrix extension correctly added.'
@@ -598,7 +602,7 @@ pro metis_l2_prep_vl_polariz 
	error_matrix = intarr(header.naxis1, header.naxis2)

	if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL')

	fits_add_checksum, extension_header, error_matrix

	mwrfits, error_matrix, out_file_name, extension_header, /no_comment, /silent

	mwrfits, error_matrix, out_file_name[3], extension_header, /no_comment, /silent



	journal, 'Error-matrix extension correctly added.'