From 67d4743f38511ab300a5a9c260d9cc13d91d9b40 Mon Sep 17 00:00:00 2001 From: yan Date: Fri, 12 Apr 2024 09:33:49 +0800 Subject: [PATCH] update ifs_sim --- .../__pycache__/csst_ifs_sim.cpython-38.pyc | Bin 0 -> 88440 bytes csst_ifs_sim/csst_ifs_sim.py | 24 ++++++++++++++---- 2 files changed, 19 insertions(+), 5 deletions(-) create mode 100644 csst_ifs_sim/__pycache__/csst_ifs_sim.cpython-38.pyc diff --git a/csst_ifs_sim/__pycache__/csst_ifs_sim.cpython-38.pyc b/csst_ifs_sim/__pycache__/csst_ifs_sim.cpython-38.pyc new file mode 100644 index 0000000000000000000000000000000000000000..3413f5b1d8cc55df760938f0976349a17a09291d GIT binary patch literal 88440 zcmd?S34B~vc{V=#rqz-x%S-I>8q13;%g*LFk@wiamL1D>5|J>BbVj3*M>F!Bk!=kl z2qZuXgt}1n0CpQ_O@Ts63xvf13QZwU%2s+y3x$ADfwn+dYL@(;=RJ4s%pJ*+2&Mn; z`+d@H&OP@%%RT#fxAUIkb#)OJ{?5OwZ7{LM>-rX7D*r-w*@0g!?Qyx3tIws}%9D2Y zx$*R-JtN*euUo@>BmTaCk*tf{?FYa4{bTvqKP2V*F zt4;exmi8?jS=P5q-s{rQk>!0W#9yCYIkKv6mG~PFZ*|{l@i(T|jI8NfGjeU;wes#s zukC9`zVp)SM%MSOw_qFkHUQR?-Z;|H*DVxNoETh&%|9iH9ldUXSyJJgNp zCOmhlZR%z`cd6~_7Cd*WTU9rnd(;lK6VJVBm)eczKD9^f#q&0`Pu+&+ezjj6!1I8* zT|EoW+topJ2+wD!9(5ScgNZ|fZgoW6@qoLpN5$1q)eGNYbxa*c=n-{7-HGQN>MnIR zo<~)mdN!WDsMS5{UexN|`#pWfRKJR;=OW#4<>_^`4-~dUb{mO!KB1!7O!U~$Saj#8 z5#78wy7{_on{V89-3`&M&aSPI$d33}ekf~fi;l-Lkx35To_4nxN@em1V|25Ezv;lf z-c%-+H^xR1nf#%wN~8ge3}ka7sR1KCp2J%polPba27KP##|~U{o574qjl`1)KxzgD z^8?w;>1_HmUwwnAJd@PuPh`_5W4p%;^`0Exoy{7`tlJwm@spQ-xH5okQE=;q2qk;BN~!>BaX<&)LnUXB@x6=8G<46B_Nx zmQ6$1k;Ep}cpkjv#`gkYpOUOGm{&Dy*X&*YBA_ap|o_DryMxf8hN zEpab#*Sn3)@H+-V4FDV`G7vKG=+2z*6g|p=&i5+sgirY<{6%#2For|Xr@YOs7MJob zb`|~fFLtTGV%J0nZYa&vN zV1Psle-_cLHFHm%iJwlOqpsfK?AyWg*x4B)&DfiJXHRLwv^s#Wbuim<0JiX)4f-Qk9Iowy$~aJ zrPArw;4Rk*$F(Tr33+{PpL>Nj;;wavJU(~5$G8n4j?pm>0InS}I^F{hr~U$25{BIz zAlnos6aWqPtj+ zbPdIZq90#sMC!&t?naG;ug+$I;woX7!GY`uH&=suhYV<=i@q&%w$kCyHm;*{J)Ilq+z6*Vgb}UG5Vm?18_1nD z1E=HZu|!S=yKxhvakLxT=-kW%^GqL`+1Q5{49ZkGk-^$+hOpgAVNJiZ61Gq@L#J^o z6L!B;pxkF9SV$JI%mgd?24S%$kb!oH8jeAG( zICBRL2$XxJ7t5=s*~O{T>{@_j&=D_~_Z0(q|8M|PBdDL8U=MjGLdDSepz>qc3?q)u zisN_0!GdbVsZoK6+Noc|#I(MO>Ek*d;`Gt`yCj~Fm`16nZD$Jk-O+M))BHE4mJSdIA&`yksaVvO$`7eP9z4>snL9D zAoZ;_%&mytV}fyGARo_k$^0{Yxl=~osAWcrm{9;|kC8}NNWx&&7>VVx+>Qo_lti#D zU}qB4cb*xv`p)!A=a_zJebX<^YX)@t5w(>D5&-&81Lz4GSYNQ1!5c{T%hLj(-37#F zZGsqg3~^IqIg|h^!B8TV970CYaH(+{MHZCO2!Yjh40_H$G+w+fDHQLB3`rgF&T zc{<>UV2>;LeD0<0HgCj>-R>;n4KQ-$A1C(?{Bm5R5-!lou0D@)C%or<*zNF7#J2R7 z^jvTc0z*ktc=8y6u!>hrb34S+WC2@+ry-|)O-0O{+8n_|t+FlEbcmNMy=@6Eq@ zsw7$nC2D1ind)K1TOc{jL+R~b(kD5z#hUSSso(y<8&v10=x9OO?TiX$BLg#^Tt> z4aK)mX~h#N#te;SfefcINi#5zA0JKROyBN32M?NoTzWvI(?Ba?Dm&JnPMEJ=^&1Cp1@$6_5=5^{M&bh!ia?{n9>YO%zHJoWAtx6d7VI^eB$ zM_fL)?`e+*&tFCSGra52rz4&gcjW1a*9Srn|9GzjNg)3SK~Bgq@^qoOyoPm@N@QGv z!Z;`f(vzf(1R=Tt54dI)mB7_eg&q!)p)d?CV~gv4(3hAaU>K;d3cuL19q5xC!np`W z!ZabEFYQ=dg?#5pE%=Lp3*c!=SVZ#XI_Nvk7(ucmNQve{KyO^b5zIUAJS*acs^SKU zA!poTP?uyO*JG@|3tgaNu*$RC6|hL2R0=2;%;=iJr!3lo*h>|COOVa1*eFN^8lTP=)8=MjQT&M zk5sAgYB;8MB!0HN#&|9L0o@Yd;Yo`e$f(3wJxe^d7pgkPR=wUV6KK1%!Yc*hI*M5jc^;?9?S<{+2?I> zw|ELGD(1&kvYj$h?A`u5_|*SdGv$JRw)w#rbDkC8`q;DL^Y-{FQ4t(u|1ZoB;|=KN z_V6t88?=f@ivT~7$(O@uy7N9NR<>3a*Ex5B~ZD2O!XLCM!;*2zS$101hYrb1h3 z%$D~(Jr^c87r@V&b}qDfYu&BbB{2S1PBZ?i^8$lkwv5svWMXq86Cf5dBeB>>R*j|U zuZ_j9^QU!)W+=3?s;{T>2%U>`$R@3D$m}x4;cQ@$iOGjNYuY1g1U3Me{pW>`;j&8| zgeMP7?|!hqPx_3td3PRoJ~>^=%e8$0L%3VHCx9s{zX}ja>@W{<{(j95;~wN8#}n&* zxbWgR@tlAW90P$oBM~Z7rE(dsT*4j8_td3(vIR9Ocse_JQag4vxMUTdYzFL#J6dR2 zeRp*ANOZM|uHGMAeJHxRH@bSTJz#pll`*}*gUELT<{5XH`Ia;-R3slkNg}x;sd~54SVpCFRzE-KP%t;LAHL; z-G;p%+-YD93Y63y58ZuiQVUKA}+O*zc}Kc432@ zB#tym&XM%UWq4e!=sS04kE=&}?t7qPr-rosWXHvTAwA$z5PF%uA)u&cU?i0p1O6I7 z`)4!Ctf5Q{s3~7A`b@vn#PsAlv8@g((+?+>tC}rlh~3kLppUYP8d;JCV1wX;5$_9O zmoq*C*wwr0C_?`yyGnt$LU%=s40E-tcG_{@VbIi% z6Z9!M<+Xv!f@WuQ4XZC*wVt0vG)Fx>0N6v4D9^O@bmwJJcnHe?${o1!!Z{5ywD|9A zkWJTWSa3gv8W3@PwRWpT9FqyJsKjb(mIUbC^!*fd8uzh3!<*?BU=mc8> zOCnmredWG6wZ!kRAs~bepgO3MQ~<9Yc%|-D)ec&M3f?~X470U)nQI*-dc)=n2!*V4+0wUc|H%UG+xakBO}UnZpeknngUbf1KJWC!$1 zDEyQZN@zwxN#8+Y$Q4@R1u{*@#Q}wU2b@U;g`kIGjCLz?aA8iq6B3n+$pro(Cdr-ft+669w)nb}A zss{Sut1Y3*MeN@bb;Y_O8KY!~B0h%Myu#b(G|Y_n>_SqA6=sg~+l(~~HPZcB@Ddmr z$*GZ?v5Jr=9MhdM-O6;2X!+=kgnt=MkFkQ@m2e8a&d!cb%_ref4)V!IG+V~ z!XW3;KZ*x1ZINXf9NzME(~A!q-{dRDeS5v>MZ=KUIM8Xjp)8U^F-~GFiKP^D29$n` z$itUUjEpkJ7$W&m`WSFOKCt%ik%Q=WXk;We1r>xSC)u_CQ1A0NYdY~+S5q5BD4_3Z zF>~}U;Gr5xBYxnCfa_R~kVd?>0B%s#XMo#c>8}MIo6LF~FF@eE7#|@lNf=!mmt6EA z?#I?chc-)ViX$B zq4QiglYMa7>$L0?8#RDut3<&kWwD=>*Yp(>q#EY=r z!Dtxj zSn(DF3RRnLI6n{K2IzK(FGwGd*E+m4O#`nNxD^YlNJHT#b6P@-*|9Ec2{Ynd#LP}7 zFho-M@mV7thUoBj$&gBZblgd~Mb< zyRI~|6?Aw9^t^^5X3}_y`7*j0CYwd;95k_b&2;3m1B{~>XAw0Q;^|naiDLthR${_J zXk#=mf?t%PGz;@4CcBo-S~x%q%9uu|lz|zvI>Gd_E@`lLL>^i3_3ts?z(9;5VlyPN zVSKcD7S|r^J!%*0dz{DosB1B#gy3pWQrH5iACx1o$1lRq=MIRnM7`|l#!6Jyp`5^d z-6k4Ycg!3&^qZ(83f5=Q)rL6@-6tp=VCI47kaz6we1{$pS#kk8o&~q#S?Pc;lAVg( zxfm4nJ#YL56l6}5`~5co^_Ut&mKia&5Rk%h@+nEnlY_ZW#%Z7y}|K+ z`S!0i9s0%c@D9wfHF&QWfKcE&+X=1XPGbWDHqzmS08N|p=ujM?T?n1R)~z+NdC?Mz zK{a_1NKeh$T~Z5`o3(CZ9Y+p+_8%u8Sq?&={t-MfPM~&jk7oj`T;g;fNw>4~yTDv3A%P3fy~AEmf`Z|Ic!1gvol?gO zVb2KO?}i(?AK$|Y?Jc}$PJ@n$moxU&jQurk+R6?M}Y#OZ)cyq&~S_ z+Oa*Pn|M1~#SD(7dW{}#u(0L~YJ_21K{Y5iqeY(0A8hW8synC!a(WN6k zk9|SzkE;dxGG1_pFy_Ja-BH~qT)*tmeh@FEJ}II<&@aNgt3l%S`MAe-hLMEc+{U-LdtZuF7Z5($X&w>XZj04ESF$Xz6?Kv35 zi{0OP-nr3#P);iP@I-&h+V-9~^&v|=W1Icf9rz47Mz?vuxP|TSY;WT>w)}oNRZE(j zKd~lK?h@Np_^g~sdKTSY-KuO+8-HZ0YMxIGT6Yn!!_t!&AoxFCY1L28sa0oQWZ0yN z`h8ZzPoq;+DLJp!rxI?mrBU?^#qct`!IJPxF+kG@8nz8h0RPwVFb+=B!Xcbpw2ld( zb=i8ZCgTuWmAh~0k0aW~tscfp9_^XXikf+1=c7a{TeYR?&5DJ%p>$xhFs;pYn*#u$u}s92Aso zg6k}!scNc*wi^(^Lxl)eM~DQhXDLpzt#PmN)k4Nn zJw-@{b4lNgU+(R2Kmvf-EaV5KheE50u!rbX0IVh)(#b=choYXM40gbS5c7cX1BMZf z;(&P!GJz)_!NA7<`OpQNTLTlD8g86)H4_kyfI%G@t|>;yoW|+4+BUGYKnCkKFqHY0 z^zihYu(Ill__uO`50fQW_|=?d3h^{BGcq%eF%G)XvC&j|*10f;S(Pdz<}{|SwqTCt z24^gWOcM)XaP63U8h;u`h@$&&z#A%UI1GEygOO+O5StH*V+p98%EBrmR9G!XFciSf z9pRa?ex#JKff3rDlk;c>=^Bl?pMbBu!3^7N3CRo_%b-9IV_8>N^M*Lb5IhvmC2#=l zBz}X?L>&aj7S-a26%v5Nfa&SUBPcgmJ%}6QnekqkAL>Z&zD(lGX3k7j&b12s7jK)V z7Ptw2J}>yIuBDi-JjjCKokPBkrGi77d#?=d=flI{&6S?VK{&WnAPLUngjTFDn$NF*pGRvW~0-3p* zBN#gJFUJD3PoMe)^6PGO8!z}jelQt%6fwE&{U2Gr+4aEV$ya^m1^X{e+}oW@U32_U z<5+ihGPW>r+kX}R=2qi0A{u8H`=<9jar&cg9q!iF72foQ$CJ-%dG~GqT6m;8`PhM4 z&-$01?oNL2Lf7S={j>X{$**5I_w_sPy?jgZxli`}-jZJo-ICnV|HaPcTYmOoW1Kk_ z=v3{fb>$W}66&1Y;7t_`4le364SpYT8{o!X`Xjw><8QCPi>RGKL4%wwNbW)PV=nU@ zG7r3i$3=AxXd93xhIvTcgCg)P=nCvo{)s?lA2G9{6?T%mvmo(39{qxJ)B~OwkCIvM z*oHu-0`~?F0o*mAS!G4Qk!{(6%sI_vQr0xmF*|#8HtLkyR4$bq`7M}gy`v+C=Ynum z9JdnbeRAXOCExqz^^tDbYM|ybgwBHQWGbWa9sil3gpn|~so>~y*sinCOh9v3YCUs3 zcjS$4p)bvVE(@6T+4R`R6jhZPjH6OE@8KvC`gA_!={RNJ33!{_TRp~u95bcyQhIU> z?Z7YBgV6FZ#Xe}LC@)O3P-WJSC)fWVo&h|t6Tjp=pSmtl8wankF z71K8^7zv``8&}^tbqu5aF0&S=*=6t$dV*$G$owzT_(7-3G<_@ym7CkJ7hH#!Dhm%CKp97v*LPrxqc(zy$ZU{+{rTFdUXrn9} z>jw^{f=iD>sa$I=V9kp0uojT!d?1gN1)TH>*TzF2g{{OrIPN&+RI?acvdt|$goBbe z61fOHvAFO&yOnmQ!EcT3&Wbr8O2jLUInA;5Z0v2bwt)Rb(|2rd?=fM*TU-Bsb|)vg zuY1m82d;D{L+@PDr*4k!NSz0PF*&*bP8V9oco&@1xtQB!jA3euN>CCYOxG9eNant^ zZ^>KL-*hM6`L3_VP9J_tw+?dsyj!R){h7y&9hYDJjTfx{$M=3rhbKP!tDSmllZ+nx zafx;IjczF3GYXIIz12{jZ>IRX@hZ+iQ+hBiO=?@;yJTM z+h&UOpG+#V0ZH?q3==9Qr%_FRA{JL;W*C%XHaCWkfNOhJP866mRz8$JLwiXhPprxUWBk*D9w7z+f;_<-9>Ii z15o9cz}P5N*ynbn3M*Y^Y@{D+ zMHRE2SxjYA5L$_r^S$y=3g0W3BJzrJCAFBk=pDT(nGHwdXB+@*kqm~}0Gbn#KJ;Zr z6B+3Qt3%ANUDT>>(Zfe)m}>=_QLC=ci+JStv%QvgFn24OgNkA+Q|K?NTocwjw5zZQ zhKWJBQ=pNGq>1tq-a?y%QH2={0Z0|$Zg<`3I@RiP<$Wh3+_3p?#pAJcN4(5hE8dHA zJjsoRR4ZcL?aFxWa@~zQ@4^pZIL;3VsT(?{JiZD}9uD$j1*UJiE^CB2fN6n|5vzgU zbzVbQLRYM|0s?6%AaeZQisKhZeMKC|MG>c=0umHRV+90dIU-wq52m7){6bmk7Xu7D`!Du2rFoZ+LiL}~)LL{P_C_~;v zTVoyiqaB#cc~*eP%^+%odcfk4P7_sG)IXEWBimFizrk)C)HfrQ);58nVAw?S@sj{j z??yuGL(8^-zbM%_X;w&K|1nbiz;aRX6KD%bb5?H$pb zneiy9$9A_qgnmNd`x8;_A7_jdb|$=zM@Nn9>6F5TVt*)ZwKorAC!(}nfX=cMd#rl% zWTG}8$A&ik2s=f|F5yug(F${O%w$vG5rw%WQpqT1Hi)(m`+QFvf+-jlfO&!hzm2^T z!vxwxCCi3b>rReCagRCFX0jRV-r#y+-`!_I!#$ptE<0q_THB^XJSN5yvFTz1tT~O? zabyn$6HZg{c$aaB&imDKmcoa|QRtYwlGv2|s07w>U zEJsez*3(|YfR&{XU?L_lz3>_DWEStDL$w4m0$6`Go2G%LtA3kb)CI5$YFDj?ZJ-PG zZWbYAjGEr5cTM1&R1*o@8l-&A{3C?4$z6h-rWX?yaziu zxZ_V#w|Rx9*4yG~!dbW`3rdPihX-&JqZBvd34m8O^II#@)dm=ovPZ}NPsOrEd5dzE z)3HmTqG0Hhb}4xBFrk3;2^<2rSpBg|c?)m02eLN9nUiG*TH(+D6uX0*81ZN_bs7`k ztOQ~K0tOf*j>o`40xGcO;T=~hnv8;!DA636;K|)wCTmJx!r@Vv5=use91*jOH*|f=9gOTVYDMLDeT^E@;r+^;im8mgz!zFD;l6BEa$-Co4*vTj;(v=9|yAsQM}@3MOllFm5q1Z<kpGb3^chIWSlP_&1H(EWy*polYfoN@RP^&4wC_?Yf<_qT2ib z|8lelm@tx4)cqLfD6U-s;|mz*rgtnq0I6U=BO_E|!nRPQ{w_+kFl3v8NgN1C3d=B= zstpPjVB!b5PNTFUXzWMMJ;njJ#_e>tAnDUF)ZUl`C9DbH$%T9++hCGfg_9(fCZ9!T zC7rc!j9<`~qH~tcOX+YupL`n}QC)Mj2S=~IwXAv$6y9Rl1x|!C|w-8m7!$4K1c(lep zS%hWXWuu2F4%3lMP{dJW?wg^JL~bZ)q?AEL_>zs|E@d3SYx^?2=*joiGu+1^XM7c2 z-Is2dQ;P%R!MP#o?(=L9o??3q9;2S6XrOb7Nesz?O8n3vj*qCg>0>8~QodOu5Ut%0 z1y6_-#>d69q8XB&hEX=7m$|UZjN{BIECD=20B!?@qZ!a+81?|H!3kugG)$|vUCEUu zqTaw|yC35295X+aThGLl_VG%X2Bc1O(41}V^WD3 z2a&u!o|X($k=86eOK{P1SqlQb)NT9s|Edr>f()TUM#)l+ENn2M>`R^5eP{>JN( zW(&?nRC6NpuafOOjSd2%RL=KmMU}Ok9b39KMCrCyP!F>$M-CO#>Ti_yckI^`RY2!X z?Bf8HO~%J^xl|mIdVwGUC$H(3N*fPwG*<26OL8OGY#s)+l6L~7-i-sA$;6)2;NVye z2=PxjJc!%4Ru;JW%$Ghf^xx1H@ps4+Ab-0N2=F(`W0(Duzp;^C%b{M8sgnS`!s8zdi0CL4^=o zt&zjkIr7yk+59BVbsB zl~-`%^fLn;WgMAG!nz$PU0EZ)gMQNTl_1XTs{Gl!Acmr60xeIYG_Jw}#zCuG$ZZb_ zH01zen`Z8?7NEfnl!@sjkdbQ(Hua6LX#tfH;NVTLZ6U6Aox9#u@2Y*eu-Y;1u9lmF zL)VJpufR{noQW}9be+dY;OK%6ygGS$_=Le?-v#oKpi89?gR(uGRm50f9B?MQI1$`e zf+-)~jzOa5JM|n&<}*GRPbRkyk{sfK2z`b9CPoN{X|9lGc5%XzqbRdIA0CFZ2U0^C zVy;Ds_bGBYs^bBRoC-zS8k~M~4dc9|YB=dR=Qj4|!(ucN5b~k4Uo~@ZC29>MsIC}|Kk6R3R?e~LGOR5Itrdr+AgPT3d%mijYfUT)-#YOv zhi|?3@Et3c4akKYKh|VpF<63Y|8Ob#|KTE7IIBW198(4KV*#CvJUy zr}yl6=Bl692faI;@*e2$X&xSL+XWB&sjaMjno%% zS1)pW;e}WKf|Ukeku>&~&UsILf!dw1&7cWgu`cPwCUKk9a!`*}+ui%v#C*08=*SkH z+$lC;p3c7?JZp7-$&YqeR-BKLJ%F!h--~HuAit&9!g+|%d=_WPD&816=+Uih=iFS8 zv*n*d&M$G~{E|}6=lrMYi{KE)wt@zB#jaB?#g6=-@#on7u7pa!(#?>%=J_NQF zj)9yseFt|QI%0aqA#Xd@3O zmu`8G%N+_0)-?|?$jV^gag*uWfB4YeLfF0;qR`OMv58a>ve^U9ZLF|db8wMxG&Q5wHhW|6xT=5NZvc4-P$Cu6(I(F%2CFj$j?`c0zP z8N*?)7~jL2$gq^W=9OA#?V+VR} zi|&l}?7cJEd*IOVgFAVVMqv@T-z9NjRJ8plQ;)7rMTM?2zJe}EzJMxQ|JkuaG&hV@ z08J`v$`-o+K(}F^LeLkj41TIRDW?#g^z0}sbrO}waJVgfM^|U(ONYI7j%x<7Tpf^t+#Og{gv>`ox+Au^!ruig!~u23-f99nJE}I}7rTDT zjhPYo*!;ss_H<6o){9soDo%X%vsx5fZI0K6&Yr zE#)*@rld*U{lhQ*J(6WxKwlTZG9^C!=5|LaT17jONCYaeX>?@P&_{@#y1^w!DmTuMIi!qK5y=l|qV zGWv>%pMR)Se%g(*#;;X;jo&3zdykh*qdjmU?--!*mm*k%8PP#7+-Zb%xH(p9! zfBo)HeDb4T=}wOQ_#2&noBhG1&j>ZJ&N*&!ShNV4HqclR5epRQZXX&)jkL%;oN+Z~0GObKu*PXu7xm`70h>`R|vK zOW%6pM;q_|D@ZR=oe)V6Ml0)C>iDdVE<5IHz>%UAq?-0KAop(Q2-~A^j z#>c*~H2;?8{T#*kkGsCxw-25B#S@SHWW4cx-N|*qm%Z%QfAGCa$$z}>ZTr?e^`lG4 zQ|lXFuy^mjb|?QhVs^jdYiOU)JNEr~;2W+T$xr{hZ_gM0>c`!aQOsZ{*N)?Nc2c-9 z4`E_!g^J_+d(DV-OVZ=Y^1S|4J0>@9_Mb*t(WsbgqK`TmldIwMv|n%h6?6G4oxi8U z2xb7xvP|}jSU4OM=8F&~QAD6V*Bpgf5;5)Pwq5?ZNOWFvZhc_ShCD|YZVfO(AXVjUFn^w7ZAiR_~*)gT-XNV zDp}E^v;adb6G3sbSLEQxX2>&ift?;51Y4(q)(D*h*gxV|RZV^#HSd7~!(tYw-7qv1#8%jF4K+0% z@(Gt9grg|TpEM*cf|xabT!&=k-{{CcB>9JMs0Ddc)XXf`-;aa>OZ9SRZLcXV_V z45GGWP!tDBqP)f^Dk>8c!)_MYY)Mj!!=nRvT!#QkA-5qqmVTA;utCbN`P?XQrRpltbRIlI?NJb(1TWNq zR$+n|SK-(dzo3uDpxkv3-!#J|y^KU*y)#>EBxT9-lw4~|xniAPJ=jPd<@j_Y#ax~#ll3C3D(742}$8i>SW zLtsK?jd9uH%^)>O68S8QdkY`jiwPi22V-s6M2FNEZgo?#gLwvV&0h|tH3p}e)+n_Y zOZ2QDcAPHM!L^X;;7wRR`Gd(857$(zolpC`JjOf~QY99z%+`3!rBc>(_nD+dGL#2&339Oszejd>7S6tV4)wI+U?N`xPq08KK2`4^IBt-_{)pLfuQ zv~?5p6Ae%jaKlJmL$UG29$00|Tk1*d&I1YRR-BYpFP2!RQuB)LLxCI0+QNwuM+9DghV?HHD|983}XS zOQ0h%cX>Hx`K*!clFT*|_hH~2ntC%CcY#V-TzPY>(D!zI zU@FdIV0AIJ;J4NE#PJ@-dmQiat)}lRPPEXCLrHvG{Nv)^m!?uTT*!n+k!BQW`XQ?U zq+kZEmZ+W_^@pI;R|-p=eJkh32Gq!fvrrJ{;zyerA@~-_I|RB^PrG-me{ny6$&jus<C3 zB)5E>>>fsW6`kK_j*WT++C{*l$ifWckOgS){@B*CcT3sZRrYQ!d5v;hBWA}nVs>03 zwpky)rnp}pdlhEO7+#^NsuwqY(c6#oa6sPoBMr%O*uiB0GlqizW>!VAH)E6E$6_A? z=9e2x7zaj@pWwQ*6&bi{Eh1d3*23?Nx}zTcYCW`Gw${DDR}Crf*jJhSoF;pRNSInA z{E0|-KRhB5mW2pr&9F?J=Mck^LxX36Re{vJi8vO-mcX%a3_#ulpCV~0=|Idj`XQU2fYArE5_I&ks{RJgEY=O~UFVIvmAjmqUwD7*_Qmh0%j(aLk|dg*`m zdKm#Lsq~U7WUy5YMW?EV>*L6N+|K*w>|Bs5F4qF!#qwWUy{R};IQom53MQ8BuRK1? z9@2t)ZXQ@h*;^1<(cmylmBQ#Z&4J4sjFie7%z%Ss0P-J)K^U8Ao<18!=7NTi>?fN* ztA)szSP0CCHTv{ay$0{-_176+E}cJAM_ES*&+%i ze@fUZ;lL0JwCz$OseEBO344dqjgv4g9bk*jrRbbsiGa3|Y#WAd`au!j~Rm zw(ID~RRF>__Twl5)ZI#(2F3wSTVTBxlCSJF^cf7CzyPr?$1&)|!3FZc2Mhp9KEx!~ zT?1|J6x*tGsKjnjmty-a1ve$9P(xR12s(&ZfnZv-1SnwO4j5+3IiFU*=`s2Vj`AHImm!!R&0rXBfCFxn%t<11vfn~6FAH0i>lxuP zfx;v688|>@2m%(;2$)=}j%%60ZPdb`h)7^WIk1i3Vm^J0#DNLIrttr;XiE6CU58_RkR5xKjsSSLes}?D#up5PoYIdd12j7s%AxbXf8T11`FEYVFP3A$8CVZ z1FzTt^)P=mHv%DrI4z(GkCa=nkc($C_Xn1+plaPC2=cI)a~;FuR?ZbPkiDaLfjH&?fMgcJF# z$k{Q%QTj6RF^Csgdo&U&Yd%(MO=iH#vY?D5qC>ygL`M<3mo+Q=M!z7LkcmkOAi6*~ zq1`NRU_z2K(W{lLQRmSkAOWODEi|Yf7b_$UZ0{-xQN$ol>+5yN-ux?Cmxly#LZ6#> zCcFD6davRI?03CD`T7TvQ}*Tg(&~C=-erIC%WQh8^t&2cL>39tDtH!FI2R5!Yjx7A zw=I{i-?nY1wA=-@+#=m_w(W$dNc^pK^Ix^yN)6-5rK{qe=fa;Hh+Xbdk>5~~v`_i=dSKpmD)P#iWQK_l=mE{357Rut|N0}Ul+LPQ;!J`AY~ z?(?{nSyXmCXix|@K^e<(EVQsk1{oM-wAnXYW942WW=DRr+{rj%?m|z?5Uwfe+-oX3 zf_rCSrS1dV2TsN-4mjfoDrH!fSr)Iqc+ZbsI3D8?WHWhu&!knyS?8HlYl>k`iYb1SQEj_Q-`Uy`CVx zAA!L%#PWy9$bM`}kA%YW# zLi|2(JWslfcSB#^ZJf*aU^w#}2;qRLX~pP8n9KMyG~?|s$iqom;|3Axm%`6))bSnn z5qNbCPMw=lXPkI+)R}tf;7f9& zb%~715*4nUp#&9($2LJ8CV72GstfEU+}|NoZY}4tU<3BVg#jDxb;cJNM2j+}w~!hw zGrX|~orrMKnYi1>qu*m>n7Ph}c_WWt-*d??GDCi}?_uXB zG*nPym+N~W@)aEnn%&?&kHb(n6G&AG;GD!XwCM>jNnlS`x={!%w;4{A>?#!&9IZ|h zCyE1pW{qZT_*J5RENPIyFeAJ79O~TC54#eEj199zzH4Wp?LU+^1T#!y&pD{3B{bth zmVOJ@@q$xQH53d=zjEw<3YJ|BwNNtS*(TTw#ed0P;DOnyrM;>Mj$X$2D{ty|;YC{1 zqL0L`ptXI=c1^**i-P~D{ReS2nT-AmSlvO{a#gisDp3s3J+pRd&w%~2(KdErcS~Uz z_;M9+=wQntQ4;v07ecz@zR5X9?8IyVdk484ELqMtMC&5uCNH=1v>3eJP4QK`Ft28F z^GZAHrqEjBbP^M-;+t+4txW7_XBLk;sUT0}a68+XyADw#-1`9Oq;O;W*!^J1o97eY z1KI=uyu4|Oj581?eb~;hvW8hE+!jr#YTY%9pN=r`7`F(iaZyPGc`?HRptZR2{Y!iXpZ?bsKqnpvt*9>zyjlp z`7CqhqoSbEmdl*GvJ4DPr7>Drz<5NoQgO?WEp`WYC9Vg0+Z9rqs4=#-3> zU!`-R0utlksepvl)fLLXjSFFonaXj=F@)F=kCqvo7eNLHwr}H0len{s>bSCN8oHUuK#!m?Q8J+*8a|I5G zUzvNaCFFH<-azN|bjq`pn2zyAI&Y%$W;zek`6ZpI>DUT1$|2m}D~P}kI44hVPS#>R zdYd3lfx`LI0e|a^zFOZj;nTnhE35ISY0BW3hFm#qp>ikQDt9aR$u+(Tr}>(amQE>- zrS)TTjKNo)635aUvg3f8JS7fpk9KfyaWdSFLk=?ewjS!`RB(Ag!%m4)r?qnIZ()C* z5~o4UbL0cQ@02+6RI?*aODT?;k0IKsF~l2TuM11FVDF;$nMeci!|59#5O?l5o`Wj1 zrCtq*BDA2Phn$SVXtQiy_9kTp3B=ANBU-;{DX}s61l0RViH#oP(;U9j@rLW{Ds-XzYa!i~D?VJc0Ss23ZMPXW5l`1HL>HD;y6bf6bGh>;~a5rg%LL;ALod3D~vch4$C&5 zWmByfeGBvgSv`jOVdJ(uhWJ(vYm^bfi2vqM!GJNX2K2!nHGhhHMP^3YBkUO?EOYMD{oB zU=0i1%G{50pbGgjZ z4%r77whLgq#bPdkAOeybI2kaC2qoh%$l)n5!Dw(ik(b_|!PFx`>d*8@ICBUghmTIl zK^uaRZ%WLmIcSrr5NvbM*rmsCm-wY8N-}fl!HCF>9cs$ATCt4!?Ms%UU`_M~bX_+p|(aSwR z<)e5SE8rA%?1nLO*a)?Vby(c2&bB`>0CPfMv5SGxu?$RSp2(Vqwsm0h9=Oh4({&akbr5Gxq95E7yv9UlZlF2>Dp?y1w>S`oKzXXUw zy~p=-8z_66^aH=#_Cs$2YqVd-X)MEWKEoNOu*8WnI6vr^7DU&{G!kyWmrxUQWd$z` zZ%{lFt<_gS$HKY~5>^XwS`0cCMXh%MzAYvJIcyNQK)*tHa3h3)r|%lRBM;u&V(sjK z8#Rsry$tAGcu#F8;@$wn-%t$Hzu0U zDF!gDsGy9viz0|?lCDW1@eTTE8Rc-v>3~lVO7KvHfIJ3hl zvB%hEm)K)GTv_6oTN#?m+wfb0mWFmZc2wH4is+?f@lNtzHN&5GVqgc<7I1kyP4rf1 zozsZJsnBgHTz|D>wwOS+H>kR0Qn?{;W%Xa-x|!11yJB^+5-*WU4Xo8eaaNA8>7DZ^ z%O|VOcQ9%7aud}2{A&n!DZ_7u&vXwP!qOM!0mu*+%zUZ?yxSE6Agb~VHJ zESL+8kko;7k5j}rOYDxg17O{Q!HbYZ2#I)>fCUkP6fXcDud^(r0hDQwf?=(m#;PDUXEMZ)!4IADvGu>a_{c?|Evi&J70!Y=8G^ZU-;~Ng#}%B|1;J& z#eo@1wdi=1F1aq4-%`MKU5Abe>jA^vTg2cTwV& zdPdOlHq=XLiJX!$2#;Z^F%R3n>G2*+N5ZHj~5( zy{1fBxymeO&k&J@;R@ceh+#Q702nMV@D=Aa#>?Kbnj&WN8j5MGNJcr55j$IG=ZZ<) zeLZ|6BQDna~+l)3z!~)zV12}g~m%uDi0MZ4x-I3E2i+8POtNgmN=y5De7^mKR5GwxAnU&E56dnvb$ZlZ$(Z@6W zBu@|$o+d81I#c3?t$ZS)ufQ<{wHDeNP_@`~EiKY5C~?JTiNLxD4L4!DK%E3jM;>B( zLD=_Q^cV>a4&3Vk`&6T_*kOf38Sa=)r8i>^5o7~H{kWzhDv>Bo)zvC*K{W1E~t*IBtkR z1GNp%zPOo`%oi*bmFaORr80e9^scNt>?4QRm1oWc1Fb3-(yUbn6KBmOGBp>gJXJ>^ zFF-SFYKdfeur)){nqj2%lcj`<*maG?ct%Bw;rt3(ER5Qcy(z~IDc6{qvm<@IT}nwW zB`oD(F{7xsG)pdxQ*sHSB)E99vLv(S0<$1iyO+umw92D%swmHNc}a_fR1+41X1#U% zimUFg;gp1UA^{-tza*X`=c$_gF9!)lOh7Ib#c|~cO$8W@FbkeRs+}em zQln`@Ke5ycW?`M<@L?MF8+Iv;A2K};6&q{07i^@joof^$PSF94K;98`;x_9aSm18`3{VXXptqgp^wngweU*t@mt--0a?*vGXT z{~6V4$!}dN=1@_qJa&DIn0c^fJ8!sFEyeP`Ong!JmWyu%h%ne6QY&%nZq=YiV5@;K zh$8JoQ?UvCx(5B)Tx{m8jW5P%RP9#0b%@sh=z4)}u%H_S-v|KSA@EHWyc6HtEctiA zw?%wg;k!IT^KEah_}!fukVZB~vLRp-39wbCmbS6<^h{BCtK%5<&TE?>LF z`r56M#;v*$w!?z$q?$3uliFoLc1!I!BHU^Z7jCuJ3fWf?a+`$gw?YnVcV(WGGTd%K z@GbQ$Dd|qe407Ik>h6T%2!15MstO9o0f}N=t-GLnG%YfqO#%<5` z*Urk<#;p{zCtRxU^CYBbg-l>WfdXEB;@5*2u8A9$Kw+63M&L7 z+YUJ)ArDz080mJ%of7g&D+D9o4!KK0e%}hgOt3@lmXLE+2xbODq$Yo0!CqAXd$k4o z!wT4IEZA$UT)FDP1~c0CbymphD?%nD$D-aoQtw}#HQpnnXgM?Hek7BU{hT|_OY@qG0* zRo`n z+k*YM#A{RU0qngN>`@JqI-+g}`3o!LeHAk#T(SO#jjif0<-31{aUW6dFZrbm+|js~ zuv{OIyc=bgE5SZ!!5*uK_aU_5hb6^F9N&BxrP8x9Y`iDGP<<48$R#W9$MNkZz&_@L zNebc#VdE*gJRi5xd;)3M*7dRnR6r#aF^LH5iI?`}zqV3+(vhkHikNMZ)2FPMe~D4Z3hdJs>@%qElTzQ$!uL7x{T+Ot7vJB*_XY8N5xy^p@5@%s?9H!Ox%`6# zV{d-dg8icfGXeW23-+}N*gspae-T){`a0r$!-D;*z}nO|0sEE(d&+`wAN{rk`;G;> z4A^%q*klFl-z?boQ2Ot~+?0jQPPFA0727DufdzNp#6p2{R`LFTh2RVZ0f9LO8&57{ zL5pCUoIDu8lPxTIadEMQBP~%s0Dggc02Y8~2d;RH`gh0sQiL5u-_{u4RsS)uO#SCX zl(7YJMg0(v9|`2o0I4-<)sF%BFAL%W2}w-np1-qKYC^Vxz%F-Ev=|wR~-5Of`WD#3BQy$8gu-W1^aac>}d@< z=TYCmZ3O1RaBXrJYd@D*10N;~~3NBC+he2pD`ts{J`6;3Q)*AIAp#rJU&g?!&y zNAZ1N{1xFHRygr~ogSD!!x3|{r0lXowqX3Z)mAIsbr_?X?f^^6q9>OZm#eqx_)70# z=dCyHU0jv9S~eI*t=737_&)bDK#=ndR?at~1ez9F0a|DUb_BF3Xd(JEEkvK*sminv zVXN&jl~;jkT8QyT3pHYXm1!YDr=x`kM~i?KqOS|SE#lh>-*w`<9;GeQLX3s-gBGH% zk``*rZIu=vw%BOAlPxfX7Gj!PrKQTW5MevYKS4#)DU6S)DU4~74;%DL|7#?)SPRv$9J=FJJt|AI?L1$Q%p|{F=Q$= zR1x2%hOC&(-9ZfzJ{>hg_zR`vWon49`zv6ih6sCcMGZ*}5%y9IleR8XLkyXo8e+(? ziZo?vh#|)-Ldw(-LrzqLl&K+xOiv9l^wbbT^-fZvkw^^@${j^hL-UQ-PN9ZY zI;f$=z&S}Htpqi+5^1ZcA;PMuA;R|Z3$GId1Qzgnd(QAdPTQ%8hVQ%8hNT4_l!kvbv_ z)DcRzV`9}+sG?Q>OH@&t@#`w8Xra+qMHMYHR#Z_%3vn0TGoy-7_i3miLYJaN+D;WQ zhRReCTWU6{h;gc^BEo)c^&;mQsUpD8;!dh)r-e0>D$?OzNmau7cUj@2igdV>D%x#@ zlPc2TVMltw=F_p3_yH@NRFMvMQbo60;iQUmIC_aag!O^cwr1L`&CJ7gq~4 zgb>-oRx46NeB<6#qK1|j&$U{mL=6FgZ{J{j`9`!1?jI*#5YRiYi>;A&?19$%wZ%2V zwZlvBc2sSH+xrmoMr_!zX<*aEJE{&$Nqc-6*mUuZs)whfJuwZeGTyb*l(`Hmm!$-& zEK3PiS(Y-aT$U27vMeQ7MOi>?UQ=9CY)$e|!dgSs?2CEZ+(yLNaf#!?8dd(_9^-Ye!jS-WeCE2Vbr992_lSAtd6 zt^}*BT^UxcNeNb2yArIjb~emW;u6$ZyAsq{yLPEv*jYPuZLe$XbD&31c-=GB1I_kB zH~4pEd(`w~aH$jQH^gvIp-{W}?$sl!)$0AL53TOi)&OZ7(5x$6kH^a;au?l5k-!x> zDp#=Y;X8&aJ8)yi>R!{^yZeCY+kN=Zkum8bSCDcE(N=BK8e`o5Vn9Mckt?L^~?cOy8~pJ9~AOha*dnrEky9qtDWS&YxaO zz&_mWrvdZ6xn+sVXvTN+yYboXC9W$?Qw^22dyW+r^5!30LoP@E_8mCZ8|AgcF#V>? z0PYMta&RwHRL~@@LPOl(b;mSH7^Pt)GqCHxu_FiWvcj({vTPyqyorv9w#)Dy>YVBC z*?HH2UVVIJKN`sl_Uygu*xsH!g_hF2hIV^)z?RwMe135i@BhdrXip>$)wgl{Z5$cd zsMN8cp=~20+j6;%TsxcM*xt!KNMJ&DZQp@|dwX^s+G{M}C~VC3!}!N(ZIwg{EfILt zOuGkthJGp3mo@G(uyAdLj_p0Tw|Dp9BZ%Qer{TiL4DN}CzB&x(=@cU!lHqC1q#FjUAm^3QMIUp!JJ_LrpZ79n7Dh&N0u^qdjqR06MDo&J$*c<%#ir z`LD9NHTrA1l!3#$Zr{86m>Gbf7TDLma*cjUV07)!(8vheZ@p6Vw3`9k7kS%WG)HBV z9oS>~5AW(du*W!wQWqMm{EzO8u7wsN5*^(c!)J^%L($l*YjjUyzzPnc&I~WqMk`8Y z4H)AdTSZEYAOw#b?mf^MyUVmZcb7bUX7I$`gNG$Zdk~~O2x6+urBs_^R;tZ0E0qK> zvIJ>QADUd|1~Zie>WbwsZ9 z#AQAyBbP^U`=JjiG2(-HTv=3=)=IW}=PsSBBQ8p${fYdU1f~@g+Iz<``ZU^TJ__A; zTwR`zmg1YC-s3$xdu}_}S=cl+<)JK8Q>g-vi}0pNdg$8{Gi(X}T> z+6&XhyL*oD`liNPTly|e0%QkLd0q}b^}CyL?HEp02`ghx`SMj`PWd)tF53&!o(8&v zio&?(s@##*bPP#^&e&#W)UFw#b~&SNnIY;H=r_k>qnUhX;Tqhol98ngE5+1u1hHxy z>G}Vv`x4-|lJmU5%wRCM2fz>n!E1Pd07zT|@B%M@I7sX+0WOH;?sAtq3*rqx3@)-~ zz$HgL)=HF9P9nv0;v|l(O0X%H^PO!aePFvRZ(2dljKd3s_4sh6;Ex}d)lgyiuRi5%XSu5 zeCnEzTK0FCV)$HF71};4TZZB=AMA zz=D{Wz_+Uge%D~~FFCiAUSG*J$_Y+fdpUXP^Z0(lVEdJ^k(aN(b~RstPuMUbipFj_ z0!maJBsd5_m~*2qQyh=VC$0`(9?RAxAK;Y6lZiPPhIzER4p-Ca5b9@}uotK4M$eob z#iF3GQy48A8KqJ9043961Sz{~Ik~!kT^_wQ>)8QF=R=UB^+13J_Iij0~@Ka_3 z#=W>QcNP_GBI69}0qMt)p8rOMUon~3U1n~KP29GoC%-u023_y{)C8=P`FDwgU5Z zYc!;Pbz*uND|aT_gp)vKU3Adspd+>n&uRH`9QdzZ%eI%4&1TkC7UyEg6h384=)_ zF+6xb@pypkcHEt%n;4u$1QN#%_#(vh1x6>I)Xt$q&GMC^ESA>SNB2dHrkN!9^3k!G>rjHv2WjBpNv#VeXL1}h?U$CVjQ0+l?i}bn zbJ7DOe!K{Dec~qGY!>Xs%xo8?437sU6|+-$p5TVBP1h`0+87jNp>eKsPfb$Rf_+VQSncLhf~B1Fru-nW8`Z-$JETePg_)@4iR zJH|bE>l(Q-JSn!vvR!xJfIhpnDyJgoAwfkSdPP>|&T(0pC3W+mYv|r&U-r1ziQ$~K zt$M}oBo`O%WVWjnSAW1|-|cy#YwKXUM^2eXcDavqYi;{XsTNn?k7wq22AJthnBTrc zw$t4eJ%?$`KYpqEaDblYWR11IHZy`FnG5S4ES>_n9o9Iu<>$j%3b`HDer#ILhqV-P zx7(}DBajVagB7%%TSsZ)J%~*RH_HVH{*br~3Aa?HT9&98r+9Zm#C%byE zCm(kfL#z*wmxm`N`&8%|kZ~I_t^&`1?6)ENmH!!#12*J<`Z+F*R<@Jo^`)KgJnPuN zW-CeA=~MrDTR66naB=nX+rqJhgzH!TYFjwAk8lI(j|Gm~fo-81BU#kG;Pey*@lUq} zL-Pn0SAVoE7+OcLe)Wgjf}wE)8&LHe1wf z04n=?PsKY=_6-cEA3=`z#Z#4xRWk@Ps9^AK&{p+J41STpFEIFX27kfee=_(p22BKx zGH7JbfB=WmS4f1qa{fE?kG`jbsFH0ogtXYit3b_SR+7Cuh76*(fF~&Ky(I-duBt=as!_dTe-f>b0Ue`OxA5{4OW) zLTaI#cXvHQs`ldi!+h8|+IQ?XzMVu>{5LQ{1&_`67Li8wM(pL=DbOO<^7pVK7fy0zM5c&08|wctWha`r&*S-Gg02D1erI@dT0n7 z2jr_bbV7^0?TNE^-l9P~4$w<%TObY3)d=GGa=5RCAC4ON{wCxB@R@}+F`71>8LPM_ zn6Bi@`Sr!XN*0?I#Pj>q7dh;3K(qpf*a+bf5(nF);4@E0Blsj$^dyz1?;!)G)RMdk zK%q9mWsEjj=0-QEgVxBPOv!{x|2`E$KTE~zXt1mZ71LB!;59jK{~(n{C%5vqUB++b+QDOR%!~Z~?yyKR#nqU>*BmP$fTLK~E)|4}>f?bPdK-Ty<4{CzO zC^2mU*}XLL19+YP#Nn0}-fJK#7@YZJW|f2?v=(^-r7n!%4( zSrDT_t7mo$-Zh_@6MO+!c=LP?jPleV+R-fSKrGOVc7A&`v%5fXEI94(o?nWLNl@&- zwZcXhv@}TsAs#Yb!!9wB;voe-?}&t*_kfcH1&& zo;WC=LNK<#Af^5*?{!x&pC{D>LiG{@9K*hMq_W?@`J5bdPjP#;F^b zpdHA)?<(qKk7E6lg5%BuI)CC)GM$E0L*fQ28P}=8D!Rl5j38L>HGizO@oPUnB1|6` zteRvirWnYr^crJg#^w{OLL+#%*scQiT}9i-`Oqrt$*!;D>kAE(VJphu5Wxn&xVkTw zmzl4V@lQBIov(FbrgkC8y7`C)JN#whA}!Xi8q@0vPNiwgi-hl0^fq4uIgy0j{Tume zJH|JeueTGl2?)ujWWH(?`D--^2{*9 z1l;$)!dtI;lS5VsvcO7${VFxyRaNL_svcLjw*l^}fpvdu`Y_amBN#4bZrq!XPbw~8 zVSK((4v1g??v?q=PX)r3vRGLZ#@6fP_)5c>aL8P0#i#J}qK%ivy<96kY!a>bz{Lg} zSa4H?qeAEnV|&0>jK2Ujl=7n7Aacu*i_MjGa5P|{iwqBTr~{A8FjoaLUihull`z^= zp8aWeofX#2HHUSTb_=`Qe{L3-Z{I=T#B>5Yed9Q0Fh;R7_c`1oMj?kbJpm@})TC#_ zQm460&M?@Go4XoDtRrMrkRQe!RgeMv^nnzTU4u_9uy5=c?bD3!pD=Eh0tyaIB|T`kr#zEU4E)1hS~Go)L_TQeK&N`t#FIan{rHs7hw*1Fx&Cx zz+72`8HWX9X6{~#eD*?OWUK=78a^g~4aWm`$%eB{9fDjm^Q$*%U`;@J)pq-tn z#Y}f<;A4yVZs6*5rR^({7;u092I7D%Ow^(MD)SUrK*CYN)M?Q0=qQHR>!CeJgK#xC z1i%M?SgHzJ(B{UGq9N7Dv9O}h(AZ%$e})GMU#dyhC!( z_JCbBTz6)tMkyrmcUileeO4{BqWtOsN;!&q^e%f8O?t1cPql39`oQ+DX+!PPk)qsP z)~=;;VEJpX%N()xTdlN!fbnb9t-4w7!aR%RcB54XWHt_9Hn!q0){5ibL0IG7fmwH8 zH{gMN_oF~=kKPR{1N09jxw+lI`Xk6|5=@73d(m>XI-q=+J$kQo_>muWsrTjf69PxJ zeX=|5MHpp`%s&1izZFNjn2uo__d6>A;{x4qJU0etw&MtXarnwL=mz%MYFqN*Iq$Gv z$T@(VTFKEe5(llLoE6gRBf1s!+vVAZBs`2=_8{O;vK#iUkUErU(}$3AR3EfD9|fT; z)vg=#Ve6Q690%hgR!mbjBSu%Y`be$~vj0_R%Td(mAdXx+>bHC9ccA`J-SL$AN72(x zjIu9v%<8hbU>2}ej1v~p0r8L;$aT7q$4elcv;%p}g#^YToKD(-Jg$%FPPG1*J}$q< z(b6z`+k@VA>7J+bwpaH87RC#?vk4}q zdjwM-F!gyb#eu0;$De|!Khv<$XmwlNANX?vy5FJth0&7Ns)Y1F05!tYpN* zfM%~ANJRAj?BvwKN-y+cus8Ik`t(3%m+qID*M~V1c;vU@*hfxa)Z_X@s^4l~3S!Nl z%naz0Tp=*)c?xzq_G9&h)Zas_A3j$RBQ|{sHBM+~@8I{O{2n6?#xOryt%1*!=T2u( z|4|v%&zanyJ`H&6Si61tjI74f)(O-$R{TyUPN4jxT|Q`+6T4G>QlH*9qwBdMaD~}n zo#L9tb*9xijZ0dqb!G|Iv&TX8`KPq*_@>sK%?#SDJC{4J&$_KUCr>}?Y28`0?yTFo z^LBZmb!QVN^;xt|yxTqw-ut|)XzXpcP@Us$hTrp<3mX?<>F$EQurWma%iKlucGx;A z*e>C_$md=(znEwEW{6cmJS{i$T#y$4Ih?6Xjp!lntN2P>6ipcB46_wZj*jLo;R`pa zFQFA9rd?>$5K^OhBsI1%ZVl^EdIg3h-b>68%L|nomvu8>#`G9qu2`ci2c8K%qQ^wf zW?b4i=JseDExOE}SOM&)VOYQzM_VpSJ80#%3una3*?;Js#a{G`=*V1%DBg3Q>!#(( zy6N!gwxZNxIb0**GRMZI>R!@3) zaSm6xSH}{4FDK&ty(i-*y5}LampvdXbxhdJgGbUvIi%!%NY-)> ztFFwZmjq#7+#rk!LS_pBcrcC~Nk3^WZrGKDVb#4dn|abqQMVXRpf2{y!d2S!C6p@2Ixgr z&lI&3Tp1LXd;=rv^Zm~HknyeX6%O=a25@xc`{*GbDi`3>NIvNtqcoq0&*7duw|*P6 ziPhu+gjeoEpJNu1Z0Xnm^G~g_si@@hrT%>ABU)xx--iHQo5#=OVhk+KMVcm;bH9Q7 zo}U|ndQuoIl-r)ky!NS~56ji-ge-VMJ`cS4Dqa{ z9a0S_f703m;jsm%w=5@cEML40nLtR*_3)*MKGTy3^qCwGdC$n1F3wS(MX862$no;)yIOb0}eLelbInr0y8<>89dD@@SI z0ud%4quwJou)ySRpfmXpq)@=QR=ks;iIfjb{2769~_i}YQ28xZARH8D)$w3=ZBS(%_ULVqG3yqlg!2s+RGe2u4} zaXf5gp;s-k0)hJA6nkVe@oF-`mQ>(Z`Vc|{9TjMRaLo&5)JG+QCGN@Zmwbf4-VeVA z>$lCgRG3DlQD8mj&6nRzKu|AVW9QoAWj$~r4^5=nDQ6QxC^o-PU4TdsB{YmEguE1a ze|#wt338o_ID%Rx+&4LEf~(2*EH%y9R$_;~foh!S_}v6QDtxW~kt7GI!ujfhU?7CR z0&4dAztt(M+B!Iw_NKyw0P~*|tSYizgC>p`n89_-DGP&@ABDg{$PVGQiZn8iGbq9k z3Xs|*)~~@Zg0FBex#<`N<#1yiE)UW{g3Trd0fdyYV3-*Eb;8VOgg8T8SrB#M{JOk_ zlru0#z_<uj0;bKDo5f27}XHy^=Y z`X?gBV}QfS$yJI69g;#w{(yUMkcDg_iLRnC?^HQRAtZX3Qh39&$szQ-1HZQ*gVYLK zh=#ra1vgu`!a8~2@kkzwNw|3+D)zNw8F5L?AktB=k$?e4xTdig|B5wg%cwc5WNHRM?3fZss>wV~F3^^d9Acv$7 z<{Z>e*H7%xF@$?1zE6XYk=UO>4tfPgA@U2Z*(1NiKxV z%6DV5kC^uivWC245CRV9L)IQ^kN7JH&Wj^AkT?#e4xz<|t-TnJeb#;HVi0(+V zfk4ucIO;+ia3MMcqEjF`39-?Ri`+5$_c-u%;IC63{-6xAwaYq~>&9OX{(3QTeV(`a z!f)SG7iK|Zmc%`IhwZ$6Pu>wbZvbzTu@F5V zv11@>9kb^`9Hl1L`(DBz8dPi?2ay9@McvO3Yhl3r$}9})6W%q1cXcvD;!Dx{)x7XW zK&m>0-|7s>F+K}#69g+7Cq0c7bQ&|}4DxHOF3gKFn2$mAgGl#?fqh@1N9vGk+znIw zLX?qfGsq{oLA0??h1~*SU534`!`uKy5S=5B0;9fLhjVAOJVT$gf*d{U5uEuT zCV>1Fy!Ulp=~;$LY)DUq-=FR3?V!^({F)bpk#{*-z$ zD?RnRX^EHCdtA2qUE7Z|-qNPgb)Mjqd9h+fIHMgw`KA%})ZH}VsOyY4>uwrlPu)$U zj=IjMv+kym_SD@p@~G>KJnL>6ZBN}zqmR17p7Kp2URqW%R+ltK3r0Y(FR<@E8Uoeq z3Cw{a#&-Fv0mS$e7zRb`lEhBR{&60)E&^)U&Mn?QMr8jOlKo@Y-ajrC?H?m0`^V@O z`^PAr$Nggrzb8Oodfz&g8wZ;rk{i#&E$kgIpj>gQ4MeERdK`Pu6@7)|LO}lLP;RGkOBwkXP{RLDZVjFY8zE9E6~2U8bP7xg4axCZwSNX8M6jEC{foC5QJ@B9+pIZ5n3 z5(xEFYDl+X7EL9F5Fge@B|e z*Ys-}V|qqkXJ0WZPT`$bt4lbmj_d2GMr+I($H_rL?%H@KH*}kI8K$gmz@YCH_S;QO zAQiM;vddmZ>PCuW&np(*9ll9$*@`d0Rs8L6_?8!nIzE`^X3}*?n_Nb z_!3U9uOs!coW4#9U z&s)|({gy^3?4P%CZv*yK!S<$Ld-JK-ZUEaYV0&BNdKxaQrwXt+O_RXAbGtWrQq@t%Z~`wn-VWWYLHvF+y`X`IHfw*Yl|V zh6g*-ZwmGWU|;ZHza!Z16k)#u?DP80Q?M^`HCS(0Z(uc~^kU|wZnqkQss1KdN}iOk zlsxIprSjfVnE&SFEFoD+Xx|Z?!Wv5zuQBWpDeMr7vilV65O{xYVf~F-*MzCGq`Nb3 zGp(0+ip5UWjrF<2dmT=A%eZUZ+IR;BqHg94F zuU9pQy#QssV#)>(;)zskQeqb9QT0F0JpmtyR(*%hTF>Xze|>wQ0M&(As;6VSOB}y{9{+ zwY_NVJ-fBNk88lNGGD_%80kJ~-R2jrcm|NEXU*wugciO7Nx2hmf&wf4Hot>Utvipp z3sQFz6y_N4;f$Un2K>&$;a3LZf5o~Z{aL{NxhVT(M)uPb=Fs{!!`FRiq5tUX{uEzkZ5$XPBA=tVhEphA+5L+TpUlb$KL6tn-K-t_E?K+}lDvp8KI z|1V#rkB9 z32}bX@bAv+O(-0^?4A->hwYaTII=LuRVN1QqanUw z@-NkKjs_XkGLhb&Y8?JZWPX2(g-2P~%b<^e%am08Ob#$OfneeDmwz`K{rvy9Xk;E| z$Q}?^nr|^iSDD$y0+*7VaJiFNSX(sAXV_CTT*5ogVD-5QuL`pkvX%_K%0M(|Mfiux1~shpZ|bz*ffQ7|@Ho`jD|7X7GT)TMPyn zpg%sfi=C$sWyKaO9qYQ+_>GB`d81i7e5f*vYlp;GiAG8+d=T(GHyv_QD z7;8c-UuRm*-9VcD{Vc9Q9rX>y{yl@Uoc;w)rZ6#wh+!`u{?8dJ;aZNf;4*^(Pq7nl z`D)lvOeW@_*PB$rdCXVM&dsf_tWy-;@C=1}D9pqX-X}yov8LFT=jDP5BT-nP!c81z zx5xxm!uUMP*t^V|WkBw3KB!<;NxAILE%~10ady0CxR>?zhU$D@^X16E`Xu)c%Ts*< z&86lpzF)AzC?6*E630cJR#&`V3wKk!f?^RG2YZ}o)o-Dl`Eo3H;=m8_cCziq7;s6c zI}8>Xq!`Q~$d^;rw^#id`j)Q{sp|wxc2X@Gzty58k&x^n!FRF>dy8_&zl5bW79?PO zK&=pAn8+o$#$;VVs+*axkajmYHv(A)^#(gDJXqPfunn5uhlbVrDvPJt<7#2ire)IR zq4x~6XA_5wI2d>!dox5qfrW?A)s|nwY^<^qy zWgHG{$P9;Jj*?-lZnHW$!D@~%i@`jD1qOE*EHX$T$OnZ5%`mQ(n7hnig~2L=H3oMX z++&~^WEiX?Fr4Bv4R64hDs$gs@PNTX23ZClWAF!DpPI28gA#^qHTiWI3bJau8o-p> z0J!gfmjxT7)?GO@Xvy%S4qptq(=Fv-x`J(L^2vboLykNgY%SYErX|>ysF zA-|5STjUrO1h9F7ED3+BDp(b&D+lWqCD7^hFm}nLjDR%@hNnMhp3gR*6^>RfTTmxh zf`Nm=wxo8N;FZ<{%bUTFji4>;kKZ5urr+=XCK<6IXtNr<^5(KVP@ zHf2Tis)EECDZ={H6tF(E3|Kv7U?sc%x9@j_Fv7v2<%kEtL!$r?N7$j4i zEOTf-BknT2!x&M#Wt%0s{tb1fk7-rxrQd zZ8PA(~ z&5fL&QDCY2Rx^}!oEYB_#?WQ-8*OM&?U8v2wI2DU-G;_ey$sJ1<|*9$cNsA<-$(G5 zZo$VIlKD{I0}0l@;bIo|d7m-&-tkxbd@1xQKzlvm_fnb5f4ib7R z*!F+`nF{&>=oz#QGQu(7@Oa(m1?m0GpF)5+XSS0Km}S+`^yuD{Nsc-x0MMph*$oThO%|t;mT= z4q3|njolWQ`Cxpobcc0NMrMzMdqGa$514&=2S`(gkcwsYQ#2tJ)BB zJEMBkg*;pW>7*UVb{FzU38a&DAUnW5If{4Ot~=znBh$7)R*u;Nj#BzN*Qq->15@qL zvpB4e0s1iJZ@a88|E&!Cl^#Zoj!Y*+K8#kxaiovQSae}5$nbZ@k>uXvy89_((WA-8 zKkmWAbeCZ21tv1~9ZY?|1QP30FvYQ}R9GDr$g;V99fwds!0HrY_%RTZJt-lEd(xXr zWrp{ssB^>_FQH8guhRXAYUon*Kzt@D?6Qz(SU|f&$7Mzy$9!r6yQ~Yl;Q-!BTn|9Q zqFadJC!k?*f+Gvk_espi7SXV{fmn}dSlqy_*$d%_0q7)9M`8fF2mEa(4#qHdcUZk3 zXrBV{eTNXVPeV-Nl;|*2fuMa__L0+8A8PFsg7zusFq{${29`S-7;#VeNqDvd4rs06 zY!CeXAY1RS25h~G`%h_8$EG$7La)NK=`6$|2HiFpO^36dHVvXpgKnD~y^BJd2BB9` zgEkF9$sTQj7{+~B!`PWYAU?;v7{88Y#d&=mS_~&JZ}}A*wob|zU`Mr1AqMRT{$hst zQkR5)j4`rxFrEd8(aH^F+Mv5|kzWOT#|DLz%$y;%!pYI~(7(8(ClV5@rDuP8)Y`RmHTP`Vx7(-}CDpd_HrK>}L|e6w zH=q3sL;`sfr=;wjVp=54_$HPqXiPPSY_pj4 z7&|~TLd0|ForXHM6*aT_ngN(n( znSYMWEXe4)lKHAY8564lCFHU1BQhap35n~)rSnDL{u6xo^NQ$~Naq(c3`DO$$W%gj z8{&AB(5gtyQw4(BC;8lMW%kwIWsUa<_Z$-YBt(g6y{#nn)dW$Axc+4!cp2iAjQ$cU zPObIkj32G<$9#^45n*BJuy~`sXqFzJ?$!%(k35i5<8m3i$-rJCRxsfWQ$T zVFM<9jLbSNpO9k-&(RgD7Gd@z z7gXRF6VBA)7#xO7_)Z*~MNa2YnOu6?Ko}_3VU3GD57&@rjNrgq2A8nW9B9cn+^yno z97pyj#;#0-B=3{Rg9RqE9mW|pW5z|VO#LFU8VG|M+6EZa1h3e>L=(#r(bXa8Kb8OO z)>~0XvQs*%qz6u2mP>s+n;0@Xbd`(~-*srC8B*^YwVjZjM_NjDA+{5Ftr+!c{MF!% zMVErOLjQ55KDB#e56aEFVQM!BAhozu)>^GCd=I;E$*RLh?zMslahbOEKog-pwVzTk z{Im8#YPjBDvi1>|k(ZI25^TC27*P{M6L~Tu+@BA=!8(uwDJ2)dAIK~ltQXeRFiR1k z0Rx026F^TuF41O2`9BEb+Q3md0J~Ve#6b`Oe2GKQCh#Q=Q-2|GggKah2kFjrl~=kl z>yWIG$5u6GJV_@4!8PLBXftuVtwTKi<9T-N)jyNEho$bHZ?0Q_lp05Xp-o_a&$9-7 za^01-a}}>b=EU>F!+hfEr{y_EQ8%atJy`uT#b9+fy~h#0exu#$l)2C$@naJ2TzyF} z92X3k0tP#Vem=Ek1S`0UZEW|Gb?qvr5J!q+UV! zuw8z^#(T-&-54>ysZrp`0B+2NyJ*9W+vOehUMf4PkR@Qo!ked?<^w`(?DXc^kUic! z1Bt8j6Fl#-c_$d1D;pCu03f}s*m&s!A4rt@Kv0Z=lxa$bsbe(nH$>9Mp?erFsddS& z_40Gns<_7)!)J|X8`@_aWk$f6cG zqw_Iy^Q8aGk3W{gG)K03&l+)gPP^wGzdlp8B~ z=!%zUbH%1vdogOSMv%7{sRM8efv>YbZpD{SknOs@wmt`Omo)qsz*8PXA?TZsUX$*s zwFU9}@!sO{vTZ2P$nKosQ@)RQ${G9Hj8!rA{fK3c58u7J{E%*Q;7)Oz7-k?Yke*Mf ziIoLh+*cRVcam1O`IkN9Li2_qW)`k)12ngmhA)`u*@x*(KnGk*x0m4>fDL&1!^EqE zF(_6F)aK&3wq~3IJqr#)w@L>>v+%?DaMMFOTzoM#3o~Q12{1ZK*7YXP!lNjJ7Xd~M7l-3+Ou|yg#b#XKsq8LayTMvB=iX(3M;ZGPV?V{%PcycM zAYW$eSD5$f41Si$|Hjx!mVSlF-(&KZ82c7uWh|&+u!F%*vEX-Du#d?yHrOQ=<;&rO zb0w{Q89*J)MwTb+Y}mFa5SHx%Oi!-8m#|pim>|a`fu#BQJu+QW@k|!!|g=7F@siiyS8m zwlGeq14Ha6lo%NVk~$oCe`i^@zbaUR_VOv@yJvTlUcf~bZ4$Dp|C^J;Z)8?!Xc1X^ z@zOeusej6)@?#7>$5~L3NLXugiG=zAW_^(p@2d>H#(=_r3RDiCiZGB3;VxtM89Zhn z%RyFxEP{W@Jem5x%GlQ!$Rd#0EVEu#!dA1sg|kW)kj(m{O!hI5Wh0BAg|k;g_hb=& zl}TCCzt31<3QDY(J^JhD!mpz+eHST+4b+C?;nwgQ;fC<7sL@IpHoHso>o?Zp>+cG*4QdK(d&fp;M* zIX}C;oZ%tr=<3>0+XktZJz{rHP{5$=mZ;3Qm@}JjN#qE}H=1vvOT+{mAI~J{b1yMH zH8qnB+upEy3_B%;m)j9GW_>sljp*(%u8LPtlnO+8Ho9Ec% zE8Wfuh{0Q!>73;GO0~Xfwt{9*+0y)G`^5e7VfYyWt;j|nXJM#n4J8F1;pu#}C!xLzVCs_yZbs?ki6WER91PkY;2+6S=AfXm3V}AV+2<6juZ+6U)4zR5B!Nzrf*IU~qx~gmM63nwf^@ zpxI?HBLo>-yM>Euiqu7js++E?<^9MaVj$hC zGQP=;l=hcc_A3nJ2==Rt$@%K*jD3v(&sXYq7>F>oh*fhFP!w@i0}M_v7-BHN;3Wn> z!5UW?dzXO-F^i<}`%H>eD4E`zVq0y)`S3Zo>x7Hv3!D$>UqBG7q$?cb42Qq4IfKui z#Oz!OxA8ye{4}!5bm$tz?c_fW32=hY67I$BJP`R>Vd?9=QskS=r9{%mdZ~q_Mu3*go literal 0 HcmV?d00001 diff --git a/csst_ifs_sim/csst_ifs_sim.py b/csst_ifs_sim/csst_ifs_sim.py index fe44fd5..d052f58 100644 --- a/csst_ifs_sim/csst_ifs_sim.py +++ b/csst_ifs_sim/csst_ifs_sim.py @@ -13,10 +13,21 @@ import os from scipy.interpolate import interp1d import astropy.coordinates as coord import ctypes -from csst_ifs_sim.support import IFSinstrumentModel -from csst_ifs_sim.support import cosmicrays -from csst_ifs_sim.support import logger as lg -from csst_ifs_sim.CTI import CTI + +# from csst_ifs_sim.support import IFSinstrumentModel +# from csst_ifs_sim.support import cosmicrays +# from csst_ifs_sim.support import logger as lg +# from csst_ifs_sim.CTI import CTI + + +sys.path.append('../') + +from CTI import CTI +from support import logger as lg +from support import cosmicrays +from support import IFSinstrumentModel + + # from optparse import OptionParser import configparser as ConfigParser import cmath @@ -355,6 +366,7 @@ class StrayLight(object): DESCRIPTION. """ + filterIndex = {'nuv': 0, 'u': 1, 'g': 2, 'r': 3, 'i': 4, 'z': 5, 'y': 6} sat = (ctypes.c_double*3)() sat[:] = self.sat @@ -392,6 +404,8 @@ class StrayLight(object): DESCRIPTION. """ + filterIndex = {'nuv': 0, 'u': 1, 'g': 2, 'r': 3, 'i': 4, 'z': 5, 'y': 6} + sat = (ctypes.c_double*3)() sat[:] = self.sat ob = (ctypes.c_double*3)() @@ -402,7 +416,7 @@ class StrayLight(object): self.slcdll.ComposeY(ob, py1, py2) earth_e1 = (ctypes.c_double*7)() - self.slcdll.EarthShine(self.jtime, sat, ob, py1, earth_e1) + self.slcdll.EarthShine(self.jtime, sat, ob, py1, earth_e1) # e[7]代表7个波段的照度 earth_e2 = (ctypes.c_double*7)() -- GitLab