; IDL Version 6.2 (linux x86 m32) ; Journal File for schmahl@hessi ; Working directory: /home/schmahl/footpt_sizes ; Date: Fri Mar 3 11:15:24 2006 time=['mar 18','nov 06','nov 19','oct 24','aug 20','feb 20','jul 23','jun 09','jun 10','oct 30','sep 08','sep 20','jun 17'] print,time ;mar 18 nov 06 nov 19 oct 24 aug 20 feb 20 jul 23 jun 09 jun 10 oct 30 sep 08 sep 20 jun 17 help,time flux1=fltarr(6,13) flux2=flux1 flux1(*,0)=[0,0,0,0,0,0] flux1(*,2)=[26.7,5.9,8.2,8.2,1.7,0]] ; % Syntax error. flux2(*,2)=[6.9,22,7.4,7.4,10.5,0] flux1(*,3)=[4.9,9.3,6.5,4.9,0,0] flux2(*,3)=[4.4,6.3,5.0,5.1,0,0] flux1(*,4)=[20,12.4,0,0,0,0] flux2(*,4)=[5.8,10.2,0,0,0,0] flux1(*,6)=[3.5,3.1,3.4,0,0,0] flux2(*,6)=[17,8.8,8.1,0,0,0] for j=0,12 do print,j,' ',time(j) ; 0 mar 18 ; 1 nov 06 ; 2 nov 19 ; 3 oct 24 ; 4 aug 20 ; 5 feb 20 ; 6 jul 23 ; 7 jun 09 ; 8 jun 10 ; 9 oct 30 ; 10 sep 08 ; 11 sep 20 ; 12 jun 17 flux1(*,7)=[16.6,0,0,0,0,0] flux2(*,7)=[11.3,0,0,0,0,0] flux1(*,8)=[6.8,4.4,13,0,0,0] flux2(*,8)=[5.5,8.1,16,0,0,0] flux1(*,9)=[6.5,5.8,4.7,0,0,0] flux2(*,9)=[20,18.1,13.8,0,0,0] flux1(*,10)=[10.7,8.2,7.4,5.9,7.7,7.8] flux2(*,10)=[4.7,3.2,3.8,5.8,3.9,0.5] flux1(*,11)=[8.9,9.8,2.1,3.9,0,0] flux2(*,11)=[36,5.0,10.8,4.0,0,0] fwhm1=flux1 fwhm2=flux2 flux1=0*flux1 flux2=0*flux2 save,time,flux1,flux2,fwhm1,fwhm2,file='fwhm_flux.sav' fwhm1(*,12)=[6.2,5.3,4.4,4.8,0,0] fwhm2(*,12)=[3.4,3.8,3.7,3.4,0,0] save,time,flux1,flux2,fwhm1,fwhm2,file='fwhm_flux.sav' plot,fwhm1,fwhm2,psym=6 flux1=[0,0,0,0,0,0] help,flux2 flux1=flux2 flux1(*,2)=[26.7,5.9,8.2,8.2,1.7,0] flux2(*,2)=[8.3,20.2,6.9,6.9,4.2,0] flux1(*,3)=[3.7,9.3,6.6,3.7,0,,0] ; % Syntax error. flux1(*,3)=[3.7,9.3,6.6,3.7,0,0] flux2(*,3)=[2.2,7.2,6.1,3.1,0,0] xdoc,'f_div ;---- Module: f_div.pro ;---- From: /usr/local/ssw/gen/idl/util ;---> Call: FUNCTION F_DIV, NUMERATOR, DENOMINATOR, DEFAULT = DEFAULT print,f_div(flux1/flux2) ; % Variable is undefined: DENOMINATOR. ; % Program caused arithmetic error: Floating illegal operand print,f_div(flux1,flux2) ; % Variable is undefined: DENOMINATOR. xdoc,'f_div ; % CHKARG: recalling f_div.pro from memory ;---- Module: f_div.pro ;---- From: /usr/local/ssw/gen/idl/util/ ;---> Call: FUNCTION F_DIV, NUMERATOR, DENOMINATOR, DEFAULT = DEFAULT help,flux1,flux2 print,f_div(flux1,flux2) ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 3.21687 0.292079 1.18841 1.18841 0.404762 0.00000 ; 1.68182 1.29167 1.08197 1.19355 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 print,f_div(fwhm1,fwhm2) ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 1.11364 1.47619 1.30000 0.960784 0.00000 0.00000 ; 3.44828 1.21569 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.205882 0.352273 0.419753 0.00000 0.00000 0.00000 ; 1.46903 0.00000 0.00000 0.00000 0.00000 0.00000 ; 1.23636 0.543210 0.812500 0.00000 0.00000 0.00000 ; 0.325000 0.320442 0.340580 0.00000 0.00000 0.00000 ; 2.27660 2.56250 1.94737 1.01724 1.97436 15.6000 ; 0.247222 1.96000 0.194444 0.975000 0.00000 0.00000 ; 1.82353 1.39474 1.18919 1.41176 0.00000 0.00000 plot,f_div(fwhm1,fwhm2),f_div(flux1,flux2),psym=6 plot,f_div(fwhm1,fwhm2),f_div(flux1,flux2),psym=6,xtit='fwhm1/fwhm2',ytit='flux1/flux2' for j=0,12 do print,j,' ',time(j) ; 0 mar 18 ; 1 nov 06 ; 2 nov 19 ; 3 oct 24 ; 4 aug 20 ; 5 feb 20 ; 6 jul 23 ; 7 jun 09 ; 8 jun 10 ; 9 oct 30 ; 10 sep 08 ; 11 sep 20 ; 12 jun 17 flux1(*,12)=[40.3,49.2,46.6,31.4,0,0] flux2(*,12)=[30.0,37.5,37.7,24.1,0,0] plot,f_div(fwhm1,fwhm2),f_div(flux1,flux2),psym=6,xtit='fwhm1/fwhm2',ytit='flux1/flux2' print,fwhm1 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 4.90000 9.30000 6.50000 4.90000 0.00000 0.00000 ; 20.0000 12.4000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 3.50000 3.10000 3.40000 0.00000 0.00000 0.00000 ; 16.6000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 6.80000 4.40000 13.0000 0.00000 0.00000 0.00000 ; 6.50000 5.80000 4.70000 0.00000 0.00000 0.00000 ; 10.7000 8.20000 7.40000 5.90000 7.70000 7.80000 ; 8.90000 9.80000 2.10000 3.90000 0.00000 0.00000 ; 6.20000 5.30000 4.40000 4.80000 0.00000 0.00000 print,fwhm2 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 6.90000 22.0000 7.40000 7.40000 10.5000 0.00000 ; 4.40000 6.30000 5.00000 5.10000 0.00000 0.00000 ; 5.80000 10.2000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 17.0000 8.80000 8.10000 0.00000 0.00000 0.00000 ; 11.3000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 5.50000 8.10000 16.0000 0.00000 0.00000 0.00000 ; 20.0000 18.1000 13.8000 0.00000 0.00000 0.00000 ; 4.70000 3.20000 3.80000 5.80000 3.90000 0.500000 ; 36.0000 5.00000 10.8000 4.00000 0.00000 0.00000 ; 3.40000 3.80000 3.70000 3.40000 0.00000 0.00000 print,f_div(fwhm1,fwhm2) ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 1.11364 1.47619 1.30000 0.960784 0.00000 0.00000 ; 3.44828 1.21569 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.205882 0.352273 0.419753 0.00000 0.00000 0.00000 ; 1.46903 0.00000 0.00000 0.00000 0.00000 0.00000 ; 1.23636 0.543210 0.812500 0.00000 0.00000 0.00000 ; 0.325000 0.320442 0.340580 0.00000 0.00000 0.00000 ; 2.27660 2.56250 1.94737 1.01724 1.97436 15.6000 ; 0.247222 1.96000 0.194444 0.975000 0.00000 0.00000 ; 1.82353 1.39474 1.18919 1.41176 0.00000 0.00000 print,fwhm1(*,10) ; 10.7000 8.20000 7.40000 5.90000 7.70000 7.80000 print,fwhm2(*,10) ; 4.70000 3.20000 3.80000 5.80000 3.90000 0.500000 fwhm2(4,10)=0 plot,f_div(fwhm1,fwhm2),f_div(flux1,flux2),psym=6,xtit='fwhm1/fwhm2',ytit='flux1/flux2' print,f_div(fwhm1,fwhm2) ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 1.11364 1.47619 1.30000 0.960784 0.00000 0.00000 ; 3.44828 1.21569 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.205882 0.352273 0.419753 0.00000 0.00000 0.00000 ; 1.46903 0.00000 0.00000 0.00000 0.00000 0.00000 ; 1.23636 0.543210 0.812500 0.00000 0.00000 0.00000 ; 0.325000 0.320442 0.340580 0.00000 0.00000 0.00000 ; 2.27660 2.56250 1.94737 1.01724 0.00000 15.6000 ; 0.247222 1.96000 0.194444 0.975000 0.00000 0.00000 ; 1.82353 1.39474 1.18919 1.41176 0.00000 0.00000 print,(f_div(fwhm1,fwhm2))[*,10] ; 2.27660 2.56250 1.94737 1.01724 0.00000 15.6000 print,fwhm2(*,10) ; 4.70000 3.20000 3.80000 5.80000 0.00000 0.500000 fwhm2(5,10)=0 $recent sav fwhm2_bak=fwhm2 restore,'fwhm_flux.sav',/ver print,fwhm2 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 6.90000 22.0000 7.40000 7.40000 10.5000 0.00000 ; 4.40000 6.30000 5.00000 5.10000 0.00000 0.00000 ; 5.80000 10.2000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 17.0000 8.80000 8.10000 0.00000 0.00000 0.00000 ; 11.3000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 5.50000 8.10000 16.0000 0.00000 0.00000 0.00000 ; 20.0000 18.1000 13.8000 0.00000 0.00000 0.00000 ; 4.70000 3.20000 3.80000 5.80000 3.90000 0.500000 ; 36.0000 5.00000 10.8000 4.00000 0.00000 0.00000 ; 3.40000 3.80000 3.70000 3.40000 0.00000 0.00000 print,fwhm2_bak ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 6.90000 22.0000 7.40000 7.40000 10.5000 0.00000 ; 4.40000 6.30000 5.00000 5.10000 0.00000 0.00000 ; 5.80000 10.2000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 17.0000 8.80000 8.10000 0.00000 0.00000 0.00000 ; 11.3000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 5.50000 8.10000 16.0000 0.00000 0.00000 0.00000 ; 20.0000 18.1000 13.8000 0.00000 0.00000 0.00000 ; 4.70000 3.20000 3.80000 5.80000 0.00000 0.00000 ; 36.0000 5.00000 10.8000 4.00000 0.00000 0.00000 ; 3.40000 3.80000 3.70000 3.40000 0.00000 0.00000 print,flux2 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 print,flux1 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 $emacs idlsave.dat& $emacs idlsave.pro& flux1(*,12)=[40.3,49.2,46.6,31.4,0,0] flux2(*,12)=[30.0,37.5,37.7,24.1,0,0] fwhm2(5,10)=0 print,fwhm2 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 6.90000 22.0000 7.40000 7.40000 10.5000 0.00000 ; 4.40000 6.30000 5.00000 5.10000 0.00000 0.00000 ; 5.80000 10.2000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 17.0000 8.80000 8.10000 0.00000 0.00000 0.00000 ; 11.3000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 5.50000 8.10000 16.0000 0.00000 0.00000 0.00000 ; 20.0000 18.1000 13.8000 0.00000 0.00000 0.00000 ; 4.70000 3.20000 3.80000 5.80000 3.90000 0.00000 ; 36.0000 5.00000 10.8000 4.00000 0.00000 0.00000 ; 3.40000 3.80000 3.70000 3.40000 0.00000 0.00000 print,fwhm1 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 4.90000 9.30000 6.50000 4.90000 0.00000 0.00000 ; 20.0000 12.4000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 3.50000 3.10000 3.40000 0.00000 0.00000 0.00000 ; 16.6000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 6.80000 4.40000 13.0000 0.00000 0.00000 0.00000 ; 6.50000 5.80000 4.70000 0.00000 0.00000 0.00000 ; 10.7000 8.20000 7.40000 5.90000 7.70000 7.80000 ; 8.90000 9.80000 2.10000 3.90000 0.00000 0.00000 ; 6.20000 5.30000 4.40000 4.80000 0.00000 0.00000 print,f_div(fwhm1,fwhm2) ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 1.11364 1.47619 1.30000 0.960784 0.00000 0.00000 ; 3.44828 1.21569 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.205882 0.352273 0.419753 0.00000 0.00000 0.00000 ; 1.46903 0.00000 0.00000 0.00000 0.00000 0.00000 ; 1.23636 0.543210 0.812500 0.00000 0.00000 0.00000 ; 0.325000 0.320442 0.340580 0.00000 0.00000 0.00000 ; 2.27660 2.56250 1.94737 1.01724 1.97436 0.00000 ; 0.247222 1.96000 0.194444 0.975000 0.00000 0.00000 ; 1.82353 1.39474 1.18919 1.41176 0.00000 0.00000 plot,f_div(fwhm1,fwhm2),f_div(flux1,flux2),psym=6,xtit='fwhm1/fwhm2',ytit='flux1/flux2' save,fwhm1,fwhm2,flux1,flux2,'fwhm_flux.sav' ; % SAVE: Expression must be named variable in this context: . save,fwhm1,fwhm2,flux1,flux2,file='fwhm_flux.sav' print,flux1 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 40.3000 49.2000 46.6000 31.4000 0.00000 0.00000 for j=0,12 do print,j,' ',time(j) ; 0 mar 18 ; 1 nov 06 ; 2 nov 19 ; 3 oct 24 ; 4 aug 20 ; 5 feb 20 ; 6 jul 23 ; 7 jun 09 ; 8 jun 10 ; 9 oct 30 ; 10 sep 08 ; 11 sep 20 ; 12 jun 17 flux1[*,2]=[26.7,5.9,8.2,8.2,1.7,0] flux2[*,2]=[8.3,20.2,6.9,6.9,4.2,0] print,flux1 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 26.7000 5.90000 8.20000 8.20000 1.70000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 40.3000 49.2000 46.6000 31.4000 0.00000 0.00000 print,fwhm1 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 4.90000 9.30000 6.50000 4.90000 0.00000 0.00000 ; 20.0000 12.4000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 3.50000 3.10000 3.40000 0.00000 0.00000 0.00000 ; 16.6000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 6.80000 4.40000 13.0000 0.00000 0.00000 0.00000 ; 6.50000 5.80000 4.70000 0.00000 0.00000 0.00000 ; 10.7000 8.20000 7.40000 5.90000 7.70000 7.80000 ; 8.90000 9.80000 2.10000 3.90000 0.00000 0.00000 ; 6.20000 5.30000 4.40000 4.80000 0.00000 0.00000 flux1(*,3)=[0,4.9,9.3,6.5,4.9,0] flux2(*,3)=[0,4.4,6.3,5.0,5.1,0] save,fwhm1,fwhm2,flux1,flux2,file='fwhm_flux.sav' flux1[*,3]=[0,3.7,9.3,6.6,3.4,0] flux2[*,3]=[0,2.2,7.2,6.1,3.1,0] save,fwhm1,fwhm2,flux1,flux2,file='fwhm_flux.sav' print,flux1 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 26.7000 5.90000 8.20000 8.20000 1.70000 0.00000 ; 0.00000 3.70000 9.30000 6.60000 3.40000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 40.3000 49.2000 46.6000 31.4000 0.00000 0.00000 print,flux2 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 8.30000 20.2000 6.90000 6.90000 4.20000 0.00000 ; 0.00000 2.20000 7.20000 6.10000 3.10000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 30.0000 37.5000 37.7000 24.1000 0.00000 0.00000 plot,f_div(fwhm1,fwhm2),f_div(flux1,flux2),psym=6,xtit='fwhm1/fwhm2',ytit='flux1/flux2' flux1[*,4]=[27,4.8,0,0,0,0] flux2[*,4]=[2.3,2.2,0,0,0,0] flux1[*,5]=[5.4,5.3,3.7,1.6,0,0] flux2[*,5]=[38,27.2,11.8,3.5,0,0] save,fwhm1,fwhm2,flux1,flux2,file='fwhm_flux.sav' plot,f_div(fwhm1,fwhm2),f_div(flux1,flux2),psym=6,xtit='fwhm1/fwhm2',ytit='flux1/flux2' print,f_div(fwhm1,fwhm2) ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 1.11364 1.47619 1.30000 0.960784 0.00000 0.00000 ; 3.44828 1.21569 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.205882 0.352273 0.419753 0.00000 0.00000 0.00000 ; 1.46903 0.00000 0.00000 0.00000 0.00000 0.00000 ; 1.23636 0.543210 0.812500 0.00000 0.00000 0.00000 ; 0.325000 0.320442 0.340580 0.00000 0.00000 0.00000 ; 2.27660 2.56250 1.94737 1.01724 1.97436 0.00000 ; 0.247222 1.96000 0.194444 0.975000 0.00000 0.00000 ; 1.82353 1.39474 1.18919 1.41176 0.00000 0.00000 print,f_div(flux1,flux2) ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 3.21687 0.292079 1.18841 1.18841 0.404762 0.00000 ; 0.00000 1.68182 1.29167 1.08197 1.09677 0.00000 ; 11.7391 2.18182 0.00000 0.00000 0.00000 0.00000 ; 0.142105 0.194853 0.313559 0.457143 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 1.34333 1.31200 1.23607 1.30290 0.00000 0.00000 print,flux1(*,4) ; 27.0000 4.80000 0.00000 0.00000 0.00000 0.00000 print,flux2(*,4) ; 2.30000 2.20000 0.00000 0.00000 0.00000 0.00000 flux1[*,4]=[0,4.8,0,0,0,0] flux2[*,4]=[0,2.2,0,0,0,0] save,fwhm1,fwhm2,flux1,flux2,file='fwhm_flux.sav' plot,f_div(fwhm1,fwhm2),f_div(flux1,flux2),psym=6,xtit='fwhm1/fwhm2',ytit='flux1/flux2' for j=0,12 do print,j,' ',time(j) ; 0 mar 18 ; 1 nov 06 ; 2 nov 19 ; 3 oct 24 ; 4 aug 20 ; 5 feb 20 ; 6 jul 23 ; 7 jun 09 ; 8 jun 10 ; 9 oct 30 ; 10 sep 08 ; 11 sep 20 ; 12 jun 17 flux1[*,6]=[0.55,0.87,0.94,0,0,0] flux2[*,6]=[0.84,0.95,1.14,0,0,0] save,fwhm1,fwhm2,flux1,flux2,file='fwhm_flux.sav' plot,f_div(fwhm1,fwhm2),f_div(flux1,flux2),psym=6,xtit='fwhm1/fwhm2',ytit='flux1/flux2' flux1[*,7]=[0.38,0.71,0,0,0,0] flux2[*,7]=[2.69,5.1,0,0,0,0] flux1[*,8]=[0,0,4.1,2.2,2.1,0] flux2[*,8]=[0,0,2.5,2.0,1.3,0] save,fwhm1,fwhm2,flux1,flux2,file='fwhm_flux.sav' plot,f_div(fwhm1,fwhm2),f_div(flux1,flux2),psym=6,xtit='fwhm1/fwhm2',ytit='flux1/flux2' flux1[*,9]=[37,36,33,0,0,0] flux2[*,9]=[62,46,31.,0,0,0] save,fwhm1,fwhm2,flux1,flux2,file='fwhm_flux.sav' plot,f_div(fwhm1,fwhm2),f_div(flux1,flux2),psym=6,xtit='fwhm1/fwhm2',ytit='flux1/flux2' flux1[*,10]=[19,15,13.7,6.5,2.4,0] flux2[*,10]=[4.6,4.7,4.9,3.6,1.2,0] save,fwhm1,fwhm2,flux1,flux2,file='fwhm_flux.sav' plot,f_div(fwhm1,fwhm2),f_div(flux1,flux2),psym=6,xtit='fwhm1/fwhm2',ytit='flux1/flux2' flux1[*,11]=[0,16.9,6.9,6.3,0,0] flux2[*,11]=[0,10.3,4.7,7.1,3.0,0,0] ; % Array subscript for FLUX2 must have same size as source expression. flux2[*,11]=[0,4.7,7.1,3.0,0,0] flux2[*,11]=[0,4.7,7.1,3.0,0] ; % Array subscript for FLUX2 must have same size as source expression. flux2[*,11]=[0,4.7,7.1,3.0,0,0] save,fwhm1,fwhm2,flux1,flux2,file='fwhm_flux.sav' plot,f_div(fwhm1,fwhm2),f_div(flux1,flux2),psym=6,xtit='fwhm1/fwhm2',ytit='flux1/flux2' print,fwhm ; % PRINT: Variable is undefined: FWHM. print,fwhm1 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 4.90000 9.30000 6.50000 4.90000 0.00000 0.00000 ; 20.0000 12.4000 0.00000 0.00000 0.00000 0.00000 ; 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 3.50000 3.10000 3.40000 0.00000 0.00000 0.00000 ; 16.6000 0.00000 0.00000 0.00000 0.00000 0.00000 ; 6.80000 4.40000 13.0000 0.00000 0.00000 0.00000 ; 6.50000 5.80000 4.70000 0.00000 0.00000 0.00000 ; 10.7000 8.20000 7.40000 5.90000 7.70000 7.80000 ; 8.90000 9.80000 2.10000 3.90000 0.00000 0.00000 ; 6.20000 5.30000 4.40000 4.80000 0.00000 0.00000 plot,f_div(fwhm1,fwhm2),f_div(flux1,flux2),psym=6,xtit='FWHM1/FWHM2',ytit='Fx1/flux2' plot,f_div(fwhm1,fwhm2),f_div(flux1,flux2),psym=6,xtit='FWHM1/FWHM2',ytit='FLUX1/FLUX2' plot,f_div(fwhm1,fwhm2),f_div(flux1,flux2),psym=6,xtit='FWHM1/FWHM2',ytit='FLUX1/FLUX2',TIT=TI ti='FLUX RATIO VS FWHM RATIO FOR RHESSI FLARES' plot,f_div(fwhm1,fwhm2),f_div(flux1,flux2),psym=6,xtit='FWHM1/FWHM2',ytit='FLUX1/FLUX2',TIT=TI TT=TVRD() write_png,'flux_ratio_vs_fwhm_ratio.png',255b-tt $xv flux_ratio_vs_fwhm_ratio.png $recent sav