C PROGRAM makegjf_g09 C C make Gaussian09 dx.gjf from REST_g09 (x=0,001-950) C C 2021.4.8 C C read data 1. field point charge ex) 0.1 C 2. basis set (A15) aug-cc-pVDZ C 3. RHF, MP2, MP3, MP4,CISD,B3LYP C 'ReUse ' (A6,A20) RHFReUse C C 4. iop(6/20=x) (a1) 8 C 5. Massage(A7), # of additional AO (i4) Massage 10 C 6. charge multiplicity 0 1 C 7. %NProcShared=3 (*) =0 is default =2 4 C C ex) makegjf_g09int.sh C C # C # write fort.3 gjf C # write fort.4 sym.dat C # write fort.7 new_d0_scf.rest C # C ln -s ../init_scf.rest fort.1 C ln -s makegjf_g09.inp fort.2 C /home/naka/source/pol_src_06/makegjf_g09int.exe C rm -fr fort.* C # C IMPLICIT REAL*8 (A-H,O-Z) DIMENSION X(200),Y(200),Z(200),NT(200),IAT(200),sx(950),sy(950), & sz(950),pot(950),isym(950),jsym(950),ksym(2,950),p(200) CHARACTER*15 BAS CHARACTER*15 NP CHARACTER*6 RHF CHARACTER*20 INT CHARACTER*8 XRHF CHARACTER*2 AT(200) CHARACTER*6 D0 CHARACTER*8 DX(950) CHARACTER*7 MAS CHARACTER*30 MASDAT(100) CHARACTER*1 iop character*25 chachg data ToAng/0.529177d0/ DATA D0/'d0.gjf'/ DATA DX/'d001.gjf','d002.gjf','d003.gjf','d004.gjf','d005.gjf', & 'd006.gjf','d007.gjf','d008.gjf','d009.gjf','d010.gjf', & 'd011.gjf','d012.gjf','d013.gjf','d014.gjf','d015.gjf', & 'd016.gjf','d017.gjf','d018.gjf','d019.gjf','d020.gjf', & 'd021.gjf','d022.gjf','d023.gjf','d024.gjf','d025.gjf', & 'd026.gjf','d027.gjf','d028.gjf','d029.gjf','d030.gjf', & 'd031.gjf','d032.gjf','d033.gjf','d034.gjf','d035.gjf', & 'd036.gjf','d037.gjf','d038.gjf','d039.gjf','d040.gjf', & 'd041.gjf','d042.gjf','d043.gjf','d044.gjf','d045.gjf', & 'd046.gjf','d047.gjf','d048.gjf','d049.gjf','d050.gjf', & 'd051.gjf','d052.gjf','d053.gjf','d054.gjf','d055.gjf', & 'd056.gjf','d057.gjf','d058.gjf','d059.gjf','d060.gjf', & 'd061.gjf','d062.gjf','d063.gjf','d064.gjf','d065.gjf', & 'd066.gjf','d067.gjf','d068.gjf','d069.gjf','d070.gjf', & 'd071.gjf','d072.gjf','d073.gjf','d074.gjf','d075.gjf', & 'd076.gjf','d077.gjf','d078.gjf','d079.gjf','d080.gjf', & 'd081.gjf','d082.gjf','d083.gjf','d084.gjf','d085.gjf', & 'd086.gjf','d087.gjf','d088.gjf','d089.gjf','d090.gjf', & 'd091.gjf','d092.gjf','d093.gjf','d094.gjf','d095.gjf', & 'd096.gjf','d097.gjf','d098.gjf','d099.gjf','d100.gjf', & 'd101.gjf','d102.gjf','d103.gjf','d104.gjf','d105.gjf', & 'd106.gjf','d107.gjf','d108.gjf','d109.gjf','d110.gjf', & 'd111.gjf','d112.gjf','d113.gjf','d114.gjf','d115.gjf', & 'd116.gjf','d117.gjf','d118.gjf','d119.gjf','d120.gjf', & 'd121.gjf','d122.gjf','d123.gjf','d124.gjf','d125.gjf', & 'd126.gjf','d127.gjf','d128.gjf','d129.gjf','d130.gjf', & 'd131.gjf','d132.gjf','d133.gjf','d134.gjf','d135.gjf', & 'd136.gjf','d137.gjf','d138.gjf','d139.gjf','d140.gjf', & 'd141.gjf','d142.gjf','d143.gjf','d144.gjf','d145.gjf', & 'd146.gjf','d147.gjf','d148.gjf','d149.gjf','d150.gjf', & 'd151.gjf','d152.gjf','d153.gjf','d154.gjf','d155.gjf', & 'd156.gjf','d157.gjf','d158.gjf','d159.gjf','d160.gjf', & 'd161.gjf','d162.gjf','d163.gjf','d164.gjf','d165.gjf', & 'd166.gjf','d167.gjf','d168.gjf','d169.gjf','d170.gjf', & 'd171.gjf','d172.gjf','d173.gjf','d174.gjf','d175.gjf', & 'd176.gjf','d177.gjf','d178.gjf','d179.gjf','d180.gjf', & 'd181.gjf','d182.gjf','d183.gjf','d184.gjf','d185.gjf', & 'd186.gjf','d187.gjf','d188.gjf','d189.gjf','d190.gjf', & 'd191.gjf','d192.gjf','d193.gjf','d194.gjf','d195.gjf', & 'd196.gjf','d197.gjf','d198.gjf','d199.gjf','d200.gjf', & 'd201.gjf','d202.gjf','d203.gjf','d204.gjf','d205.gjf', & 'd206.gjf','d207.gjf','d208.gjf','d209.gjf','d210.gjf', & 'd211.gjf','d212.gjf','d213.gjf','d214.gjf','d215.gjf', & 'd216.gjf','d217.gjf','d218.gjf','d219.gjf','d220.gjf', & 'd221.gjf','d222.gjf','d223.gjf','d224.gjf','d225.gjf', & 'd226.gjf','d227.gjf','d228.gjf','d229.gjf','d230.gjf', & 'd231.gjf','d232.gjf','d233.gjf','d234.gjf','d235.gjf', & 'd236.gjf','d237.gjf','d238.gjf','d239.gjf','d240.gjf', & 'd241.gjf','d242.gjf','d243.gjf','d244.gjf','d245.gjf', & 'd246.gjf','d247.gjf','d248.gjf','d249.gjf','d250.gjf', & 'd251.gjf','d252.gjf','d253.gjf','d254.gjf','d255.gjf', & 'd256.gjf','d257.gjf','d258.gjf','d259.gjf','d260.gjf', & 'd261.gjf','d262.gjf','d263.gjf','d264.gjf','d265.gjf', & 'd266.gjf','d267.gjf','d268.gjf','d269.gjf','d270.gjf', & 'd271.gjf','d272.gjf','d273.gjf','d274.gjf','d275.gjf', & 'd276.gjf','d277.gjf','d278.gjf','d279.gjf','d280.gjf', & 'd281.gjf','d282.gjf','d283.gjf','d284.gjf','d285.gjf', & 'd286.gjf','d287.gjf','d288.gjf','d289.gjf','d290.gjf', & 'd291.gjf','d292.gjf','d293.gjf','d294.gjf','d295.gjf', & 'd296.gjf','d297.gjf','d298.gjf','d299.gjf','d300.gjf', & 'd301.gjf','d302.gjf','d303.gjf','d304.gjf','d305.gjf', & 'd306.gjf','d307.gjf','d308.gjf','d309.gjf','d310.gjf', & 'd311.gjf','d312.gjf','d313.gjf','d314.gjf','d315.gjf', & 'd316.gjf','d317.gjf','d318.gjf','d319.gjf','d320.gjf', & 'd321.gjf','d322.gjf','d323.gjf','d324.gjf','d325.gjf', & 'd326.gjf','d327.gjf','d328.gjf','d329.gjf','d330.gjf', & 'd331.gjf','d332.gjf','d333.gjf','d334.gjf','d335.gjf', & 'd336.gjf','d337.gjf','d338.gjf','d339.gjf','d340.gjf', & 'd341.gjf','d342.gjf','d343.gjf','d344.gjf','d345.gjf', & 'd346.gjf','d347.gjf','d348.gjf','d349.gjf','d350.gjf', & 'd351.gjf','d352.gjf','d353.gjf','d354.gjf','d355.gjf', & 'd356.gjf','d357.gjf','d358.gjf','d359.gjf','d360.gjf', & 'd361.gjf','d362.gjf','d363.gjf','d364.gjf','d365.gjf', & 'd366.gjf','d367.gjf','d368.gjf','d369.gjf','d370.gjf', & 'd371.gjf','d372.gjf','d373.gjf','d374.gjf','d375.gjf', & 'd376.gjf','d377.gjf','d378.gjf','d379.gjf','d380.gjf', & 'd381.gjf','d382.gjf','d383.gjf','d384.gjf','d385.gjf', & 'd386.gjf','d387.gjf','d388.gjf','d389.gjf','d390.gjf', & 'd391.gjf','d392.gjf','d393.gjf','d394.gjf','d395.gjf', & 'd396.gjf','d397.gjf','d398.gjf','d399.gjf','d400.gjf', & 'd401.gjf','d402.gjf','d403.gjf','d404.gjf','d405.gjf', & 'd406.gjf','d407.gjf','d408.gjf','d409.gjf','d410.gjf', & 'd411.gjf','d412.gjf','d413.gjf','d414.gjf','d415.gjf', & 'd416.gjf','d417.gjf','d418.gjf','d419.gjf','d420.gjf', & 'd421.gjf','d422.gjf','d423.gjf','d424.gjf','d425.gjf', & 'd426.gjf','d427.gjf','d428.gjf','d429.gjf','d430.gjf', & 'd431.gjf','d432.gjf','d433.gjf','d434.gjf','d435.gjf', & 'd436.gjf','d437.gjf','d438.gjf','d439.gjf','d440.gjf', & 'd441.gjf','d442.gjf','d443.gjf','d444.gjf','d445.gjf', & 'd446.gjf','d447.gjf','d448.gjf','d449.gjf','d450.gjf', & 'd451.gjf','d452.gjf','d453.gjf','d454.gjf','d455.gjf', & 'd456.gjf','d457.gjf','d458.gjf','d459.gjf','d460.gjf', & 'd461.gjf','d462.gjf','d463.gjf','d464.gjf','d465.gjf', & 'd466.gjf','d467.gjf','d468.gjf','d469.gjf','d470.gjf', & 'd471.gjf','d472.gjf','d473.gjf','d474.gjf','d475.gjf', & 'd476.gjf','d477.gjf','d478.gjf','d479.gjf','d480.gjf', & 'd481.gjf','d482.gjf','d483.gjf','d484.gjf','d485.gjf', & 'd486.gjf','d487.gjf','d488.gjf','d489.gjf','d490.gjf', & 'd491.gjf','d492.gjf','d493.gjf','d494.gjf','d495.gjf', & 'd496.gjf','d497.gjf','d498.gjf','d499.gjf','d500.gjf', & 'd501.gjf','d502.gjf','d503.gjf','d504.gjf','d505.gjf', & 'd506.gjf','d507.gjf','d508.gjf','d509.gjf','d510.gjf', & 'd511.gjf','d512.gjf','d513.gjf','d514.gjf','d515.gjf', & 'd516.gjf','d517.gjf','d518.gjf','d519.gjf','d520.gjf', & 'd521.gjf','d522.gjf','d523.gjf','d524.gjf','d525.gjf', & 'd526.gjf','d527.gjf','d528.gjf','d529.gjf','d530.gjf', & 'd531.gjf','d532.gjf','d533.gjf','d534.gjf','d535.gjf', & 'd536.gjf','d537.gjf','d538.gjf','d539.gjf','d540.gjf', & 'd541.gjf','d542.gjf','d543.gjf','d544.gjf','d545.gjf', & 'd546.gjf','d547.gjf','d548.gjf','d549.gjf','d550.gjf', & 'd551.gjf','d552.gjf','d553.gjf','d554.gjf','d555.gjf', & 'd556.gjf','d557.gjf','d558.gjf','d559.gjf','d560.gjf', & 'd561.gjf','d562.gjf','d563.gjf','d564.gjf','d565.gjf', & 'd566.gjf','d567.gjf','d568.gjf','d569.gjf','d570.gjf', & 'd571.gjf','d572.gjf','d573.gjf','d574.gjf','d575.gjf', & 'd576.gjf','d577.gjf','d578.gjf','d579.gjf','d580.gjf', & 'd581.gjf','d582.gjf','d583.gjf','d584.gjf','d585.gjf', & 'd586.gjf','d587.gjf','d588.gjf','d589.gjf','d590.gjf', & 'd591.gjf','d592.gjf','d593.gjf','d594.gjf','d595.gjf', & 'd596.gjf','d597.gjf','d598.gjf','d599.gjf','d600.gjf', & 'd601.gjf','d602.gjf','d603.gjf','d604.gjf','d605.gjf', & 'd606.gjf','d607.gjf','d608.gjf','d609.gjf','d610.gjf', & 'd611.gjf','d612.gjf','d613.gjf','d614.gjf','d615.gjf', & 'd616.gjf','d617.gjf','d618.gjf','d619.gjf','d620.gjf', & 'd621.gjf','d622.gjf','d623.gjf','d624.gjf','d625.gjf', & 'd626.gjf','d627.gjf','d628.gjf','d629.gjf','d630.gjf', & 'd631.gjf','d632.gjf','d633.gjf','d634.gjf','d635.gjf', & 'd636.gjf','d637.gjf','d638.gjf','d639.gjf','d640.gjf', & 'd641.gjf','d642.gjf','d643.gjf','d644.gjf','d645.gjf', & 'd646.gjf','d647.gjf','d648.gjf','d649.gjf','d650.gjf', & 'd651.gjf','d652.gjf','d653.gjf','d654.gjf','d655.gjf', & 'd656.gjf','d657.gjf','d658.gjf','d659.gjf','d660.gjf', & 'd661.gjf','d662.gjf','d663.gjf','d664.gjf','d665.gjf', & 'd666.gjf','d667.gjf','d668.gjf','d669.gjf','d670.gjf', & 'd671.gjf','d672.gjf','d673.gjf','d674.gjf','d675.gjf', & 'd676.gjf','d677.gjf','d678.gjf','d679.gjf','d680.gjf', & 'd681.gjf','d682.gjf','d683.gjf','d684.gjf','d685.gjf', & 'd686.gjf','d687.gjf','d688.gjf','d689.gjf','d690.gjf', & 'd691.gjf','d692.gjf','d693.gjf','d694.gjf','d695.gjf', & 'd696.gjf','d697.gjf','d698.gjf','d699.gjf','d700.gjf', & 'd701.gjf','d702.gjf','d703.gjf','d704.gjf','d705.gjf', & 'd706.gjf','d707.gjf','d708.gjf','d709.gjf','d710.gjf', & 'd711.gjf','d712.gjf','d713.gjf','d714.gjf','d715.gjf', & 'd716.gjf','d717.gjf','d718.gjf','d719.gjf','d720.gjf', & 'd721.gjf','d722.gjf','d723.gjf','d724.gjf','d725.gjf', & 'd726.gjf','d727.gjf','d728.gjf','d729.gjf','d730.gjf', & 'd731.gjf','d732.gjf','d733.gjf','d734.gjf','d735.gjf', & 'd736.gjf','d737.gjf','d738.gjf','d739.gjf','d740.gjf', & 'd741.gjf','d742.gjf','d743.gjf','d744.gjf','d745.gjf', & 'd746.gjf','d747.gjf','d748.gjf','d749.gjf','d750.gjf', & 'd751.gjf','d752.gjf','d753.gjf','d754.gjf','d755.gjf', & 'd756.gjf','d757.gjf','d758.gjf','d759.gjf','d760.gjf', & 'd761.gjf','d762.gjf','d763.gjf','d764.gjf','d765.gjf', & 'd766.gjf','d767.gjf','d768.gjf','d769.gjf','d770.gjf', & 'd771.gjf','d772.gjf','d773.gjf','d774.gjf','d775.gjf', & 'd776.gjf','d777.gjf','d778.gjf','d779.gjf','d780.gjf', & 'd781.gjf','d782.gjf','d783.gjf','d784.gjf','d785.gjf', & 'd786.gjf','d787.gjf','d788.gjf','d789.gjf','d790.gjf', & 'd791.gjf','d792.gjf','d793.gjf','d794.gjf','d795.gjf', & 'd796.gjf','d797.gjf','d798.gjf','d799.gjf','d800.gjf', & 'd801.gjf','d802.gjf','d803.gjf','d804.gjf','d805.gjf', & 'd806.gjf','d807.gjf','d808.gjf','d809.gjf','d810.gjf', & 'd811.gjf','d812.gjf','d813.gjf','d814.gjf','d815.gjf', & 'd816.gjf','d817.gjf','d818.gjf','d819.gjf','d820.gjf', & 'd821.gjf','d822.gjf','d823.gjf','d824.gjf','d825.gjf', & 'd826.gjf','d827.gjf','d828.gjf','d829.gjf','d830.gjf', & 'd831.gjf','d832.gjf','d833.gjf','d834.gjf','d835.gjf', & 'd836.gjf','d837.gjf','d838.gjf','d839.gjf','d840.gjf', & 'd841.gjf','d842.gjf','d843.gjf','d844.gjf','d845.gjf', & 'd846.gjf','d847.gjf','d848.gjf','d849.gjf','d850.gjf', & 'd851.gjf','d852.gjf','d853.gjf','d854.gjf','d855.gjf', & 'd856.gjf','d857.gjf','d858.gjf','d859.gjf','d860.gjf', & 'd861.gjf','d862.gjf','d863.gjf','d864.gjf','d865.gjf', & 'd866.gjf','d867.gjf','d868.gjf','d869.gjf','d870.gjf', & 'd871.gjf','d872.gjf','d873.gjf','d874.gjf','d875.gjf', & 'd876.gjf','d877.gjf','d878.gjf','d879.gjf','d880.gjf', & 'd881.gjf','d882.gjf','d883.gjf','d884.gjf','d885.gjf', & 'd886.gjf','d887.gjf','d888.gjf','d889.gjf','d890.gjf', & 'd891.gjf','d892.gjf','d893.gjf','d894.gjf','d895.gjf', & 'd896.gjf','d897.gjf','d898.gjf','d899.gjf','d900.gjf', & 'd901.gjf','d902.gjf','d903.gjf','d904.gjf','d905.gjf', & 'd906.gjf','d907.gjf','d908.gjf','d909.gjf','d910.gjf', & 'd911.gjf','d912.gjf','d913.gjf','d914.gjf','d915.gjf', & 'd916.gjf','d917.gjf','d918.gjf','d919.gjf','d920.gjf', & 'd921.gjf','d922.gjf','d923.gjf','d924.gjf','d925.gjf', & 'd926.gjf','d927.gjf','d928.gjf','d929.gjf','d930.gjf', & 'd931.gjf','d932.gjf','d933.gjf','d934.gjf','d935.gjf', & 'd936.gjf','d937.gjf','d938.gjf','d939.gjf','d940.gjf', & 'd941.gjf','d942.gjf','d943.gjf','d944.gjf','d945.gjf', & 'd946.gjf','d947.gjf','d948.gjf','d949.gjf','d950.gjf'/ C WRITE(6,*) '% start makeinp_g09 %' C 10 FORMAT(9X,I4) C 11 FORMAT(I2,F16.10,2(2X,F16.10)) C IU1=1 IU2=2 IU3=3 IU4=4 IU7=7 KK=0 KKK=1 OPEN(IU1,FORM='FORMATTED',ACCESS='SEQUENTIAL') OPEN(IU2,FORM='FORMATTED',ACCESS='SEQUENTIAL') OPEN(IU4,FILE='sym.dat',FORM='FORMATTED',ACCESS='SEQUENTIAL') OPEN(IU7,FILE='new_d0_scf.rest',FORM='FORMATTED', & ACCESS='SEQUENTIAL') READ(IU1,'(20x,i5,a25)') NAT,chachg C READ(IU1,10) NAT write(6,*) '% NAT',nat DO I=1,NAT READ(IU1,'(i2,4(f16.10,2x))') NT(I),X(I),Y(I),Z(I),p(i) C write(6,*) NT(I),X(I),Y(I),Z(I),p(i) ENDDO C READ(IU2,*) CHARG READ(IU2,'(A15)') BAS READ(IU2,'(A6,a20)') RHF,INT if(INT.eq.' ')then INT='(ultrafine,Acc2E=12)' endif read(IU2,'(A1)') iop READ(IU2,'(A7,i4)') MAS,NMAS read(IU2,*) icharge,imulti read(IU2,*) nproc C if (nproc.eq.1)then NP='%NProcShared=1 ' else if (nproc.eq.2)then NP='%NProcShared=2 ' else if (nproc.eq.3)then NP='%NProcShared=3 ' else if (nproc.eq.4)then NP='%NProcShared=4 ' else if (nproc.eq.5)then NP='%NProcShared=5 ' else if (nproc.eq.6)then NP='%NProcShared=6 ' else if (nproc.eq.7)then NP='%NProcShared=7 ' else if (nproc.eq.8)then NP='%NProcShared=8 ' else if (nproc.eq.9)then NP='%NProcShared=9 ' else if (nproc.eq.10)then NP='%NProcShared=10' else if (nproc.eq.11)then NP='%NProcShared=11' else if (nproc.eq.12)then NP='%NProcShared=12' else if (nproc.eq.13)then NP='%NProcShared=13' else if (nproc.eq.14)then NP='%NProcShared=14' else if (nproc.eq.15)then NP='%NProcShared=15' else if (nproc.eq.16)then NP='%NProcShared=16' else NP='%NProcShared=2 ' endif C if(MAS.eq.'Massage')then do i=1,nmas read(iu2,'(a30)') MASDAT(i) enddo endif if(RHF.eq.' RHF')then XRHF='HF ' else if (RHF.eq.' MP2')then XRHF='(HF,MP2)' else if (RHF.eq.' MP3')then XRHF='HF ' else if (RHF.eq.' MP4')then XRHF='HF ' else if (RHF.eq.' CISD')then XRHF='(HF,CI) ' else XRHF='HF ' endif C OPEN(IU3,FILE=D0,FORM='FORMATTED',ACCESS='SEQUENTIAL') WRITE(IU3,'(a10)') '$ RunGauss' write(iu3,'(a7)') '%chk=d0' if(INT.eq.'ReUse '.or.INT.eq.'REUSE '.or.INT.eq.'reuse ')then write(iu3,'(a7)') '%int=d0' endif write(iu3,'(a15)') NP write(iu3,'(a12)') '%mem=30000MB' if(INT.eq.'ReUse '.or.INT.eq.'REUSE '.or.INT.eq.'reuse ')then WRITE(IU3,'(a1,a6,a1,a15,a17,a1,a16,a8,1x,a7,a6)') & '#',RHF,'/',BAS,' pop=MK iop(6/20=',iop,') NoSym Density=', & XRHF,MAS,' polar' else WRITE(IU3,'(a1,a6,a1,a15,a17,a1,a16,a8,1x,a7,a26,a6)') & '#',RHF,'/',BAS,' pop=MK iop(6/20=',iop,') NoSym Density=', & XRHF,MAS,' Int=(ultrafine,Acc2E=12) ',' polar' endif WRITE(6,'(a1,a6,a1,a15,a17,a1,a16,a8,1x,a7,a26,a6)') & '#',RHF,'/',BAS,' pop=MK iop(6/20=',iop,') NoSym Density=', & XRHF,MAS,' Int=(ultrafine,Acc2E=12) ',' polar' WRITE(IU3,*) ' ' WRITE(IU3,*) 'Gaussion09 data 0' WRITE(IU3,*) ' ' WRITE(IU3,'(2i2)') icharge,imulti DO I=1,NAT C IF(NT(I).eq.1) then C AT(I)='H ' C else if(NT(I).eq.6) then C AT(I)='C ' C else if(NT(I).eq.7) then C AT(I)='N ' C else if(NT(I).eq.8) then C AT(I)='O ' C else if(NT(I).eq.15) then C AT(I)='P ' C else C AT(I)='X ' C endif WRITE(IU3,'(i2,3F16.10)') NT(I),X(I)*ToAng,Y(I)*ToAng,Z(I)*ToAng ENDDO WRITE(IU3,*) ' ' if(MAS.eq.'Massage')then do i=1,nmas write(iu3,'(a30)') masdat(i) enddo WRITE(IU3,*) ' ' endif CLOSE(UNIT=IU3) C READ(IU1,*) NSF write(6,*) ' # NSF',NSF do i=1,nsf READ(IU1,'(4f15.9)') SX(i),SY(i),SZ(i),pot(i) C write(6,*) i,SX(i),SY(i),SZ(i),pot(i) isym(i)=0 jsym(i)=0 enddo C C find same ESPotential C do i=1,nsf if(isym(i).eq.0)then do j=i,nsf if(pot(i).eq.pot(j))then isym(j)=i jsym(j)=i endif enddo endif enddo kk=0 do i=1,nsf k=0 do j=i,nsf if(jsym(j).ne.0)then if(isym(i).eq.jsym(j))then k=k+1 jsym(j)=0 endif endif enddo if(k.gt.0)then kk=kk+1 ksym(1,kk)=i ksym(2,kk)=k write(IU4,'(3i5)') kk,ksym(1,kk),ksym(2,kk) endif enddo write(6,*) 'isym' C write(6,'(20I3)') (isym(i),i=1,nsf) C write(6,*) 'jsym' C write(6,'(20I3)') (jsym(i),i=1,nsf) C write(iu7,'(20x,i5,a25)') nat,chachg do i=1,nat write(iu7,'(i2,4(F16.10,a2))') & NT(I),X(I),' 1',Y(I),' 1',Z(I),' 1',p(i),' 1' enddo write(iu7,'(20x,i5)') kk DO J=1,kk OPEN(IU3,FILE=DX(J),FORM='FORMATTED', & ACCESS='SEQUENTIAL') WRITE(IU3,'(a10)') '$ RunGauss' write(iu3,'(a9)') '%chk=temp' if(INT.eq.'ReUse '.or.INT.eq.'REUSE '.or.INT.eq.'reuse ')then write(iu3,'(a9)') '%int=temp' endif write(iu3,'(a15)') NP write(iu3,'(a12)') '%mem=30000MB' WRITE(IU3,'(a1,a6,a1,a15,a17,a1,a6,a20,a12)') & '#',RHF,'/',BAS,' pop=MK iop(6/20=',iop,') Int=',INT, & ' Guess=Read ' write(iu3,'(a25,a8)') & 'NoSym Massage Density=',XRHF WRITE(IU3,*) ' ' WRITE(IU3,*) 'Gaussion09 data',J WRITE(IU3,*) ' ' WRITE(IU3,'(2i2)') icharge,imulti DO I=1,NAT WRITE(IU3,'(i2,3F16.10)') & NT(I),X(I)*ToAng,Y(I)*ToAng,Z(I)*ToAng ENDDO WRITE(IU3,*) ' ' WRITE(IU3,'(i3,a4,f5.2,3f16.10)') & NAT+1,' Ch ',CHARG,SX(ksym(1,j))*ToAng,SY(ksym(1,j))*ToAng, & SZ(ksym(1,j))*ToAng write(iu7,'(4f15.9)') SX(ksym(1,j)),SY(ksym(1,j)), & SZ(ksym(1,j)),pot(ksym(1,j)) if(MAS.eq.'Massage')then do i=1,nmas write(iu3,'(a30)') masdat(i) enddo endif WRITE(IU3,*) ' ' CLOSE(UNIT=IU3) ENDDO C CLOSE(UNIT=IU1) CLOSE(UNIT=IU4) WRITE(6,*) '% stop makeinp_g09 %' STOP END