From bf63ae852851da4aea222a3caec6e06f90744fa6 Mon Sep 17 00:00:00 2001 From: Emmanuel Giner Date: Tue, 26 Mar 2019 14:03:27 +0100 Subject: [PATCH] added the g09_zmat_to_qp_create.py --- Manuscript/G2-srDFT.pdf | Bin 476640 -> 476646 bytes Manuscript/G2-srDFT.tex | 11 ++--- Ref/Molecules/g09_zmat_to_qp_create.py | 57 +++++++++++++++++++++++++ 3 files changed, 63 insertions(+), 5 deletions(-) create mode 100755 Ref/Molecules/g09_zmat_to_qp_create.py diff --git a/Manuscript/G2-srDFT.pdf b/Manuscript/G2-srDFT.pdf index c0e444df6d8396860b7da67c465b3effb8364509..b794092dfeb63f16097a5426f56cad32e5dfc5c0 100644 GIT binary patch delta 7702 zcmai&RZJWVpsaCsUEHO(ySuwP#g+w%yDU=NrMNp3*W&I_U~#873q^~!aQ}Of^L#S- zl9^=QCz<_tij#PXCQDm5et~prTW0W_iK^$808amf;f(=EhjtJyexhdSc~O~BBX4<0 zc8Q$a)Q7K~VJ$8F&z`UJ4Qyj7WY}l zg{c?sa9}nh1t4e7FW#-PTP5_*?Hx26-`VkjgrBJfJXW`zHeqCE%t2tx#{dibteUQn z9PkU~nS=3vO;hJR5)&_MG?;O;hs?Z;1@xPJ{rE}L8l49odFD*oZq%GP^CO{2Dak~j z=>(h2(lF-=)72Nr(q$+!X6pee=b(=5pQhRtl=np0Ym!;2Boc9k^T(3d5T$U0<_}BIqn*MU;gxD+z zRGbF)2OFRci{`{?3E`r@ytsTvu}WGD3R*a2N*vdW*OE6!afJgaLg4hDjFxxCzNj8E zKV(Vw7zZ?WzYJ(voJ%wsZU9_L-7*qw6PzklJd&u0<6jO=qV@9(yEMuALl@LYH;Itd zSX;4iwMJFY-*8gdgj5_|KD|gGN;xck@uxPZa9Ckf3JR-JZ_n1o zpinT7`%Y-Jx)ORe3qVd0gGtfWqR~m4a>p6cl9NdV#^lmh_=3}IyQOrF=`Yb7@zVBx zRZd3v9>bbPJKgB>4C;&tW;FR1;+(NKC=|}O=VOj5TcDp@V92- z$=%Af8lXUgVPh-#>-CSDT_KGKUAncf6(+IYPJb*;iS@S;Q#J7!^aK_6O7Z7nQ=;!z zV$WMpX@-GnG{C*+XoC>qgR(ose?M;ITMfFVRcvxNvD?QG7k9%0zqC8dY%h+36)U*3M{{FLNx*hRWC|{@<%)x=EfoTT*)gVl zJy4n7hf-8V-eH^y?Im$vAH!v(PMu22dJCsj8*Z9eD$jD*iU_|E)Rv}++}Y_l2~}%B z69}+$bm0*eIhHAQ&Qm3uG?H5yO4VF6>8Am)Ejw88m8Lo;{A)^IU}Kpk`@THjJ5t8- zDQH)g!Gk9pJO@(tqHB|Z@rh-oif3+Q=IQ>pSWT(fr+RlcRoT}p2P|>S335UATxoPo zyOv9-Nk^`c^@@<#vtA-?NNOeFb=mg|ZEH4unU)46;XcgpYl>Jd5IUt@jfR0uO)8I9 zE==3PVrOYw@|K7(ChFo@KA!+Fxl3(Jpnxwh4ZJnUO5tzUKcTVIbz;?4ps-1|+rMN! z23v+Ue`I`v`2uL0k(C>}yio=t7uA*f=rfVT#<9hRhZ!0#)9bLp9GclFpE5UHejKsH z!{8^*R+r9&EQQo#>S3Hbse?X()8&*(9mO}T_PzDthV_CGL+gG9hiW04Ch>J(^jJut zfE8^BEv19D14&r<+2d@^Bjljm3>a6n_gri>G$J67+R8`d?b5{kkb3(bnQ2&!6XF0> zDV*+=dV~?ESt?1PuAEYNu()#5CM%Uns@CWjaQw6(bglf*7a zvJe0CNl2+7XejKvk(s6YDHr6DC;FL8xC^T6U9_(#98~(IG+;CJ4P0&3?rpD{gSVqr z%gwarI26kA$)Z-<^>+(eI^MB{kns!(tCddF&PwwN-Q**26OLlt=v(gwJA7NG7Zz8A zu7ST&E`r$um6o8fiQHYm*`)O7cp*%RbQ{#V2cb!=%LDrfiEVC6B*NMIH2NvOLs2%* z^^;zFJkI8s?s^DakX^6|hsM8Jf_YPML1mGD)=pSnTX@yV50{@u_P`O>Z&S125%cot zo(<@&kGWUb8n$j8$MtZuZ>y8R;|ybmbX5)>L$G8!r%luAaPA5(Xf%g0I9=Huigj$$ z!@!^pQ`gVBUmyuzhVqPkOsM;aX7|O`)62uo+U5T{*Uum{Zb4pZF6s!7Q#8m4_CMhM z1O7h{{sXZSvG@b|o&xLZ(1h^h3;~{w9tC z09u}3f?tBY1R+N^Tymerk2<$F8|HYeTyY97-1WIgTD{7_L!b?=GcGv&lMGwh(1=wW z0XOm&LlFXPV!~HP4bx26r{EIr+Q+t$NrS_V?@d!;mlh(Ie!u*5ZsB@Rpryx6Kgj5O ze<%@Z1Q}=Rq@$Yy)CQ_<*VplUOD9l2CBZULlbE9)vMIrskIT>Evs%$c{L+vQNIOQQ zWH4NZ+#vXo6g_&;S}}luu%L!fmjalaMn4_G23l5q)IXDGG%BQt(sm1^E-WyotQbaP zw3&Czgo`Ad4qVLZaK9LCSY3u=S<@oFVrZ#`{mDd45owAqBs?38{Zt!Mu&A+5h~Q=~ z-A{vV=Vm<*nAKw}24OspPzmv!mDI~)M^K93tdCX^7U7V{!tjCpd6RT%s1@jB2A4NO zRkmBj?gOxkrDvpICYodbO-)229t0OOG#+2zZP5Vj>cx>Ln`oa)sVNj#nCTsgtCv+B zY&KO(DXBvCw&TAW8 zFF^Afb{i5O&OJCL^egn74p>F`ije2#Yz8Z;%FmI=AxEPM!XY%^OI1Oa{vG`Bukf|R z4;?f59|bY=A#V@e>3bOb^%OSs{>9eR(atXTolZTuaxMPg?Hb*?ynvdqAPK1c&6Pun zX@82i^3Q0ABTNF3R5`$AA5DJ@3X3lwnP_Xl~d<&H5D=0`jIRG53 zwl2Q-`3e2g>t&#WMI&In$6CU0IAaaB>k06aL1%C*dGL}6z@fRIY?%37N;Y@{41&Ku zu0$dwR0EUC0j~O7bjbk$ss%_3|K7F&YBR_Nm9AZhrC1U+pdR~HVox0_cKlqqDn&>? z^_~PfM6loxC{^962bR}=gjMaK`B-JO`|M0Im}zq@+_Y@AxEo`^nGG5w2A?A{hJ;~M zJ%pC98$>{Y9QRoU0hrFb7-zekrEoO1QdG?asNg5(R;MuYKntGbh7Dh~Ca29{?53gI zfmJP8BktvWA09n#w_x!xEmw4f0LZY^{2^o+@B1$!98+TXm zsx5ASy$kB1)^^w*55zW)4vX>#eMTW}KxU9pi$MaZ$6?8lH=;cN19;e}=fLbyRq2aT z;A_jBpll7|Y+pXWdmeQ#SzIv@Qj208EFw{+1!hV+wH*q)VWi6^k>5IUgcu>Bifb;& z+!%EcRjqFPDE?0ZNFI2))=+Ne6@oIPfr@uu{eB0GyoV1CxZrv?#h%hZZE!u@n;H-A ztZVHcQUfwt{w1~0R6DW~v5Xz^)J9+%c+7b@(^WR5@t|0yZvjV-v0CddOL(=!qH|UM zjt^;wgj5vSO6Qi)MY-a_;zjs%ZY8PucW~Xn7U6J#Ho(+`5hTYYG77x`Sbd=cBhaMF z_r!&yce@_bTNUqT-y86xke8GpZIe|_O{$4y`=)wRqe!`gFBfV6;CY;+l_EF%q6`dp z2$_cYx2R)nr{qq5b2TIW;z8PEN~AmbjXO_nyI(+{#FIG$5Ws)O$__{j4n1@gh=)j> zG$d2ra~tmCJ`lHZ2;XfoZkncmojS549{h)S9StS3-f9l=xszY#-y^DZ5$vRyx(?0j zu8)X5@f$Xm-=&EnAK@kNlE(S_t0kkA<1#uqG6}5(!8k;CC_4t(5okR=Bp7mR;)8ek zDSHZJ?PuTUjas}Wb%EI-%c3xy!V#+!X^qFKJhN%U4QR285tY>-eK;c%TFP#$ztKhwCzf zSw?Nq+Ecl!c)+UcPz~k43h$D}^^6I#Z^AUmAb0*M#67DXeG?7FvLiS*QGw(apry=e zPDQ}%B6m;(L7sP}AsG?1NS|#9x5MT@h=a<`E`r6*t`}%4+`tCzcqpR0Exs+{Fn~nK zEew!7e~~Sk8)M05dL`O^GGme}zJo1}8BnyMO|1?RXC-;GCcpJ>IDTn$aQd|0V!APu zFvlm$PjV<@pkL9iSw|+&tK`ahp{UO|4l-_Kb0Xd&7(n-}E{dXc@y#}}f^D5< zqp_BYpZ&-T^Y;oo`$V_w_juEO_qcFZU#{-`Ha_z>?70AcQ8L|HroLK9atx*>q?r2@ zdFdHXswm>Tgi5qAG~?W?5@R=ntjrDf(`Fgou!>BZDM?58C!O~EjKq<<5kIw-90OdJH7ddtM!0>OoSrFPKM+EF04kPQ_#cSk<#} z0p{3b1Ewa4%GK`bET!{D%&_`*Ug-&w1=A2F{er+X$5H_F#4ql@_yBGR#o7QVP^g@? z>Tm-)(gG=)=v*~T4bEl@a|=SGn49y&oQgN0be3ljCsO1h?H&2KGI2aho-MjSr=7K3 zg&*9%KRCb%M3X6IFNPug-UXt(wp07kghM5C0Z3Y(6~XLLexs;g;&iHfC+uDlYini& z&l1>{sj%>mD3O&kMEZ|)*tiedKJ8Jc1Djk-y`L>p(Db6+>r`3#Fn26dxidu?Li~%_ z>^(_;!~M55-fTghk_F|aj7Pvkz%8;5R@so=q+Q2G0;a#J^7yYnr7x&mxIr-l zP*`F7_`tY{5)G`WcyW?i;s^ry&%nOWTpvSHZDHZc?dVdJx)E{1RL|UB#+p(I2CHMi z5~C0S7R3cp{SKfk((6(lseHfr4X@$_aBb99Pk*j)vihv$QXAyhxm z3R#TJ3ya*l$O>&V#R0p2g-9*UI8X))54?gpuzlOBC^@6jZ#RY>#%%R60--CPy&hOI zK5O5b9I6$4fBns7{pw0o_ z{CcOLa*mbb)4Z9*cGB;$-8$ekP;KQJ9lx$;yHkKT(0-4P29q(1=~IgV;?~ zBo_3(uPF266(o|&KE77=_IpV%E}Gdl+H7O2zE`$>8^*ZKN9zX@|Cjl&?HQ9>ja|#t z=_iwa8Z+iFKbN+qY}>#%^7(3wqEKf1n#H2d5#>aIjwhRMg4=-uT2$j;X@aFKV@te` zC6;UcrY?Ncy%BWM*L-u?wswPD@D-xt@F%gFOZl#&u>=nodYpVdLzu|B>!@^`DgJ^f$@?0_WWK+irdszjutK2sZ+_wfZ19hs<{jj=GK6p2xs;ui^S}EA>Up9Q zr~ksJYy)GyVk2=`{Dzx{MB$!VRMB$YgGN>($FgOaX8hvAOOV=h;6Ahf zZsq%FqehY?-1edRx`+ak@yw}8H)jv@JOG;0Eu3;P(hfe00?~C-0QHHNH?S`?B>s-LIIY;_D@YF5Fd$o7@K#bt!+a?bd?e#s(F7Q=KO@i zh%(hgff)Gjy@O)3W?|6N=jW$C8*c9Y$s7#tD*U&RGi!<$(|f0*nQs>-r<+d{QkvCM&Mffeofwy6rxP)g%=hE1ty|x*(-*54 zaj?Oh915)>E|-i^HeV0F|JuNSxZ=rFM!1AXbd#QLdopY}-S>Xm{&{&MO*QcTaCFKP z%r*6Z#N|lYarB1vPPuoEm0fA6YuJ}7a;8p%r;DuK+?>Q~xB1iTR8Gf}+H@^^N9OX;- zLuBv*4ed}zd>_OiwhszIR~dC_h2h8R>{l-hJsi5C#~1_sRzvL+t{qSVjf*s#(Hk?H z`d?2@-qI^7#h!P%2oq@@h|EXSPBV*^r0vWr*G{R>vyZlp&$jzX*sLNh^;bu+USyZV zYP<`XLy{Im$Xy?u#VL#G2H$8kq~lDe)Lr(IQ^z0Q;>JygYHFTCkj?H<%-i~?3U%U# zWZ@jKS7Kz;MUaRy2S(Iim!SS1`7^`@kd81N2}ct$^eZ~2kNJm%Ip99#JG@o_{=O=^b6A5Uixvz%VK1Vc<~ zWlsFe=^fOqw!$_wW0Oj@bxP~zh-_nNwQG&_%E%&y;XlOA5Ffod57w_6JE*1>^2Ut7 zF~AI#xV~a5#PDnp!ZLEV{T=+cnz!x4h~1;vZ}y9Gp7QF@tkvF?;g&GB+k}{NexS6h z&dUT2BkS7aabUR9naw=wjEj}UG0-CUn}%xiRzXr=?UbF0ffsj$I<2zLSV8uVV^0;O zS>@agRoiFcif49UPgSo;rHg#J>^$a01Ec=d!uWi|$zO@H%$eXo=oI*w994goO5psG zH}1q{=Ol07WTfg9J*^_+UmM|`dL!o)2uK#=$hyrK_t-nS#Z| zh+@1AnAhZ!z3qt!TZy=y9cTqD*;(>fTDi;Zpj7)n?48W!LCeBRN5ace5l4=M?gSYURZ{US4Ns!j+dKXfmdFTmxqg+M^1uT z?EkEoqj9R)eR1#tQu7N7Q2)>ONRct+p`lF>IWLhGh>Itugb-^H!x+35!-0ZF`wA?b z9vqyUoTW3Wh-11TnEuLQBYSv{g`}72s@`uBxc5=mGr#Nf@u}huXH|x7dF0hsk>(^w zC|b&pvIm(~RMBJF{gQ%KnC{VE6HI~drK6xb#wLnem%_72K>VpZ@??aQ862_SFve$7dZhRacZKScaPIDl> z+Kom78}0fIV}nG+hlULhEeocJy28Ut-in^yo@iZV+3Wop+5nINdY=3>QiAx^Y`s&a z6@89_ujgi7%cy%elW?_`D|8qJjZ8n=%K8xfp4~-Tlf9z)bhqk**(6~vt5Y7C8daDN Mje$X4QvvP&0rfS-C;$Ke delta 7686 zcmV+h9{J(siyh#L9k5?y0y8<2VPq(O)mlq)+cprs>sRnrX;LD<0*l9~)A-RzCP|yr z(?io76kB#QwI$b*^Z55&5PU!)C6g4LrYA!*2z>kP?uT7SJh@7E^2_Mi+34x1P%3fE zY0$~pWrAp!(vTQLW;H#!NZyTc^8W0_i05x@{|!r6T0$$oM5?T1(qfWIV8Y3Nu*TU$ zx$)1>UY<-+Kt6uXCMj}0KAx{u%Smd+=kw`$zFdU=&K9$~tFK}F7f0^*r)_k+jS$FC zOM;NCf1J-DD4IDVt)D|9lu0U>R2Zhss%h0yl4+rnXcp}$!tV1N(ckCO>!&x%>+orL z8UC(6&cfiRX?J6owrXR0Xx`p`C=qY}V_5XprbPv3$X)+MU+-8{S@H0q$efk&qSM7i zSoXAOS>TMfA)vM_5k1R-;{JNJ3LBg+Z*K3`+0%T>&-3-iF#0sTn=kM4i8r(L$1I{Q z@&_fx57X5=f3}*fS8OAyPtHdFj06pNA`(=I5fsEIlJlF(g= zO@iDq?vTdU$?u~#+wYTAy36ZjHbzbqj|0oa6h&V!pWgVbx?+X48f(AT*PP z+~R;~W%eNgidbo3A0mB_F5&=vCw`55%hRtI`RR)H+N5J5O&XsK$RB|2ZLD zpef}OdPR-T!ebddmTLTehMsESqSuk07-h`(q_(lT38oPXbrX_ODW7M&UP~qO3o}zphTo3lfxi_zKDb((GV4XNlNb~$+Ur_ib>cO z(kLXs-Bq3qZ;J!E3fEw(WPpW=Zzy5?JSS!?r}JFUn?lE-`{YzWoA6aYzDbp)j1>xl zpXyml38MI+St|tToPb2rp4O`3*f--Z1i4kJ0i{}jm1nro$t=FK7ay=RR=aGVyg)UW z4IT^NrVUSC9#4ROYN=Md3ktx(7UH9{1^qa;MYBe^N{hu(Umj?&1BEo>^wadZ6&e{a zi%gYhv>%sfcr&wfAOj4>EZJoEHnYDyJ85NyW}>u1W(*gsy3!70q7SVeLUx$9Ce>s_ z1~36wom31uR3jA~$--bMht_3qZ#>I@VPQd0v2kr%3t?)1nXSf@v6TT!OAClY0hMA>eKMcphZvzlulvw*}C$|s6h(MohAScbabt#;^VeTRzpHGI43pdf+MFYyT`eIeUscqR_Vwp^DR70=bD4TJdNZd z-t6hN;($?wXR<7KL_T_XxsRxr9X9E2V3U>9uIeNuQISPa(DN#bnnY>^V3?~c=^O=? z{Z6j32r$vejLO^W*Ui`3K*eBo-=y~6H%Wvfn%1EW23}PqAWYd;0Lb#(2MBer*C;$mYS3@Vpc`WnzD-x<37xG>!|&I{#CSA<&&pVM%2pP;t!X()^1;7?TJOgPsgqstnu~ss3;g< zN1k;I%whUoA0Q?ICnIpU4X#eSq)tr^Mm+7~xErMJ_1bD9ann0U*N2ndU3YESZ*J>< zE_Q_0STncEx9axahX9rsRM~DnW+w!k1bV3Z5-K;DXh_)ZX&P&GVB+yA;7OFc(Y(Q{U20&&bZ(?C0 zXpyl5dgOo*1&Pkl6ixyQiGI#TOc;IP%k(M~-Y9#Z%Chht!cryhD>aA$5J5vS3aAAK zb`@BP(vZuth??{QzUV=-k|x2chyd<1$XCRqguo7z0$HFI&lFI@KLDgKLTW&&6h=r# zxdgtAis)jjBHw?K5Gzyi4)iLoi4f*ffs!JoDp+*@(4b`$KIYI0G?7oOF)d(}nhF&u zaR*TnI31KOxPXq34j=+&z&cVPlm+ipH0|R-Na_IsNQrUUfMwix%mb*QJ)l2O{3Vwr zQC8^ywy4;kq@bo)4arj|)@WN$iL8nGz)|Tf)C86QYXE;j$viy9Qf;gUfz(4A(7~`X z1W-_aU=XTE3RwaTMTJcYNTpR}Jox_m1Mx%g78(xwJ1xZPH@_EJ$AJOw!45#+?frQC z?%+TFNnKASxAWrr?+bAvu^tdhgHC8o(KZd}HQ>n&ML>Mp%@K>H0mx(1fRTnpH$$ZA z26QYeW{ZE6S2M;_LUn_psc#KQPYqBlUQ8z!uc~?Rwh%9VI4Q(g^=V$T44r+xtMFoY zRUL>2V6M8I&j>n)ABfXxHo2c(RI}t5(xYe9dRdfY|MWtKkzWzt$l+ZD0)j72GqE%eSv{7 zW&k7I-p_$l7_*k;z(~KnuPoz0p_5W!U|>aCJ8S6&4g&WjIFw)NtjvJ^%{e{L%UZi` z<`M?t6g0eyJ9aM1frAO|)ERzYp!^reFj5W}+I$?8>Pr(x2V;aOYd4^K--!x@kbzfj zM-zYU*=E-QyY{?fs9m$~SYX$*`yOC}HGi$4I=5n|&g^BV;^_8vf~npXjP#_1fu7f; zrBlBH@#_A3p4v}Gw|^Xn5RWEjlLFHG*L-+Bt_FYd@aFErWHQIP)nNEs{BzltF0k5{|tZ?~U>qj!+uM*Zz}&Pr(!b>zC}yQZr|*J%H# z{_d)!Z?$32Ku9vE2^hw9I8>_Ea3ImOP&@-Vcp#?_OkfHIO^@QhJ4%dYn<+f;MJJ7N zPV>EF5rWj4c`yvtBQDsOuMXyBK@fjZS@er^oYZ}tZX70Zs7DN92`5Td=%@ns3Iw3U z!Nf?4Jqtk)fZ{|09LplHWZHt53bV0{tt@3C7ebqTkVRUcPw#E8um<624)kZiWFc>h zfMy^AYf^=j=2DnLXTskBQ!8#N0jOcr=W0mif|UD?)w_U3<5oBh+6 z`aa-Y>Q}!{&qL&zKF|8k{=1}39+Wv5(i3ZZ{UrAYp=tS&%1iEyW1YithGJfx)1m}v zgvnysB%EtLfU`n9&bA4}2KeKeIdNEr$k-$uxSh9#W?$52=?1 z&GafgMO{LVfEofmOR+~V)B*I!J?IHnASCcbqI=PmR%n!626P9(-F3P z;5;T@>9Q47N~ig}O8b9tWwEROw1@|d|Kz?I`;Z9wcLm?A=(Ex0eV}i_bIXlF1S;zQ zO7iyzOe7MeuVpF2F^nTNueH2>_G~Q$WtmBC8tf(+s>zrJ!!d@wPpmsOyrj__6~RCv z27|b1S=KJLzu5Bdgt5smHP1WGAbm2@EZ%YEkdaS$LS7@EDItH9E-7M{|A+o}!&~?w zZO*wb&ghNb`kp~;(>tnF?%^GSehSBKVH3#?luGS3jI$WD$Dv^dt1rB7`l)(duGIurj^&tiH0DDGs!)8q^()SA;anGYyPSJQviqvTj~YihY}E-C*}ksh@GKfnYD1EGB_B3i{OKpoW}ViW|fczBqkji zrI;~gx@ra$ykp5@X?2E)Y%VAaEhgqf*0EsV-a85mMl;QjUNW+|R3=&0uOQA~Ed$@K zZAcZ5RVjZ7Nb~WlmH}ZbG`huwkzl|ug8IY&QUfH4G9~O%^GlEf%SJM4q*MrfSrlo4 zLLNA`;h|G? z%aQu?P(IDWPdGeCOZI|3VV3R@lnvXc99eWL@NIv%HLF=FGfg(h)02L~0?eg!5Mx-N zu^$UK4>{l+3lU?XnJFna^fT>!>>6d@pCKk)#*zQiczd#$pgPdX%?w09eghK_>q#|<)UMK zNpgQx%!lUpp(H*g10i1(mO?JG%b~Q(M+j*iev0D3&zSHbsZe*wnONLYJO*{L3pUx4 zqIxk%Q6!4i3#e&#rAjK7t|5=1C>-l>%WdSGX7D2ylR=#}#%;cKMZ)!BHXECgPtQ*e z*?W*zelTUN53^H_aFxkKY2y?bA=jSuIWT{rDZ*?FKXX!K>62*A|5oaGA5(rvwahsq z_W*1T&KyUxO&c68`4*`xBZmoF&Ch1@OR9kR_)S%FU0ItSRgX1LndCw<_p+?buS`d4 zeyUFB_m#CO>m&5dkFur-()7Gf>k?n$=)nh*comr^m#H6d3&AC~N%}3nx*M7IQj&k0 zEYrt$x~bT>KYosC{^iT#J9}&v)0= ztr&^1xDmHvA||($n2MR0i}}ZixEG)Qo%uN#fd-jBZuR%PH>33O$FtKnk9sH_Z->ce zn2@jZaFRWobnM|2HsFN3p#zkP%RoI4Kh$7xNB*h!Sv(h~39?_rui`iH2GE@owid1r z;;IYTA1!Ei;;x#GCYL?rK8b(N;_ohcs+9Yn_xSA1Pj8;|(R;oLJ{U=cPq8-N2R`?Y z;Cqyydj{|tZ6jnEtPk&(0PE`5GGy1&s;YteOME1FedzQ)(Ea%2;qwPS61u0O^J+RD z-B#P}(>98SD7Qea?BZ4~P%Ae;=Hng;_oT+{Igoost9p z^Q4I5$!$j+(?X*9^L{w4(Z=z9B=Ui0@8{yf=r4$3azE`$;(uD@y#+Q4k?=d*xcV^9 zchj^$EvLKj{Y?BR{=A>etNenJn^k{d!rADPn2m?CYfAF{T4{>E#oyI*(if`TN2m{; zKRtbX)EDYjccN8j+=+iQ3#lS`>$_dsvj@4+9=d*4>#lVFysu{S5ros+yl-y${Q7

iZ9#7`BLMm17q7Uwa%AYSK_bg)4aO9oQ+6%4_vudpnCQE)!FGwfP1HZnMVEpu0_xOg+}B6zXczePvgk)cx~&(p5sdJ|=yo*!EN1t3wBpa@UjF>m zn`bYNPFn5`Gg0fA^kmX=Ux%A`9d~Ym+T@FrUrOZ;V{0y(SUs=iaM0Bl=jv`YVt#RR zFYc~K;%fV;j5crpXt3E@7dWezdriuyI)R!doe`qaXFeiJv^R_FE_+(%dlP6 zJJeBk&7=%HGx5DR<~C-rN_VWp_pIc#nD#hRdwIR1$7iQME^&4`xfyQY3{_}STjb{4 z=IrL@()JxdPpCm{Dg#Zo4D7|{$B4)%^G4H)`tUai%^OeMNr}tS|rYj;Bw7TL@$O#jN*dGW+wT^3tdK8B4uW|St~vwFE+Ji@izt? zt>dpiYPEW|69$+r<&K3i&0a;(;L?;X8JmB?#6BkQaBE7epOC~Zs;tEDqP|C(EZ!uJ zJ0)vMPrr*%mPVtB5x*6{nax*%#LkU9A&oMB11 z0u*wPwdLq9jaroA_q9lM^nX%Yzw4ie_G?Q@r?9$}3SyObfcPa&T6J6)V8*!zE`LSTDq>gG_;&r=sF{b|A}o{y%pIbCBR;X0oVn}))a zzl|>E*R%8!tJJrFGX17$Z7kSE)`ow*!)3bEGF`^>=u#cJ!ndGn)2`eyVP=Ie+aoM2 z2=nn<5Vm2@Ynd*$OqcEnl0`)em|LbR+W$SeBzGnkgt65S7UVKvq1>OaP5Vj9l!ax= zobD0kJBX=Qxn#}WO5!Utb+s;DUVlYLS?RO{C~n`i?0x+aR&+6hO(F?;#hs}AY1>i{F4JbV2QA-4$$tylwv<=P zw5@2e>$zds)ugx!ZTxxAjK6;my61M^fp{{yL_+2BFPmic8AbDJ-R4Js%`czoUyv-N zG|%tnh`gLNky+GulFH^jFDKh^qd5w)=9fuF$j0(r{llZ?7j`eE)nELva+irT(#2+E z^ZUG|kw5b%Ctj@?Y1@&iKeB5p&7zrydY+nzL$x%cY-Y5fH)EHHo-coHsZ6uG9z=&r z7+K%WZhEsfGp6ZHxyq?+(Y1YaP1zl}23^ylYuuKRU8wCEs7(vCQER73Bi9npLaM`N zrd5xCuA#KFC4g;Te7e7uZUNSPfK^w1dJTWN1z4AxN2<`FQoGe@RExH9Yp2O0*ASpu zfR)+;QxEN)maDFo*R%-rVEQjSshEtHQ6B~qmr^|ki7GWQIW;*nH!CnPConK4DGD!5 zZ)8MabY&nYL^?7sF*z_ZGBq(dH90glw*@{2>;zIqH8L?VK}13^Lo_i)I7USESAd_KaD1RN4OLI+85QTSd$+;m2B7`6zL=cag_gjRJ5JW`0UXQ&pi16(8b z6txrV0(Z$hLA8S2;2ya(R2$d>2=y4X7qkP8(j(M9u%D3X^^mXw93Z4fJs>;?4iQqR z?h|%`EoFa@T8kK(&|EK>Ui7kM7cfD2$2 zTm+W@g=emi{YkEnsUUNYOnuqqQd;H;vU8Rq)qneVE|^`_DsZ1Y08hpv^DI4hyTtav zSK#f6t-SV6aJF!^@U!r8ji-5uT~>bbuL?tMm6)vYAmX?2CMO? z-u$(ih)Ug7lTk~3{A@K9)wd?A>8Sqvw8Af1joS2k!a|=5eJ<=>> text = "Allowed Hello Hollow" + >>> tuple(findall('ll', text)) + (1, 10, 16) + """ + index = 0 - len(sub) + try: + while True: + index = string.index(sub, index + len(sub)) + yield index + except ValueError: + pass + +filepath = sys.argv[1] +basis = sys.argv[2] +icount = 0 +zmat = [] +with open(filepath, "r") as fp: + cipsi = [] + for line in fp: + if (icount == 0 ): + a=line.split(',') + charge = int(a[0]) + spin = int(a[1]) + else: + a=line.split() + if(len(a)>0): + n=a[0].count(',') + if(n > 6): + listfind=list(findall(',', a[0])) + char=a[0] + a[0]=char[0:listfind[n-1]] + + b=a[0].replace(',', ' ') + b=b.replace('=', ' ') + zmat.append(b) + else: + b=" " + zmat.append(b) + + icount += 1 + +c = filepath.split('.') +file_origin = c[0] +ezfio = file_origin+".ezfio" +zmat_file = file_origin+".zmt" +file_zmat = open(zmat_file,"w+") +for l in zmat: + file_zmat.write(l +'\n') + +print("qp_create_ezfio -m %s -c %s -b %s -o %s %s" % (spin,charge,basis,ezfio,zmat_file))