; IDL Version 8.1, Mac OS X (darwin x86_64 m64)
; Journal File for Bob@steinr.pa.msu.edu
; Working directory: /Volumes/d2/MHDV/mhd12v/Stokes
; Date: Sat Nov  5 00:52:41 2011
 
journ
f=read_prf('mhd12_6km_342.6.prf')
;key:          -1
;nregion:           1
;nw:         200
;nprf, nx, ny =         200        1008        1008
v=reform(f.prf[*,3,*,*])
vhin=v
vatst=v
for n=0,199 do vhin[n,*,*]=seeing(reform(v[n,*,*],size=12,ap=50,ora=0.344)
; % Syntax error.
for n=0,199 do vhin[n,*,*]=seeing(reform(v[n,*,*]),size=12,ap=50,ora=0.344)
;  0.000374446  0.000287484
;  0.000397101  0.000304719
;  0.000422971  0.000324367
;  0.000449781  0.000344701
;  0.000480452  0.000367940
;  0.000512375  0.000392106
;  0.000549069  0.000419827
;  0.000587409  0.000448749
;  0.000631689  0.000482100
;  0.000678183  0.000517067
;  0.000732141  0.000557549
;  0.000789110  0.000600208
;  0.000855616  0.000649883
;  0.000926281  0.000702509
;   0.00100936  0.000764177
;   0.00109830  0.000829940
;   0.00120392  0.000907527
;   0.00131804  0.000990770
;   0.00145513   0.00108968
;   0.00160543   0.00119663
;   0.00178768   0.00132455
;   0.00198779   0.00146330
;   0.00222963   0.00162908
;   0.00249492   0.00180938
;   0.00282447   0.00202959
;   0.00319725   0.00227339
;   0.00364133   0.00255803
;   0.00421704   0.00291380
;   0.00490971   0.00332069
;   0.00585855   0.00384624
;   0.00708923   0.00448029
;   0.00887604   0.00534699
;    0.0112202   0.00645188
;    0.0145606   0.00803734
;    0.0186824    0.0101055
;    0.0240744    0.0130480
;    0.0303400    0.0167824
;    0.0382431    0.0218832
;    0.0469745    0.0279210
;    0.0569301    0.0352287
;    0.0662986    0.0424767
;    0.0749127    0.0494385
;    0.0809516    0.0544898
;    0.0842901    0.0572912
;    0.0845420    0.0572057
;    0.0819591    0.0543569
;    0.0769774    0.0494113
;    0.0694206    0.0421872
;    0.0606263    0.0334723
;    0.0512330    0.0229443
;    0.0453981    0.0149468
;    0.0456098    0.0165263
;    0.0511320    0.0250690
;    0.0583605    0.0334469
;    0.0658130    0.0409165
;    0.0716061    0.0463248
;    0.0755479    0.0498474
;    0.0769627    0.0510419
;    0.0760849    0.0502019
;    0.0733123    0.0477234
;    0.0688163    0.0437880
;    0.0635014    0.0392163
;    0.0573789    0.0340572
;    0.0514704    0.0292392
;    0.0455020    0.0246040
;    0.0401948    0.0207350
;    0.0351013    0.0172725
;    0.0307393    0.0145060
;    0.0266469    0.0120821
;    0.0231641    0.0101545
;    0.0199109   0.00845938
;    0.0171774   0.00710712
;    0.0146503   0.00591772
;    0.0125228   0.00497113
;    0.0105491   0.00414580
;   0.00889981   0.00349731
;   0.00746163   0.00296369
;   0.00616302   0.00250632
;   0.00510675   0.00214974
;   0.00416571   0.00184199
;   0.00341792   0.00160220
;   0.00278273   0.00139637
;   0.00230925   0.00123626
;   0.00192349   0.00109745
;   0.00163523  0.000985942
;   0.00139681  0.000885132
;   0.00121698  0.000801299
;   0.00106381  0.000723507
;  0.000943290  0.000657454
;  0.000838017  0.000594988
;  0.000753261  0.000541115
;  0.000676256  0.000489578
;  0.000611613  0.000444597
;  0.000550720  0.000401068
;  0.000498047  0.000362673
;  0.000447259  0.000325097
;  0.000402466  0.000291553
;  0.000358619  0.000258311
;  0.000319520  0.000228284
;  0.000280869  0.000198157
;  0.000246210  0.000170580
;  0.000213571  0.000143894
;  0.000181651  0.000116668
;  0.000154043  9.14763e-05
;  0.000129288  6.60774e-05
;  0.000112140  4.43098e-05
;  0.000104577  3.07647e-05
;  0.000110003  3.86674e-05
;  0.000128332  6.18204e-05
;  0.000155037  8.90466e-05
;  0.000190228  0.000120711
;  0.000229305  0.000153486
;  0.000275828  0.000190867
;  0.000325639  0.000229737
;  0.000384582  0.000274608
;  0.000448769  0.000321978
;  0.000526797  0.000377612
;  0.000613157  0.000437378
;  0.000718095  0.000508821
;  0.000834783  0.000587043
;  0.000979263  0.000682447
;   0.00114313  0.000789212
;   0.00135244  0.000922702
;   0.00160195   0.00107602
;   0.00193787   0.00127233
;   0.00235056   0.00150326
;   0.00288749   0.00179041
;   0.00364938   0.00217892
;   0.00466005   0.00266778
;   0.00615940   0.00336730
;   0.00814352   0.00430077
;    0.0110139   0.00569793
;    0.0146860   0.00759055
;    0.0196579    0.0103636
;    0.0255475    0.0139365
;    0.0330173    0.0188254
;    0.0413124    0.0246032
;    0.0509753    0.0316813
;    0.0605685    0.0390212
;    0.0703327    0.0468143
;    0.0786298    0.0537168
;    0.0855025    0.0596476
;    0.0895906    0.0632032
;    0.0907678    0.0638896
;    0.0887624    0.0611962
;    0.0836895    0.0550972
;    0.0767372    0.0470504
;    0.0681834    0.0373373
;    0.0600770    0.0277600
;    0.0536776    0.0189877
;    0.0517830    0.0159108
;    0.0547501    0.0210080
;    0.0616197    0.0305954
;    0.0694351    0.0401378
;    0.0767677    0.0485792
;    0.0816145    0.0540743
;    0.0839402    0.0569055
;    0.0835641    0.0570150
;    0.0808020    0.0548835
;    0.0763956    0.0512416
;    0.0704428    0.0462577
;    0.0639845    0.0408422
;    0.0568509    0.0348897
;    0.0501331    0.0293709
;    0.0434980    0.0240948
;    0.0377498    0.0197520
;    0.0323606    0.0159446
;    0.0278240    0.0129739
;    0.0236458    0.0104402
;    0.0201728   0.00848634
;    0.0170074   0.00682437
;    0.0144139   0.00554649
;    0.0120902   0.00446626
;    0.0102008   0.00364327
;   0.00850303   0.00295264
;   0.00711936   0.00242726
;   0.00593282   0.00200679
;   0.00487154   0.00165425
;   0.00401601   0.00138601
;   0.00326129   0.00115984
;   0.00265983  0.000986036
;   0.00213326  0.000838272
;   0.00172305  0.000725507
;   0.00138034  0.000631592
;   0.00112573  0.000560528
;  0.000918991  0.000500065
;  0.000768957  0.000452426
;  0.000649562  0.000410027
;  0.000562123  0.000375272
;  0.000491268  0.000343432
;  0.000439477  0.000316733
;  0.000397136  0.000291856
;  0.000364326  0.000270721
;  0.000335029  0.000250815
;  0.000310681  0.000233741
;  0.000288027  0.000217536
;  0.000268776  0.000203531
;  0.000250591  0.000190152
;  0.000234917  0.000178532
;  0.000219924  0.000167384
for n=0,199 do vatst[n,*,*]=seeing(reform(v[n,*,*]),size=12,ap=400)
;  0.000374446  0.000357281
;  0.000397101  0.000378842
;  0.000422971  0.000403428
;  0.000449781  0.000428897
;  0.000480452  0.000458023
;  0.000512375  0.000488335
;  0.000549069  0.000523138
;  0.000587409  0.000559484
;  0.000631689  0.000601440
;  0.000678183  0.000645469
;  0.000732141  0.000696516
;  0.000789110  0.000750375
;  0.000855616  0.000813180
;  0.000926281  0.000879840
;   0.00100936  0.000958087
;   0.00109830   0.00104170
;   0.00120392   0.00114062
;   0.00131804   0.00124705
;   0.00145513   0.00137399
;   0.00160543   0.00151191
;   0.00178768   0.00167796
;   0.00198779   0.00185934
;   0.00222963   0.00207809
;   0.00249492   0.00231798
;   0.00282447   0.00261276
;   0.00319725   0.00294294
;   0.00364133   0.00333411
;   0.00421704   0.00383292
;   0.00490971   0.00441954
;   0.00585855   0.00520360
;   0.00708923   0.00619726
;   0.00887604   0.00763437
;    0.0112202   0.00954812
;    0.0145606    0.0123306
;    0.0186824    0.0158458
;    0.0240744    0.0205461
;    0.0303400    0.0261247
;    0.0382431    0.0333174
;    0.0469745    0.0414434
;    0.0569301    0.0508910
;    0.0662986    0.0599204
;    0.0749127    0.0683256
;    0.0809516    0.0742753
;    0.0842901    0.0775855
;    0.0845420    0.0777876
;    0.0819591    0.0750923
;    0.0769774    0.0699589
;    0.0694206    0.0621693
;    0.0606263    0.0529587
;    0.0512330    0.0427997
;    0.0453981    0.0362817
;    0.0456098    0.0367784
;    0.0511320    0.0432737
;    0.0583605    0.0512281
;    0.0658130    0.0590881
;    0.0716061    0.0650342
;    0.0755479    0.0690047
;    0.0769627    0.0703977
;    0.0760849    0.0694790
;    0.0733123    0.0666663
;    0.0688163    0.0621457
;    0.0635014    0.0568419
;    0.0573789    0.0507834
;    0.0514704    0.0450019
;    0.0455020    0.0392429
;    0.0401948    0.0342124
;    0.0351013    0.0294850
;    0.0307393    0.0255219
;    0.0266469    0.0218810
;    0.0231641    0.0188460
;    0.0199109    0.0160552
;    0.0171774    0.0137312
;    0.0146503    0.0115961
;    0.0125228   0.00981894
;    0.0105491   0.00820010
;   0.00889981   0.00687596
;   0.00746163   0.00574418
;   0.00616302   0.00473899
;   0.00510675   0.00393121
;   0.00416571   0.00321761
;   0.00341792   0.00265455
;   0.00278273   0.00218097
;   0.00230925   0.00183486
;   0.00192349   0.00156054
;   0.00163523   0.00135903
;   0.00139681   0.00119079
;   0.00121698   0.00106043
;   0.00106381  0.000945969
;  0.000943290  0.000852686
;  0.000838017  0.000767160
;  0.000753261  0.000695203
;  0.000676256  0.000627565
;  0.000611613  0.000569273
;  0.000550720  0.000513372
;  0.000498047  0.000464435
;  0.000447259  0.000416848
;  0.000402466  0.000374630
;  0.000358619  0.000333096
;  0.000319520  0.000295871
;  0.000280869  0.000258869
;  0.000246210  0.000225433
;  0.000213571  0.000193652
;  0.000181651  0.000162114
;  0.000154043  0.000134256
;  0.000129288  0.000108466
;  0.000112140  8.97443e-05
;  0.000104577  8.10252e-05
;  0.000110003  8.70114e-05
;  0.000128332  0.000106878
;  0.000155037  0.000134605
;  0.000190228  0.000169945
;  0.000229305  0.000208348
;  0.000275828  0.000253392
;  0.000325639  0.000301067
;  0.000384582  0.000356811
;  0.000448769  0.000416315
;  0.000526797  0.000486915
;  0.000613157  0.000563587
;  0.000718095  0.000656219
;  0.000834783  0.000758865
;  0.000979263  0.000885576
;   0.00114313   0.00102899
;   0.00135244   0.00121051
;   0.00160195   0.00142308
;   0.00193787   0.00170307
;   0.00235056   0.00204384
;   0.00288749   0.00248152
;   0.00364938   0.00309179
;   0.00466005   0.00388876
;   0.00615940   0.00507379
;   0.00814352   0.00668443
;    0.0110139   0.00908269
;    0.0146860    0.0122320
;    0.0196579    0.0166016
;    0.0255475    0.0218919
;    0.0330173    0.0287353
;    0.0413124    0.0364740
;    0.0509753    0.0456381
;    0.0605685    0.0548664
;    0.0703327    0.0643757
;    0.0786298    0.0725443
;    0.0855025    0.0793695
;    0.0895906    0.0834402
;    0.0907678    0.0845379
;    0.0887624    0.0823159
;    0.0836895    0.0768261
;    0.0767372    0.0693069
;    0.0681834    0.0599806
;    0.0600770    0.0509368
;    0.0536776    0.0435041
;    0.0517830    0.0411713
;    0.0547501    0.0446245
;    0.0616197    0.0524936
;    0.0694351    0.0611944
;    0.0767677    0.0692124
;    0.0816145    0.0745057
;    0.0839402    0.0771574
;    0.0835641    0.0770005
;    0.0808020    0.0744044
;    0.0763956    0.0701163
;    0.0704428    0.0642579
;    0.0639845    0.0578728
;    0.0568509    0.0508217
;    0.0501331    0.0442164
;    0.0434980    0.0377606
;    0.0377498    0.0322525
;    0.0323606    0.0271899
;    0.0278240    0.0230224
;    0.0236458    0.0192661
;    0.0201728    0.0162058
;    0.0170074    0.0134653
;    0.0144139    0.0112533
;    0.0120902   0.00929755
;    0.0102008   0.00773625
;   0.00850303   0.00636420
;   0.00711936   0.00527130
;   0.00593282   0.00435494
;   0.00487154   0.00355325
;   0.00401601   0.00292003
;   0.00326129   0.00236976
;   0.00265983   0.00193422
;   0.00213326   0.00155305
;   0.00172305   0.00125772
;   0.00138034   0.00101704
;   0.00112573  0.000846148
;  0.000918991  0.000713243
;  0.000768957  0.000618959
;  0.000649562  0.000543889
;  0.000562123  0.000488005
;  0.000491268  0.000440383
;  0.000439477  0.000402587
;  0.000397136  0.000368803
;  0.000364326  0.000340869
;  0.000335029  0.000314991
;  0.000310681  0.000293003
;  0.000288027  0.000272265
;  0.000268776  0.000254427
;  0.000250591  0.000237446
;  0.000234917  0.000222739
;  0.000219924  0.000208656
print,ib[1]
;(      17.0000,      129.000)
i=1
help,w
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
; % Syntax error.
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=5
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
ps
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
psend
x
$mv idl.ps V_sim_atst_hinode_12v_17-129.ps
i=211
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
ps
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
psend
x
$mv idl.ps V_sim_atst_hinode_12v_909-422.ps
print,i,ib[i]
;     211(      891.000,      679.000)
help,ib
for i=0,258 do begin plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1  & print,i,ib[i] & wait,2
;       0(      14.0000,      131.000)
;       1(      17.0000,      129.000)
;       2(      26.0000,      213.000)
;       3(      26.0000,      214.000)
;       4(      42.0000,      241.000)
;       5(      43.0000,      241.000)
;       6(      44.0000,      344.000)
;       7(      45.0000,      345.000)
;       8(      45.0000,      351.000)
for i=0,258 do begin plot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1  & print,i,ib[i] & wait,1
;       0(      14.0000,      131.000)
;       1(      17.0000,      129.000)
;       2(      26.0000,      213.000)
;       3(      26.0000,      214.000)
;       4(      42.0000,      241.000)
;       5(      43.0000,      241.000)
;       6(      44.0000,      344.000)
;       7(      45.0000,      345.000)
;       8(      45.0000,      351.000)
;       9(      46.0000,      346.000)
;      10(      64.0000,      309.000)
;      11(      68.0000,      780.000)
;      12(      94.0000,      34.0000)
;      13(      94.0000,      35.0000)
;      14(      102.000,      761.000)
;      15(      103.000,      760.000)
;      16(      110.000,      761.000)
;      17(      110.000,      762.000)
;      18(      111.000,      760.000)
;      19(      111.000,      761.000)
;      20(      114.000,      576.000)
;      21(      115.000,      576.000)
;      22(      116.000,      576.000)
;      23(      120.000,      407.000)
;      24(      127.000,      341.000)
;      25(      129.000,      345.000)
;      26(      129.000,      346.000)
;      27(      129.000,      348.000)
;      28(      133.000,      626.000)
;      29(      134.000,      625.000)
;      30(      135.000,      361.000)
;      31(      141.000,      624.000)
;      32(      186.000,      983.000)
;      33(      188.000,      259.000)
;      34(      189.000,      259.000)
;      35(      192.000,      973.000)
;      36(      223.000,      932.000)
;      37(      225.000,      931.000)
;      38(      226.000,      930.000)
;      39(      271.000,      599.000)
;      40(      272.000,      599.000)
;      41(      286.000,      612.000)
;      42(      299.000,      291.000)
;      43(      307.000,      275.000)
;      44(      311.000,      973.000)
;      45(      311.000,      974.000)
;      46(      312.000,      973.000)
;      47(      312.000,      974.000)
;      48(      318.000,      949.000)
;      49(      318.000,      950.000)
;      50(      323.000,      849.000)
;      51(      324.000,      848.000)
;      52(      324.000,      849.000)
;      53(      338.000,      272.000)
;      54(      354.000,      154.000)
;      55(      354.000,      161.000)
;      56(      356.000,      147.000)
;      57(      365.000,      280.000)
;      58(      383.000,      69.0000)
;      59(      384.000,      71.0000)
;      60(      384.000,      72.0000)
;      61(      384.000,      73.0000)
;      62(      385.000,      210.000)
;      63(      393.000,      833.000)
;      64(      394.000,      834.000)
;      65(      398.000,      191.000)
;      66(      421.000,      759.000)
;      67(      421.000,      760.000)
;      68(      428.000,      725.000)
;      69(      430.000,      724.000)
;      70(      431.000,      724.000)
;      71(      432.000,      724.000)
flush,!journal
i=4
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=5
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=18
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=50
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
for i=50,258 do begin plot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1  & print,i,ib[i] & wait,1
;      50(      323.000,      849.000)
;      51(      324.000,      848.000)
;      52(      324.000,      849.000)
;      53(      338.000,      272.000)
;      54(      354.000,      154.000)
;      55(      354.000,      161.000)
;      56(      356.000,      147.000)
;      57(      365.000,      280.000)
;      58(      383.000,      69.0000)
;      59(      384.000,      71.0000)
;      60(      384.000,      72.0000)
;      61(      384.000,      73.0000)
;      62(      385.000,      210.000)
;      63(      393.000,      833.000)
;      64(      394.000,      834.000)
;      65(      398.000,      191.000)
;      66(      421.000,      759.000)
;      67(      421.000,      760.000)
;      68(      428.000,      725.000)
;      69(      430.000,      724.000)
;      70(      431.000,      724.000)
;      71(      432.000,      724.000)
;      72(      439.000,      707.000)
;      73(      442.000,      304.000)
;      74(      445.000,      711.000)
;      75(      448.000,      246.000)
;      76(      449.000,      246.000)
;      77(      461.000,      32.0000)
;      78(      462.000,      32.0000)
;      79(      463.000,      32.0000)
;      80(      464.000,      31.0000)
;      81(      468.000,      579.000)
;      82(      489.000,      29.0000)
;      83(      490.000,      696.000)
;      84(      497.000,      1003.00)
;      85(      500.000,      1003.00)
;      86(      527.000,      1004.00)
;      87(      527.000,      1005.00)
;      88(      528.000,      483.000)
;      89(      528.000,      1003.00)
;      90(      529.000,      1003.00)
;      91(      529.000,      1004.00)
;      92(      530.000,      1001.00)
;      93(      530.000,      1002.00)
;      94(      546.000,      892.000)
;      95(      551.000,      543.000)
;      96(      552.000,      538.000)
;      97(      552.000,      539.000)
;      98(      552.000,      540.000)
;      99(      552.000,      541.000)
;     100(      555.000,      908.000)
;     101(      556.000,      535.000)
;     102(      556.000,      536.000)
;     103(      556.000,      909.000)
;     104(      557.000,      536.000)
;     105(      557.000,      537.000)
;     106(      557.000,      538.000)
;     107(      559.000,      534.000)
;     108(      559.000,      535.000)
;     109(      560.000,      533.000)
;     110(      560.000,      534.000)
;     111(      560.000,      535.000)
i=104
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=105
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=103
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=95
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=94
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=75
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=52
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
for i=100,258 do begin plot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1  & print,i,ib[i] & wait,1.5
;     100(      555.000,      908.000)
;     101(      556.000,      535.000)
;     102(      556.000,      536.000)
;     103(      556.000,      909.000)
;     104(      557.000,      536.000)
;     105(      557.000,      537.000)
;     106(      557.000,      538.000)
;     107(      559.000,      534.000)
;     108(      559.000,      535.000)
;     109(      560.000,      533.000)
;     110(      560.000,      534.000)
;     111(      560.000,      535.000)
;     112(      560.000,      536.000)
;     113(      561.000,      534.000)
;     114(      561.000,      535.000)
;     115(      561.000,      536.000)
;     116(      561.000,      537.000)
;     117(      562.000,      534.000)
;     118(      562.000,      535.000)
;     119(      562.000,      536.000)
;     120(      564.000,      898.000)
;     121(      565.000,      530.000)
;     122(      565.000,      531.000)
;     123(      566.000,      529.000)
;     124(      566.000,      530.000)
;     125(      569.000,      529.000)
;     126(      569.000,      530.000)
;     127(      573.000,      914.000)
;     128(      575.000,      262.000)
;     129(      577.000,      261.000)
;     130(      577.000,      919.000)
;     131(      577.000,      920.000)
;     132(      578.000,      260.000)
;     133(      578.000,      920.000)
;     134(      578.000,      921.000)
;     135(      579.000,      260.000)
;     136(      580.000,      922.000)
;     137(      581.000,      921.000)
;     138(      581.000,      922.000)
;     139(      582.000,      921.000)
;     140(      582.000,      922.000)
;     141(      583.000,      920.000)
;     142(      583.000,      921.000)
;     143(      583.000,      922.000)
;     144(      584.000,      919.000)
;     145(      584.000,      920.000)
;     146(      585.000,      919.000)
;     147(      586.000,      919.000)
;     148(      588.000,      186.000)
;     149(      590.000,      918.000)
;     150(      590.000,      919.000)
;     151(      615.000,      904.000)
;     152(      615.000,      905.000)
;     153(      616.000,      906.000)
;     154(      620.000,      187.000)
;     155(      626.000,      153.000)
;     156(      627.000,      153.000)
;     157(      627.000,      154.000)
;     158(      628.000,      154.000)
;     159(      628.000,      155.000)
;     160(      629.000,      155.000)
;     161(      669.000,      868.000)
;     162(      671.000,      323.000)
i=110
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=150
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=151
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=152
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=160
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
flush,!journal
for i=150,258 do begin plot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1  & print,i,ib[i] & wait,1.5
;     150(      590.000,      919.000)
;     151(      615.000,      904.000)
;     152(      615.000,      905.000)
;     153(      616.000,      906.000)
;     154(      620.000,      187.000)
;     155(      626.000,      153.000)
;     156(      627.000,      153.000)
;     157(      627.000,      154.000)
;     158(      628.000,      154.000)
;     159(      628.000,      155.000)
;     160(      629.000,      155.000)
;     161(      669.000,      868.000)
;     162(      671.000,      323.000)
;     163(      672.000,      322.000)
print,i
;     164
i=161
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
for i=150,258 do begin plot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1  & print,i,ib[i] & wait,1.5
;     150(      590.000,      919.000)
;     151(      615.000,      904.000)
;     152(      615.000,      905.000)
;     153(      616.000,      906.000)
;     154(      620.000,      187.000)
;     155(      626.000,      153.000)
;     156(      627.000,      153.000)
;     157(      627.000,      154.000)
;     158(      628.000,      154.000)
;     159(      628.000,      155.000)
;     160(      629.000,      155.000)
;     161(      669.000,      868.000)
;     162(      671.000,      323.000)
;     163(      672.000,      322.000)
print,i
;     164
i=161
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
for i=160,258 do begin plot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.4,.4],yst=1  & print,i,ib[i] & wait,1.5
;     160(      629.000,      155.000)
;     161(      669.000,      868.000)
;     162(      671.000,      323.000)
;     163(      672.000,      322.000)
;     164(      672.000,      323.000)
;     165(      672.000,      888.000)
;     166(      673.000,      323.000)
;     167(      673.000,      888.000)
;     168(      675.000,      399.000)
;     169(      676.000,      399.000)
;     170(      676.000,      653.000)
;     171(      677.000,      399.000)
;     172(      726.000,      866.000)
;     173(      727.000,      866.000)
;     174(      733.000,      880.000)
;     175(      734.000,      880.000)
;     176(      735.000,      879.000)
;     177(      735.000,      880.000)
;     178(      768.000,      891.000)
;     179(      768.000,      892.000)
;     180(      769.000,      891.000)
i=171
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=174
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
for i=170,258 do begin plot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.4,.4],yst=1  & print,i,ib[i] & wait,1.5
;     170(      676.000,      653.000)
;     171(      677.000,      399.000)
;     172(      726.000,      866.000)
;     173(      727.000,      866.000)
;     174(      733.000,      880.000)
;     175(      734.000,      880.000)
;     176(      735.000,      879.000)
;     177(      735.000,      880.000)
;     178(      768.000,      891.000)
;     179(      768.000,      892.000)
;     180(      769.000,      891.000)
;     181(      770.000,      917.000)
;     182(      771.000,      885.000)
;     183(      772.000,      883.000)
;     184(      773.000,      880.000)
;     185(      773.000,      881.000)
;     186(      773.000,      921.000)
;     187(      774.000,      836.000)
;     188(      774.000,      880.000)
;     189(      775.000,      834.000)
i=182
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=178
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
; % No commands containing '288' exist in the command buffer.
; % No commands containing '287' exist in the command buffer.
for i=170,258 do begin plot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.4,.4],yst=1  & print,i,ib[i] & wait,1.5
;     170(      676.000,      653.000)
;     171(      677.000,      399.000)
;     172(      726.000,      866.000)
;     173(      727.000,      866.000)
;     174(      733.000,      880.000)
;     175(      734.000,      880.000)
;     176(      735.000,      879.000)
;     177(      735.000,      880.000)
;     178(      768.000,      891.000)
;     179(      768.000,      892.000)
;     180(      769.000,      891.000)
;     181(      770.000,      917.000)
;     182(      771.000,      885.000)
;     183(      772.000,      883.000)
;     184(      773.000,      880.000)
;     185(      773.000,      881.000)
;     186(      773.000,      921.000)
;     187(      774.000,      836.000)
;     188(      774.000,      880.000)
;     189(      775.000,      834.000)
;     190(      775.000,      835.000)
;     191(      790.000,      835.000)
;     192(      800.000,      684.000)
;     193(      801.000,      683.000)
;     194(      801.000,      685.000)
;     195(      802.000,      838.000)
;     196(      803.000,      689.000)
;     197(      804.000,      687.000)
;     198(      804.000,      688.000)
;     199(      807.000,      860.000)
;     200(      826.000,      638.000)
;     201(      826.000,      639.000)
;     202(      835.000,      655.000)
;     203(      835.000,      656.000)
;     204(      845.000,      900.000)
;     205(      852.000,      892.000)
;     206(      855.000,      668.000)
;     207(      863.000,      891.000)
;     208(      864.000,      890.000)
;     209(      874.000,      876.000)
;     210(      884.000,      699.000)
;     211(      891.000,      679.000)
;     212(      909.000,      422.000)
;     213(      910.000,      420.000)
;     214(      910.000,      421.000)
;     215(      910.000,      422.000)
;     216(      910.000,      431.000)
;     217(      910.000,      433.000)
;     218(      911.000,      421.000)
;     219(      911.000,      422.000)
;     220(      912.000,      423.000)
;     221(      912.000,      424.000)
;     222(      912.000,      436.000)
;     223(      912.000,      437.000)
;     224(      913.000,      423.000)
;     225(      913.000,      439.000)
;     226(      913.000,      441.000)
;     227(      914.000,      441.000)
i=209
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=212
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
for i=210,258 do begin plot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.4,.4],yst=1  & print,i,ib[i] & wait,1.5
;     210(      884.000,      699.000)
;     211(      891.000,      679.000)
;     212(      909.000,      422.000)
;     213(      910.000,      420.000)
;     214(      910.000,      421.000)
;     215(      910.000,      422.000)
;     216(      910.000,      431.000)
;     217(      910.000,      433.000)
;     218(      911.000,      421.000)
;     219(      911.000,      422.000)
;     220(      912.000,      423.000)
;     221(      912.000,      424.000)
;     222(      912.000,      436.000)
;     223(      912.000,      437.000)
;     224(      913.000,      423.000)
;     225(      913.000,      439.000)
;     226(      913.000,      441.000)
;     227(      914.000,      441.000)
;     228(      924.000,      206.000)
;     229(      924.000,      207.000)
;     230(      924.000,      208.000)
;     231(      924.000,      209.000)
;     232(      925.000,      209.000)
;     233(      940.000,      415.000)
;     234(      940.000,      683.000)
;     235(      940.000,      684.000)
;     236(      941.000,      683.000)
;     237(      942.000,      683.000)
;     238(      944.000,      471.000)
;     239(      944.000,      472.000)
;     240(      944.000,      473.000)
;     241(      945.000,      472.000)
;     242(      945.000,      473.000)
;     243(      951.000,      737.000)
;     244(      952.000,      170.000)
;     245(      953.000,      170.000)
;     246(      954.000,      170.000)
;     247(      957.000,      168.000)
;     248(      962.000,      488.000)
;     249(      962.000,      489.000)
;     250(      962.000,      490.000)
;     251(      962.000,      491.000)
;     252(      976.000,      252.000)
;     253(      993.000,      411.000)
;     254(      999.000,      148.000)
;     255(      999.000,      150.000)
;     256(      999.000,      151.000)
;     257(      1003.00,      429.000)
;     258(      1004.00,      481.000)
i=258
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=237
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=213
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
; % No commands containing 'afda' exist in the command buffer.
help,ir
for i=0,100 do begin plot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.4,.4],yst=1  & print,i,ir[i] & wait,1.5
;       0(      52.0000,      58.0000)
;       1(      60.0000,      313.000)
;       2(      221.000,      520.000)
;       3(      221.000,      521.000)
;       4(      222.000,      520.000)
;       5(      222.000,      521.000)
;       6(      303.000,      617.000)
;       7(      343.000,      938.000)
;       8(      345.000,      934.000)
i=0
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
flush,!journal
for i=0,100 do begin plot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.4,.4],yst=1  & print,i,ir[i] & wait,1.5
;       0(      52.0000,      58.0000)
;       1(      60.0000,      313.000)
;       2(      221.000,      520.000)
;       3(      221.000,      521.000)
;       4(      222.000,      520.000)
;       5(      222.000,      521.000)
;       6(      303.000,      617.000)
i=2
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=3
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
for i=0,100 do begin plot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.4,.4],yst=1  & print,i,ir[i] & wait,1.5
;       0(      52.0000,      58.0000)
;       1(      60.0000,      313.000)
;       2(      221.000,      520.000)
;       3(      221.000,      521.000)
;       4(      222.000,      520.000)
;       5(      222.000,      521.000)
;       6(      303.000,      617.000)
;       7(      343.000,      938.000)
;       8(      345.000,      934.000)
;       9(      345.000,      935.000)
;      10(      346.000,      934.000)
;      11(      346.000,      935.000)
;      12(      347.000,      934.000)
;      13(      348.000,      438.000)
;      14(      348.000,      933.000)
;      15(      348.000,      934.000)
;      16(      349.000,      933.000)
;      17(      359.000,      430.000)
;      18(      360.000,      430.000)
;      19(      360.000,      434.000)
;      20(      360.000,      435.000)
;      21(      360.000,      436.000)
;      22(      361.000,      430.000)
;      23(      361.000,      431.000)
;      24(      361.000,      435.000)
;      25(      362.000,      116.000)
;      26(      362.000,      117.000)
;      27(      362.000,      431.000)
;      28(      362.000,      435.000)
;      29(      363.000,      112.000)
;      30(      363.000,      113.000)
i=14
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=15
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
for i=10,100 do begin plot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.4,.4],yst=1  & print,i,ir[i] & wait,1.5
;      10(      346.000,      934.000)
;      11(      346.000,      935.000)
;      12(      347.000,      934.000)
;      13(      348.000,      438.000)
;      14(      348.000,      933.000)
;      15(      348.000,      934.000)
;      16(      349.000,      933.000)
;      17(      359.000,      430.000)
;      18(      360.000,      430.000)
;      19(      360.000,      434.000)
;      20(      360.000,      435.000)
;      21(      360.000,      436.000)
;      22(      361.000,      430.000)
i=17
IDL> label,lin=0,'ATST
; % Syntax error.
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=18
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=20
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=21
IDL> label,x=.6,y=.3,lin=2,'Sim'
; % Syntax error.
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
for i=20,100 do begin plot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.4,.4],yst=1  & print,i,ir[i] & wait,1.5
;      20(      360.000,      435.000)
;      21(      360.000,      436.000)
;      22(      361.000,      430.000)
;      23(      361.000,      431.000)
;      24(      361.000,      435.000)
;      25(      362.000,      116.000)
;      26(      362.000,      117.000)
;      27(      362.000,      431.000)
;      28(      362.000,      435.000)
;      29(      363.000,      112.000)
;      30(      363.000,      113.000)
;      31(      363.000,      114.000)
;      32(      363.000,      115.000)
;      33(      363.000,      116.000)
;      34(      363.000,      432.000)
;      35(      474.000,      41.0000)
;      36(      482.000,      274.000)
;      37(      502.000,      933.000)
;      38(      503.000,      928.000)
;      39(      503.000,      929.000)
;      40(      531.000,      732.000)
;      41(      532.000,      731.000)
;      42(      532.000,      732.000)
;      43(      533.000,      731.000)
;      44(      533.000,      732.000)
;      45(      534.000,      731.000)
;      46(      534.000,      732.000)
;      47(      535.000,      731.000)
;      48(      535.000,      732.000)
;      49(      536.000,      731.000)
;      50(      536.000,      732.000)
;      51(      536.000,      733.000)
;      52(      536.000,      736.000)
;      53(      537.000,      732.000)
;      54(      537.000,      733.000)
;      55(      538.000,      732.000)
;      56(      538.000,      734.000)
;      57(      538.000,      735.000)
;      58(      539.000,      733.000)
;      59(      539.000,      734.000)
i=57
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=44
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
for i=50,100 do begin plot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.4,.4],yst=1  & print,i,ir[i] & wait,1.5
;      50(      536.000,      732.000)
;      51(      536.000,      733.000)
;      52(      536.000,      736.000)
;      53(      537.000,      732.000)
;      54(      537.000,      733.000)
;      55(      538.000,      732.000)
;      56(      538.000,      734.000)
;      57(      538.000,      735.000)
;      58(      539.000,      733.000)
;      59(      539.000,      734.000)
;      60(      539.000,      737.000)
;      61(      540.000,      734.000)
;      62(      540.000,      736.000)
;      63(      540.000,      737.000)
;      64(      574.000,      198.000)
;      65(      636.000,      353.000)
;      66(      636.000,      354.000)
;      67(      636.000,      355.000)
;      68(      638.000,      355.000)
;      69(      782.000,      613.000)
;      70(      811.000,      610.000)
;      71(      812.000,      610.000)
;      72(      812.000,      611.000)
;      73(      813.000,      609.000)
;      74(      813.000,      610.000)
;      75(      813.000,      611.000)
;      76(      814.000,      610.000)
;      77(      823.000,      637.000)
;      78(      838.000,      928.000)
;      79(      840.000,      250.000)
i=52
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=60
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=65
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
for i=60,100 do begin plot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.4,.4],yst=1  & print,i,ir[i] & wait,1.5
;      60(      539.000,      737.000)
;      61(      540.000,      734.000)
;      62(      540.000,      736.000)
;      63(      540.000,      737.000)
;      64(      574.000,      198.000)
;      65(      636.000,      353.000)
;      66(      636.000,      354.000)
;      67(      636.000,      355.000)
;      68(      638.000,      355.000)
;      69(      782.000,      613.000)
;      70(      811.000,      610.000)
;      71(      812.000,      610.000)
;      72(      812.000,      611.000)
;      73(      813.000,      609.000)
;      74(      813.000,      610.000)
;      75(      813.000,      611.000)
;      76(      814.000,      610.000)
;      77(      823.000,      637.000)
;      78(      838.000,      928.000)
;      79(      840.000,      250.000)
i=69
IDL> label,lin=4,'Hinode'
; % Syntax error.
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
for i=70,100 do begin plot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.4,.4],yst=1  & print,i,ir[i] & wait,1.5
;      70(      811.000,      610.000)
;      71(      812.000,      610.000)
;      72(      812.000,      611.000)
;      73(      813.000,      609.000)
;      74(      813.000,      610.000)
;      75(      813.000,      611.000)
;      76(      814.000,      610.000)
;      77(      823.000,      637.000)
;      78(      838.000,      928.000)
;      79(      840.000,      250.000)
;      80(      840.000,      251.000)
;      81(      878.000,      182.000)
i=75
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
i=79
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
for i=80,100 do begin plot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.4,.4],yst=1  & print,i,ir[i] & wait,1.5
;      80(      840.000,      251.000)
;      81(      878.000,      182.000)
;      82(      879.000,      182.000)
;      83(      879.000,      183.000)
;      84(      879.000,      184.000)
;      85(      895.000,      973.000)
;      86(      899.000,      220.000)
;      87(      911.000,      13.0000)
;      88(      912.000,      12.0000)
;      89(      912.000,      13.0000)
;      90(      939.000,      995.000)
;      91(      940.000,      994.000)
;      92(      940.000,      995.000)
;      93(      941.000,      993.000)
;      94(      941.000,      994.000)
;      95(      942.000,      994.000)
;      96(      945.000,      994.000)
;      97(      946.000,      993.000)
;      98(      946.000,      994.000)
;      99(      946.000,      995.000)
; % Attempt to subscript IR with I is out of range.
i=88
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.35],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.30],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
ps
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.30],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=4
label,x=.6,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
psend
x
print,ir[i]
;(      912.000,      12.0000)
$mv idl.ps V_sim_atst_hinode_12vr_912-12.ps
flush,!journal
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.30],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=4
label,x=.9,y=.3,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
psend
; % Keyword CLOSE not allowed in call to: DEVICE
x
ps
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.30],yst=1
; % Variable is undefined: W.
for ii=-3,3,2 do b
; % Attempt to call undefined procedure/function: 'B'.
psend
x
ps
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.30],yst=1
; % Variable is undefined: W.
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
; % Variable is undefined: W.
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=0
; % Variable is undefined: W.
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=4
; % Variable is undefined: W.
label,x=.6,y=.9,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
psend
x
retall
ps
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.30],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=4
label,x=.6,y=.9,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
psend
x
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.30],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=4
label,x=.6,y=.9,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
$mv idl.ps V_sim_atst_hinode_12vr_912-12.ps
ps
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.30],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=4
label,x=.5,y=.8,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
psend
x
$mv idl.ps V_sim_atst_hinode_12vr_912-12.ps
ps
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.30],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=4
label,x=.55,y=.8,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
psend
x
$mv idl.ps V_sim_atst_hinode_12vr_912-12.ps
print,i
;      88
i=0
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.30],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=4
label,x=.55,y=.8,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
ps
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.30],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=4
label,x=.55,y=.8,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
psend
x
print,ir[i]
;(      52.0000,      58.0000)
$mv idl.ps V_sim_atst_hinode_12vr_52-58.ps
i=3
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.30],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=4
label,x=.55,y=.8,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.15,.40],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=4
label,x=.55,y=.8,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
ps
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.15,.40],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=4
label,x=.55,y=.8,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
psend
x
print,ir[i]
;(      221.000,      521.000)
$mv idl.ps V_sim_atst_hinode_12vr_221-521.ps
print,i
;       3
print,i,ir[i]
;       3(      221.000,      521.000)
ps
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.15,.40],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=4
label,x=.6,y=.8,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
psend
x
$mv idl.ps V_sim_atst_hinode_12vr_221-521.ps
i=15
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.3,.30],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=4,lin=4
label,x=.55,y=.8,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
ps
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=4,yr=[-.3,.25],yst=1
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vatst[*,real_part(ir[i]),imaginary(ir[i])],thick=3,lin=0
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=5,lin=4
label,x=.55,y=.8,lin=2,'Sim'
label,lin=0,'ATST
label,lin=4,'Hinode'
psend
x
print,i,ir[i]
;      15(      348.000,      934.000)
$mv idl.ps V_sim_atst_hinode_12vr_348-934.ps