4Misc_Start ROGI'XXp>><4!*8 WINSPOOLHP LaserJet 4MLPT1:HP LaserJet 4M@g XXLetter 8 1/2 x 11 in@MSUDNHP LaserJet 4M;d ROGI'XXp>><4!*8 WINSPOOLHP LaserJet 4MLPT1:HP LaserJet 4M@g XXLetter 8 1/2 x 11 in@MSUDNHP LaserJet 4M;d ROGI'XXp>><4!*8 WINSPOOLHP LaserJet 4MLPT1:HP LaserJet 4M@g XXLetter 8 1/2 x 11 in@MSUDNHP LaserJet 4M;d ROGI'XXp>><4!*8 WINSPOOLHP LaserJet 4MLPT1:HP LaserJet 4M@g XXLetter 8 1/2 x 11 in@MSUDNHP LaserJet 4M;d oSummary_Layout_2D;ROGI'XX tREP/winspoolHP LaserJet 4V/4MV PostScripthp4mvHP LaserJet 4V/4MV PostScriptG 4dXA4PRIV'''' oSummary_Layout_2D_2;ROGI'XX tREP/winspoolHP LaserJet 4V/4MV PostScripthp4mvHP LaserJet 4V/4MV PostScriptG 4dXA4PRIV'''' oRValues_contour_2D_layout;ROGI'XX tREP/winspoolHP LaserJet 4V/4MV PostScripthp4mvHP LaserJet 4V/4MV PostScriptG 4dXA4PRIV'''' oRValues_3D_layout;ROGI'XX tREP/winspoolHP LaserJet 4V/4MV PostScripthp4mvHP LaserJet 4V/4MV PostScriptG 4dXA4PRIV'''' oLayout0;ROGI'XX tREP/winspoolHP LaserJet 4V/4MV PostScripthp4mvHP LaserJet 4V/4MV PostScriptG 4dXA4PRIV'''' osummary_layout_3D;ROGI'XX tREP/winspoolHP LaserJet 4V/4MV PostScripthp4mvHP LaserJet 4V/4MV PostScriptG 4dXA4PRIV'''' oSummary_layout_3D_2;ROGI'XX tREP/winspoolHP LaserJet 4V/4MV PostScripthp4mvHP LaserJet 4V/4MV PostScriptG 4dXA4PRIV'''' oSummary_layout_3D_1;ROGI'XX tREP/winspoolHP LaserJet 4V/4MV PostScripthp4mvHP LaserJet 4V/4MV PostScriptG 4dXA4PRIV'''' ROGI'XXp>><4!*8 WINSPOOLHP LaserJet 4MLPT1:HP LaserJet 4M@g XXLetter 8 1/2 x 11 in@MSUDNHP LaserJet 4M;d  ROGI'XXp>><4!*8 WINSPOOLHP LaserJet 4MLPT1:HP LaserJet 4M@g XXLetter 8 1/2 x 11 in@MSUDNHP LaserJet 4M;d ^Graph*@@??_RWDashSettings#  !CommandHNormal@ Marlett<HHHH$$HistoryHNormal@ Times New Roman<HHHH$$4History4History4History4History4History4HistoryProcedure GNormal@ Arial<HHHH$$4Procedure 4Procedure 4Procedure 4Procedure 4Procedure 4Procedure homeYQdZxC:Documents and Settings:Ian:Desktop:NewWeb:Files:C:Documents and Settings:Ian:Desktop:NewWeb:FilesIan\Desktop\NewWebFilesH!`3E(9wXH wwklwwB_{temppathG?dB6Q:meis:dysi:dysisimuls:dy100t48:Q:meis:dysi:dysisimuls:dy100t48dysi\dysisimulsdy100t48top\NewWebFilesH!`3E(9wXH wwklwwB_ 4Misc_End tXOPState_Start@ Data Browserdysisimuls\Surface Plottermuls:dy100t48:4XOPState_End "J#C҉h>f'&&,pc@韮eb@Xz`pa@[Ha@=yL`@(>\`@x6n`@$pJK(|a@=ʹ`@7`@h%r/L`@3`@[x`@:p%`@o bd@uc@b@Ovj6b@yyGea@3z.a@r`@=tĉ`@~u``@P`@[`@'cctqd@bnc@oac@ -!uib@T+a@oa@:a@d N `@ X`@h!Xԗ`@fD`@=wd@u~c@"5c@uEb@i~2b@В*a@ʙ-U}a@gmMa@E2)"2a@dS&a@'ׁ-3a@΍|d@Dd@>7Y{c@}b@色b@);b@4:a@{a@岴a@]a@se, a@.sim2t?$+Nx&?HzG?w/?"~j?B`"?MbX9?(\?A`"?X9v?MbX?q= ףp?/$?p= ף?MbX?Mb?Dl?"~j?Cl?&1?Q?i|?5?(\?sh|??y&1?+?(\? r?Gz?l?On?5^I ?(\?#~j?1Zd?v? rh?q= ףp?/$?On?p= ף?NbX9?"~? r?@5^I ?> ףp=?kt?> ףp=?^I +?X9v?y&1?MbX9?bX9?Q?1Zd?ʡE?K7?Mb?l?(\?i|?5?~jt?K7A?HzG?"~?"~j?(\?5^I ?Gz?-?Dl? rh?m?x&?}?5^I?S㥛?$C? ףp= ?(\?m?MbX?"~j?)\(?x&1?Q?S?jt?S?Zd;O?x&1?y&1?E?(\?|?5^?m?`"?bX9?S?x&1?{Gz???oʡ?|?5^?Gz?zG?zG?7A`?-?rh|?Cl?x&1?Zd;O?rh|?Zd;?rh|?333333?.sim0t?$+ND@D@333333D@LD@gfffffD@D@D@33333D@D@gffffD@E@E@333333E@LE@gfffffE@E@E@33333E@E@gffffE@F@F@333333F@LF@gfffffF@F@F@33333F@F@gffffF@G@G@333333G@LG@gfffffG@G@G@33333G@G@gffffG@H@H@333333H@LH@gfffffH@H@H@33333H@H@gffffH@I@I@333333I@LI@gfffffI@I@I@33333I@I@gffffI@J@J@333333J@LJ@gfffffJ@J@J@33333J@J@gffffJ@K@K@333333K@LK@gfffffK@K@K@33333K@K@gffffK@L@L@333333L@LL@gfffffL@L@L@33333L@L@gffffL@M@M@333333M@LM@gfffffM@M@M@33333M@M@gffffM@N@N@333333N@LN@gfffffN@N@N@33333N@N@gffffN@O@O@333333O@LO@gfffffO@O@.exp1t?+N+N?'?F ?[4A ?槙? )?yoE?%oکD?9ҭ?OM? xڌ?0T?Ωu$2?~?ծ e?y|E?vv?(4?CFu?*t?}?~?v?.&?*XEp?ja?T=? K|?E1X? ?j?tO]G?J6?S?S_~5?mBl?vH ?s@ r@r@5Cr3CV1CH4CO@ECfWCdCE{CE(CzCԢC͞C"CVCqvC~CaC>_CJC*\CsΘCC5|C3hCgCuCXCCCSCC}C +uCoC:|bCV[C,x^C@K]C=\CRiC:pCl8zC{CC;ɊC:CCۖCCYǙCC?A?\m??/ ?^ z63? s=X?Jx}?h*U?TJ??J ?(g.?P?xdq-q?R&?|חײ?h$̀x?ޙ?F=?l .?| ӒK?3h?3? ??H>?$Y?h ?,#?8e?aٲBD?_I?M?)Q? dvU? iX?h4Z?D6_w\?`^?d9/^?> _?LjpƼ^?pe]?s=\? ,Z?!X? q U?b>ȔQ?˦~M?DHH?C?Xx=?,~7?4J0?ʽA)?W!?x(Q?x \?.1eexp_fitt?+N+Nh۸?! ?0N?1-`K?lY?D8?vD?Dk?>9T?6.?#0?PXR̿?qw|?J a?ܳZ?&g!_r?(:!?,V?V?+Ys?q'?fJ)?5j?Vz?cw?(H_?8Fm?HT-p}?xy?p_v?7ils?+Tp?uiZm?J?ԔDmJ? ]0:J?"g&J? 0J?zZZJ?ߢJ?J K?b)!K?ˏ5L?L?ƻM?XN?fO?kd?Q?^*R?ԣT?8DU? ;sW?}`oMY?#F[?M^]? s_?!,go?0vr?Q u? x?]|?U?Ll?)t[?35B?H?.n?-^?kt?&q?;FS8?2?֮?|_Գ?x? <,?]}?*?ė?(@N?y$?Ԛd?m-?Q_?5!?,zz!?򉿰?E^?eC+ ?@"?"?*KL"?QW_coef?+N+NYb$`@PB: קc?F6fit_P1_RValuenw z|?+N+N#CC0BCCCUC C[CyCr4CCCoC0CCXCցC0KCfCyCgC2CUC\*CCC CClC|LC.ChCPBBxB!cBC"CWCCCCS?C}C;CyCDCC\C CgCC@C_UCZC1CSCtCC(cCLCL! C( C CuL C C2 C[ C` CAd C CG C CVfit_P3_RValueZp \?11|l>)|>`> >=>>? ?g???IW ? ?n ? ?/?.n?l??+??Da???R??0?~??,????;?;??c??3?z?6????!?Y??~ ?t! ? ??#?/N?b?a?M>8>">>!>>>K>>9>~T>C>>F> >I>]>E>>>=>(`z>0zj>NI>c8>u'>>>x4=&==8b==<F6fit_P2_RValuenw z|?+N+Nb%CK$C#C#CQ"C!Ct Cq C{aCCCBC)CvC6C7CC.=CC[CVCCCSC2CSCC -C5CWCCC rCdChC:CgCF Cn C C;~ C C̕ C$ C CH C Ct C CC$FCkCC!*CC wC C,CyC)CAC CDCCCrC0CCCNwC=CCBCC'kC;CACCCMCmCIC(C CBBqBDB~BBBUB'BBwB=jByaB\BB\B_BvgB8sBB BBIBBBn8BeBRvalue0d?$$CNC8CNC+ CuC{7C00RC0Cy&aD)DrCCh9nCҜCC CY'C@*{Ew@DIf@: X@e^uK U@dǜY@Af[@@@l[@^dc@t~kk@<(w@E@T"jhu@q'N`@+Q@NSjeT@#}\@ f^@jxPBeX@3]@(|f@x'ht@@<~JQs@Jgݟ\@w1"vP@ڱV@my`!a@)Ca@0uR@#pb@"ko@ΫJ8x@ߩnmk@@lX@Q@^l`@隸f@}d@V "rW@R*!BKQ@ ~!a@,A3o@P8kXw@gh@2eW@H hU@ǜwnc@wʷh@? d@b`y2Y@|_R@V+`@6$n@"v@Ɔ$h@xۄY@CUY@nIf@Jj@3Хf@Hq]@ʊOV@رlsa@Fn@ΣI0{]s>[MF?0È?q>2>Ahx?B>`^?>7>T?p?K??'>,>?p6?iI0>Dr?#??΋> =:>>9?/>o>r?Uv?q>/+S=R=I>>*?V>J?']?a>,=U=1Y>TW ?)u?pHY?=??o:/?@?D_>8D>OBr>%??1Ӕ?zK???@?D>>L?̶?1!?P?+?c?tKh?4&?zc?@ń8!11QMatrix318 ????h?BX6?^>4>"G>dž=R6=U>>>E>v1>!?7>=W~=k<>`p>9?/>o>r?Uv?q>/+S=R=I>\d@?e+?+>@i?z?p?O>FZ=Q>(?*D?M?fD?ɝ?I??m>=6:?:q?^?%?HK?IP?^j?*>|c>ۻ ?4|?T???i?#?f ?r?>E#3?w|?6g?l_W?}?2q??VD?(ò>os??q? e7?B???Л?1>yfP3_RValue ?119?/>o>r?Uv?q>/+S=R=I>sim3?+N+N?@Nsim4?+N+N,@P@sim5?+N+Nףp= <@Pd@sim6?+N+N?@+sim7?+N+N??Hsim8?+N+NQ?5=sim9?+N+N @8sim10?+N+N??BMsim11?+N+NQ?Q?sim12?+N+N@@sim13?+N+N??KOsim14?+N+NQ?sim15?+N+N#~j@^I +@wsim16?+N+NQ?Q?@;@=@<@sim61?+N+N????Lsim62?+N+NQ?Q??sim63?+N+N#~j@ @^I +@Asim64?+N+Nw/=5@Gzt4@w/=5@Gzt4@fsim65?+N+N(\µ?(\µ?(\µ?(\µ?Csim66?+N+N?@A@@@@@sim67?+N+N????Lsim68?+N+NQ?Q?'Ksim69?+N+NZd@|?5^? @Jsim70?+N+Nףp= 7@9v_8@9v_8@ףp= 7@jsim71?+N+N(\µ?(\µ?(\µ?(\µ?CKsim72?+N+NB@B@A@C@sim73?+N+N????Ksim74?+N+NQ?Q?sim75?+N+N#~j@Zd@|?5^?^I +@Hsim76?+N+NM;@:@:@M;@dsim77?+N+N(\µ?(\µ?(\µ?(\µ?9Ksim78?+N+NE@C@D@D@sim79?+N+N????Jsim80?+N+NQ?Q?Asim81?+N+N#~j@ @^I +@Psim82?+N+N/ݤ>@\(=@/ݤ>@\(=@hsim83?+N+N(\µ?(\µ?(\µ?(\µ?.cVeg1t?+N+NA`"??np?Eh?&1?"?d;O?Zd;߇??Gz\?ʡE[?Zd3?{G?Zd;߿?(\?Gz^?Gz4?K?l?^I ?(\?Pns?Q?A`"[?!rhm?\(\?J +?ҝ?43333?,w?sh|/?ʡE6?Dl?Mb8?!rhM?&1?%C?J +?V-B?S?(\r?/$?Zd;?V?v?@5^I|?zG?Zd;? r?)\?#~j ?V6CО?[yfg?. Veg2t?+N+N= ףp=?~jt?Mbp?MbH?|?5^?v?+^?jth?n? rh?x&??ʡ%?m? rhA?-?K7A? b?Q? r?ˡE?n?*\(?GzN?n ?9v_?l?L7A`?B`"?Zd;O?n?Cl?|?5^z?Mb?sh|?Zd?-?|?5^J?QF?$~j?~jt?FԨ?lr?h|?E?Q?+?ʡE?zG?x?x&1`?"?/$1?SC?m2?/$?Zd;?&1?tf?Zd;O?o!?s?n!? +9?|?5^?ףp= ?"~?"~?Hz?i|?5?Zd;?Q?gffffV?(\?S+?"~z?R?(\E?m?h|??7A`?K7A@?ʡ?|?5?T㥛P?y&1?Ƌ?K79?`"? +?&1\?S[?/$?"~?J +?tf?0$?z&1?^I ?R1?Mb?(\/? ףp=?'1Z?\(?cX98?w/?v??i|??,?SC?Cl[?Qn?(\u?.cVeg3t?+N+Nx&Q?I +&?"~j?1Z$?Cl[?+Η?Cl[?-?L7A`E?K79?/$??I +7?B`"?Dl? ףp=Z?_I {?V?Zd;/?x?p= C?nz?̬?jt?Gz?"rh?V-?x&?K?^I ?S[?Zd ?q= ף?Q?Cl? ףp=z?y&1?S3?t.?43333?"~?Hz?MbX?K7)?n?ףp= ?Mb?S㥛?S㥛|? J?Pn?S?I +?,?`"?V-]?T㥛?rhL?43333?t?QX?֣p= ?v?8A`?#~j|?^I ?n?$s?+?&1L?(\U??5^I ?V-?ʡE?(\?V-/?Q?~jt3?ʡEf?-?"?`"9?fffff?Gz?Q?OnB?t?ʡE?h|?e?`"I?jtT?S㥛?Zd;?#?&C{?v?KW?y&1? rha? ?Zd[?uV?+N?Qn?Gzn?S?A`B?R?~jt?L7A`?E?x)?l:?2ZL?/$Y? rh]?.8Veg4t?+N+Nna?n0?n`?K7A?~jt?lR?NbX9d?~jt?Pn??䥛 0?:v/?e;OG?/$?p= ף?&1 ? r?\( ףp?̡E? ףp=?M?nx?bX9@?Mb?x?X9Ⱦ?|?5^?K7?li?x6?zG?Zd;O?\(?{G?Zd;ߟ?Q;?J +?(\5?v?@5^I?(\??"?,?(1?$c?Zd{?J +?{G:?Q?/$?C ?㥛 ?+.?Zd;?Mb?Q?tVn?Q?S?1Z4?43333?NbX9?L7A`E?On?`"I?S㥛?n?(1z?vJ?|?5>?NbX9T?"~j?V-? r8?ʡE?"?"~z?K7?Q?n ?p= ׃??Gzn?%C?7A`? ףp=?K?M?(\?Mb?"?7A` ?X9v.?S㥛 ףp]?Zd;O%?jt?v?NbX9?43333?El?Zd;OU?F ?|?5?zG?uV?Cl?1Zt?K?}?5^?C?n?n?ʡE#?kt?rh?̡Em?2ZD?J +V?Q?n?K7Ap?zG?Q??K7a? +9?6^I ?A`Ђ?ʡES?la?)\X?|?5?{Gz?S?6^I ?Zd;?SS?F?Mb?#~jL?xV?v?ʡ?l?Dl?{Gz?Clg?"?V-r?Dl?w/? ?䥛 ?(\2?8A` ?n?S?x?$?/$??Q}?rhD?@5^I ?Dl?p= ף?Cl?= ףpm?zGY?nA?&1 ?tf?FԈ?jt?K7?Zd;OM?(\?}?~jt?ʡE??/$?-臨?K7A?lQ?/$&?M2?i?ʡ?+?䥛 p?"n?"~?q= ף?J +?{G?/$?m?d;O?E$?^I k?-?ʡ5?i|??E$?sh|?)\?$c?E4?t&?Mb8?Ra?L7A`? rh?}?5^i?33333?n@?Cl?C?;O?/T??Pn3?w/?|?5^? ?8A`?MbH?~jt?Cl?"~?Q?/$?i|??/?v"?S㥛$?.)Veg7t?+N+NV-o?bX9(?X9v?](? ?l?&1?n?A`?~j??Q?GzN?"^?x&1?Dl??ʡE?x?~j?K7A?6^I b?}?5^)?V-?p= ׻?v?Zd;o?CT? B?n*?+?p= 3?d;OW?Onb?/$Q?V-"?o?"^?S?e;O?S㥛?x?I +?T㥛İ?R1?S?F?X96?Cl?l?\(?S㥛?@5^I?+?Zd;_?Y9v?~j?sh|?ʡEv?MbX?,?$Ck?x&?"~Z?Q?)\?bX9h?Dl9?I +?x&1?K79? rhm?Cl?Mb?~j?K7?:v_?"~?6^I b??䥛 p?S? +I?x/?6^I "?~jt?|?5?A`B?x?F?0$?Q?l??S㥛?/$?(\?.Veg8t?+N+NZd;Om? ?E?-罹?B`"{?(\?FԸ?K7?^I +?,?OnB?l?/$E?p= 3??Q{?^I ?8A`0?Pn?p= ?m?#~j?On?^I  ?ʡE?T㥛?X9v^?^I ? r(?S㥛?Cl?@5^I ףp? r?Q??/$?X9?R? ףp=?-阮?ˡE?,Χ?.V?zGQ?ʡE?ʡES?Sk?Cl?ʡ?sh|? ףp=?~jt?EH?$s?:v_?+?](\?gffff?Dl?L7A`?Zd;?(\•?ʡE?/$A?Mb?zGA?(\?8A`p?;?tV?Zd;?Pn?/$?"~jL?zG?jt?CL?V-?L7A`5?$?"~?(\?Cl'?Mb?ʡE?> ףpm?A`?n*?S?fffff?rh|/?fffffv?x? ?Zd?(\?I +?@5^I?bX9?x&1?.xVeg10t?+N+N{GZ?T㥛?-?~j?rh?lq?Zd;O ?2Z??K7A?B`"?Cl?/$&?Q?M?K?On?E?7A`?Gz>?'1Z?{G?䥛 0?n`?rh?Gz? rh?L7A`%?QE?z&1l?"~?Zd;??&1?'1*?rht?I +?/T?${?Gzn?Q>?R ?NbX9?K7?rhd??x!?S㥛 ?ˡE?%C?jtX?MbX?x&?ʡ?8A`?ҝ?zGa?V-?Mb?S㥛4?ˡE-?䥛 ? ףp ?X9v?Q?Q?i?jt ?ˡE-?̡E-?a" ?J +?-6?B`"?ʡE?V?&C?Zd;Om?> ףp=?S㥛t?M?%C ?%CK?p= ?K7a?0$?y&1L?uVN?\(?tf?M?33333? +?2Z?~jt?(\E?B`"k?uV?K7?+>?l?"~*?L7A`?C|?S㥛?EԸ?n?h|?u?Hz?J +?\(?ʡe?7A`?~j4?kt?T㥛0?Mb?n?x&?ˡE?n?X9V??#~j ףp}?43333?q= ף?'1?+َ? ?V-o?-'?v?(\R?jt?R!?x&1?~?5^?E$?zGY?KO?}?5^!?䥛 ?t?y?SM?E$?+?B`"?M?Gz?X?V-?Q?33333?sh|?V-?X9??2Z?{G?X9?I +W?fffff&?|?5^? rh?S?Mb0?;O?V-/?}?5^??5^I ?v/ ?Cl?Q?la?S㥛?A`Т?uVN??%C?p= ?Zd;O=?|?5?x&1?V-O?> ףp?&1?|?5^:?uV?On?Dli?A`"k?(\•?{G?x&?fffffv?S?CL?Zd;?7A`P?ˡE?2ZD?Dl?RA?ˡE?uVn?Zd;?B`"K?vV?m ?> ףp?+&?V-7?ףp= O?Qf?Mbx?+~?.Veg14t?+N+Njt8? ?HzW?,?Mb?sh|??MbX?%C+?(\?K7A?33333?Qe?X9&?Q?Pn? Zd?S?p= #?~jt? +? +)?Zd;O ?S?jt?_I K?Zd;?(\u?x? rh?|?5?S㥛?0?1ZD??jt?&C?PnC? ףp ?RA?RA?ʡE?w/=?Qn?^I ?|?5^? rh1?K7?sh|?A`В?/$F?"~?%C?p?J +V?la?M?!rh?(\2?S㥋?zG? ףp=j?uV?(\r?33333?Clk?V-?Ra?Mb?p= ׃? rh?bX9X?+?fffff?bX9 ?8A`(?V-5?Zd;G?q= ףX?Zd;Oe?`"i?.OVeg15t?+N+NClG?p= ?QH? rh?/D??t?Zd;?A`Ђ?C\?(\?Gz?أp= w?> ףp?@5^I|?v?/$a?;O^? rh?B`"K?@5^I|?R? rh? ףp=j?l?n?|?5^:?$?/$a??8A` ?33333?la?`?kt?̌?x&1?x/-?Cl?e;O?Gz?~j|? ףp?8A`??x&1h?&1 ?+?On?"~Z??p= C?i|?e?kt?Q?V-]?6^I ?V--??o1?Zd?K7A@? +?i|?%?Q?> ףp?x&1?㥛 ?Zd;O]?"?p= ?43333?MbX?!?&1,?rh|7?X9>?x&A?.Veg17t?+N+NʡEV?~jt?v?h|?U? +y?V-O?Zd;?zGZ?Gz>?v?E?Sk?L7A`?v?zG ?p= ?+^?bX9(? rh=?E$?1Z? r8?F?B`"[?kt?"?Cl?(\?CL?Zd;?/ݔ? Q>?K?Gz?/$?/$q?S㥛? ?|?5^?x&1?I +?X9vF?"?x?أp= ?K7?np?FX?&1?|?5?On"?-F?V}?Hz? r(?Q?R?䥛 p?{G?Onr?\(?n?ʡE?Y?X9v?h|?5? +?E? rX?̬?|G??p ?kt ?#~j ?i|? ?.?Zd;?.QVeg19t?+N+NHzG?I +?K7A??~j?r= ף?Cl?t?|?5?~jt?Hz?Q?Hzw?Gz?Onr?R?#~j ףp?䥛 ?I +?jt(?|?5^?֣p= ?bX9? rh?V-?|?5^Z?y&1 ףpE?rh|'?Q??S?gffffF?!rh?I +?FԘ?zGq?Zd;??? r?^I +?-?/$?$c?V-2?bX9h?~j?}?5^?(\?ʡE ?x&1 ףp?C|?L7A`?fffff?U㥛?ClW?o?,?Mb?Mb?V-?zGA?T㥛?@5^Il?t?uV~?/$?gffff&?x?9v?PnC?fffff?~jt?L7A`?(\?ʡE?/$?+?ʡ?~jd?-?@5^I?X9v?/$?kt?tN?V-m?*\X?Zd#?S㥛?K7A?Hz?Zd;OM?K7)?i|? ?rh|?tV?x&?/?"^? ףp=z?*\h?:v??M?E?xv?v ?2Z? rh?Gz>?|?5^ ?X9v>?Q?'1?> ףp?䥛 ?m?> ףp?Zd;/?S㥛?+َ?Hz?MbX?T㥛?-?Clg?{G?K'?oʑ??@5^I?SS?p= ??ˡE?+?Dl?Q(?T㥛Ā?.?Mb0?@5^I?ʡ?%C?I +?-g?x?'1j?Q?SK?K?Cl?T㥛0?"~?L7A`?K?33333S?V-}?-?l?~j?ʡu?v/]?SS? rhQ?.,Veg23t?+N+NF?)\?L7A`U?Gz??~?5^9?z&1?|GZ?(\b?q= ף?L7A`?$?Hzg?n?\(|?x?l ?K7?ClI?E?"~?rh|?GzN?X9V?K7!?J +?uVN?ˡE?~j?zGQ?zG?Zd?"~j,?x&Q?(\?;O?\(D?fffffv?xi?^I ;? rh?+?Zd;O?Zdk? rhA?Zd;? rh?(\?Gz?S[?:v?fffff?A`?S㥫?ClK?Mb?V-/?Gz?C?$?2ZT?L7A`E?ƛ?/$V?MbXi?vV?(\B?6^I ?x&q?v?t?ET?Q?Mb?C?^I ?l?6^I ?S㥛D?0$e?A`В?Cl?ʡ5?|?5?K7)?A`?~?5^y?E?Gz?}?5^?p= ף?ʡE?6^I ?fffff&? ףp?|?5^?Gz?oa?uV?i|??|?5^?0$?m?> ףp}?q= ף?jt?4?Ra?a"Q?l!?K7A?(\?1Z|?lQ? r(?K7?`"?t?V?FH?-?V-?v?E?Q(?|?5^?Dl?S㥛?rh|?x&?:v??Mb0?Mb?jtD?-V?/$?Zd;?8A`?L7A`5?X9v~?V-r? ?&1\?ˡE?S?^I ?䥛 P?{G?_I ?vV?p1?T㥛Ā?|?5? ףp=J??x&1?X96?X9v?x?x&q?33333s?ʡE?Pn?Zd;/?(\?䥛 ?MbX?Zd?Zd;O]?{G?nQ?x?!rhM? +?'1z??xY?ˡE?^I ?w/%?ʡ-?x/=?y&1T?@5^Il?v?/$?.Veg25t?+N+N:v?On?SC?M2?x ?ʡE?Mb?Q?K7A?v?/$v?ʡ?)\?(\•?Pn3?`?Dl)?_I  ?ʡE?K7?$?? rh?+?%C+? ?(\?Q?A`Ђ?/$?w/]?(\?NbX9?8A`?Q?,ί?33333#?GzL? :?x& ?ףp= ?ʡ?xa?S5?= ףp ?n?y&1?M?xq?@5^I4?V-?R?C?KW?Zd;?Mb?ʡE?S㥋?$c?I +?t&?fffff?Gzn?v*?B`";?p= ׃?R?.v?q= ף?~j4?ףp= '?"?&1 ?}?5^9?Y?DlI?K?&1?+Η?NbX9?'1j?A`Т?Gz??x&a?S?2Z?jtx? ?K7a?K?S?p= ?NbX9t?6^I ?$c?n ?rh?T㥛Đ? Zd[?43333C?^I K?> ףpm?V?ʡ? rh?\(?-6?R?&1:?Gz?B`"K?S?A`2?S?^I +?FԨ?B`" ?I +g?MbX?l ??Mb ?Mb(?"~2?V=?S㥛D?I +G?.{wVeg28t?+N+NCl;?x?Cl?+W?Zd;?ʡE?K7?(\R?Dl ?R?B`"K?zGa?Dl?xY?&1?ʡEF?1Z?jt8?Zd;?d;O?(\E?(\?vz?@5^I?On?sh|?o?Mb?v?:v?J +?On?,?'1?~?5^?W-??Pn?C?Mb?(\?+G?|?5^?n?zG?Cl?jtp?SU?rh4?Mb??5^I?Pn#?\( ?Cl?"~?Cl'?n?(\E?x?+G?(\?X9?(\?Zd?/?B`"?i?-? "?/$U? rh=?|?5^? ףp= ?/$&?R1?(\?q= ף?E?Zd;/?'1j?}?5^?K7!?fffff?Zd;?̡E?K?$c?n:? 2?(1J?e;Ow?V?`"?M??K7Y?i|??X9v^?}?5^?sh|o?ףp= ? ףp}? ?v?ʡE?S㥛d?d;O? rh=? rh?Gz?~jtc? ףp=??B`" ?bX9?S?\(?E?S?/$?.Veg30t?+N+N+?jt?t?On?e;O? rhQ?*\?K?E?bX9?+?V?[d;O?gffff?v?Gz^?Zd;O}?&1?jt?㥛 ? +y?Clk?/$!?vV?n ?̌?Q?Zd;O?S?~?5^?Zd{?jt?(\e?d;O?֣p= ?֣p= _?I +?{G?{G?ʡ}?(\E?Hz?K7A?ףp= ?Zd;?E\?أp= ??V-'?S?Hz?(\?Mb?&1?/$v?zG1?$~j?6^I r?Q ?Cl?I +?S㥛?NbX9?^I ?L7A`?)\?l?B`"?q= ף@?h|?? r?S?%C ?0$E?(1Z? Zd[?Y9v.?"~?R?HzW?~jt?~j?V}??A`Т?+g?RA??&1 ?p= ף ?QU?ʡ?L7A`? ףp=J?ƻ?#~j,?S㥛?(\?V?ףp= '?B`"?Q?C|? +?KG?S㥛?C ?V]?ҭ?*\?(\?(\?uV?Dl?V-??Zd?.Veg31t?+N+NCl?(\? rh?q= ף0?ˡE?R1? r?m"?L7A`%?l?C ?ʡE?7A`?(\?"~?K?|?5^?|?5^J?C ףp]?-?x&1(?~j?|?5^:?"~z?E?\(?OnJ?(\•?Q?I +? r`?p= ף(?$?C?~jt?Ed?43333C?;O&?K7A?rh?cX9X?A`В?&1?"~?-G?d;O? +?SS??#~j?S? r?Ra?Mb?7A`p?T㥛p? ?V-?Cl?ʡ5?䥛 P?X9&?+Χ?Zd;?6^I ?sh|?v?"~J?L7A`?V-?,F?8A`?(\2?S?z&1l?"~:?> ףp?R?@5^I?M"?ʡEc?Q?X9v?Se?^I ?QN?Gz?ʡ5?Dl?lA?$?NbX94?Zd;O?~jt?QK?jt?I +?K7AP?I +? ףp=?~jt?\(?J +?-?V?x?)\?.Veg32t?+N+Nl? ?~jt?{G?GzN?|?5^?"~:?Qn?p= ׃?+7?/$F??uV?V-?V-??K7i?x&Q?S㥛?Mb?p= ?Q?jt?S?PnC?ƫ?MbX?Hz?/$?xF?ʡES?"~?.V?Mb?A`?/$?Pn3?y?^I ?T㥛x?~jt@?x& ?V-?ףp= ?ʡm?QH? rh)?Zd;?M7A`? Q?Cl'?Sc?:vo?/$Q?F?B`"?R?Mb0?Zd?Hz?Zd;߿?Zd;OM? ףp?+'?23333?ʡe?\(?K?!rh-?xI?Zdk?"?Cl?q= ף`?^I ?K7?Mb ?fffff??"~jl??v*?Zd;ߏ?rh ?S?ZdK?|?5?/$?x&1?"?!rh?V-/?v? rh?"~:?h|??Gz?1Z?5^I ?S㥛t?$?Qu?"?SU?V-?X9?;?/$?K7?ʡE#?㥛 `?$C?ʡ?r= ף?Qn?2Z?vj?X9v^?Zd[?.Veg34t?+N+NZd;?V-?ʡE?x&?Gz4?S{?\( ?v_?x&1?On?xF?p= ף0?uV~?rh,?:v??I +?Zd;?Cl?A`"?x&1?\(,?}?5^?MbX?8A`?"~Z?Gz?$?/?!rh?ClK?J +?x&1?v?n*?bX9x?"~?@5^IL?V-?^I ?jt8? rh?@5^I? rx?$C{?/$?%C?Cl7?e;O?Q?K7A`?$?/$e?Pn?DlY?v?(\U?zG?M?E?la?/$?v?zG)?n1?lA? +Y?lq?NbX9?^I ?. Veg36t?+N+NS㥛?x?䥛 0?'1:?rh\?Gz?Hz?d;O?n?|?5~?jtD?T㥛?zG?p= ?MbX?x?Rq?K7?V-?|?5^Z?_I K?i|?u?ˡE?zGa?`"۹?K?ʡE?Cl'?@5^I,?ʡE?V-?Mb?d;O??Cl ?~j?zG1?z&1T?kt?(1?-?|?5?jt?zG1?B`"?+g?t? r(?Zd;ߟ?Gzw?(\?}?5^?q= ףp?Gz?xf?Mb?S? rh ?xF?> ףp]?~jth?nJ?Pn?Dl9?Sk?nʑ?Hz?^I ?x&?e;O?ҍ?:v??Gz?|?5^?J +v?Pnc?oq?Gz?n?/D?̜?R?Q{?K7?Q?n ?~jt?L7A`?i|?u? ףp= ?ףp= ?fffff?uVn??Mb?|?5^*?A`2?x&1@?}?5^Q?(\b?xn?~jts?.Veg37t?+N+N(\?Dl?> ףp?C ?.?8A`?fffff??ʡE&?d;O?Ƌ?ˡE ?On?-?Mb?|?/$?&1|?rh ףp?ʡ?Zd;_?T㥛?Clw?K?l ?'1?fffff? ףp="?أp= ?䥛 ?J +?Q`?> ףp-?K7?x?(\­?B`"?`"a?\( ףp ?S??NbX9D?jt?la?8A`?Q? rh1?Q?/$?K7i?RA?ʡe?Pn?ʡE?|?Dl? + ?MbX?i?bX9?ʡE?(\?K7y?x ?|?5^?h|??Mb?= ףp ?x/m?Zd;O?/$a?,? r?#~j?Y?QE?K7AP?S㥛t?2Z?Cl ?Gzn?On?Gz>?Dl?MB?&1?SS?Gz?_I ;?Gz?S3?nʱ?h|??Rq?~jt?/$?(\?\($?\(,?.6?K7A?lI?CL?.9Veg39t?+N+N&1?Cl?T㥛?zGa? ףp= ? Zd?䥛 ?L7A`E?Pn?Cl?MbX?V-_?Pn?A`?MbX ?K7?On? rhA?n?6^I ?Zd;?V-?0$?|G:?Pn#?/ݤ?Q ?(\o?~jt?K7?n?MbH?l? "?x& ?@5^I?K7?Dl ?w/?Zd;O?{?GzF?6^I ?#~j??5^I?B`"?Gzl?ʡM?K'?K7?|G?S?~?5^?s?ˡE?n?'1J?"~?v?(\?33333?x?K7AP?Q(?QH?"?(\?PnC?/$?X9v?~?5^?hffff?SC?B`"K?NbX9D?S㥛?T㥛Ġ?1Z?w/-?S㥛d??S?[d;O?n?\(?/$?p= c?|?5>?K79?5^I R?:v?h|??5^I "?v?~j?Sc??jth?A`?jtx?x&?MbXY?h|??33333C?S㥻?`"?/$q?C?~jt??X9v?C$?Zd+?|?5^2?rh|7?K79?.HVVeg40t?+N+Np= ף?Gz?&1?Q??K7?(\?(\?"~*?n?ʡU?_I {?$~j?*\?K7?n:?|GJ?x?PnS?Pn?ʡE?@5^IL?&1?Q>?~jtS?ʡE?vo?sh|?p= ?"~:?Zd;?Gzt?v?Q%?Q?K?X9v?Q?+?ҕ?E\?x&1(???J +?Kw?^I S?|?56?K7A?nʑ?v?$?Dl?h|?E?A`?E?x)?](?K7A?gffff?p= ף?I +?|?5.?/? "?ʡEc?Zd;?V-?p= C?2Zd?S;?L7A`?+? ?Hz?NbX9?ˡE-?On?n?Mb?QU??GzD?S?(\•?Se? ףp=:?Mb ?Zd+?OnS?v?L7A`?/$6?$?K7?n?\(?ʡ?n?Cl? ?W-o?E?\(L?+پ? +?Cl?(\?8A`?Zd;??Cl?K7A? +?Dl? ףp=?.Veg41t?+N+N^I  ?x?y&1??|?5^:?Gz?(\?33333?/݄?K79?ףp= ? ףp=?K?(\o?K7i?/?/$?R?~jt?|?5?Q(?Hz? Q~?Cl[?+Η?+N?X9?z&1L?x9?|?5n?n?|?5?%C?F(?|?5?#~jt?/$?S?\(?Ry?T㥛@?y&1 ?"?Mb?`"ہ?.$]?1Z?䥛 ?أp= G?rh?y&1L?fffff??Gzt?0$?gffff?5^I ?$?Cl?kt?(\?2Z?w/M?Q?A`?6^I B?K?Hz'?Zd?x&1h?F8?\(? rh??ʡE&?}?5^i?m?J +?Gzn?43333?xV?ʡE?=?V-?`"I??#~j ףp?x?q= ףP?sh|/?Zd;o??h|?%?ʡ?Clg?vz?M?Gz?zG? ףp=? r?K?~jt#?,G?~jtC?-'?Q?Cl?Zd;?-'?MbX?S? ףp=?SE?S ?jt? ףp=?"?sh|?M2?Gz?/$?Zd;??n?2Z$?ƛ?'1 ?|??Zd;O}?I +?|?5^?B`"?/$?KG?V-?V?ˡE-? +i?$?ҭ?FԸ?%C?C?A`r?-f?$c?.Veg45t?+N+No?Pn?1Z4?w/?na?/ݴ?v/?/T? r?S?/݄?nz?Cl?/$v? r?ʡ?I +O?Zd#?x?rh ףp?NbX9?أp= W?p= ף?Fh?Gzn?MbXy?rh?p= ףP?Pn?x&1?bX9(?Y9v^?K7q?p?"~?$C?bX9?K7!?X9?L7A`M?m?%CS?t?Gz?Gz?X9v~?ʡEV?jt,?"?v?5^I ?!rhe?NbX9?/$V?zGA?Gz?v?L?rh?e;OG?Dl?)\?$S?/$?nZ?䥛 ?:v?!rh?v/?S㥛?oʡ?t?|?5^J?ʡE&?zG? ףp=?X9v?v?ʡE?x&A? rh?Zd;?V-?@5^I?xi?V-?OnB?h|??l?zG1?y&1?x?d;Ow?A`"{?ʡE?\(?(1:?](?"rh?Ed?+?xi?J +?rh\? ?%C[? r?x?"~?L7A`e?`"?`"?rh,?(\5?L7A`E?NbX9\?Gzt?K?uV?.Veg47t?+N+NPn?ʡE?n?V-B?p= ף?n:?أp= ?Mb?cX9?cX98?S㥛?x&?Cl ?Zd;o?Q8?MB?&C?Q?x&?MbX ?|?5^?kt?䥛 ?E?zGZ? rh?On?n?Gz?S㥛?}?5^Y?K7A?C|?K7!?p= #??Dl9?"V?n:?\(?M??Qc?"~:? +?ʡE?MbX?Q~?mR?/$ ?Q.?*\?S㥛?K7?K'?uV?8A`0?K?kt?v:? ףp= ?7A`@?-?Mb?A`?^I k?K?Hzg?R?l?B`"?(\R?Gzw?Qx?(\o?jtD?jt?R1?tf?%C?V-??1Z$?R?Hz?Q%?#~j?S?(\?/$^?~j,?x?n?Dl?ʡE{?X9vV?%C+?|?5?K7?43333?y&1?x&18?/$?$Ck?y&1? r?S?-?jt?Q?+ٞ?n? +?/$&?V?(\?xF?B`"k?̡E=?!rh??n?&1?gffffv?M?MbX?V-?9v?/$?rh|_?(\?}?5^Y? rh?R?I +?xV?QE?A`R?K7y??n?lq?L7A`?lA?rh|?xI?p= ?EX?l?ˡE=?l?-7?I +?v?-v?Hz?x&1?Zd;? ףp]?)\8?fffff6?x&Q?V-?h|??5^I "?v?/?jtd?6^I ?y&1l?Cl?ʡE}?J +?ˡE]?F?/$F?Zd;?(\?x&1x?ʡE?$C ? rh?"?E$?$~j,?z&14?x9?%C;?.Veg51t?+N+Nv?ʡE?ʡE?T㥛?M?ʡE3??Ed?~j?V-?~j$?/$?Zd;O-??rh|/??~?5^y?r= ף?6^I "?K?(\/?y&1?Y9v?Mb?x&?L7A`?p= ׃?zGA?~?5^Y?\(?rh|O?5^I ?v*?L?bX9?rh|?*?Mb?Q?K7?nX?Q%?{G?V-?w/?zGr?zGJ?`")?rh?Mb?Zd;?Dli?K7A0?|?5?i|??Q(??p= s?Clk??`"ۙ?Zd;O?Sc?p= S?S?jt?#?|?5n?%C?/$?Q?E?^I  ?(1?#~j?̜?S#?|?5^z?{G?a"?B`"K?Cl?jt8?R?> ףp?? ? ףp=?"~?$?:v?(\?J +F?p= C?I +?hffff?$~j|?S%? r?Gz?Cl{?T㥛Ā?vZ?Zd;?1Z?Mb?Cl?Q>?&1Z?vO?5^I ?jtl? ?K7?NbX9?/$U?C?Pn?bX9? ?ktd??y&1??On?/?S%?2Z?K?m?T㥛0?T㥛?Hz?^I ;?:v/?> ףp? +Y?MbX?^I ?Mb?/4?+w?S㥛?.?rh ?x&?Q?K7A0?lB?Gz.?jt? rh?tF?M?ˡE?S;?> ףp?y&1?V-?lb?/4?X9?^I ?jt?~jt#?Fh?A`в?J +?Gzn?43333?أp= W?I +?n@?-?tVN?;O?Gz>?d;O?uV?J +V?Zd;߿?~j?+^?~j?t?X9v?I +?K7?ʡE?n?)\? +?.7IVeg54t?+N+NHz?MbXi?"rhM?v?S?E?J +V?&1,?Zd;?K7?|?5?$?Zd;O ?Gzn?q= ף ?B`"?FԨ?Zd;O?= ףpM?Zd;OM?,?e;O?K7A?M"??Q?Cl? ףp=?h|??y&1 ףp?F(?K7A?rh ?X9v~?ˡE?q= ף?B`"?Mb?v?1Z?MbXI?Pn?ʡE?Zd;/?k?ʡE?n?Q?Zd;O?|?5?(\u?Clk?Dli?.޴Veg56t?+N+Nx?"~?Q%?a"۩?"rh?K7?|?5^J?(\B?T㥛?GzT?Zd;_?ףp= ?> ףp?B`"?+?fffffV?Qm?V-"?K7A?NbX9?33333?jth?V-=?Q?d;O?/$A?}?5^?l1?Gz?d;O?#~j?1ZD?I +?B`"?M?d;OG?Zd?Q??/?&C#?ʡ?%C?J +?Mb?~jt?^I +?MbH?y&1l? ?gffff?S㥋??E?> ףpm?X9v?uV?Q?Mb?C?\(?zG1?y?E?vJ??K7I?p= ף?S3?~?5^?.6?`"? rhq?bX9?FX?/$?t?p= ף(?5^I 2?/D?/$a?I +??&C? ףp=?K?Gz?NbX9$?33333?#~j?$?"~:? +?Q?e;Ow?~?5^?`"?n!?K7A(?K7A0? :?i|?E?L?(\O?.Veg61t?+N+NB`"?Mb?Cl?S{?B`"?;O?x?On?V-O?I +?v/?䥛 p?"~z?T㥛Ġ?ʡE?ʡEs?v??9v/?5^I ?A`?~jt8?Zd;O?-阮?n!?Cl[?V-O?A`B?sh|??x&1x?p= ?-淪?X9v>?ʡEs?tn?q= ף8?/$?(1?Zd;?jt?"~?K7Ap?V=?L7A` ?~jt?m?S?Pn[? 2?K7?#~j\?k?a"I?2Z?Q?V-??#~j?Zd{? +9?QH?X9?ʡ? Zd ? ףp?V-?"~?MbX?3?Hz7?m?43333?Mb?y&1?h|?u?S?7A`?= ףp=?Q?+?Mb?Q?I +W?Q.?)\?+'?lR?Zd;O?`"?Q;? +?أp= ?.阮? ?{G?~j? +?+?V-r?أp= ?M?n?X9v?p= s?v?rh?m?-? ףp=?Zd??t?uV?.aVeg63t?+N+N%C? ףp=z?fffff? ?x?mR?x&1?ʡE?E?n?+7?ƫ?m?SE?nP?)\x?~?5^ ? ףpu?vV&?)\?Q?x&?v?(\m?333333?MbX?K7?fffff?Zd;?Zd;OU?|G*?(\?Hz?+?l?M?? +I?^I ?E?Zd;??M?(\?v?I +w?}?5^?|?5^z?7A`?p= ?`"?~jtc?n?K7?kt?\(|?zG?V-?S㥫? rx?K7)?&C?t?y&1L?㥛 ?K?[?jt??MbX?K71?(\?~j?S#?y&1\?ҝ?tV?2ZT?.?Cl7?K? "?"~?3?Hz?+?|G?cX9?d;OW?E??&1l??5^I?X9?n ?I +?fffff?n ?t?$~j??.Veg64t?+N+Nd;O?T㥛`?"^?_I K?2Zd?2Z$?L7A`E?On?%CK?Mb?xy?\(?L7A`e?R?n?/$?Vm?~jt8?ʡE?Qx?V-"? rhq?M"?Clg?K?J +V?y?ҍ?sh|?~j?S+?|?5^?M?t?5^I ? z?S?E? ? rP?Gz?v/?V?x?A`b?l9?6^I ?X9v?vz?q= ף?i|??x?Cl{?|G?zG?lq?x&!?~j?T㥛?Clk?(\X?Dl?d;OW?Zd[?^I ?n?|?5^*?QU?Zd;o?Kg?jt?B`"k?Sk?$C?Zd;O ?ƻ?=?Zd?K?~jt3?~jt?zG?Mb?ClY?&1,?sh|?ʡ?(\? "?xi?/ݴ?Dl ? r?jt?#~j\??ktD?x?/$Q?zG?㥛 @?K7?x&?^I [?1Z?d;O?p= ף`?d;O?v?p= ?C?x?Q?> ףp?`"?MbX?.GVeg65t?+N+Nʡ? ףp=J? ףp ?Ed? +?MbX?rht? r0?Q?x&1?FԘ?S㥛l?X9vF?sh|?jt?V-?`"I?Zd;?X9ȶ?v?GzN?Cl?ʡ?#~jL?8A`?Q?rh|?"~J?-7?X9Ȇ?T㥛0?o1?Gzn?Cl?C?|?5^ ? +? ?t?ʡ?i|??](?(\?v/M?;O?0?(\?X9?Mb@?)\?{Gj?bX9(?Cl?jt? ? rh?-'?^I {?|?5?x&?nʑ?"~ ?p= ף?"?Sc?h|??^I k??0$U?*\?,?I +W?x&1?fffff?^I K?Q?`"۹??uV?Gz?ҽ?V?d;O?x?._Veg66t?+N+N ?p= 3?C?Zd?n?&C?Zd;?%CK?I +?K7A?%C?? B?|?5^?t&? ?,?rh?d;O?nʑ?K7A?C?p= #?V-?/$?p!?(\O?Zd[?> ףp}?"~?-v?T㥛? +?Cl?rh?rhL?Gz~?ʡE?㥛 X??d;O?l?'1z?M?MbX)?v/?jt?T㥛Ā?K7?X9vn?x?vo?K7A ?rh?'1j?K'?5^I ? r?9v?fffff&?T㥛?w/]?(\?S?+>?8A`?Q?= ףp?zG? +?+N?zG? ףp=?zGq?S5?X9v?n`?/$?Zd;?> ףp}? rh?S㥛?R1?@5^I?Cl?|?5?Q?Dl?n0?ҍ?ʡE?T㥛@?> ףp?*\(?/$?+?x&?q= ף?S??NbX9d??h|??ZdK?fffff?R?jt4?np?S㥛?t?/$?~jt?~j?^I {?+n?k?.]Veg67t?+N+N~?5^?l?E?-? ףp=?S?Mb`?&1*?L7A`?t??K?$c?xY?x&1?Cl?ƃ?B`"?x&?j|??X9?/$?R?x&?C?V-2?"?"~?Zd;??A`?+ٮ?V]?Dl?33333;?x&A?kt?`"i?~j?c?Zd;O-?ʡE?l?I +?T㥛p?(\B?^I ?tV?Q?ClK??rh?;O? rhq?֣p= ?@5^I?V- ?+Η?MbX?$?Gz'?Cl;?(\µ??Gz?+?Ƌ?Gz$?S?rh?v??M?^I 3?x&18?x&1(??"~?x&?-?(\"?+>?K7A`?ʡE?V ?fffff?,?W-?أp= g?K7 ?Q?ҍ?-?&1?V?jt4?+~?Zd?Zd;O? ףpM?^I ?SE? rX?{G2?-?Q?~jt?-^?+7?\( ?`"?~jt?v_?x&?@5^I?rh?Zd;Om?v?t?lA?v?.f?̡E-?QX?Dl?ʡE ?K?ʡU? rhq??x&Q?ʡE?I +G?|?5^?zG?7A`p?Dl?^I ?~j?I +w?K71?&1? +?T㥛P?NbX9t?nʡ?7A`?Zd;O]?Q?_I {?"rh-?kt?9v?A`r?S㥛d?ʡEv?`"۩?S?|?5^J?M?Cl?l?Mb?V-?X9v?I +?{G?C|?V-?X9v?&1 ?s?Pn?kt?gffff.?x6?@5^ID?U?-f?^I s?d;Ow?.Veg70t?+N+NV-?+~?(\?+'?V?Ra?Pn#?S㥋?ʡEf?~j??+.?"?23333C?bX9?Gz?K7A?la?bX9?J +?+ٮ?33333?Qe?^I ?rh?㥛 ?Zd?rh?nZ?~j?!rh?+7?Ss?Gzv?أp= O?Zd;? rh1?-???Q?(\?K7Ap?SC??x&? +?[d;O}?A`J?Pn?~jt?.v?ClG?Hz?X9v?Gz?M?䥛 P?|?5^?NbX9D??(\?!rhm?#~j ףp?Q?{?V-M?-?jt?+w? rh?x&1?&1?ףp= G?&1?@5^I|?\(,?M?jt?X9?rh|o?tV~?K7?\(?ʡE?F(?x?S㥛?p= ?X9v?0$?o?x&1?Cl?|?5?y?|?5?&1?`"?OnB?ʡE?%C?X9vN?C?ʡE?:vO? ףp=*?Dl ?R?(\"?&1\?:v?8A`?gffffV?Q? ףp=:?Cl?-'?v?Cl7?}?5^?C,?_I ?F?QX?x?rh|?(\o?"?x?n ?d;O?J +?V- ?ףp= ?Gz?S?.uVeg75t?+N+N{G?xF?6^I R?Mb?*\?"~j?MbXI?~jt?i|??&C?Zd;O?zG?m?̡E}?B`"[?ˡE=?ʡE?Q?0$u?ʡE#?NbX9?d;O?jtx?\(?Cl?'1:?L7A`?Pn?Dly?h|??i|??-?K7A ?#~j?Q+?Q?+η?Zd;O?gffff?K7I?I +?Hz?MbX?ףp= ?`?Cl3?v?S㥻?I?R?v? rhq?Cl?V-?h|?U?"~ ?ʡ? +?&1? rhQ?rh|_?Zd;?X9Ȗ?+ٮ?v?xY??S㥻? ףp=?-?Gz^?6^I ?8A`?/$V?;O?K7?x&1?l?"? +)?/$?A`?Zd?7A`P?ʡE&?Dl ?/$?Pn?o!?MbXi?/ݴ?Dl ? r?jt?#~j\??ktD??RQ??S㥛D?V-?p= ?Cl[?Gz?d;O?p= ף`?d;O?v?p= ?C?x?Q?> ףp?`"?MbX?.\Veg76t?+N+N/݄?@5^I,?Cl+?/$? ףp?ʡE?K7A?S㥛?E?1Z?Fx?V-?w/?Kw?l?Mbp?fffff?Hz?v?8A`?(\o?rh?NbX9$?+?> ףpM?0?'1?ʡE?h|?m?h|??Dl?K7q?Dl)?+?jt?43333?MbXi?B`"C? r?$C?\(?v?mb?Gzn?SE?> ףp?R?x&1(?S?$?MbXy?}?5^?v/?{G:? +?Zd;o?X9Ȇ?S?(\%?$C[?K7q?Qu?&1Z?|?5^?Dl9?~jt(?ˡE?ʡE?lB?R?fffff&?^I {?uV?(\5?V?䥛 `?Mb ?Mb?> ףp?> ףp?n?ףp= '?{Gz?|?5?x&?nʑ?"~ ?p= ף?"?Sc?L7A`?k?Q?MbX?|G?jt?KW?x&1?J +?K?> ףp?"~j?S?v?S㥛?ҽ?V?d;O?x?.Veg77t?+N+N(\o?p= ?'1 ?Cl?jt?M"?"~z? ?zG?^I ;?t?~jt?V-?T㥛? r?E$?Q>?Mb? rH?w/=?Pn3?J +6?Zd;?x&1H?I +g?S?&1|?NbX9?`"?SK?V-?Pnc?(\5?p= ׋?Q?X9vV?Q?d;O?T? ?㥛 ?Cl?-v?ClK?X9v&?@5^I?"~?Zd[?Cl?K7A0?@? ףp=?V-?mR?C?= ףp?A`Ђ? ףp=Z?)\x?rh ?-?43333?ED?K7Y?V-?V-?x&1?R!?v?x?Zd;ߏ?Zd;O?̼?l?ףp= 7?NbX9?p= S?C?Q?Qu?Cl?Dly?jt(?h|??n?F?n?K?V-/?V-?ʡE?T㥛@?!rh?~?5^)?Pn?V-?E??/$?(\?ʡEf?(\?K?rhL? r?Gz?t6?zGq?ƛ?t?/$?~jt?~j?^I {?+n?k?.xVeg78t?+N+Nr= ףp?j|?U?|?5^ ?B`"?Dl?L? rhq?\( ?ףp= W?~jt?(\%?$?p= ׳??Y9v?Mb?w/?I +?{G:?S#?%C?%C?x&1?!rh?Vm? rh1?"?"~*?(\?Zd;o?䥛 ?أp= ? rh?(1Z?S?.? Qn?p= ׃?uV^?gffff&?n?(\?&1?ʡEn?Gz>?5^I ?+ί?)\p?K71?t?̜?mb?Dl ?fffff?}?5^?On?{G:?V- ?MbP?A`?$CK?Cl?"? ףp=?-g?rh?E?v?Mbx?:v?> ףp%?֣p= O?23333K?n1?~j?ʡE?m?~j?l?".?㥛 P?x?/$?X9v~?6^I ? ףp=?Onb?J +?v?V-?-?C?\(? rh1?|?5^z?x?x/M?V-?DlI?8A`?NbX94?#~j?B`";?uV?ktt?}?5^?}?5^Y?m?(\?Dl)?Zd3?|?5F?Dla?d;O?.淪?(\?.Veg79t?+N+N$~j|?Zd[?Q? +?S㥛?Pn?8A`P?On"?/ݔ?~j$?%Ck?'1?NbX9? rh?x&1?#~j? ףp?&1?V-?nq? rx?/$?ˡE? ףp-?zG?/$?K?@5^I\? ףp?I?Q?p= ף?Mb?v/?(\?Zd;O?-阮?h|??+?x&?x?Onr?+ٮ?K?|?5^j?x&?^I ?C\?-?Mb?43333C?:v?&C?ʡE?X96?Cl?rh|?K7?V--?d;O?J +?Qn?RQ?~jtS?7A`p?@5^I?v?Zd[?%C?gffff&?kt?#~j,?p= ׳?fffff6?+ٮ??&1?"?-v?\(L?n@?lQ? +y??n?oq?/$?$C?5^I ?ZdK?~j?}?5^Y?ʡE?zGA?ʡ?y&1?K7A?(\?43333c?Zd;O?ʡE?&1|?QK?@5^I,?䥛 0?T㥛P?S?"~?Cl'?䥛 ? ףp=?Fh?ʡE?p?n??Cl?"^?Q?vJ?ʡE?6^I "?K7y?~j?&1 ?y&1?V-"?Hz'?|?5.?/$5?43333;?w/=?.Veg84t?+N+NK?~jtS?z&1?2ZD?ףp= ?^I  ?`")?Zd;?/$?X9v?p= ׳?(\?ʡE?$S?ʡE? ?ףp= ?MbP?|Gj?Q?S?Zd;?V ?x)?x&?X9v~??Gz?Zd;o?l?Gz?|?5N?jtH?"?FP?&1?/$?x?/$?S}?jtD?@5^I?䥛 ?䥛 ?ʡE?a?On*?l?43333s?&1?V-?xi?Cl ?n?S㥛4?/?2Z?-?S?`"y?"~?Zd;??-&?QU?Mb?I +'?rh|?7A`?v?rh? rha?\(?[d;O}?K7?Q?E?(\?"~Z?|?5?bX9? "?V?E?+η?x&1x?|?5^J?l"?7A`?jt$?RQ?v/?`"?Q;?K7?,?jt?h|??̌?/$?K7?+?V-r?أp= ?M?n?X9v?p= s?v?E?S?8A`?Q?2Z?w/?Gz?uV?.dVeg85t?+N+NQ?zGA?~?5^i?Zd;?GzN?|?5^?fffff6?I +'?OnR?8A`?jt?/$?uV?6^I "? r(?(\? ףp=:?x??i|?E?t?GzT?x?\( ףp?NbX9T?X9?x&18?MbX?t&?tV?K7?}?5^?C,?_I ?F?QX?x?(\?;On?V?x?n ?I +?X9v?V- ?ףp= ?Gz?S?.FVeg86t?+N+N|?5~?V-?QE? +y? ףp=?Cl?J +F?`?V-?K7Q?NbX94?m?e;O?w/?t?X9Ȇ?)\?7A`?K71?Dl?ʡEc?"?J +?rh\?m?Q+?K7A?x&1??^I k?7A` ?䥛 ?B`"{?On?Cl;?S㥋?$?E?~?5^?V-B?K?)\?E?&1?&1\?Cl+?6^I ?NbX9?/$?8A`P?CL?K?.臨?/D?I +? ?Sk?6^I R?v?tV>?Vm?ʡE?v?MbX ?/$f?-?Cl?㥛 ?l!?Pn?#~j?x&?n?-g?43333?n?Gz$?Gz?m?Zd;O?MbXy?L7A`?(\?I +G?ˡE? ?S?K7A?7A` ?K7i?/ݴ?Dl ? r?jt?#~j\??ktD?x?/$Q?zG?㥛 @?K7?x&?^I [?1Z?d;O?p= ף`?d;O?v?p= ?C?x?Q?> ףp?`"?MbX?.5Veg87t?+N+NDli?Q? "?rh?䥛 ?:v?V-]?tV?A`?Zd;?Mbp?p= ף ?A`"?zG?V]? ףp=?ףp= ?^I [?@5^I?KW?n?zG?l?$C?Q?Zd;ߏ?K?On?Gz?~?5^?~?5^I?ҽ?Q?I +?ʡE.?Zd;Ou?/$?l?rhl?^I #?-?Pn?rh|?+f?@?6^I ??-g?Cl?8A` ?(\?S?x&1?~j? ףp=?= ףp}?x&1L?$3?`"i?Zd?QH?fffff?ʡ?V-? ףp=:?L7A`?"~?K?Gz?K?v:?Sk?MbXI? ףp?"~j?S?v?S㥛?ҽ?V?d;O?x?.Veg88t?+N+N~jtS??A`? r?rh|??Mb?Mbx?V?PnC?w/?Zd?)\8?ʡ?Qk?%C?ʡEf?Zd;O?Pn?x&1?+?S?n?{G?jt?/$?K?B`"?X9v^?!rh=?V-?Zdk?L7A`?q= ף?nʩ?Q?!rh]?ʡ?Zd;O?X9vN?M?-?= ףp?{Gr?RI?43333#?V?K7A?S5?L7A`?r= ף?~jt?V-?V-]?x?Mb?fffffV?x&1(?zG?xF?+?#?(\?Q?p= ׳? ?Zd;_?MbX?jt?ʡEs?/D? ?T㥛?uV?ʡE?K79?Q?fffffF?v?S ?x&1h?Dl?w/m?V-?\(?x?T㥛?K7?R?,?̌?/$?T㥛@?> ףp?*\(?/$?+?x&?q= ף?S??NbX9d??L7A`?](L?Q?Gz?t6?zGq?ƛ?t?/$?~jt?~j?^I {?+n?k?.OVeg89t?+N+NGzT?MbX9?K7?x?R?Gzt?I +w?kt?$?̡E=?Pn?/$?K?!rh=?T㥛p?/$1?M?V?Cl?n?S;?jt4?/$A?{G:?أp= 7?#~j ףpE?Clg?"~Z?Zd;7?l?Mb?U㥛P?'1?#~j?Cl?Cl;?M?V-?Rq?bX9?6^I ?^I [?Mb?S? ףp=?(\…?^I ?'1?+.? rx?t?@5^IL??ClG?+پ?(\2??.6?`"?zGq?-?xV?X9v?V?K7A(?Cl1?/D?/$a?I +? ףp ?w/?䥛 p?1Z?HzW?I +? ףp=?zG?#~j,?{G?Gz?I +?jtd?NbX9?J +?&1?8A`p?x&1x?L7A`?y&1? ףp=:?rh?rh|?Clg?)\?Qh?jt?^I [?R?B`"[? r?x?"~?L7A`e?~?5^?x?+?Gz4?#~jD?p= [?#~jt?K?uV?.2Veg91t?+N+NDli?tF??xy?!rh}?Q?n@?K?K7?S㥛4?\(?K7A@?$C?cX98?ʡE?Dl?ʡE?Pn?I +?S?+N?v?ʡEv?Cl?+.?> ףp}?p= ף?y?h|??Zdk?S+?~jt?ʡE?|?5?Sk?K7?{GJ?GzT?x&?MbX?M? ?)\X?v2?|?5^?ʡE?Zds?Pn3?^I ?rh ?S㥛?+ٮ?K?~j?/$a?Zd;?(\µ?? rh??Q1?X9v?T㥛?-? rh?"rh-?S㥫?_I  ?(\E? rH?ʡE?Q?أp= ?V-?T㥛p?v?|?5??> ףp-?/T?5^I ?5^I ?ʡEC?^I ? ףp=j?Gz?J +?~j?}?5^i?@5^I\?|?5n?On?:v?ʡE?+ٞ?h|??䥛 ?t?^I ? ?Gz?x&1?vz? rh?tV?C ?33333s?x&?n?E,?w/5?DlA?QS?rhd?q?/$v?.LVeg92t?+N+N +y?|?5N?x&?FH?أp= 7?q= ף?zG1?Q;?x?x&?X9v?"n?xV?+'?MbX???Gz?Zd;?On?K7?y&1?Cl{?\(?ʡE?@?㥛 ?ˡE=?th|?6^I ?Y9v?nq?"~j ?V-?S?y&1l??(\5?S㥛 ףp?"?!rh? ףp?v??$c?|?5?33333?Gz?'1?'1z?:v_?S㥛t?(\r?(\??㥛 ? +Y?Hz?l?K?㥛 P?8A`?S?@5^IL?Zd;O?%C?vj?rh? ףpM?nʑ?J +?rh?Zd;O?ʡE?bX9X?h|??\(?ktD?"?@5^I?]?|?5?bX9?1Z$?Zd;ߟ?2Z4?~jt?/ݤ?Clk?5^I B?Hz7?'1J?mr?/ݴ?V ?V-o?Gz?T㥛@?> ףp?/$F??L7A`U?+?B`";?|?5?jt4?p= ׳?jt?Su?e;O?䥛 ?sh|?ʡ%?S-?)\8?C?PnK?> ףpM?.ZfVeg94t?+N+Nx?!rhM?x?Gz? Zd{?S㥛T?fffff&?+َ?x&1h?ˡE?x&a?l?Vm?T㥛?~?5^I? ףp=Z?fffff?Hz?m?(\?v?ףp= w?w/?%C?n? rh1?h|??K7?n?3?jt?FԈ?!rh}?X9v?K7Ah?+?l ?v ?e;O?$~j?Q[?NbX9,? +?n?[d;O?MbXq?F0?"?~jd?l?)\x?3?V-?fffffV? rh?p= ף?|?5^i?V-_?ʡE?tV?p= ?Cl?ƛ?x?x&1H?V-?Zd;?/$U?|?5n?!rh]?y&1?+?|?5?~?5^?+N?Qn?~?5^Y?L7A`? ףp=?5^I 2?2Zt?(\?fffffV?l?bX9?Rq?V-B?i|?%?|?5^*?vJ?&1|?2Z?R!?n?/$?h|?e? ?$Ck?/$?%C{?Pn?'1Z?.$?jtD?X9v?Zd;?x&1x?ʡE?$C ? rh?x?S㥛$?~jt+?n2? r8?"~:?.Veg95t?+N+NB`"{?B`";?m?M7A`e?:v?|?5^*?/$&?%C?(\?,g?\(?$?MbXi?8A`?Cl?E?/$?\(,?V-?Clg?> ףp?zG?^I k?(\?y&1?Fx?ʡES?l1?#~j?sh|o?n0?,η?ƛ?K7?]? rh?x? Q?|?5^?v/u?;?L7A` ?v?33333?V-?|?5V?7A`? ?Hz7?$c?&1L?MbX ? ?y&1,?I +?K7A?X9vN?SE?n?䥛 p?(\?I +w?Qu?(\?:v?kt?t?> ףp ? +?T㥛?%C?Mb?+Η?n@?M?~jts?bX9?h|?E?= ףp?$C?ʡ?W-?*\?Gz?:vo?V-B?{G?x ?̡E?VM?FԈ?Q?/$6?ʡE?On?43333?p= ף?&1?S㥛? +?+?V-r?أp= ?M?l? rh?/$q?L7A`?}?5^?K7A?Gz?Mb?cX9?K7?m?Q?.=Veg96t?+N+Np= ףp?'1*?n`?{G?Gz?bX9?K7A0?/$?> ףp?J +?v?gffff?ףp= g?V-?ˡE?Zd;OM?43333C?l?(\?jt?S㥛T?Cl?xI?p= ?Gz?㥛 ?V-?S㥛?t?K?Qe?/?(1?S?ʡU?ʡE?S?NbX9?J +?xV?#~j?,?C? rh?A`"s?S;?Q?S?d;O?3?䥛 ??Qu??v?"~j?֣p= ? ףp=?ʡEV?(1?n?:vO?V-?S㥛4?ʡ?n ?I +?0?~jt?x&?/݄?V-?V-R? rh?Dl?`"i?x?43333?أp= ?Q ?On ?/$?S?m?Pn?.GVeg97t?+N+NMb?t?V-=?Q?Q+?rh?V-=?h|?%?MbH?33333? r?X9&?|?5^?zGZ? ףp=*?L7A`?p= ף?S?`"??{G?ףp= '?K7A?Q?/$?:v?q= ף?MbX?:v_?x/?$?K7I?Pn?ʡE?Zd;/?k?ʡE?n?Q?Zd;O?|?5?(\u?Clk?Dli?.Veg100t?+N+Na"9?V-?B`"?E? +9?e;O?tv?ʡ?a"ۙ?%C?t???5^I ?v??On?.$?,罹?'1?1Zl?p= ף?jtp?L7A`}?Zd;O??T㥛?$?Gz?J +f?F0? rh?Dl?ʡE>?"?x&?z&1?|?5?I +'?ˡE?F8?d;O?zG?C?Zd;?!rh?#~j|?S[?~jts?{G?9?-?Clg?ʡE? r?1Z ?;O^?|?5^z?ʡEf?~jt;?Gz ?+ٮ?7A`0?rh?t? ףp= ?Zd+?6^I r?T㥛?S㥛d?S?S㥛?|?5N?/$?̬?R?V-}?$?p= ?K7)?V-r?Q?ED?/$?nA?jt?B`"+?K7?X9v.?l?K7i?n?q= ףP?F?R?C$?ˡE-?T㥛@?rh\?Pn{?"~?$C?.3Veg101t?+N+NfffffF?jt(?o?Cl?#~j?+N?KW?V?|?5^ ?n`?VM?|?5^:?jt?T㥛İ?n?ʡE&?ףp= g?p= 3?5^I ?-? rhQ?x&? +9??zG?I +g?J +?Zd?7A`?zG?KG?Cl?zG? ףp=?43333s?-?.^?-g?jt8?Q?ʡE?d;O?Clo?oI?fffff?M?ʡEv?L7A`-?W-?T㥛?x?Dli?o?x&1?h|??}?5^?&1z?Ex?$C?Q?/$f?S㥛D?V]?/$?jt$?^I ?K7A@?oʡ?v?I +?\(?K?+G?Mb?Q?w/m?Zd;ߟ?&1?\(|?K?x&1 ?|?5?xI?n!?@5^I?Zd;O?Mb?Et?2Z?v?&1 ףp?-?n?6^I ?E$?= ףp-?x&1?0$u?z&1?RI?:vO??](|? rh?C?K7?h|??+n?(1?Pn?(\U?tV?v?cX9? rh?Zd;o??+?+n?i?> ףp?xv?I +?/4?Zd;7?v/?@5^I?2Z?V]?3?gffff?S?v?zGI?V?zGq?K7?%Ck?@5^I?V-?)\8?= ףp?:v?|?5N?SS?33333?^I ?1Z4?v?K'?(\r?V-?Mbx?.?Zd+?CL?|?5>?B`"?{G? ףp?`"ۙ?7A` ?+Χ?{G*?S㥛?ʡ?l?zG!?x&1?x&?V-r?K7?Zd?Q#?V-*?2Z4?V-B?P? ףp=Z?ˡE]?.,RVeg104t?+N+NQe?X96?Pn?Cl'?cX9?MbX?(\?+?> ףp?C?|?5?Gz?"?Q~?K?}?5^ ?+?a"I? + ? ףp= ?L7A`?Mb0?S;?ʡE?|?5^)?ףp= ?K7A?ʡE?[d;O?&C[?8A` ?&C?K?jt?\(t?Mb?|?5??n?w/?rhl?+>?|?5? r?8A`?&1|?(12?r= ף?bX9H?-g?lA?On?-阮?Zd;O ?(\¥?+^?$3?S;?^I ?7A`?-?q= ף?bX9?~jT?Zd;O?V-O?6^I ?`"?zG?V?1Z? r?l?n ?ʡE?Zd;/?2Z?ʡE?(1J?Zd;O}?Hz?l?Mb?,'?Q?Dl?oa?K79?䥛 0?C?Gzn?n?x&1?i?ˡE?Zd;?~?5^?OnB?"~?GzN?X9?l2?t?Cl+?S㥫?Mb?V-m?(\?x&1 ?Zd?K7!?o)?/4?Hz??"F?MbXI?.XVeg105t?+N+N@5^Il?I +7?W-?E?> ףp?uVn??ףp= W?S?`"?+N?uV??pa?x/??jt?ʡ?@5^Il?X9Ȇ?2Z??&1?|?5^:?}?5^I?S;?$C?W-O?th|?v?ʡEc?,?/?S?np?䥛 ??Gz?(\??Zd;O?x!?p= ף?S?S?v/e?/$?ƻ?n!?Q>???v_?(\?V-?RA?](?L7A`%?(\?֣p= ?Gz?(\?y&1?X96?w/?}?5^)?x&?{G?ףp= ?V?x&1?MR??X9v? B?`"?@5^I ףp?+?Cl?.@Veg107t?+N+N`"Y?ףp= ?%C[?rh?㥛 ?K7!?#?~j?-罹?:vo?K?T㥛0?8A`? ףp?$?@5^IL?t?}?5^ ?V-R?-?|?5^?~jt?@5^I?K?`"?V-?Q?"~?.NVeg109t?+N+NʡE6?_I ?أp= ?ˡE]?V-?K7?h|?E?S#?">?|?5^ ?h|?E?$s?\(?ʡE?x&1?Qh?v?/?l? +?~jtS?jtD?Zd;ߏ?Mb?S;?Gz?:v?z&1 ?uV?X9v?%C+??V-?Q?C?S㥛|?K7?n?Q]?ˡE ?A`?MbX? ?S㥛\?1?R?jtx?Gz?(1j?-?w/m?/$?M?أp= 7?Zd;?h|??ʡE?i|???/?Zd;O?-G?Gz?F8?,W?x&A?n?J +?v ?ʡE?8A`?Q?x?K7?rh?x&1X?V?~j$?B`"k?X9?M2?OnR?9? + ? ?cX98?I?Mb?K7A? ףp=*?K?^I ?FԈ?Zd;O?vVN?r= ף?$S?t?K?.$?C?fffff?M?K7A ?M?#~j?J +6? ?+?Mb?J +?Fx?d;O? rx?x&?d;OW?E?Zd?/$A?|G?(\?Hz'?$c?Gz?Q?S?Q?t? rhm?ʡEc?na?.HVeg111t?+N+NMbX?R??S㥛?n`? r?nq?ʡE?ʡE?K7i?NbX9?/4?V-?"~?*7?J +n?n?K?Zd;/?m?x&Q?+?"~ ?d;Ow?Zd;O?vo?d;O?V?Zd;O ?NbX9?ʡEV?Zd;O?Q?/ݜ?"~r?sh|?2Zl?J +v?J +F?n ?,?S?S㥛?X9v^?Gz&?q= ף?|?5^r?jt ?ʡ?-?tv? ?ʡE?K79?|?5^?Qx?ףp= G?)\X?NbX9?n?6^I ?p= ף?C?V-?K7?@?B`"?Zd;O-?q= ף`?V-*?+v?Q?ʡEn?`"9? ףp=?"?K7A?ʡE?ʡE?MbX? r?Gz^?;O?L7A`U?Y9v?x&1?6^I B?Dl? ?}?5^y?Kw?Gz?X9v?Gz$?;On?^I ?sh|??|?5?}?5^9?Mb?m"?MbX?ʡE&?x?8A``?d;O?J +F?Gz?/ ?Mb ?o)?&1 ףp}?NbX9?.!Veg113t?+N+N~jt3?S㥛?'1?)\?K?-'?NbX94?M?"?.阮?Zd?(\??}?5^?Cl?+?x&1 ?}?5^Y?Q?S?%Ck?ʡ?Zd;?V-?th|?K7?(\•?(\?{G?fffff6?t?|?Zd;o?v?@5^I|?Q?ˡEE?jtH?\(??V-?jtp?RI?V%?n?cX9?䥛 H?/?y&1L?̡E]?S+?p= ?Gzd?n?1Z?䥛 @?%C?M2?S㥻??Mr?(\o?Zd;O?zG?x&q?zG?~jt?%C?Dl?/?\(,?Cl;?Cl?ʡE?p= S?MbX?V-R?&C?/$?l1?\?%C?8A` ?Zd;?PnS?x&1?l?(\?&C[?x&Q?Ed?MbX?0$?(1:? ?F?Pns?Hz?Zd{?Q?ktt?,? ףp=j?/$?ʡE?\(?A`b?T㥛?(\?$?E,?o9?SK?~j\? rhi?ʡEm?.שVeg114t?+N+NktD?+?V-?i?䥛 ?_I ?l!?o?9v??y&1 ?](?o?Gz?K7 ?V-?n?Gz?NbX9?ZdK?ʡE&?7A`?`"?uV?d;O??`"?ʡ?Cl?"~j?rh|?A`B?.?Zd? r?Zd;O}?V-?V-/?:v/?|?5^? r?K7?QS?#~j,?Zd;?Q?K7A?l2?Q?L7A`%?Gz4?x&?'1?Q;?A`?(\_?/$!?S?X9v?X9Ȧ?Cl?ClW?V-R?sh|?Gz?$S?Gz?.F?Zd;?2Z?Dly?#~j?8A`? r?Qh?S㥛?e;Ow?&1?~jtS?:v?K7A?o?(\U?Zd;?~jtc?Zd;O ?Hz?Hz?1ZT?V=?RA?Mb`?Cl?Q?^I K?n?h|??l?bX9?V-? "?#~j?> ףp ?x? +??)\?Dli?/$?Zd;?ףp= ?Q%?".?Q=?QK?GzV?`"Y?./Veg115t?+N+NRQ?(\%?uV?"~:?Zd;O=?0Z?Gz?/?(\?(\U?(1 ?~?5^?v?Gz?0$U?@5^I?ClG?-'?jt??V}?(1 ?A`В?R?p= ׃?F?X96?!rhm?I +?T㥛?|?5~?.?-?"~?NbX9T?|?5.?ףp= '?S;?Se?Q?? rha?(\?433333?K7? ףp=:?M?SE? ףp?$~j?Hz?`?T㥛p?Zd;? ףp=?S1?8A`?6^I "?L7A`?33333?ʡE ?Zd;O?K7?K7? r?On?\(\?v/-?M?Dl? ףp ?X9v>?B`"{?Cl? r(?/$?~j?.$u?l?S{?E?FԈ?X9v?oa?.?> ףp=?zG?v/ ?p= ף`?E?Pn?Mb?S ?x&1?n?MbX?n?S?.gsVeg118t?+N+N/D?X9?"~Z?-G?أp= ?a"9?~j?0$u?K7AP?Gz.?@5^I?$c?n ?"~jl?ClG?A`в?Gz>?"~?ʡE?Dl?43333?:v?(\?ʡ?sh|?y&1?ˡE]?Q?M?&1|?"~:?r= ף?M"?&1?|Gb?/$?cX9?x?S㥛?/$A?-?33333?C?Hz?a?(\?rh?ʡ%?Gz~?V?n`?@5^I ?ףp= ?S?NbX9?!rh?Hz? rh?V-B?Zd;O?x&?:v?gffff?ZdK? ?V-2?@5^Il?Zd;Om?lR?v/?{G?Q?RQ?_I ??33333?E4?q= ף?Zd;?L7A`?(\?bX9H?Q? rh1?~jt?;O?Qh?/$?Q-?~j?ClW?zG?sh|?e;O?/$q?"~"?ףp= ?1Z?Gz? +q?ClC?NbX9?V}?x?nQ?7A``?X96?~jt?\(l?㥛 ?#~j?{? ףp=z?S㥛Ġ?0Z$?K71? ?&1?? ?Q?rh|?n0?(\%?:v?䥛 ?S5?Cl;?(\?A`r?Gz?Sk?"~?C?q= ף?o?VM?jt?Clg?ʡ%?Q?A`?&1?1Z?x&1?43333S?8A`?-?+^?1Z?bX9H?MbX?A`2?x&1?Q>?|?5?rh,?i|??Zd?jtD?ҭ?K7?DlI?:v?V-?$C?jt?v?l?Cl?S?M?._Veg120t?+N+Nx&?q= ף?%C?&C?+N?6^I ?ˡE-?y&1? rh?NbX9?Dl9? +?Zd;O ? rh?x/m?l?S?"?v/?l2?Q?bX9?~j?"~*?K7 ?L7A`?(\5?n?;O.?n?On? ?(\2?7A`?OnJ?Gz|?J +?L7A`?Zd;OU?E?l?x&1?Mz?v/U? r(?jt?\(L?n?m"?A`2?@5^I ?Zd?PnC?gffff?L7A`u?MbXY?Q[?On?$? ףp]?Qx?Hz?K7? ?@5^I?d;O?zGJ?"~?/d?/$e?A`?w/?Cl?K7A?K7?Ss??sh|??Mb?MB?-?Q?@5^I?NbX9?|?5?'1?Mb?'1z?NbX9?tV.?"~?-?Zd;ߏ?tV?q= ףp?(\?K7Ap?)\?(\O?̜?43333?F8? ? rh?rh|?Zd[?S?䥛 ?B`"?rh?}?L7A`e?Q[?Y?d$ȏ+N+Nymatrix128 ????O.C+=&C"]CMCC Ca CN4Cd`CCCP+C!#CdCCCB ChC#CqCCaC(C C>DC۫C COCCECClCC~&CC&C5CJ CJ!CCCBBCw%C]C yCbHC CMC7CCBBLB $CC{@CCC/ C!C5C[7CBfB Cg#CȆCq9CCD CVCФC!`CCC+CX#C,CCUC* C,u CNC!NCCCC#CwcC$CKC]C CCCCdC5C#$CCC/bC4CVqC Co C C6 C1 C %C]6 CCCCICCC C CIhCF6{diff_p2_rvaluenw z|?+N+N?>6?s?s?s?s? s?ks?s?s?s?s?s??s? s?s?s?s? s??s?s?s?s? s?s?s? s?s?s??s?s?Rt?s??s? s?s?s?s?s? s?s?s? s?s?s?s?s? s?s?s?s?s?s?s??s?s?Rt?s??s?s?s?s?s?s?s??s?s?Rt?s??s?s?s?s?s?s?s?s??s?s?s?s?s??s?s?s?s?s?s?s?s?s?s?s?s?s?s?s?s?s?s?s?s?s?s?s?s?s?s?"6?>F6diff_p1_rvaluenw z|?+N+N*`>%?\?C]?\?\?$]?$]?\?\?C]?\?C]?C]?\?\?\?C]?C]?\?C]?8]?\?\?\?C]?C]?\?C]?C]?\?\?\?l\?H]?C]?g\?\?l\?g\?g\?g\?\?C]?C]?\?h\?h\?h\?i\?\?C]?C]?\?g\?g\?\?l\?=]?a\?g\?\?\?C]?$]?C]?\?\?$]?\?\?N]?$]?C]?\?\?C]?\?C]?C]?\?C]?C]?\?\?$]?C]?\?\?\?\?C]?C]?\?C]?C]?\?\?C]?C]?\?u\?u\?C]?C]?\?C]?C]?%?r> Packages l3X0c sbcu27s l3X0cN > Asim0x?%Ɗ>N@N@333333N@LN@ffffffN@N@N@33333N@N@fffffN@O@O@333333O@LO@ffffffO@O@O@33333O@O@fffffO@P@ P@P@fffff&P@333333P@@P@LP@YP@ffffffP@33333sP@P@̌P@P@fffffP@33333P@P@P@P@fffffP@33333P@Q@ Q@Q@fffff&Q@333333Q@@Q@LQ@YQ@ffffffQ@33333sQ@Q@̌Q@Q@fffffQ@33333Q@Q@Q@Q@fffffQ@33333Q@R@ R@R@fffff&R@333333R@@R@LR@YR@ffffffR@33333sR@R@̌R@R@fffffR@33333R@R@R@R@fffffR@33333R@S@ S@S@fffff&S@333333S@@S@LS@YS@ffffffS@33333sS@S@̌S@S@fffffS@33333S@S@S@S@fffffS@33333S@T@ T@T@fffff&T@333333T@@T@LT@YT@ffffffT@33333sT@T@̌T@T@fffffT@33333T@T@T@T@fffffT@33333T@U@ U@U@fffff&U@333333U@@U@LU@YU@ffffffU@33333sU@U@̌U@U@fffffU@33333U@U@U@U@fffffU@33333U@V@ V@V@fffff&V@333333V@@V@LV@YV@ffffffV@33333sV@V@̌V@V@fffffV@33333V@V@V@V@fffffV@33333V@W@ W@W@fffff&W@333333W@@W@LW@YW@ffffffW@33333sW@W@̌W@W@fffffW@33333W@W@W@W@fffffW@33333W@X@ X@X@fffff&X@333333X@@X@LX@YX@ffffffX@33333sX@X@̌X@X@fffffX@33333X@X@X@X@fffffX@33333X@Y@ Y@Y@fffff&Y@333333Y@@Y@LY@YY@ffffffY@33333sY@Y@̌Y@Y@fffffY@33333Y@Y@Y@Y@fffffY@33333Y@Z@ Z@Z@fffff&Z@333333Z@@Z@LZ@YZ@ffffffZ@33333sZ@Z@̌Z@Z@fffffZ@33333Z@Z@Z@Z@fffffZ@33333Z@[@ [@[@fffff&[@333333[@@[@L[@Y[@ffffff[@33333s[@[@̌[@[@fffff[@33333[@[@[@[@fffff[@33333[@\@ \@\@fffff&\@333333\@@\@L\@Y\@ffffff\@33333s\@\@̌\@\@fffff\@33333\@\@\@\@fffff\@33333\@]@ ]@]@fffff&]@333333]@@]@L]@Y]@ffffff]@33333s]@]@̌]@]@fffff]@33333]@]@]@]@fffff]@33333]@^@ ^@^@fffff&^@333333^@@^@L^@Y^@ffffff^@33333s^@^@̌^@^@fffff^@33333^@^@^@^@fffff^@33333^@_@ _@_@fffff&_@333333_@@_@L_@Y_@ffffff_@33333s_@_@̌_@_@fffff_@33333_@_@_@_@fffff_@33333_@`@fffff`@ `@33333`@`@ `@fffff&`@,`@333333`@9`@@`@fffffF`@L`@33333S`@Y`@``@ffffff`@l`@33333s`@y`@`@fffff`@̌`@33333`@`@`@fffff`@̬`@33333`@`@`@fffff`@`@33333`@`@`@ [| Platform=WindowsNT, IGORVersion=3.100 NewPath/Z temppath "Q:meis:dysi:dysisimuls:dy100t48:" Simulation_summary() MoveWindow/C 5.25,280.25,584.25,468.5 dummy_graph() r_value() comparison() P2_RValue_Graph() MoveWindow/P 226.5,71,730.5,517.25 P1_RValue_Graph() Window P1_RValue_Graph() : Graph PauseUpdate; Silent 1 // building window... Display /W=(470.25,37.25,680.25,167.75) P1_RValue,fit_P1_RValue ModifyGraph mode(P1_RValue)=3 ModifyGraph marker(P1_RValue)=19 ModifyGraph lSize(fit_P1_RValue)=2 ModifyGraph rgb(fit_P1_RValue)=(0,0,0) Label left "Parameter one R-value" Label bottom "Parameter one" Textbox/N=text0/F=0/A=MC/X=3.00/Y=40.00 "\\JCmin. value@ \\{p1_min}\rsigma = \\{p1_sigma}" Tag/N=text1/F=0/X=-4.04/Y=50.00 fit_P1_RValue, 2.66055045871559636, "- sigma" Tag/N=text2/F=0/X=-0.92/Y=50.00 fit_P1_RValue, 5.77981651376146832, "+ sigma" Cursor A fit_P1_RValue 2.66055;Cursor B fit_P1_RValue 5.77982 EndMacro Window P2_RValue_Graph() : Graph PauseUpdate; Silent 1 // building window... Display /W=(374.25,281.75,584.25,412.25) P2_RValue,fit_P2_RValue ModifyGraph mode(P2_RValue)=3 ModifyGraph marker(P2_RValue)=19 ModifyGraph lSize(fit_P2_RValue)=2 ModifyGraph rgb(fit_P2_RValue)=(0,0,0) Label left "Parameter two R-value" Label bottom "Parameter two" Textbox/N=text0/F=0/A=MC/X=3.00/Y=40.00 "\\JCmin. value@ \\{p2_min}\rsigma = \\{p2_sigma}" Tag/N=text1/F=0/X=0.82/Y=50.00 fit_P2_RValue, 7.52293577981651395, "- sigma" Tag/N=text2/F=0/X=3.30/Y=50.00 fit_P2_RValue, 10, "+ sigma" Cursor A fit_P2_RValue 7.52294;Cursor B fit_P2_RValue 10.4587 EndMacro Window comparison() : Graph PauseUpdate; Silent 1 // building window... Display /W=(0,37.25,459.75,237.5) sim1 vs sim0 as "Comparison of experiment and simulation" AppendToGraph exp1 vs sim0 ModifyGraph mode(exp1)=4 ModifyGraph marker(exp1)=19 ModifyGraph lSize(sim1)=2 ModifyGraph rgb(exp1)=(0,0,0) ModifyGraph nticks(bottom)=10 Label left "scattering yield [layers]" Label bottom "scattering angle [degrees]" Legend/J/N=text0/F=0/A=MC/X=35.00/Y=-48.00 "\\JL\\Z09\\s(sim1) simulation \\s(exp1) experiment" Cursor A exp1 2;Cursor B exp1 113 Button Restart,pos={69,30},size={55,17},proc=Vegas,title="Restart" SetVariable Browse,pos={129,7},size={110,18},proc=Select_current_sim,title="Get sim." SetVariable Browse,limits={1,Inf,1},value= sim_num Button start,pos={69,51},size={55,17},proc=batch_start,title="Run" Button reset,pos={69,8},size={55,17},proc=reset,title="Reset" EndMacro Window r_value() : Graph PauseUpdate; Silent 1 // building window... Display /W=(201.75,250.25,459.75,370.25) RValue as "progress chart" ModifyGraph zero(bottom)=1 ModifyGraph nticks(bottom)=10 ModifyGraph minor(bottom)=1 Label left "R-value" Label bottom "number of files" Textbox/N=text0/F=0/A=MC/X=-6.27/Y=39.39 "min. R-value = \\{Rv}" Cursor A RValue 54 ShowInfo ControlBar 17 Button button0,pos={5,1},size={91,16},proc=save_rvalue,title="save R-value" Button button1,pos={122,1},size={90,16},proc=load_Rvalue,title="load R-value" Button button2,pos={244,1},size={52,16},proc=revert_Rvalue,title="revert" EndMacro Window dummy_graph() : Graph PauseUpdate; Silent 1 // building window... Display /W=(4.5,42.5,400.5,249.5) expt_bin1 as "dummy" EndMacro Window Simulation_summary() : Panel PauseUpdate; Silent 1 // building window... NewPanel /W=(0,43.25,641.25,484.25) as "Simulation analysis" SetDrawLayer UserBack SetDrawEnv fsize= 14,fstyle= 1 DrawText 33,31,"This program produces R-values from Vegas simulations" SetDrawEnv fstyle= 1 DrawText 30,80,"1." SetDrawEnv fstyle= 1 DrawText 30,110,"2. " SetDrawEnv fstyle= 1 DrawText 30,140,"3. " SetDrawEnv fstyle= 1 DrawText 30,170,"4. " SetDrawEnv fstyle= 1 DrawText 30,200,"5. " SetDrawEnv fstyle= 1 DrawText 30,230,"6. " SetDrawEnv fstyle= 1 DrawText 30,260,"7." SetDrawEnv fstyle= 1 DrawText 30,290,"8." SetDrawEnv fstyle= 1 DrawText 30,320,"9." SetDrawEnv fstyle= 1 DrawText 25,350,"10." SetDrawEnv fstyle= 1 DrawText 25,380,"11." SetDrawEnv fstyle= 1 DrawText 25,410,"12." SetDrawEnv fstyle= 1 DrawText 315,410,"13." SetVariable setvar0,pos={100,66},size={288,18},title="Experimental data; " SetVariable setvar0,frame=0,limits={-Inf,Inf,1},value= expdata SetVariable setvar1,pos={102,95},size={474,18},title="Simulation folder; " SetVariable setvar1,frame=0,limits={-Inf,Inf,1},value= simfolder Button button0,pos={470,307},size={91,53},proc=run,title="Run" CheckBox check0l,pos={50,123},size={185,20},proc=multi,title="Multiple simulation files",value=1 CheckBox check1,pos={50,306},size={147,20},proc=printbest,title="Print best result?",value=0 CheckBox check2,pos={50,243},size={147,20},proc=smoothset,title="Smooth simulation",value=1 CheckBox check3,pos={50,273},size={147,20},proc=spikeset,title="De-spike simulation",value=0 CheckBox check4,pos={50,212},size={147,20},proc=skewset,title="Skew data",value=1 CheckBox check4_1,pos={50,152},size={135,20},proc=element_1set,title="Use first element",value=0 CheckBox check4_2,pos={185,152},size={147,20},proc=element_2set,title="Use second element",value=1 CheckBox calibflag,pos={50,182},size={122,20},proc=calibset,title="Use calibration",value=0 SetVariable setvar4,pos={50,333},size={225,18},title="Multi parameter number 1 ;" SetVariable setvar4,limits={1,100,1},value= K11 SetVariable setvar5,pos={50,363},size={225,18},title="Multi parameter number 2 ;" SetVariable setvar5,limits={1,100,1},value= K12 SetVariable setvar6,pos={50,393},size={225,18},title="Multi parameter number 3 ;" SetVariable setvar6,limits={1,100,1},value= K13 Button find_expt,pos={45,61},size={50,20},proc=find_expt,title="find" Button find_sim,pos={46,94},size={50,20},proc=find_sim,title="find" SetVariable setcalib,pos={175,182},size={215,18},proc=SetCalibration,title="Value: counts per layer ;" SetVariable setcalib,limits={-Inf,Inf,1},value= calibration SetVariable set_calib_range,pos={399,182},size={94,18},title="% range" SetVariable set_calib_range,limits={-Inf,Inf,1},value= calibration_range SetVariable setvar7,pos={199,273},size={207,18},title="set spike tolerance level" SetVariable setvar7,limits={0,2,0.05},value= spikelevel SetVariable setvar8,pos={199,244},size={127,18},title="set amount " SetVariable setvar8,limits={3,Inf,2},value= smooth_amount SetVariable setvar2,pos={337,393},size={238,18},title="Parameter polynomial fit order" SetVariable setvar2,limits={3,9,2},value= poly_amount EndMacro т |version 2 - 6/3/96 - expt and sim matched to average between 60 and 135 degrees |version 3 - 12/3/96 - low angle intensity fall-off corrected for |version 4 - 15/3/96 - new exptl. data |version 5 - 18/3/96 - overall curvature messed about with |version 6 - 19/3/96 - smoothing added to simulation |version 7 - 19/3/96 - IGOR pro |version 8 - 21/3/96 - batch processing and smooth curve correction for each file |version 9 - 8/9/96 - working towards all singing, all dancing versuion |version 10 - 12/1/97 - ??? |version 11 - 15/1/97 - Contours added |version 12 - 21/1/97 - R-value error fixed |version 59 - 28/1/98 - files loaded by name (for Dave Spence, York) |version 60 - 11/2/98 - oops! test Xsq expression was being used (only minor effect) |version 61 - 19/3/98 - Rvalue parameter fit now uses 5th order polygon |version 62 - 20/3/98 - Rvalue scaling error fixed (sim_avg) |version 63 - 20/3/98 - first attempt at including absolute calibration |version 68 - 23/3/98 - refined |version 70 - 25/3/98 - 2D output changed |version 71 - 1/6/98 - matching section removed |version 72 - 3/7/98 - calibration bug in batch_end fixed |version 73 - 3/8/98 - modified for the new element specific output from vegas |version 74 - 5/8/98 - check included for matching file length between data and simulation |version 75 - 27/8/98 - bug fix in 2D layout macro ("error cursor not on top graph") - not true!!!!! |version 76 - 10/9/98 - order of parameters changed to match the order in Vegas |version 77 - 24/9/98 - moved from mac to pc platform |version 78 - 23/10/98 - fixed (hopefully) "error cursor not on top graph" |version 79 - 9/12/98 - contour axes corrected + few minor mods |version 80 - 23/3/99 - "error cursor not on top graph" fixed for 3D simulations (still don't understand it) |version 81 - 23/3/99 - de-spiking option included |version 82 - 23/3/99 - smoothing of simulation can be user-defined |version 83 - 14/4/99 - 3D summary prinout changed |version 84 - 15/4/99 - 3D contours tweaked |version 85 - 21/4/99 - improved de-spiking |version 86 - 22/4/99 - few minor changes |version 87 - 30/4/99 - changed contour levels to auto |version 88 - 19/5/99 - moved skew and calib into subroutines |version 89 - 20/5/99 - STILL found things to change! |version 90 - 2/6/99 - Element specific options introduced. |version 91 - 4/6/99 - Bugs introduced in versions 88 onwards fixed! |version 92 - 7/6/99 - Minor bugs fixed. |version 93 - 8/6/99 - Control buttons moved to LHS of comparison graph. |version 94 - 10/6/99 - Very subtle bug destroyed! |version 95 - 22/7/99 - Modified to accept just one multimodel parameter. #pragma rtglobals=0 //sod this silly new global parameter stuff!! macro set_variables() K11=5; K12=5; K13=5 string/G S_filename string/G tempath string/G rawdata="no_filename" string/G expdata="no_filename" string/G simdata="no_filename" string/G simfolder="no_foldername" string/G str1 variable/G E_angle1=45 //start angle for expt. data variable/G E_angle2=E_angle1+24 //end angle for expt. data variable/G S_angle1, S_angle2 //start and end angles for simul. variable/G delta_angle=0.5 //angle step for simul. variable/G Rv, index, sim_num variable/G printflag=0 variable/G multiflag=1 variable/G element_1flag, element_2flag=1 variable/G smoothflag=1, smooth_amount=3 variable/G spikeflag=0, spikedetect=0 variable/G spikelevel=0.5 variable/G calibflag=0 variable/G skewflag=1 variable/G Ne1, Ne2, N_over=0 variable/G calib_angle, offset, expand //original calibration angle, calculated offset, calculated expansion variable/G exp_avg, sim_avg, Ang_avg=1 //average y-value of expt. data and sim data variable/G P1_sigma, P2_sigma, P3_sigma variable/G p1_min, P2_min, P3_min variable/G avg_y, scale, calibration, calibration_range //scale is the emperical calibration variable/G poly_amount=3 Make/D/O/N=100 sim0, sim1, sim2, sim3, sim4, sim5, sim_smooth, sim_fit Make/D/O/N=100 exp1, expt0, expt1, expt_bin0, expt_bin1, exp_orig, exp_fit Make/D/O/N=100 RValue, Rvalue0, Rvalue2, Rvalue_orig, W_coef KillWaves/A/Z/F DoWindow/B Contour_12 DoWindow/B Contour_23 DoWindow/B Contour_31 DoWindow/B Comparison DoWindow/B R_Value DoWindow/B Exp_Show DoWindow/B Simulation_summary DoWindow/F Match_Data tidy_now() end Menu "Tidy" "-" "tidy_now/0" end Macro batch_set_up() silent 1; pauseupdate variable M string pathname, filename if (k12==0) k12=1 endif if (k13==0) k13=1 endif pathname="temppath" filename=indexedfile($pathname, 0, "????") Loadwave/Q/G/D/O/W/P=temppath/N=sim/L={0,0,0,0,6} filename duplicate/O sim1 exp1, exp_fit, sim_fit, sim_smooth dowindow/k comparison2 dowindow/k comparison dowindow/k comparison2_1 dowindow/k comparison_1 execute "comparison()" if (numpnts(expt_bin1)!=numpnts(sim1)) tidy_now() beep abort "the data and simulation files are not the same length" endif if ((element_1flag==0) %& (element_2flag==0)) beep abort "you have to select a least one element!" endif if ((element_1flag==1) %& (element_2flag==1)) sim1=sim1+sim2 //using both elements endif if ((element_1flag==0) %& (element_2flag==1)) sim1=sim2 //using second element only endif exp1=expt_bin1 duplicate/O exp1, exp_orig Ne1=2 Ne2=numpnts(exp1)-3 if (spikeflag>0) spike() endif if (smoothflag>0) Smooth smooth_amount, sim1 endif WaveStats/Q/R=[Ne1, Ne2] exp1 exp_avg=V_avg WaveStats/Q/R=[Ne1, Ne2] sim1 Sim_Avg=V_avg exp1*=Sim_avg/exp_avg WaveStats/Q/R=[Ne1, Ne2] sim0 Ang_Avg=V_avg Legend/K/N=text0 Legend/J/N=text0/F=0/A=MC/X=18.13/Y=-43.81 "\\JL\\Z09\\s(sim1) simulation \\s(exp1) experiment" if (calibflag>0) AppendText "scaling factor = \\{scale/calibration}, calibration range = \\{calibration_range}%" endif cursor/C=(0,0,65000)/P A, exp1, Ne1 cursor/C=(0,65000,0)/P B, exp1, Ne2 N_over=ne2-ne1+1 if (skewflag>0) skew() endif if (calibflag>0) calib() endif DoWindow/B simulation_summary DoWindow/B match_data DoWindow/F R_Value DoWindow/F Comparison KillWaves/A/Z/F end macro batch_start(ctrlName) : ButtonControl String ctrlName, filename, pathname, str1 pauseupdate silent 1 spikedetect=0 Ne1=pcsr(A) Ne2=pcsr(B) N_over=Ne2-Ne1+1 WaveStats/Q/R=[Ne1, Ne2] sim0 Ang_Avg=V_avg DoWindow/K R_Value DoWindow/K Comparison R_Value2() comparison2() DoWindow/B contour12 DoWindow/B contour23 DoWindow/B contour13 print "Started on "+date()+" at "+time() batchprocess() end Function Calc_RValue_is() variable N Rv=0 do Rv+=((sim1[N+Ne1]-exp1[N+Ne1])/exp1[N+Ne1])^2 N+=1 while (N_over>(N-1)) RValue[index]=sqrt(Rv)*100/N_over end Function Calc_RValue_Xsq() variable N, R N=Ne1-1 do N+=1 R+=(sim1[N]-exp1[N])^2/exp1[N]/(sin(sim0[N]/57.3/2)^4) while (N0) spike() endif if (smoothflag>0) Smooth smooth_amount, sim1 endif sprintf sim_name, "Veg" +"%g" index duplicate/O sim1 $sim_name //creates wave from modified simulation WaveStats/Q/R=[Ne1, Ne2] sim1 sim_avg=v_avg WaveStats/Q/R=[Ne1, Ne2] exp_orig exp_avg=v_avg scale=exp_avg/sim_avg //scale is the emperical calibration exp1=exp_orig/scale //scales experimental data - simulation is left unchanged if (skewflag>0) skew() endif if (calibflag>0) calib() endif Calc_RValue_Xsq() index+=1 doupdate while (index<(filenum+1)) simdata=S_Filename Execute "batch_end()" End Macro batch_end() print "Finished on "+date()+" at "+time() if (spikeflag>0) print "% of spikes detected per file= ", (spikedetect/index/(Ne2-Ne1)*100) endif tidy_now() silent 1 variable N pauseupdate wavestats/Q RValue sim_num=V_minloc //sim_num points to simulation with lowest R Value String sim_name sprintf sim_name, "Veg"+"%g" sim_num sim1=$sim_name WaveStats/Q/R=[Ne1, Ne2] sim1 sim_avg=v_avg WaveStats/Q/R=[NE1,Ne2] exp_orig exp_avg=v_avg scale=exp_avg/sim_avg //scale is the emperical calibration exp1=exp_orig/scale //scales experimental data - simulation is left unchanged DoWindow/F Comparison if (skewflag>0) skew() endif if (calibflag>0) calib() endif DoWindow/F R_Value if ((K13>=2) %& (K12>=2)) contour_3D() RValue_Graphs_3D() endif if ((K13==1) %& (K12>=2)) contour_2D() RValue_Graphs_2D() endif if ((K13==1) %& (K12==1)) RValue_Graphs_1D() endif Rv=RValue[sim_num-1] if (printflag==1) printout() endif tidy_now() end Macro single() silent 1 pauseupdate string pathname, filename pathname="temppath" Loadwave/Q/G/D/O/W/P=temppath/N=sim/L={0,0,0,0,6} filename duplicate/O sim1 exp1, exp_fit, sim_fit exp1=expt_bin1 duplicate/O exp1, exp_orig Ne1=2 Ne2=numpnts(exp1)-3 WaveStats/Q/R=[Ne1, Ne2] exp1 exp_avg=V_avg WaveStats/Q/R=[Ne1, Ne2] sim1 Sim_Avg=V_avg exp1*=Sim_avg/exp_avg WaveStats/Q/R=[Ne1, Ne2] sim0 Ang_Avg=V_avg cursor/C=(0,0,65000)/P A, exp1, Ne1 cursor/C=(0,65000,0)/P B, exp1, Ne2 N_over=ne2-ne1+1 DoWindow/B simulation_summary DoWindow/B match_data DoWindow/F R_Value DoWindow/F Comparison KillWaves/A/Z/F if ((element_1flag==0) %& (element_2flag==0)) beep abort "you have to select a least one element!" endif if ((element_1flag==1) %& (element_2flag==1)) sim1=sim1+sim2 //using both elements endif if ((element_1flag==0) %& (element_2flag==1)) sim1=sim2 //using second element only endif if (smoothflag>0) Smooth smooth_amount, sim1 endif WaveStats/Q/R=[Ne1, Ne2] sim1 exp1=exp1*V_Avg/exp_avg //scales experimental data - simulation is left unchanged if (skewflag>0) skew() endif if (calibflag>0) calib() endif Calc_RValue_Xsq() doupdate PRINT "the rvalue is", Rv end Macro Printout() silent 1 Summary_Layout_3D() printlayout Summary_Layout_3D DoWindow /K Summary_Layout_3D Return Window r_value() : Graph PauseUpdate; Silent 1 Display /W=(202,250,460,370) RValue as "progress chart" ModifyGraph zero(bottom)=1 ModifyGraph nticks(bottom)=10 ModifyGraph minor(bottom)=1 Label left "R-value" Label bottom "number of files" Wavestats/Q RValue Textbox/N=text0/F=0/A=MC/X=-6.27/Y=39.39 "min. R-value = \\{Rv}" ShowInfo ControlBar 17 Button button0,pos={5,1},size={91,16},proc=save_rvalue,title="save R-value" Button button1,pos={122,1},size={90,16},proc=load_Rvalue,title="load R-value" Button button2,pos={244,1},size={52,16},proc=revert_Rvalue,title="revert" EndMacro Window comparison() : Graph PauseUpdate; Silent 1 Display /W=(0,30,460,230) sim1 vs sim0 as "Comparison of experiment and simulation" Append exp1 vs sim0 ModifyGraph mode(exp1)=4 ModifyGraph marker(exp1)=19 ModifyGraph lSize(sim1)=2 ModifyGraph rgb(exp1)=(0,0,0) ModifyGraph nticks(bottom)=10 Label left "scattering yield [layers]" Label bottom "scattering angle [degrees]" Legend/J/N=text0/F=0/A=MC/X=35/Y=-48 "\\JL\\Z09\\s(sim1) simulation \\s(exp1) experiment" if (calibflag>0) AppendText "scaling factor = \\{scale/calibration}, calibration range = \\{calibration_range}%" endif cursor/C=(0,0,65000)/P A, exp1, Ne1 cursor/C=(0,65000,0)/P B, exp1, Ne2 Button Restart,pos={69,30},size={55,17},proc=Vegas,title="Restart" SetVariable Browse,pos={129,7},size={110,18},proc=Select_current_sim,title="Get sim." SetVariable Browse,limits={1,Inf,1},value= sim_num Button start,pos={69,51},size={55,17},proc=batch_start,title="Run" Button reset,pos={69,8},size={55,17},proc=reset,title="Reset" EndMacro Function reset(ctrlName) : ButtonControl String ctrlName execute "batch_end()" End Window Simulation_summary() : Panel PauseUpdate; Silent 1 // building window... NewPanel /W=(0,42.8,640.2,483.2) as "Simulation analysis" SetDrawLayer UserBack SetDrawEnv fsize= 14,fstyle= 1 DrawText 33,31,"This program produces R-values from Vegas simulations" SetDrawEnv fstyle= 1 DrawText 30,80,"1." SetDrawEnv fstyle= 1 DrawText 30,110,"2. " SetDrawEnv fstyle= 1 DrawText 30,140,"3. " SetDrawEnv fstyle= 1 DrawText 30,170,"4. " SetDrawEnv fstyle= 1 DrawText 30,200,"5. " SetDrawEnv fstyle= 1 DrawText 30,230,"6. " SetDrawEnv fstyle= 1 DrawText 30,260,"7." SetDrawEnv fstyle= 1 DrawText 30,290,"8." SetDrawEnv fstyle= 1 DrawText 30,320,"9." SetDrawEnv fstyle= 1 DrawText 25,350,"10." SetDrawEnv fstyle= 1 DrawText 25,380,"11." SetDrawEnv fstyle= 1 DrawText 25,410,"12." SetDrawEnv fstyle= 1 DrawText 315,410,"13." SetVariable setvar0,pos={100,66},size={288,18},title="Experimental data; " SetVariable setvar0,frame=0,limits={-Inf,Inf,1},value= expdata SetVariable setvar1,pos={102,95},size={474,18},title="Simulation folder; " SetVariable setvar1,frame=0,limits={-Inf,Inf,1},value= simfolder Button button0,pos={470,307},size={91,53},proc=run,title="Run" CheckBox check0l,pos={50,123},size={185,20},proc=multi,title="Multiple simulation files",value=1 CheckBox check1,pos={50,306},size={147,20},proc=printbest,title="Print best result?",value=0 CheckBox check2,pos={50,243},size={147,20},proc=smoothset,title="Smooth simulation",value=1 CheckBox check3,pos={50,273},size={147,20},proc=spikeset,title="De-spike simulation",value=1 CheckBox check4,pos={50,212},size={147,20},proc=skewset,title="Skew data",value=1 CheckBox check4_1,pos={50,152},size={135,20},proc=element_1set,title="Use first element",value=0 CheckBox check4_2,pos={185,152},size={147,20},proc=element_2set,title="Use second element",value=1 CheckBox calibflag,pos={50,182},size={122,20},proc=calibset,title="Use calibration",value=0 SetVariable setvar4,pos={50,333},size={225,18},title="Multi parameter number 1 ;" SetVariable setvar4,limits={1,100,1},value= K11 SetVariable setvar5,pos={50,363},size={225,18},title="Multi parameter number 2 ;" SetVariable setvar5,limits={1,100,1},value= K12 SetVariable setvar6,pos={50,393},size={225,18},title="Multi parameter number 3 ;" SetVariable setvar6,limits={1,100,1},value= K13 Button find_expt,pos={45,61},size={50,20},proc=find_expt,title="find" Button find_sim,pos={46,94},size={50,20},proc=find_sim,title="find" SetVariable setcalib,pos={175,182},size={215,18},proc=SetCalibration,title="Value: counts per layer ;" SetVariable setcalib,limits={-Inf,Inf,1},value= calibration SetVariable set_calib_range,pos={399,182},size={94,18},title="% range" SetVariable set_calib_range,limits={-Inf,Inf,1},value= calibration_range SetVariable setvar7,pos={199,273},size={207,18},title="set spike tolerance level" SetVariable setvar7,limits={0,2,0.05},value= spikelevel SetVariable setvar8,pos={199,244},size={127,18},title="set amount " SetVariable setvar8,limits={3,Inf,2},value= smooth_amount SetVariable setvar2,pos={337,393},size={238,18},title="Parameter polynomial fit order" SetVariable setvar2,limits={3,9,2},value= poly_amount EndMacro Macro element_2set(ctrlName, flag) : checkboxControl String ctrlName variable flag element_2flag=flag EndMacro Macro element_1set(ctrlName, flag) : checkboxControl String ctrlName variable flag element_1flag=flag EndMacro macro find_expt(ctrlName) : ButtonControl String ctrlName loadwave /G/O/A/Q/N=expt_bin expdata=S_filename exp_show() End macro find_sim(ctrlName) : ButtonControl String ctrlName newpath/O temppath pathinfo temppath simfolder=S_path End Macro multi(ctrlName, flag) : checkboxControl String ctrlName variable flag multiflag=flag EndMacro Macro printbest(ctrlName, flag) : checkboxControl String ctrlName variable flag printflag=flag EndMacro Macro smoothset(ctrlName, flag) : checkboxControl String ctrlName variable flag smoothflag=flag EndMacro Macro spikeset(ctrlName, flag) : checkboxControl String ctrlName variable flag spikeflag=flag EndMacro Macro calibset(ctrlName, flag) : checkboxControl String ctrlName variable flag calibflag=flag EndMacro Macro skewset(ctrlName, flag) : checkboxControl String ctrlName variable flag skewflag=flag EndMacro Macro run(ctrlName) : ButtonControl String ctrlName DoWindow/B Simulation_summary DoWindow/B Match_data if (multiflag==1) batch_set_up() else single() endif EndMacro Function ButtonRestart(ctrlName) : ButtonControl String ctrlName dowindow/f Match_data End //________________________________________________________________ // START OF CONTOUR STUFF //________________________________________________________________ //K11 is number of values of multi-model paramter 1 //K12 is number of values of multi-model paramter 2 //K13 is number of values of multi-model paramter 3 macro test() K16=trunc(k19/K11/k12) K15=trunc((k19-K16*K11*K12)/K11) K14=k19-K15*K11-K16*K12*K11 print k19 print k14 print k15 print k16 end function contour_2D() variable N, M wavestats/q rvalue V_minloc-=1 K14=trunc(V_minloc/K12) K15=V_minloc-K14*K12 Make/N=(K12, K11)/O matrix12 N=0 M=0 do | picks K11 and K12 do matrix12[N] [M]=rvalue[(N+M*K12)] N+=1 while (N0) skew() endif if (calibflag>0) calib() endif End Macro Vegas(ctrlName) : ButtonControl String ctrlName DoWindow/K exp_show DoWindow/F Simulation_summary End Window comparison2() : Graph PauseUpdate; Silent 1 Display /W=(0,30,680,320) sim1 vs sim0 as "Comparison of experiment and simulation" Append exp1 vs sim0 ModifyGraph mode(exp1)=4 ModifyGraph marker(exp1)=19 ModifyGraph lSize(sim1)=2 ModifyGraph rgb(exp1)=(0,0,0) ModifyGraph nticks(bottom)=10 ModifyGraph lblMargin(bottom)=3 ModifyGraph lblLatPos(bottom)=-1 Textbox/N=text0/S=3/A=MC/X=32.14/Y=-44.25 "simulation number = \\{index}" cursor/C=(0,0,65000)/P A, exp1, Ne1 cursor/C=(0,65000,0)/P B, exp1, Ne2 EndMacro Window R_Value2() : Graph PauseUpdate; Silent 1 Display /W=(170,330,510,480) RValue as "progress chart" ModifyGraph zero(bottom)=1 ModifyGraph nticks(bottom)=10 ModifyGraph minor(bottom)=1 SetAxis/A/E=1 bottom EndMacro function tidy_now() PauseUpdate; Silent 1 dowindow/k comparison2 dowindow/k r_value2 dowindow/k comparison dowindow/k r_value dowindow/k comparison2_1 dowindow/k r_value2_1 dowindow/k comparison_1 dowindow/k r_value_1 dowindow/k Graph0 dowindow/k Graph1 dowindow/k table0 dowindow/k table1 wavestats/Q RValue execute "r_value()" cursor/p A, RValue, (sim_num-1) execute "comparison()" end function chi_square_precision_3D() variable N, sigma execute "CurveFit/Q/L=(numpnts(p3_rvalue)*10) poly poly_amount, P3_RValue /D " wavestats/Q fit_p3_rvalue p3_min=V_minloc Duplicate/O fit_p3_rvalue diff_p3_rvalue Differentiate diff_p3_rvalue Differentiate diff_p3_rvalue sigma=(abs(2/ diff_p3_rvalue[9+10*trunc(V_minloc)]))^0.5 p3_sigma=sigma dowindow/F P3_Rvalue_graph cursor A, fit_p3_rvalue, V_minloc-sigma cursor B, fit_p3_rvalue, V_minloc+sigma Textbox/N=text0/F=0/A=MC/X=10/Y=40 "\\JCmin. value@ \\{p3_min}\rsigma = \\{p3_sigma}" Tag/N=text1/F=0/X=(xcsr(A)-6.7)/Y=50 fit_P3_RValue, xcsr(A), "- sigma" Tag/N=text2/F=0/X=(xcsr(B)-6.7)/Y=50 fit_P3_RValue, xcsr(B), "+ sigma" execute "CurveFit/Q/L=(numpnts(p1_rvalue)*10) poly poly_amount, P1_RValue /D " wavestats/Q fit_p1_rvalue p1_min=V_minloc Duplicate/O fit_p1_rvalue diff_p1_rvalue Differentiate diff_p1_rvalue Differentiate diff_p1_rvalue sigma=(abs(2/ diff_p1_rvalue[9+10*trunc(V_minloc)]))^0.5 p1_sigma=sigma dowindow/F P1_Rvalue_graph cursor A, fit_p1_rvalue, V_minloc-sigma cursor B, fit_p1_rvalue, V_minloc+sigma Textbox/N=text0/F=0/A=MC/X=3/Y=40 "\\JCmin. value@ \\{p1_min}\rsigma = \\{p1_sigma}" Tag/N=text1/F=0/X=(xcsr(A)-6.7)/Y=50 fit_P1_RValue, xcsr(A), "- sigma" Tag/N=text2/F=0/X=(xcsr(B)-6.7)/Y=50 fit_P1_RValue, xcsr(B), "+ sigma" execute "CurveFit/Q/L=(numpnts(p2_rvalue)*10) poly poly_amount, P2_RValue /D " wavestats/Q fit_p2_rvalue p2_min=V_minloc Duplicate/O fit_p2_rvalue diff_p2_rvalue Differentiate diff_p2_rvalue Differentiate diff_p2_rvalue sigma=(abs(2/ diff_p2_rvalue[9+10*trunc(V_minloc)]))^0.5 p2_sigma=sigma dowindow/F P2_Rvalue_graph cursor A, fit_p2_rvalue, V_minloc-sigma cursor B, fit_p2_rvalue, V_minloc+sigma Textbox/N=text0/F=0/A=MC/X=10/Y=40 "\\JCmin. value@ \\{p2_min}\rsigma = \\{p2_sigma}" Tag/N=text1/F=0/X=(xcsr(A)-6.7)/Y=50 fit_P2_RValue, xcsr(A), "- sigma" Tag/N=text2/F=0/X=(xcsr(B)-6.7)/Y=50 fit_P2_RValue, xcsr(B), "+ sigma" end function chi_square_precision_2D() variable N, sigma execute "CurveFit/Q/L=(numpnts(p2_rvalue)*10) poly poly_amount, P2_RValue /D " wavestats/Q fit_p2_rvalue p2_min=V_minloc Duplicate/O fit_p2_rvalue diff_p2_rvalue Differentiate diff_p2_rvalue Differentiate diff_p2_rvalue sigma=(abs(2/ diff_p2_rvalue[9+10*trunc(V_minloc)]))^0.5 p2_sigma=sigma dowindow/F P2_Rvalue_graph cursor A, fit_p2_rvalue, V_minloc-sigma cursor B, fit_p2_rvalue, V_minloc+sigma Textbox/N=text0/F=0/A=MC/X=3/Y=40 "\\JCmin. value@ \\{p2_min}\rsigma = \\{p2_sigma}" Tag/N=text1/F=0/X=(xcsr(A)-6.7)/Y=50 fit_P2_RValue, xcsr(A), "- sigma" Tag/N=text2/F=0/X=(xcsr(B)-6.7)/Y=50 fit_P2_RValue, xcsr(B), "+ sigma" execute "CurveFit/Q/L=(numpnts(p1_rvalue)*10) poly poly_amount, P1_RValue /D " wavestats/Q fit_p1_rvalue p1_min=V_minloc Duplicate/O fit_p1_rvalue diff_p1_rvalue Differentiate diff_p1_rvalue Differentiate diff_p1_rvalue sigma=(abs(2/ diff_p1_rvalue[9+10*trunc(V_minloc)]))^0.5 p1_sigma=sigma dowindow/F P1_Rvalue_graph cursor A, fit_p1_rvalue, V_minloc-sigma cursor B, fit_p1_rvalue, V_minloc+sigma Textbox/N=text0/F=0/A=MC/X=3/Y=40 "\\JCmin. value@ \\{p1_min}\rsigma = \\{p1_sigma}" Tag/N=text1/F=0/X=(xcsr(A)-6.7)/Y=50 fit_P1_RValue, xcsr(A), "- sigma" Tag/N=text2/F=0/X=(xcsr(B)-6.7)/Y=50 fit_P1_RValue, xcsr(B), "+ sigma" end function chi_square_precision_1D() variable sigma p1_rvalue=rvalue execute "CurveFit/Q/L=(numpnts(p1_rvalue)*10) poly poly_amount, P1_RValue /D " wavestats/Q fit_p1_rvalue p1_min=V_minloc Duplicate/O fit_p1_rvalue diff_p1_rvalue Differentiate diff_p1_rvalue Differentiate diff_p1_rvalue sigma=(abs(2/ diff_p1_rvalue[9+10*trunc(V_minloc)]))^0.5 p1_sigma=sigma dowindow/F P1_Rvalue_graph cursor A, fit_p1_rvalue, V_minloc-sigma cursor B, fit_p1_rvalue, V_minloc+sigma Textbox/N=text0/F=0/A=MC/X=3/Y=40 "\\JCmin. value@ \\{p1_min}\rsigma = \\{p1_sigma}" Tag/N=text1/F=0/X=(xcsr(A)-6.7)/Y=50 fit_P1_RValue, xcsr(A), "- sigma" Tag/N=text2/F=0/X=(xcsr(B)-6.7)/Y=50 fit_P1_RValue, xcsr(B), "+ sigma" end Function save_rvalue(ctrlName) : ButtonControl String ctrlName Save/G RValue as "RValue" End Function load_Rvalue(ctrlName) : ButtonControl String ctrlName ctrlName = S_Filename Duplicate/O RValue Rvalue_orig loadwave /G/O/A/Q/N=Rvalue Rvalue+=Rvalue0 Rvalue/=2 S_Filename = ctrlName execute "batch_end()" End Function revert_Rvalue(ctrlName) : ButtonControl String ctrlName Rvalue=Rvalue_orig execute "batch_end()" End Function calc_rvalue(ctrlName) : ButtonControl String ctrlName Ne1=pcsr(A) Ne2=pcsr(B) N_over=Ne2-Ne1+1 Calc_RValue_Xsq() execute "print rvalue(index)" End Function SetCalibration(ctrlName,varNum,varStr,varName) : SetVariableControl String ctrlName Variable varNum String varStr String varName End Window dummy_graph() : Graph PauseUpdate; Silent 1 Display /W=(5,42,400,250) expt_bin1 as "dummy" EndMacro Window Summary_layout_3D() : Layout PauseUpdate; Silent 1 // building window... Layout/C=1/W=(4.8,42.2,573.6,445.4) Contour_12(601.8,226.2,768.6,369.6)/O=1,comparison(72.6,72.6,414.6,216)/O=1 Append r_value(72.6,226.2,415.8,369.6)/O=1,P1_RValue_Graph(425.4,72.6,592.2,216)/O=1 Append P2_RValue_Graph(425.4,226.2,592.2,369.6)/O=1,P3_RValue_Graph(425.4,378.6,592.2,522)/O=1 Append Contour_31(601.8,72.6,768.6,216)/O=1,Contour_23(601.8,378.6,768.6,522)/O=1 Textbox/N=text0/F=0/A=LB/X=0.60/Y=0.8 "\\Z12Printed on "+ date() + "@" + time() AppendText "simulation file; "+simdata+".\{sim_num}" AppendText "exptl. data file; "+expdata sprintf str1 "The R-value is %g", RValue[k8-1] AppendText str1 sprintf str1 "For %g data points", (ne2-ne1) //sprintf str1 "The R-value-2 is %g", ((RValue[k8-1]-1)/avg_y*1e3) AppendText str1 sprintf str1 "from \{hcsr(A)} to \{hcsr(B)} degrees" AppendText str1 if (calibflag>0) AppendText "Using absolute calibration" else AppendText "Not using absolute calibration" endif if (spikeflag>0) sprintf str1 "percentage of spikes in file=\{spikedetect/index/(Ne2-Ne1)*100}" AppendText str1 sprintf str1 "spike tolerance range=\{spikelevel}" AppendText str1 else AppendText "Simulation not de-spiked" endif if (smoothflag>0) sprintf str1 "Simulation smoothed, value = \{smooth_amount}" AppendText str1 else AppendText "Simulation not smoothed" endif if (skewflag>0) sprintf str1 "Simulation skewed" AppendText str1 else AppendText "Simulation not skewed" endif EndMacro Function spike() variable M=0 do sim_smooth=sim1 smooth 15, sim_smooth if (sim1[M+1]>(sim_smooth[M]+spikelevel)) // checks for a greater than 'spikelevel' increase.... sim1[M+1]=sim_smooth[M+1] // between adjacent points spikedetect+=1 endif M+=1 while ((M+1)(1+calibration_range/100))) |sets the upper limit to the range of calibration allowed exp1=exp1*scale/(calibration*(1+calibration_range/100)) endif if (scale/calibration<(1-calibration_range/100))) |sets the lower limit to the range of calibration allowed exp1=exp1*scale/(calibration*(1-calibration_range/100)) endif end Function skew() CurveFit/Q poly 3, exp1 (xcsr(A),xcsr(B)) /X=sim0 /D=exp_fit exp_fit= W_coef[0]+ W_coef[1]*sim0+ W_coef[2]*(sim0^2) CurveFit/Q poly 3, sim1 (xcsr(A),xcsr(B)) /X=sim0 /D=sim_fit sim_fit= W_coef[0]+ W_coef[1]*sim0+ W_coef[2]*(sim0^2) exp1*=sim_fit/exp_fit end Window exp_show() : Graph PauseUpdate; Silent 1 Display /W=(0,40,640,480) expt_bin1 vs expt_bin0 as "Experimental data" ModifyGraph mode=4 ModifyGraph marker=19 ModifyGraph lStyle=1 ModifyGraph nticks(bottom)=10 Label left "scattering yield corrected for sin\\S4\\M\\F'Symbol'(q/2)" Label bottom "scattering angle [degrees]" Textbox/N=text0/F=0/S=2/A=MC/X=-7.93/Y=47.96 "Data from file - " + expdata Button exp_OK,pos={523,209},size={70,20},proc=Vegas,title="Continue" EndMacro