Документ взят из кэша поисковой машины. Адрес оригинального документа : http://hea-www.harvard.edu/~saku/modmapcookbookmod.pdf Дата изменения: Tue Apr 10 23:21:56 2012 Дата индексирования: Mon Oct 1 21:54:40 2012 Кодировка: Поисковые слова: п п п п п п п п п п п п п п п п п п п п п п п п п п п п п п п п п п п п п п

Steps for running Modmap at the CfA
In the following computer output is in small type, prompts end with ":" or ">", comments are in boldface (with command names in quotes), and that which the types is in italics: Either of the following should bring up Modmap: /home/dsteeghs/Software/bin/modmap /data/saku10/software/modmap/bin/modmap Following which you should see this start-up text:
======================================== MODMAP running on keiko Tue Oct 5 15:48:24 2010 Version 1.0 Danny Steeghs ======================================== modmap> Commands: C ... Conventional Doppler mapping V ... Setup velocity grid/rest F ... Fake data nw: 100 nphi: 80 D ... Load molly Data O ... Generate optimal starting image P ... Plot fit R ... Annotate plots with lobe M ... MEM iterations S ... Save maps/data to disk file L ... Load maps from disk file Q ... Quit

To load a molly file (.mol file) into modmap we use the command "d" and enter the filename. Generally the file should consist of a set of spectra trimmed to capture the line of interest (and including a reasonable amount of continuum on each side).
modmap> d MOLLY file to load ? [empty]test.mol Rest wavelength:[0= 6562.76] 6562.76 MOLLY File opened : test.mol Error reading MOLLY headers 43 spectra read. DEBUG: 36.461834604179366 251 0. 0. 6529.491 Data is uniform in velocity Covering 6431.364 - 6629.919 Angstrom Data covers phases 0.73045653 - 1.5451138 *** Warning, there were 1320 ( 29.803568 %) negative pixels! Readmolly done Found 103 pixels between - 1868.9 and 1868.9 km/s Velocity dispersion of data : 36.3573 km/s


Velocity resolution of tomogram : 34.32323 FWHM of local profile : 60. Filled array OBSN with real data... 103 x 43 Filled array SIGMA with real errors...oc

The command "v" is used to setup the parameters of the system. If unknown, use your best guess for (radial velocity); you can check this using different values in modmap. (250 pixels at 0.8angstroms/pixel; better to oversample)
modmap> v Rest wavelength (A) : 6562.75788 Gamma velocity : -44 Number of pixels in velocity grid : 100 Maximum velocity in tomogram: 2000 FWHM of local line profile: 100 Velocity grid set to 100 x 100 pixels Spanning +/- 2000. @ 40.40404 km/s/pixel Origin corresponds to 6562.758 A + -44. km/s The rest wavelength of the feature to be used The system's radial velocity

Next, we use the command "o" to generate the optimal starting image:
modmap> o Performing optical scaling of starting image =============== op0 === 5 0. pregauss: 50000 0.00008499622 scale factor = 0.15342808

This optimises Chi squared starting on the mean map (sets pixel level). The command "m" gets us into 'MEMSYS', to fit the data.
modmap> m Ready for MEMSYS iterations Velocity grid = 100 x 100 Data grid = 103 x 43 Using 5 images to map modulations image pixels= 50000 data values= 4429

MEMCORE - configuring MEM data and image buffers. Image dimensions (X,Y,Z)= 50000 1 1 Data dimension = 4429 9 images in 9 files. 1 50001 100001 150001 200001 250001 -9999999 -9999999 -9999999 -9999999 -999999 9 -9999999 -9999999 -9999999 -9999999 -9999999 -9999999 400001 300001 350001 8 data vectors in 8 files.


450001 454430 458859 463288 467717 472146 476575 481004 -9999999 -9999999 -9999 999 -9999999 -9999999 -9999999 -9999999 -9999999 -9999999 -9999999 -9999999 -999 9999 Image pixels: 450000 = 92.700935 % Data values: 35432 = 7.2990656 % Total storage: 485432 = 17.978964 % of buffer: 2700000 Divisor for weights: 0.9 values: 0.0001 0.15342808 defaults: 0.001 0.16342808 weights: 0.5 1. =============== op0 === 5 0. Chi^2/ 4429 = 24.33169 =============== fixdef === ------------------------------------------------Iter= 0 CPU= 0. CPU/Iter= 0. Elapsed realtime= 0. Mode= 1030 Level= 20 Ndat= 4429 Nimg= 50000 SumF= 0. Caim= 1. Chi^2/N= 24.33169 Sold= 0. Cold= 0. Snew= 0. Cnew= 0. Alp= 0. -2AlpS= -0. Stru= -0.371579 Nfit= 0. Rmax= 0.2 Step= 0. TEST= 0. Take 1 iterations on each GO. Update default every -1 iterations. MEMITER >

This subroutine uses iterations to reduce the value of chi squared to extremely unlikely that system will ever get a reduced Chi square of say 30 (in the example above it is 24.33), start by aiming for 20 - if it steadily reduce to 10 and so on. To set final Chi square value to aim ask for an initial guess and number of iterations:
MEMITER > c == Auto C-aim stepper === test version Enter final C-aim 8 Enter test value to go ahead with 2.0

try and find a good fit. It is 1; if the initial Chi square is looks like it is improving, for, use "c". This will also


Enter number of iterations per step 2 Enter number of steps to reach final CAIM 2 New C-aim ; 8. 1

