;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

; .r aia_image_psf_correction_test_v1.pro

; path to the raw data AIA FITS file
path0='~/data/sdo/aia/data/FITS/'
; raw data AIA FITS filename
aia_raw_data_fits='aia.lev1.131A_2014-09-10T16_20_08.63Z.image_lev1.fits'

; channel
channel=131

; reading the raw data AIA FITS file -> index_raw,data_raw
read_sdo,path0+aia_raw_data_fits,index_raw,data_raw

; data_raw size
data_raw_size=size(data_raw)
; required PSF data size
data_size=max([data_raw_size(1),data_raw_size(2)])

; creation of the FFT PSF for the selected channel and data size
aia_psf_fft_creation_jr, channel, data_size, aia_psf_fft_channel, error=error
if error ne 0 then stop ; due to wrong action of the procedure 'aia_psf_fft_creation_jr'

; raw AIA data correction for the PSF influence
aia_psf_correction_jr, data_raw, aia_psf_fft_channel, data_cor, error=error
if error ne 0 then stop         ; due to wrong action of the procedure 'aia_psf_correction_jr'

end