subroutine fit3d(r12,r13,r23,e,der) ************************************************************* * PES of 1 ^1A'' state of OH+F * correlating to the O(^1D) + HF asymptote * * S. Gomez-Carrasco, A. Aguado, M. Paniagua and O. Roncero * J. Photochemistry and Photobiology A, 190, (2007), 145 * * INPUT/OUTPUT * r12 --> OH * r13 --> OF * r23 --> HF * e : energy * der(3): derivatives * * All units arr in atomic units ************************************************************ implicit real * 8 (a-h,o-z) dimension der(3) call diat12(r12,e12,d12) call diat13(r13,e13,d13) call diat23(r23,e23,d23) call triabc(r12,r13,r23,e123,der) e=e12+e13+e23+e123 der(1)=d12+der(1) der(2)=d13+der(2) der(3)=d23+der(3) return end ************************************************************************ subroutine diat12(r,ener,der) ************************************************************************ * This subroutine computes the energies of a diatomic potential * fitted to 99 points * rms = 0.03663237 kcal/mol * emax = 0.15496720 kcal/mol ************************************************************************ implicit real*8 (a-h,o-z) dimension cf( 7) data cf( 1)/0.166099112348E+02/ data cf( 2)/-.234610020110E-01/ data cf( 3)/-.903833124354E+00/ data cf( 4)/-.520290893216E+02/ data cf( 5)/0.360070446947E+03/ data cf( 6)/-.975763508620E+03/ data cf( 7)/0.102060054092E+04/ e0= 0.0000000E+00 der=0.d0 vex1= 0.1038258E+01 vex2= 0.3230380E+01 aux = 1.d0/r bux = dexp(-vex2*r)*aux cux = dexp(-vex1*r) ener=e0+cf(1)*bux dux=1.d0 eux=r*cux do 1 i=2, 7 der=der+(i-1)*cf(i)*dux dux=dux*eux ener=ener+cf(i)*dux 1 continue der=der*(1.d0-vex1*r)*cux der=der-cf(1)*(vex2+aux)*bux return end ************************************************************************ subroutine diat23(r,ener,der) ************************************************************************ * This subroutine computes the energies of a diatomic potential * fitted to 55 points * rms = 0.16634385 kcal/mol * emax = 0.53446400 kcal/mol ************************************************************************ implicit real*8 (a-h,o-z) dimension cf( 8) data cf( 1)/0.330295692487E+01/ data cf( 2)/-.979572803700E-01/ data cf( 3)/-.713361014661E+00/ data cf( 4)/0.253860503519E+02/ data cf( 5)/-.514434028330E+03/ data cf( 6)/0.230586585626E+04/ data cf( 7)/-.409176348526E+04/ data cf( 8)/0.227494651541E+04/ e0= 0.0000000E+00 der=0.d0 vex1= 0.1000000E+01 vex2= 0.1000000E+01 aux = 1.d0/r bux = dexp(-vex2*r)*aux cux = dexp(-vex1*r) ener=e0+cf(1)*bux dux=1.d0 eux=r*cux do 1 i=2, 8 der=der+(i-1)*cf(i)*dux dux=dux*eux ener=ener+cf(i)*dux 1 continue der=der*(1.d0-vex1*r)*cux der=der-cf(1)*(vex2+aux)*bux return end ************************************************************************ subroutine diat13(r,ener,der) ************************************************************************ * This subroutine computes the energies of a diatomic potential * fitted to 54 points * rms = 0.14129135 kcal/mol * emax = 0.46773777 kcal/mol ************************************************************************ implicit real*8 (a-h,o-z) dimension cf( 7) data cf( 1)/0.328464606941E+03/ data cf( 2)/-.751930069879E-01/ data cf( 3)/0.389150989866E+01/ data cf( 4)/-.193735988491E+03/ data cf( 5)/0.160255154498E+04/ data cf( 6)/-.523667432040E+04/ data cf( 7)/0.742190031796E+04/ e0= 0.0000000E+00 der=0.d0 vex1= 0.1094583E+01 vex2= 0.4382704E+01 aux = 1.d0/r bux = dexp(-vex2*r)*aux cux = dexp(-vex1*r) ener=e0+cf(1)*bux dux=1.d0 eux=r*cux do 1 i=2, 7 der=der+(i-1)*cf(i)*dux dux=dux*eux ener=ener+cf(i)*dux 1 continue der=der*(1.d0-vex1*r)*cux der=der-cf(1)*(vex2+aux)*bux return end ************************************************************* subroutine triabc(r12,r13,r23,ener,der) ************************************************************* * This subroutine computes the energies of a 3D PES * for the ABC system class fitted to 9367 points * rms = 0.96407541 kcal/mol * emax = 6.29268373 kcal/mol ************************************************************* implicit real*8(a-h,o-z) dimension i1( 192),i2( 192),i3( 192),cf( 192) dimension f12(0: 8),f13(0: 8),f23(0: 8) dimension der(3) data cf( 1)/-.6268587080551566E+01/ data i1( 1)/ 0/,i2( 1)/ 1/,i3( 1)/ 1/ data cf( 2)/0.2893295878141408E+01/ data i1( 2)/ 1/,i2( 2)/ 0/,i3( 2)/ 1/ data cf( 3)/-.2539207287068814E+01/ data i1( 3)/ 1/,i2( 3)/ 1/,i3( 3)/ 0/ data cf( 4)/0.3957308256487705E+03/ data i1( 4)/ 0/,i2( 4)/ 1/,i3( 4)/ 2/ data cf( 5)/0.4329056010299598E+02/ data i1( 5)/ 0/,i2( 5)/ 2/,i3( 5)/ 1/ data cf( 6)/-.1447675682802629E+03/ data i1( 6)/ 1/,i2( 6)/ 0/,i3( 6)/ 2/ data cf( 7)/-.3965004690959770E+02/ data i1( 7)/ 1/,i2( 7)/ 1/,i3( 7)/ 1/ data cf( 8)/0.8379513725549847E+01/ data i1( 8)/ 1/,i2( 8)/ 2/,i3( 8)/ 0/ data cf( 9)/-.7549994486219265E+02/ data i1( 9)/ 2/,i2( 9)/ 0/,i3( 9)/ 1/ data cf( 10)/0.1764231683235486E+03/ data i1( 10)/ 2/,i2( 10)/ 1/,i3( 10)/ 0/ data cf( 11)/-.3906896146690058E+04/ data i1( 11)/ 0/,i2( 11)/ 1/,i3( 11)/ 3/ data cf( 12)/-.5947287276921391E+04/ data i1( 12)/ 0/,i2( 12)/ 2/,i3( 12)/ 2/ data cf( 13)/0.5892893844729035E+03/ data i1( 13)/ 0/,i2( 13)/ 3/,i3( 13)/ 1/ data cf( 14)/0.3998939421129372E+04/ data i1( 14)/ 1/,i2( 14)/ 0/,i3( 14)/ 3/ data cf( 15)/-.3250104860171731E+04/ data i1( 15)/ 1/,i2( 15)/ 1/,i3( 15)/ 2/ data cf( 16)/0.3147809553675454E+04/ data i1( 16)/ 1/,i2( 16)/ 2/,i3( 16)/ 1/ data cf( 17)/0.8463116419330050E+03/ data i1( 17)/ 1/,i2( 17)/ 3/,i3( 17)/ 0/ data cf( 18)/-.2103533981054228E+03/ data i1( 18)/ 2/,i2( 18)/ 0/,i3( 18)/ 2/ data cf( 19)/0.1372949729647830E+04/ data i1( 19)/ 2/,i2( 19)/ 1/,i3( 19)/ 1/ data cf( 20)/-.2015372231439722E+04/ data i1( 20)/ 2/,i2( 20)/ 2/,i3( 20)/ 0/ data cf( 21)/0.9271896083478023E+03/ data i1( 21)/ 3/,i2( 21)/ 0/,i3( 21)/ 1/ data cf( 22)/-.3095151524315159E+04/ data i1( 22)/ 3/,i2( 22)/ 1/,i3( 22)/ 0/ data cf( 23)/0.6862156017023400E+04/ data i1( 23)/ 0/,i2( 23)/ 1/,i3( 23)/ 4/ data cf( 24)/0.7090480142383749E+05/ data i1( 24)/ 0/,i2( 24)/ 2/,i3( 24)/ 3/ data cf( 25)/0.3053818135078462E+05/ data i1( 25)/ 0/,i2( 25)/ 3/,i3( 25)/ 2/ data cf( 26)/-.9685585755449309E+04/ data i1( 26)/ 0/,i2( 26)/ 4/,i3( 26)/ 1/ data cf( 27)/-.4541359666383811E+05/ data i1( 27)/ 1/,i2( 27)/ 0/,i3( 27)/ 4/ data cf( 28)/0.5665212865989852E+05/ data i1( 28)/ 1/,i2( 28)/ 1/,i3( 28)/ 3/ data cf( 29)/-.2040104048765717E+05/ data i1( 29)/ 1/,i2( 29)/ 2/,i3( 29)/ 2/ data cf( 30)/-.1531726690554016E+05/ data i1( 30)/ 1/,i2( 30)/ 3/,i3( 30)/ 1/ data cf( 31)/-.8721674163691947E+04/ data i1( 31)/ 1/,i2( 31)/ 4/,i3( 31)/ 0/ data cf( 32)/0.1144451078644147E+05/ data i1( 32)/ 2/,i2( 32)/ 0/,i3( 32)/ 3/ data cf( 33)/-.1070949724496771E+05/ data i1( 33)/ 2/,i2( 33)/ 1/,i3( 33)/ 2/ data cf( 34)/-.2772964755516135E+05/ data i1( 34)/ 2/,i2( 34)/ 2/,i3( 34)/ 1/ data cf( 35)/0.1825326194937988E+05/ data i1( 35)/ 2/,i2( 35)/ 3/,i3( 35)/ 0/ data cf( 36)/0.4437230998597901E+04/ data i1( 36)/ 3/,i2( 36)/ 0/,i3( 36)/ 2/ data cf( 37)/-.8129348038845644E+04/ data i1( 37)/ 3/,i2( 37)/ 1/,i3( 37)/ 1/ data cf( 38)/0.2064501589592981E+05/ data i1( 38)/ 3/,i2( 38)/ 2/,i3( 38)/ 0/ data cf( 39)/-.8141890868526445E+04/ data i1( 39)/ 4/,i2( 39)/ 0/,i3( 39)/ 1/ data cf( 40)/0.2465877373495054E+05/ data i1( 40)/ 4/,i2( 40)/ 1/,i3( 40)/ 0/ data cf( 41)/0.7020820932894455E+05/ data i1( 41)/ 0/,i2( 41)/ 1/,i3( 41)/ 5/ data cf( 42)/-.3483230869457474E+06/ data i1( 42)/ 0/,i2( 42)/ 2/,i3( 42)/ 4/ data cf( 43)/-.2560150940746429E+06/ data i1( 43)/ 0/,i2( 43)/ 3/,i3( 43)/ 3/ data cf( 44)/-.1104563376343637E+06/ data i1( 44)/ 0/,i2( 44)/ 4/,i3( 44)/ 2/ data cf( 45)/0.5571788535367280E+05/ data i1( 45)/ 0/,i2( 45)/ 5/,i3( 45)/ 1/ data cf( 46)/0.2435160232331912E+06/ data i1( 46)/ 1/,i2( 46)/ 0/,i3( 46)/ 5/ data cf( 47)/-.3459880381340098E+06/ data i1( 47)/ 1/,i2( 47)/ 1/,i3( 47)/ 4/ data cf( 48)/-.1444418532628224E+06/ data i1( 48)/ 1/,i2( 48)/ 2/,i3( 48)/ 3/ data cf( 49)/0.2545583077669026E+06/ data i1( 49)/ 1/,i2( 49)/ 3/,i3( 49)/ 2/ data cf( 50)/0.4245164234937626E+05/ data i1( 50)/ 1/,i2( 50)/ 4/,i3( 50)/ 1/ data cf( 51)/0.7704825548186528E+04/ data i1( 51)/ 1/,i2( 51)/ 5/,i3( 51)/ 0/ data cf( 52)/0.5151133355839894E+05/ data i1( 52)/ 2/,i2( 52)/ 0/,i3( 52)/ 4/ data cf( 53)/0.7176630685738183E+05/ data i1( 53)/ 2/,i2( 53)/ 1/,i3( 53)/ 3/ data cf( 54)/-.1072159585759103E+06/ data i1( 54)/ 2/,i2( 54)/ 2/,i3( 54)/ 2/ data cf( 55)/0.2282724117083552E+06/ data i1( 55)/ 2/,i2( 55)/ 3/,i3( 55)/ 1/ data cf( 56)/-.1077844806727200E+05/ data i1( 56)/ 2/,i2( 56)/ 4/,i3( 56)/ 0/ data cf( 57)/-.2858420311261140E+06/ data i1( 57)/ 3/,i2( 57)/ 0/,i3( 57)/ 3/ data cf( 58)/0.3992752211980041E+06/ data i1( 58)/ 3/,i2( 58)/ 1/,i3( 58)/ 2/ data cf( 59)/-.3229835273444669E+05/ data i1( 59)/ 3/,i2( 59)/ 2/,i3( 59)/ 1/ data cf( 60)/-.1648291469176453E+06/ data i1( 60)/ 3/,i2( 60)/ 3/,i3( 60)/ 0/ data cf( 61)/-.4629417896419612E+05/ data i1( 61)/ 4/,i2( 61)/ 0/,i3( 61)/ 2/ data cf( 62)/-.1949906333894955E+05/ data i1( 62)/ 4/,i2( 62)/ 1/,i3( 62)/ 1/ data cf( 63)/-.1439432580568045E+06/ data i1( 63)/ 4/,i2( 63)/ 2/,i3( 63)/ 0/ data cf( 64)/0.7484429702524940E+05/ data i1( 64)/ 5/,i2( 64)/ 0/,i3( 64)/ 1/ data cf( 65)/-.7719278656714057E+05/ data i1( 65)/ 5/,i2( 65)/ 1/,i3( 65)/ 0/ data cf( 66)/-.3744526655754551E+06/ data i1( 66)/ 0/,i2( 66)/ 1/,i3( 66)/ 6/ data cf( 67)/0.6940047143166066E+06/ data i1( 67)/ 0/,i2( 67)/ 2/,i3( 67)/ 5/ data cf( 68)/0.1125856661021373E+07/ data i1( 68)/ 0/,i2( 68)/ 3/,i3( 68)/ 4/ data cf( 69)/0.5044504216449235E+06/ data i1( 69)/ 0/,i2( 69)/ 4/,i3( 69)/ 3/ data cf( 70)/0.2244209651728667E+06/ data i1( 70)/ 0/,i2( 70)/ 5/,i3( 70)/ 2/ data cf( 71)/-.1217975830886483E+06/ data i1( 71)/ 0/,i2( 71)/ 6/,i3( 71)/ 1/ data cf( 72)/-.6804566677946799E+06/ data i1( 72)/ 1/,i2( 72)/ 0/,i3( 72)/ 6/ data cf( 73)/0.1136882243269914E+07/ data i1( 73)/ 1/,i2( 73)/ 1/,i3( 73)/ 5/ data cf( 74)/0.7049788104425615E+06/ data i1( 74)/ 1/,i2( 74)/ 2/,i3( 74)/ 4/ data cf( 75)/-.4409117017507916E+06/ data i1( 75)/ 1/,i2( 75)/ 3/,i3( 75)/ 3/ data cf( 76)/-.5280169489320560E+05/ data i1( 76)/ 1/,i2( 76)/ 4/,i3( 76)/ 2/ data cf( 77)/-.8336635347757976E+06/ data i1( 77)/ 1/,i2( 77)/ 5/,i3( 77)/ 1/ data cf( 78)/0.1984568740921859E+06/ data i1( 78)/ 1/,i2( 78)/ 6/,i3( 78)/ 0/ data cf( 79)/-.7497257743490884E+06/ data i1( 79)/ 2/,i2( 79)/ 0/,i3( 79)/ 5/ data cf( 80)/-.3954274625463749E+06/ data i1( 80)/ 2/,i2( 80)/ 1/,i3( 80)/ 4/ data cf( 81)/0.1735477063688864E+07/ data i1( 81)/ 2/,i2( 81)/ 2/,i3( 81)/ 3/ data cf( 82)/-.2266522918162597E+07/ data i1( 82)/ 2/,i2( 82)/ 3/,i3( 82)/ 2/ data cf( 83)/0.1092758282747234E+07/ data i1( 83)/ 2/,i2( 83)/ 4/,i3( 83)/ 1/ data cf( 84)/-.4261937781432071E+06/ data i1( 84)/ 2/,i2( 84)/ 5/,i3( 84)/ 0/ data cf( 85)/0.1733939018651450E+07/ data i1( 85)/ 3/,i2( 85)/ 0/,i3( 85)/ 4/ data cf( 86)/-.2409516331917309E+07/ data i1( 86)/ 3/,i2( 86)/ 1/,i3( 86)/ 3/ data cf( 87)/0.2094101577695866E+07/ data i1( 87)/ 3/,i2( 87)/ 2/,i3( 87)/ 2/ data cf( 88)/-.1928759228409048E+07/ data i1( 88)/ 3/,i2( 88)/ 3/,i3( 88)/ 1/ data cf( 89)/0.9662626178734894E+06/ data i1( 89)/ 3/,i2( 89)/ 4/,i3( 89)/ 0/ data cf( 90)/0.9456925231332388E+06/ data i1( 90)/ 4/,i2( 90)/ 0/,i3( 90)/ 3/ data cf( 91)/-.2403344032870927E+07/ data i1( 91)/ 4/,i2( 91)/ 1/,i3( 91)/ 2/ data cf( 92)/0.7405595729720789E+06/ data i1( 92)/ 4/,i2( 92)/ 2/,i3( 92)/ 1/ data cf( 93)/0.8928531609523881E+05/ data i1( 93)/ 4/,i2( 93)/ 3/,i3( 93)/ 0/ data cf( 94)/0.3080162027038944E+06/ data i1( 94)/ 5/,i2( 94)/ 0/,i3( 94)/ 2/ data cf( 95)/0.7007141618658493E+06/ data i1( 95)/ 5/,i2( 95)/ 1/,i3( 95)/ 1/ data cf( 96)/0.7937287896389972E+06/ data i1( 96)/ 5/,i2( 96)/ 2/,i3( 96)/ 0/ data cf( 97)/-.4672476467936527E+06/ data i1( 97)/ 6/,i2( 97)/ 0/,i3( 97)/ 1/ data cf( 98)/-.5805417271424522E+05/ data i1( 98)/ 6/,i2( 98)/ 1/,i3( 98)/ 0/ data cf( 99)/0.6489387115067258E+06/ data i1( 99)/ 0/,i2( 99)/ 1/,i3( 99)/ 7/ data cf(100)/-.1627785320094919E+06/ data i1(100)/ 0/,i2(100)/ 2/,i3(100)/ 6/ data cf(101)/-.2834282925094131E+07/ data i1(101)/ 0/,i2(101)/ 3/,i3(101)/ 5/ data cf(102)/-.2861275413984854E+06/ data i1(102)/ 0/,i2(102)/ 4/,i3(102)/ 4/ data cf(103)/-.1346663726194692E+07/ data i1(103)/ 0/,i2(103)/ 5/,i3(103)/ 3/ data cf(104)/0.3831936336563133E+05/ data i1(104)/ 0/,i2(104)/ 6/,i3(104)/ 2/ data cf(105)/0.3417595027165637E+04/ data i1(105)/ 0/,i2(105)/ 7/,i3(105)/ 1/ data cf(106)/0.9804458883866143E+06/ data i1(106)/ 1/,i2(106)/ 0/,i3(106)/ 7/ data cf(107)/-.1776523060030818E+07/ data i1(107)/ 1/,i2(107)/ 1/,i3(107)/ 6/ data cf(108)/-.1047445825964430E+07/ data i1(108)/ 1/,i2(108)/ 2/,i3(108)/ 5/ data cf(109)/-.1899993424208760E+07/ data i1(109)/ 1/,i2(109)/ 3/,i3(109)/ 4/ data cf(110)/0.5540983863194000E+07/ data i1(110)/ 1/,i2(110)/ 4/,i3(110)/ 3/ data cf(111)/-.5294993995812363E+07/ data i1(111)/ 1/,i2(111)/ 5/,i3(111)/ 2/ data cf(112)/0.4787233471979634E+07/ data i1(112)/ 1/,i2(112)/ 6/,i3(112)/ 1/ data cf(113)/-.7945274639432294E+06/ data i1(113)/ 1/,i2(113)/ 7/,i3(113)/ 0/ data cf(114)/0.1780942504451237E+07/ data i1(114)/ 2/,i2(114)/ 0/,i3(114)/ 6/ data cf(115)/0.2275388627692387E+07/ data i1(115)/ 2/,i2(115)/ 1/,i3(115)/ 5/ data cf(116)/-.4765279427427372E+07/ data i1(116)/ 2/,i2(116)/ 2/,i3(116)/ 4/ data cf(117)/-.4250071900857664E+07/ data i1(117)/ 2/,i2(117)/ 3/,i3(117)/ 3/ data cf(118)/0.1224637520047149E+08/ data i1(118)/ 2/,i2(118)/ 4/,i3(118)/ 2/ data cf(119)/-.9232261524821052E+07/ data i1(119)/ 2/,i2(119)/ 5/,i3(119)/ 1/ data cf(120)/0.1196359203235616E+07/ data i1(120)/ 2/,i2(120)/ 6/,i3(120)/ 0/ data cf(121)/-.3731288162860165E+07/ data i1(121)/ 3/,i2(121)/ 0/,i3(121)/ 5/ data cf(122)/0.5371502419543261E+07/ data i1(122)/ 3/,i2(122)/ 1/,i3(122)/ 4/ data cf(123)/0.1347970308997842E+07/ data i1(123)/ 3/,i2(123)/ 2/,i3(123)/ 3/ data cf(124)/-.8155675067383258E+07/ data i1(124)/ 3/,i2(124)/ 3/,i3(124)/ 2/ data cf(125)/0.1080335413411182E+08/ data i1(125)/ 3/,i2(125)/ 4/,i3(125)/ 1/ data cf(126)/-.1348761133674709E+07/ data i1(126)/ 3/,i2(126)/ 5/,i3(126)/ 0/ data cf(127)/-.3188148419329537E+07/ data i1(127)/ 4/,i2(127)/ 0/,i3(127)/ 4/ data cf(128)/0.4018222745713055E+07/ data i1(128)/ 4/,i2(128)/ 1/,i3(128)/ 3/ data cf(129)/-.7011194209757750E+05/ data i1(129)/ 4/,i2(129)/ 2/,i3(129)/ 2/ data cf(130)/-.4130333117906914E+07/ data i1(130)/ 4/,i2(130)/ 3/,i3(130)/ 1/ data cf(131)/-.1455728941197444E+07/ data i1(131)/ 4/,i2(131)/ 4/,i3(131)/ 0/ data cf(132)/-.2446856546956068E+07/ data i1(132)/ 5/,i2(132)/ 0/,i3(132)/ 3/ data cf(133)/0.5766183817478556E+07/ data i1(133)/ 5/,i2(133)/ 1/,i3(133)/ 2/ data cf(134)/0.1250350874281330E+07/ data i1(134)/ 5/,i2(134)/ 2/,i3(134)/ 1/ data cf(135)/0.1100068791346752E+06/ data i1(135)/ 5/,i2(135)/ 3/,i3(135)/ 0/ data cf(136)/-.2523017037819096E+06/ data i1(136)/ 6/,i2(136)/ 0/,i3(136)/ 2/ data cf(137)/-.4249693707316679E+07/ data i1(137)/ 6/,i2(137)/ 1/,i3(137)/ 1/ data cf(138)/-.1860583782114977E+07/ data i1(138)/ 6/,i2(138)/ 2/,i3(138)/ 0/ data cf(139)/0.1420723106272898E+07/ data i1(139)/ 7/,i2(139)/ 0/,i3(139)/ 1/ data cf(140)/0.7758988478826688E+06/ data i1(140)/ 7/,i2(140)/ 1/,i3(140)/ 0/ data cf(141)/-.3554283837907269E+06/ data i1(141)/ 0/,i2(141)/ 1/,i3(141)/ 8/ data cf(142)/-.6487552981978842E+06/ data i1(142)/ 0/,i2(142)/ 2/,i3(142)/ 7/ data cf(143)/0.2491522280362633E+07/ data i1(143)/ 0/,i2(143)/ 3/,i3(143)/ 6/ data cf(144)/0.6031835973819369E+06/ data i1(144)/ 0/,i2(144)/ 4/,i3(144)/ 5/ data cf(145)/-.1115377620875402E+06/ data i1(145)/ 0/,i2(145)/ 5/,i3(145)/ 4/ data cf(146)/0.1468004250586346E+07/ data i1(146)/ 0/,i2(146)/ 6/,i3(146)/ 3/ data cf(147)/-.4472801421157531E+06/ data i1(147)/ 0/,i2(147)/ 7/,i3(147)/ 2/ data cf(148)/0.2066868720061168E+06/ data i1(148)/ 0/,i2(148)/ 8/,i3(148)/ 1/ data cf(149)/-.5746272966868690E+06/ data i1(149)/ 1/,i2(149)/ 0/,i3(149)/ 8/ data cf(150)/0.2184228100282770E+06/ data i1(150)/ 1/,i2(150)/ 1/,i3(150)/ 7/ data cf(151)/0.5843324361347137E+07/ data i1(151)/ 1/,i2(151)/ 2/,i3(151)/ 6/ data cf(152)/-.1430420386347256E+08/ data i1(152)/ 1/,i2(152)/ 3/,i3(152)/ 5/ data cf(153)/0.2258847246261448E+08/ data i1(153)/ 1/,i2(153)/ 4/,i3(153)/ 4/ data cf(154)/-.2371994927017400E+08/ data i1(154)/ 1/,i2(154)/ 5/,i3(154)/ 3/ data cf(155)/0.1566583154472574E+08/ data i1(155)/ 1/,i2(155)/ 6/,i3(155)/ 2/ data cf(156)/-.8032173904579133E+07/ data i1(156)/ 1/,i2(156)/ 7/,i3(156)/ 1/ data cf(157)/0.8948961709987468E+06/ data i1(157)/ 1/,i2(157)/ 8/,i3(157)/ 0/ data cf(158)/-.5009245447333019E+06/ data i1(158)/ 2/,i2(158)/ 0/,i3(158)/ 7/ data cf(159)/-.8439805242504861E+07/ data i1(159)/ 2/,i2(159)/ 1/,i3(159)/ 6/ data cf(160)/0.2531893660514360E+08/ data i1(160)/ 2/,i2(160)/ 2/,i3(160)/ 5/ data cf(161)/-.2861383458688283E+08/ data i1(161)/ 2/,i2(161)/ 3/,i3(161)/ 4/ data cf(162)/0.2605584713716009E+08/ data i1(162)/ 2/,i2(162)/ 4/,i3(162)/ 3/ data cf(163)/-.2485705978672431E+08/ data i1(163)/ 2/,i2(163)/ 5/,i3(163)/ 2/ data cf(164)/0.1507865865974372E+08/ data i1(164)/ 2/,i2(164)/ 6/,i3(164)/ 1/ data cf(165)/-.4588341856661596E+06/ data i1(165)/ 2/,i2(165)/ 7/,i3(165)/ 0/ data cf(166)/0.3148312418399331E+07/ data i1(166)/ 3/,i2(166)/ 0/,i3(166)/ 6/ data cf(167)/-.1380815265315640E+08/ data i1(167)/ 3/,i2(167)/ 1/,i3(167)/ 5/ data cf(168)/0.1388263410042255E+08/ data i1(168)/ 3/,i2(168)/ 2/,i3(168)/ 4/ data cf(169)/-.1103457579168680E+08/ data i1(169)/ 3/,i2(169)/ 3/,i3(169)/ 3/ data cf(170)/0.1868214678011471E+08/ data i1(170)/ 3/,i2(170)/ 4/,i3(170)/ 2/ data cf(171)/-.1571250584348770E+08/ data i1(171)/ 3/,i2(171)/ 5/,i3(171)/ 1/ data cf(172)/-.1363690987410023E+07/ data i1(172)/ 3/,i2(172)/ 6/,i3(172)/ 0/ data cf(173)/0.6620780484777940E+07/ data i1(173)/ 4/,i2(173)/ 0/,i3(173)/ 5/ data cf(174)/-.4353977315870030E+07/ data i1(174)/ 4/,i2(174)/ 1/,i3(174)/ 4/ data cf(175)/0.1602896687881461E+07/ data i1(175)/ 4/,i2(175)/ 2/,i3(175)/ 3/ data cf(176)/-.1069815858666470E+08/ data i1(176)/ 4/,i2(176)/ 3/,i3(176)/ 2/ data cf(177)/0.4056929401768219E+07/ data i1(177)/ 4/,i2(177)/ 4/,i3(177)/ 1/ data cf(178)/0.5339833599125876E+07/ data i1(178)/ 4/,i2(178)/ 5/,i3(178)/ 0/ data cf(179)/-.3557603906883401E+06/ data i1(179)/ 5/,i2(179)/ 0/,i3(179)/ 4/ data cf(180)/-.6311925222558991E+07/ data i1(180)/ 5/,i2(180)/ 1/,i3(180)/ 3/ data cf(181)/0.1029311660562216E+08/ data i1(181)/ 5/,i2(181)/ 2/,i3(181)/ 2/ data cf(182)/0.1821731385786316E+07/ data i1(182)/ 5/,i2(182)/ 3/,i3(182)/ 1/ data cf(183)/-.3591212751388259E+07/ data i1(183)/ 5/,i2(183)/ 4/,i3(183)/ 0/ data cf(184)/0.5280972922289867E+07/ data i1(184)/ 6/,i2(184)/ 0/,i3(184)/ 3/ data cf(185)/-.8966372205266234E+07/ data i1(185)/ 6/,i2(185)/ 1/,i3(185)/ 2/ data cf(186)/-.3694292188135857E+07/ data i1(186)/ 6/,i2(186)/ 2/,i3(186)/ 1/ data cf(187)/0.2664209057240756E+07/ data i1(187)/ 6/,i2(187)/ 3/,i3(187)/ 0/ data cf(188)/-.8616904198637077E+06/ data i1(188)/ 7/,i2(188)/ 0/,i3(188)/ 2/ data cf(189)/0.7246137051030328E+07/ data i1(189)/ 7/,i2(189)/ 1/,i3(189)/ 1/ data cf(190)/0.7733352852316707E+06/ data i1(190)/ 7/,i2(190)/ 2/,i3(190)/ 0/ data cf(191)/-.1604165845600394E+07/ data i1(191)/ 8/,i2(191)/ 0/,i3(191)/ 1/ data cf(192)/-.1066428237112068E+07/ data i1(192)/ 8/,i2(192)/ 1/,i3(192)/ 0/ vex1=0.1199945730000000E+01 vex2=0.9500538860000000E+00 vex3=0.1089432710000000E+01 f12(0)=1.d0 f13(0)=1.d0 f23(0)=1.d0 bux12=r12*dexp(-vex1*r12) bux13=r13*dexp(-vex2*r13) bux23=r23*dexp(-vex3*r23) do 1 i=1, 8 f12(i)=f12(i-1)*bux12 f13(i)=f13(i-1)*bux13 f23(i)=f23(i-1)*bux23 1 continue ener = 0.d0 der12 = 0.d0 der13 = 0.d0 der23 = 0.d0 do 2 l=1,192 aux=f12(i1(l))*f13(i2(l))*f23(i3(l)) dux12=i1(l)*f12(i1(l)-1)*f13(i2(l))*f23(i3(l)) dux13=i2(l)*f12(i1(l))*f13(i2(l)-1)*f23(i3(l)) dux23=i3(l)*f12(i1(l))*f13(i2(l))*f23(i3(l)-1) ener=ener+cf(l)*aux der12=der12+cf(l)*dux12 der13=der13+cf(l)*dux13 der23=der23+cf(l)*dux23 2 continue der(1)=der12*(1.d0-vex1*r12)*dexp(-vex1*r12) der(2)=der13*(1.d0-vex2*r13)*dexp(-vex2*r13) der(3)=der23*(1.d0-vex3*r23)*dexp(-vex3*r23) return end