** MEM iteration: 3 ******************* TO GO 2 =============== fixdef === CAIM <- 8. =============== op0 === 5 0. C === 2.08283E+01 =============== tr1 === 5 0. S === 1.62363E-04 TEST === 0.47519 =============== op0 === 5 0. =============== op0 === 5 0. SNEW === -2.09627E-02 dist === 0.20000 CNEW === 1.87210E+01 =============== tr1 === 5 0. =============== op0 === 5 0. SNEW === -2.09614E-02 dist === 0.20000 CNEW === 1.87203E+01 RMAX <- 0.2 =============== op0 === 5 0. Chi^2/ 4429 = 18.71996 GLOBAL |C/CAIM-1|= 1.3399949 < 0.01 FAILED. LOCAL |C/CAIM-1|= 0. < 1. OK. TEST = 0.4751914 < 0.00007 FAILED.

** MEM iteration: 4 ******************* TO GO 1 =============== fixdef === CAIM <- 8. =============== op0 === 5 0. C === 1.87200E+01 =============== tr1 === 5 0. S === 9.66787E-05 TEST === 0.47423 =============== op0 === 5 0. =============== op0 === 5 0. SNEW === -2.16063E-02 dist === 0.20000 CNEW === 1.65007E+01 =============== tr1 === 5 0. =============== op0 === 5 0.


SNEW === -2.16031E-02 dist === 0.20000 CNEW === 1.64995E+01 RMAX <- 0.2 =============== op0 === 5 0. Chi^2/ 4429 = 16.499258 GLOBAL |C/CAIM-1|= 1.0624073 < 0.01 FAILED. LOCAL |C/CAIM-1|= 0. < 1. OK. TEST = 0.47422916 < 0.00007 FAILED. ** MEM iterations completed ********* Updating CAIM @@@@@@@@ Sorry, lost it!! ------------------------------------------------Iter= 4 CPU= 26.059038 CPU/Iter= 6.5147595 Elapsed realtime= 26.048044 Mode= 1030 Level= 20 Ndat= 4429 Nimg= 50000 SumF= 952.6445 Caim= 8. Chi^2/N= 16.499258 Sold= 0.000096678734 Cold= 18.719967 Snew= -0.021603143 Cnew= 16.49951 Alp= 0.053443063 -2AlpS= 0.0023090763 Stru= -0.371579 Nfit= 0.01985832 Rmax= 0.2 Step= 0.2 TEST= 0.47422916 Take 2 iterations on each GO. Update default every -1 iterations. MEMITER > ------------------------------------------------Iter= 14 CPU= 89.01747 CPU/Iter= 6.358391 Elapsed realtime= 88.9855 Mode= 1030 Level= 20 Ndat= 4429 Nimg= 50000 SumF= 996.9127 Caim= 8. Chi^2/N= 8.5955305 Sold= -0.037704766 Cold= 8.700225 Snew= -0.06027075 Cnew= 8.595371 Alp= 0.0025477018 -2AlpS= 0.0003071038 Stru= -0.371579 Nfit= 0.00094667246 Rmax= 0.2 Step= 0.2 TEST= 0.21966195 Take 5 iterations on each GO. Update default every -1 iterations. MEMITER > -------------------------------------------------


"g" does one iteration. You can also use "i 5", for example, to set "g" do (5) iterations. The "g" command will then start the iterations. "p" shows the Chi square convergence; lower values should be approached in a straight line; if values go up and down sporadically, you are trying for too low a Chi square value. In this case, you will have to restart modmap and stay within a more reasonable limit. Once you have reached the desired reduced Chi square, you will still need to iterate to minimize Maximum Entropy. When plotting, modmap will ask you which plotting device you would like to use. To open the graph in a window, choose /xs (xserve). To open multiple windows simultaneously, use successive numbers followed by /xs, e.g. 1/xs, 2/xs etc, each new number opening a new window. To save a plot, type the filename followed by /ps, e.g. file.ps/ps.
MEMITER > p Iterations: 14 1 14 Graphics device/type (? to see list, default /NULL): /xs ------------------------------------------------Iter= 14 CPU= 89.01747 CPU/Iter= 6.358391 Elapsed realtime= 88.9855 Mode= 1030 Level= 20 Ndat= 4429 Nimg= 50000 SumF= 996.9127 Caim= 8. Chi^2/N= 8.5955305 Sold= -0.037704766 Cold= 8.700225 Snew= -0.06027075 Cnew= 8.595371 Alp= 0.0025477018 -2AlpS= 0.0003071038 Stru= -0.371579 Nfit= 0.00094667246 Rmax= 0.2 Step= 0.2 TEST= 0.21966195 Take 5 iterations on each GO. Update default every -1 iterations.p

The command"f" plots the modulated maps, along with the trail and the predicted trail and the sine cosine variations:
MEMITER > f =============== plotfit === =============== op0 === 5 0. Graphics device/type (? to see list, default /NULL): /xs ------------------------------------------------Iter= 14 CPU= 89.01747 CPU/Iter= 6.358391 Elapsed realtime= 88.9855 Mode= 1030 Level= 20 Ndat= 4429 Nimg= 50000


SumF= 996.9127 Caim= 8. Chi^2/N= 8.5955305 Sold= -0.037704766 Cold= 8.700225 Snew= -0.06027075 Cnew= 8.595371 Alp= 0.0025477018 -2AlpS= 0.0003071038 Stru= -0.371579 Nfit= 0.00094667246 Rmax= 0.2 Step= 0.2 TEST= 0.21966195 Take 5 iterations on each GO. Update default every -1 iterations. MEMITER >

To save the data produced, use the command "s". This will ask you whether you wish to save the data or the maps. Choose using LOWER CASE 'd' or 'm' to save the data or maps respectively, then enter the filename to save: modmap> s Dump (D)ata or (M)aps? d =============== dumpdata === File to SAVE data in : file.dat data dumped...

Type name.ps/ps to save images onto a file name.ps.