; OSPEX script created Tue Jul 18 10:38:17 2006 by OSPEX writescript method. ; ; Call this script with the keyword argument, obj=obj to return the ; OSPEX object reference for use at the command line as well as in the GUI. ; For example: ; ospex_script_18_Jul_2006, obj=obj ; ; Note that this script simply sets parameters in the OSPEX object as they ; were when you wrote the script, and optionally restores fit results. ; To make OSPEX do anything in this script, you need to add some action commands. ; For instance, the command ; obj -> dofit, /all ; would tell OSPEX to do fits in all your fit time intervals. ; See the OSPEX methods section in the OSPEX documentation at ; http://hesperia.gsfc.nasa.gov/ssw/packages/spex/doc/ospex_explanation.htm ; for a complete list of methods and their arguments. ; pro ospex_script_18_Jul_2006, obj=obj if not is_class(obj,'SPEX',/quiet) then obj = ospex() obj-> set, spex_specfile= 'D:/analysis/working/hsi_spectrum_20031113_041638d4.fits' obj-> set, spex_drmfile= 'D:/analysis/working/hsi_srm_20031113_041638.fits' obj-> set, spex_source_angle= 90.0000 obj-> set, spex_source_xy= [-968.000, 8.00000] obj-> set, spex_erange= [3.0000000D, 10.000000D] obj-> set, spex_fit_time_interval= [['13-Nov-2003 04:43:48.000', $ '13-Nov-2003 04:44:48.000'], ['13-Nov-2003 04:44:48.000', '13-Nov-2003 04:45:48.000'], $ ['13-Nov-2003 04:45:48.000', '13-Nov-2003 04:46:48.000'], ['13-Nov-2003 04:46:48.000', $ '13-Nov-2003 04:47:48.000'], ['13-Nov-2003 04:47:48.000', '13-Nov-2003 04:48:48.000'], $ ['13-Nov-2003 04:48:48.000', '13-Nov-2003 04:49:48.000'], ['13-Nov-2003 04:49:48.000', $ '13-Nov-2003 04:50:48.000'], ['13-Nov-2003 04:50:48.000', '13-Nov-2003 04:51:48.000'], $ ['13-Nov-2003 04:51:48.000', '13-Nov-2003 04:52:48.000'], ['13-Nov-2003 04:52:48.000', $ '13-Nov-2003 04:53:48.000'], ['13-Nov-2003 04:53:48.000', '13-Nov-2003 04:54:48.000'], $ ['13-Nov-2003 04:54:48.000', '13-Nov-2003 04:55:48.000'], ['13-Nov-2003 04:55:48.000', $ '13-Nov-2003 04:56:48.000'], ['13-Nov-2003 04:56:48.000', '13-Nov-2003 04:57:48.000'], $ ['13-Nov-2003 04:58:04.000', '13-Nov-2003 04:58:48.000'], ['13-Nov-2003 04:58:48.000', $ '13-Nov-2003 04:59:48.000'], ['13-Nov-2003 04:59:48.000', '13-Nov-2003 05:00:48.000'], $ ['13-Nov-2003 05:00:48.000', '13-Nov-2003 05:01:48.000'], ['13-Nov-2003 05:01:48.000', $ '13-Nov-2003 05:02:48.000'], ['13-Nov-2003 05:02:48.000', '13-Nov-2003 05:03:48.000'], $ ['13-Nov-2003 05:03:48.000', '13-Nov-2003 05:04:28.000'], ['13-Nov-2003 05:05:08.000', $ '13-Nov-2003 05:05:48.000'], ['13-Nov-2003 05:05:48.000', '13-Nov-2003 05:06:48.000'], $ ['13-Nov-2003 05:06:48.000', '13-Nov-2003 05:07:48.000'], ['13-Nov-2003 05:07:48.000', $ '13-Nov-2003 05:08:52.000'], ['13-Nov-2003 05:09:12.000', '13-Nov-2003 05:09:48.000'], $ ['13-Nov-2003 05:09:48.000', '13-Nov-2003 05:10:48.000'], ['13-Nov-2003 05:10:48.000', $ '13-Nov-2003 05:11:48.000'], ['13-Nov-2003 05:11:48.000', '13-Nov-2003 05:12:48.000'], $ ['13-Nov-2003 05:12:48.000', '13-Nov-2003 05:13:48.000'], ['13-Nov-2003 05:13:48.000', $ '13-Nov-2003 05:14:48.000'], ['13-Nov-2003 05:14:48.000', '13-Nov-2003 05:15:48.000'], $ ['13-Nov-2003 05:15:48.000', '13-Nov-2003 05:16:48.000'], ['13-Nov-2003 05:16:48.000', $ '13-Nov-2003 05:17:48.000'], ['13-Nov-2003 05:17:48.000', '13-Nov-2003 05:18:48.000'], $ ['13-Nov-2003 05:18:48.000', '13-Nov-2003 05:19:48.000'], ['13-Nov-2003 05:19:48.000', $ '13-Nov-2003 05:20:48.000'], ['13-Nov-2003 05:20:48.000', '13-Nov-2003 05:21:48.000'], $ ['13-Nov-2003 05:21:48.000', '13-Nov-2003 05:22:48.000']] obj-> set, spex_bk_sep= 1 obj-> set, spex_bk_order=[0, 0, 3, 3, 3, 3, 0] obj-> set, spex_bk_eband=[[3.0000000D, 6.0000000D], [6.0000000D, 12.000000D], [12.000000D, $ 25.000000D], [25.000000D, 50.000000D], [50.000000D, 100.00000D], [100.00000D, $ 300.00000D]] obj-> set, this_band=0, this_time=['13-Nov-2003 05:25:00.000', '13-Nov-2003 05:27:40.000'] obj-> set, this_band=1, this_time=['13-Nov-2003 05:26:00.000', '13-Nov-2003 05:27:40.000'] obj-> set, this_band=2, this_time=[['13-Nov-2003 04:45:04.000', $ '13-Nov-2003 04:52:08.000'], ['13-Nov-2003 05:16:08.000', '13-Nov-2003 05:28:20.000']] obj-> set, this_band=3, this_time=[['13-Nov-2003 04:48:44.000', $ '13-Nov-2003 04:54:04.000'], ['13-Nov-2003 05:08:52.000', '13-Nov-2003 05:27:48.000']] obj-> set, this_band=4, this_time=[['13-Nov-2003 04:47:20.000', $ '13-Nov-2003 04:55:04.000'], ['13-Nov-2003 05:10:08.000', '13-Nov-2003 05:26:56.000']] obj-> set, this_band=5, this_time=['13-Nov-2003 04:45:04.000', '13-Nov-2003 05:26:40.000'] obj-> set, spex_uncert= 0.000000 obj-> set, fit_function= 'vth+multi_therm_pow+3pow+drm_mod' obj-> set, fit_comp_params= [0.000100000, 2.00000, 1.00000, 0.000000, 0.500000, 4.00000, $ 1.00000, 1.00000, 0.000000, 3.00000, 30.0000, 4.00000, 400.000, 2.00000, 1.00000, $ 0.000000, 1.00000] obj-> set, fit_comp_minima= [1.00000e-020, 0.500000, 0.0100000, 1.00000e-010, 0.500000, $ 1.00000, 0.0100000, 0.0100000, 1.00000e-020, 0.500000, 1.00000, 0.500000, 10.0000, $ 0.500000, 0.500000, -1.00000, 0.0500000] obj-> set, fit_comp_maxima= [1.00000e+020, 50.0000, 10.0000, 100.000, 4.00000, 10.0000, $ 100.000, 10.0000, 1.00000e+020, 12.0000, 20000.0, 12.0000, 20000.0, 12.0000, 1.50000, $ 1.00000, 1.50000] obj-> set, fit_comp_free_mask= [1B, 1B, 0B, 0B, 0B, 0B, 0B, 0B, 0B, 0B, 0B, 0B, 0B, 0B, 0B, $ 0B, 0B] obj-> set, fit_comp_spectrum= ['full', 'full', '', ''] obj-> set, fit_comp_model= ['chianti', 'chianti', '', ''] obj-> set, spex_autoplot_units= 'Flux' obj-> set, spex_eband= [[3.00000, 6.00000], [6.00000, 12.0000], [12.0000, 25.0000], $ [25.0000, 50.0000], [50.0000, 100.000], [100.000, 300.000], [300.000, 800.000]] ;obj -> restorefit, file='D:/analysis/working/ospex_results_18_Jul_2006.fits' end