1) Calculation routine in in fsi_cr.py --> only for single image
2) CR_timeseries.py -looping routine that execute fsi_cr.py

3) Start code:
conda activate sunpy_test 
ipython --pylab 
plot_CRhists(Ith=500,save_figure=False, figx=9,figy=7, fontsize=16, ticklabelsize=14)


; SAVING RESULTS
CR=plot_CRhists(Ith=500,save_figure=False, figx=9,figy=7, fontsize=16, ticklabelsize=14)




Code generation:
conda activate sunpy_test 
ipython --pylab 
%run CR_timeseries_v2.py
CR=plot_CRhists(Ith=500,save_figure=True, figx=9,figy=7, fontsize=16, ticklabelsize=14, data_save=True)
res=plot_CR_timeseries(filename_txt=output_path+'/all.txt', filename_save=output_path+'/final_result.png')



Code generation:
conda activate sunpy_test 
ipython --pylab 
%run CR_timeseries_v2.py
exi
res=plot_CR_timeseries(filename_txt=output_path+'/all.txt', filename_save=output_path+'/final_result.png')