; This script runs the firdem code on synthetic AIA data made from test DEMS composed ; of the sum of several random log-normal distributions, whose statistical properties are ; described by npeaks, emrange, sigrange, and trange, below. Comment out the 'delvar,seed' ; lines to obtain a different set of random peaks with each run of the code. delvar,seed ; Ensures that we obtain the same set of random peaks each time... .compile ../estimate_quantile.pro .compile ../firdem.pro .compile aia_firdem_wrapper.pro .compile ../firdem_get_pixel.pro .compile aia_firdem_test_4.pro set_plot,'ps' !p.font=0 device,file='plots/aia_dem_multimod_test_plots.eps',/encapsulated,/times device,/inches,xsize=8.0,ysize=4.5 ;set_plot,'x' ;ni=1000 ;npeaks=1 ;emrange=[5.0e31,5.0e31] ;sigrange=[4.0,4.0] ;trange=[6.5,6.5] npeaks=5 emrange=[5.0e27,5.0e28] sigrange=[0.05,0.15] trange=[5.75,7.0] ;window,0 tempsin = [5e5,1e6,5e6,1e7] npx=3 npy=3 !p.multi = [0,npx,npy] for i=0,npx*npy-1 do aia_firdem_test_4, npeaks, emrange, sigrange, trange, temps_in=temps_in, $ dem_in=dem_in, temps_out=temps_out, dem_out=dem_out, niter=ni, tr_struct=tr_struct, $ its=its, chi2s=chi2s, seed=seed device,/close set_plot,'x' delvar,seed