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.166099112432E+02/ data cf( 2)/-.234610012185E-01/ data cf( 3)/-.903833157870E+00/ data cf( 4)/-.520290888529E+02/ data cf( 5)/0.360070444062E+03/ data cf( 6)/-.975763500496E+03/ data cf( 7)/0.102060053232E+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.53446401 kcal/mol ************************************************************************ implicit real*8 (a-h,o-z) dimension cf( 8) data cf( 1)/0.330295692403E+01/ data cf( 2)/-.979572789592E-01/ data cf( 3)/-.713361129504E+00/ data cf( 4)/0.253860529707E+02/ data cf( 5)/-.514434052377E+03/ data cf( 6)/0.230586596180E+04/ data cf( 7)/-.409176370741E+04/ data cf( 8)/0.227494669622E+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.328464606788E+03/ data cf( 2)/-.751930083968E-01/ data cf( 3)/0.389150997186E+01/ data cf( 4)/-.193735989756E+03/ data cf( 5)/0.160255155466E+04/ data cf( 6)/-.523667435431E+04/ data cf( 7)/0.742190036246E+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 **** points * rms = 0.61240154 kcal/mol * emax = 6.89876087 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)/-.1451173891367916E+00/ data i1( 1)/ 0/,i2( 1)/ 1/,i3( 1)/ 1/ data cf( 2)/-.2746708279836622E+01/ data i1( 2)/ 1/,i2( 2)/ 0/,i3( 2)/ 1/ data cf( 3)/-.1039483125132600E+01/ data i1( 3)/ 1/,i2( 3)/ 1/,i3( 3)/ 0/ data cf( 4)/-.5446907233365197E+02/ data i1( 4)/ 0/,i2( 4)/ 1/,i3( 4)/ 2/ data cf( 5)/0.1984600473186336E+02/ data i1( 5)/ 0/,i2( 5)/ 2/,i3( 5)/ 1/ data cf( 6)/0.7274514186503724E+02/ data i1( 6)/ 1/,i2( 6)/ 0/,i3( 6)/ 2/ data cf( 7)/0.1239164746177456E+02/ data i1( 7)/ 1/,i2( 7)/ 1/,i3( 7)/ 1/ data cf( 8)/-.2624649196008197E+02/ data i1( 8)/ 1/,i2( 8)/ 2/,i3( 8)/ 0/ data cf( 9)/0.5817878880256682E+02/ data i1( 9)/ 2/,i2( 9)/ 0/,i3( 9)/ 1/ data cf( 10)/0.1114838277506441E+03/ data i1( 10)/ 2/,i2( 10)/ 1/,i3( 10)/ 0/ data cf( 11)/0.1004207816239764E+04/ data i1( 11)/ 0/,i2( 11)/ 1/,i3( 11)/ 3/ data cf( 12)/0.6844360495886535E+02/ data i1( 12)/ 0/,i2( 12)/ 2/,i3( 12)/ 2/ data cf( 13)/0.9686542839185411E+02/ data i1( 13)/ 0/,i2( 13)/ 3/,i3( 13)/ 1/ data cf( 14)/-.1332286225313906E+04/ data i1( 14)/ 1/,i2( 14)/ 0/,i3( 14)/ 3/ data cf( 15)/0.7574526396812965E+03/ data i1( 15)/ 1/,i2( 15)/ 1/,i3( 15)/ 2/ data cf( 16)/-.8193678366360208E+03/ data i1( 16)/ 1/,i2( 16)/ 2/,i3( 16)/ 1/ data cf( 17)/0.3461957169900290E+03/ data i1( 17)/ 1/,i2( 17)/ 3/,i3( 17)/ 0/ data cf( 18)/-.6977075706158557E+03/ data i1( 18)/ 2/,i2( 18)/ 0/,i3( 18)/ 2/ data cf( 19)/0.2891986345044372E+03/ data i1( 19)/ 2/,i2( 19)/ 1/,i3( 19)/ 1/ data cf( 20)/-.5493968984591789E+03/ data i1( 20)/ 2/,i2( 20)/ 2/,i3( 20)/ 0/ data cf( 21)/-.1114102466250037E+04/ data i1( 21)/ 3/,i2( 21)/ 0/,i3( 21)/ 1/ data cf( 22)/-.2204827552603339E+04/ data i1( 22)/ 3/,i2( 22)/ 1/,i3( 22)/ 0/ data cf( 23)/-.1924679891093269E+04/ data i1( 23)/ 0/,i2( 23)/ 1/,i3( 23)/ 4/ data cf( 24)/-.1358895926855879E+05/ data i1( 24)/ 0/,i2( 24)/ 2/,i3( 24)/ 3/ data cf( 25)/0.6287552061803304E+04/ data i1( 25)/ 0/,i2( 25)/ 3/,i3( 25)/ 2/ data cf( 26)/-.4884612405531008E+04/ data i1( 26)/ 0/,i2( 26)/ 4/,i3( 26)/ 1/ data cf( 27)/0.1855200797994766E+05/ data i1( 27)/ 1/,i2( 27)/ 0/,i3( 27)/ 4/ data cf( 28)/-.3141944006864828E+05/ data i1( 28)/ 1/,i2( 28)/ 1/,i3( 28)/ 3/ data cf( 29)/0.2471397494584910E+05/ data i1( 29)/ 1/,i2( 29)/ 2/,i3( 29)/ 2/ data cf( 30)/-.2601214041815175E+04/ data i1( 30)/ 1/,i2( 30)/ 3/,i3( 30)/ 1/ data cf( 31)/-.9213001714837226E+03/ data i1( 31)/ 1/,i2( 31)/ 4/,i3( 31)/ 0/ data cf( 32)/-.9848047179725659E+03/ data i1( 32)/ 2/,i2( 32)/ 0/,i3( 32)/ 3/ data cf( 33)/0.1566255089262676E+05/ data i1( 33)/ 2/,i2( 33)/ 1/,i3( 33)/ 2/ data cf( 34)/-.9959570978049873E+04/ data i1( 34)/ 2/,i2( 34)/ 2/,i3( 34)/ 1/ data cf( 35)/0.1778230091563868E+04/ data i1( 35)/ 2/,i2( 35)/ 3/,i3( 35)/ 0/ data cf( 36)/0.5031559988934252E+04/ data i1( 36)/ 3/,i2( 36)/ 0/,i3( 36)/ 2/ data cf( 37)/-.1179161812726421E+05/ data i1( 37)/ 3/,i2( 37)/ 1/,i3( 37)/ 1/ data cf( 38)/0.1100896528909468E+05/ data i1( 38)/ 3/,i2( 38)/ 2/,i3( 38)/ 0/ data cf( 39)/0.1572119125685420E+05/ data i1( 39)/ 4/,i2( 39)/ 0/,i3( 39)/ 1/ data cf( 40)/0.2204795414668714E+05/ data i1( 40)/ 4/,i2( 40)/ 1/,i3( 40)/ 0/ data cf( 41)/-.3050737898469690E+05/ data i1( 41)/ 0/,i2( 41)/ 1/,i3( 41)/ 5/ data cf( 42)/0.1067279159830858E+06/ data i1( 42)/ 0/,i2( 42)/ 2/,i3( 42)/ 4/ data cf( 43)/-.8087212043749463E+04/ data i1( 43)/ 0/,i2( 43)/ 3/,i3( 43)/ 3/ data cf( 44)/-.1209963534370361E+05/ data i1( 44)/ 0/,i2( 44)/ 4/,i3( 44)/ 2/ data cf( 45)/0.3087438790465187E+05/ data i1( 45)/ 0/,i2( 45)/ 5/,i3( 45)/ 1/ data cf( 46)/-.1374653729373905E+06/ data i1( 46)/ 1/,i2( 46)/ 0/,i3( 46)/ 5/ data cf( 47)/0.3494728544446885E+06/ data i1( 47)/ 1/,i2( 47)/ 1/,i3( 47)/ 4/ data cf( 48)/-.3017832940234218E+06/ data i1( 48)/ 1/,i2( 48)/ 2/,i3( 48)/ 3/ data cf( 49)/0.4080391153719974E+05/ data i1( 49)/ 1/,i2( 49)/ 3/,i3( 49)/ 2/ data cf( 50)/0.7063517441542405E+04/ data i1( 50)/ 1/,i2( 50)/ 4/,i3( 50)/ 1/ data cf( 51)/0.6768090586093130E+04/ data i1( 51)/ 1/,i2( 51)/ 5/,i3( 51)/ 0/ data cf( 52)/0.2154349658201849E+05/ data i1( 52)/ 2/,i2( 52)/ 0/,i3( 52)/ 4/ data cf( 53)/-.1659031757132856E+06/ data i1( 53)/ 2/,i2( 53)/ 1/,i3( 53)/ 3/ data cf( 54)/0.2276706517250694E+06/ data i1( 54)/ 2/,i2( 54)/ 2/,i3( 54)/ 2/ data cf( 55)/-.7482144498531142E+05/ data i1( 55)/ 2/,i2( 55)/ 3/,i3( 55)/ 1/ data cf( 56)/-.4936026359455009E+04/ data i1( 56)/ 2/,i2( 56)/ 4/,i3( 56)/ 0/ data cf( 57)/0.2357607178184999E+05/ data i1( 57)/ 3/,i2( 57)/ 0/,i3( 57)/ 3/ data cf( 58)/-.8004667236860760E+05/ data i1( 58)/ 3/,i2( 58)/ 1/,i3( 58)/ 2/ data cf( 59)/-.4194157689618540E+05/ data i1( 59)/ 3/,i2( 59)/ 2/,i3( 59)/ 1/ data cf( 60)/0.1988951786328808E+05/ data i1( 60)/ 3/,i2( 60)/ 3/,i3( 60)/ 0/ data cf( 61)/-.5191131519488317E+05/ data i1( 61)/ 4/,i2( 61)/ 0/,i3( 61)/ 2/ data cf( 62)/0.1822246641337865E+06/ data i1( 62)/ 4/,i2( 62)/ 1/,i3( 62)/ 1/ data cf( 63)/-.9127844879886131E+05/ data i1( 63)/ 4/,i2( 63)/ 2/,i3( 63)/ 0/ data cf( 64)/-.1195291925623979E+06/ data i1( 64)/ 5/,i2( 64)/ 0/,i3( 64)/ 1/ data cf( 65)/-.1360502176637492E+06/ data i1( 65)/ 5/,i2( 65)/ 1/,i3( 65)/ 0/ data cf( 66)/0.1725772030259211E+06/ data i1( 66)/ 0/,i2( 66)/ 1/,i3( 66)/ 6/ data cf( 67)/-.2897809681706423E+06/ data i1( 67)/ 0/,i2( 67)/ 2/,i3( 67)/ 5/ data cf( 68)/-.1744137309422895E+06/ data i1( 68)/ 0/,i2( 68)/ 3/,i3( 68)/ 4/ data cf( 69)/0.1954196391706023E+06/ data i1( 69)/ 0/,i2( 69)/ 4/,i3( 69)/ 3/ data cf( 70)/-.1263895549558077E+06/ data i1( 70)/ 0/,i2( 70)/ 5/,i3( 70)/ 2/ data cf( 71)/-.6350788823841917E+05/ data i1( 71)/ 0/,i2( 71)/ 6/,i3( 71)/ 1/ data cf( 72)/0.5177055780609085E+06/ data i1( 72)/ 1/,i2( 72)/ 0/,i3( 72)/ 6/ data cf( 73)/-.1664469452119266E+07/ data i1( 73)/ 1/,i2( 73)/ 1/,i3( 73)/ 5/ data cf( 74)/0.1569519768083164E+07/ data i1( 74)/ 1/,i2( 74)/ 2/,i3( 74)/ 4/ data cf( 75)/-.3768473494529913E+06/ data i1( 75)/ 1/,i2( 75)/ 3/,i3( 75)/ 3/ data cf( 76)/-.8444256450507481E+04/ data i1( 76)/ 1/,i2( 76)/ 4/,i3( 76)/ 2/ data cf( 77)/0.8084783758145325E+05/ data i1( 77)/ 1/,i2( 77)/ 5/,i3( 77)/ 1/ data cf( 78)/-.7660512009708767E+05/ data i1( 78)/ 1/,i2( 78)/ 6/,i3( 78)/ 0/ data cf( 79)/0.3796380092747029E+05/ data i1( 79)/ 2/,i2( 79)/ 0/,i3( 79)/ 5/ data cf( 80)/0.5743748245516301E+06/ data i1( 80)/ 2/,i2( 80)/ 1/,i3( 80)/ 4/ data cf( 81)/-.1096579450613331E+07/ data i1( 81)/ 2/,i2( 81)/ 2/,i3( 81)/ 3/ data cf( 82)/0.7698468000319052E+06/ data i1( 82)/ 2/,i2( 82)/ 3/,i3( 82)/ 2/ data cf( 83)/-.2438051278436467E+06/ data i1( 83)/ 2/,i2( 83)/ 4/,i3( 83)/ 1/ data cf( 84)/0.9165280168627475E+05/ data i1( 84)/ 2/,i2( 84)/ 5/,i3( 84)/ 0/ data cf( 85)/-.3040476154502116E+06/ data i1( 85)/ 3/,i2( 85)/ 0/,i3( 85)/ 4/ data cf( 86)/0.8117880581033842E+06/ data i1( 86)/ 3/,i2( 86)/ 1/,i3( 86)/ 3/ data cf( 87)/-.8404844456801022E+06/ data i1( 87)/ 3/,i2( 87)/ 2/,i3( 87)/ 2/ data cf( 88)/0.5720314773407144E+05/ data i1( 88)/ 3/,i2( 88)/ 3/,i3( 88)/ 1/ data cf( 89)/0.2672576924130543E+05/ data i1( 89)/ 3/,i2( 89)/ 4/,i3( 89)/ 0/ data cf( 90)/-.6557366493845494E+04/ data i1( 90)/ 4/,i2( 90)/ 0/,i3( 90)/ 3/ data cf( 91)/-.1955888886004594E+05/ data i1( 91)/ 4/,i2( 91)/ 1/,i3( 91)/ 2/ data cf( 92)/0.8075756614475244E+06/ data i1( 92)/ 4/,i2( 92)/ 2/,i3( 92)/ 1/ data cf( 93)/-.2133565651866055E+06/ data i1( 93)/ 4/,i2( 93)/ 3/,i3( 93)/ 0/ data cf( 94)/0.3393675611236237E+06/ data i1( 94)/ 5/,i2( 94)/ 0/,i3( 94)/ 2/ data cf( 95)/-.1170692649637329E+07/ data i1( 95)/ 5/,i2( 95)/ 1/,i3( 95)/ 1/ data cf( 96)/0.3995760458370105E+06/ data i1( 96)/ 5/,i2( 96)/ 2/,i3( 96)/ 0/ data cf( 97)/0.4776863980662093E+06/ data i1( 97)/ 6/,i2( 97)/ 0/,i3( 97)/ 1/ data cf( 98)/0.5192287125777223E+06/ data i1( 98)/ 6/,i2( 98)/ 1/,i3( 98)/ 0/ data cf( 99)/-.3319067783514287E+06/ data i1( 99)/ 0/,i2( 99)/ 1/,i3( 99)/ 7/ data cf(100)/0.2563664771934466E+06/ data i1(100)/ 0/,i2(100)/ 2/,i3(100)/ 6/ data cf(101)/0.5935498088614707E+06/ data i1(101)/ 0/,i2(101)/ 3/,i3(101)/ 5/ data cf(102)/-.2349642680047289E+06/ data i1(102)/ 0/,i2(102)/ 4/,i3(102)/ 4/ data cf(103)/-.2818499844600155E+06/ data i1(103)/ 0/,i2(103)/ 5/,i3(103)/ 3/ data cf(104)/0.4654128240007957E+06/ data i1(104)/ 0/,i2(104)/ 6/,i3(104)/ 2/ data cf(105)/0.1500378920159702E+05/ data i1(105)/ 0/,i2(105)/ 7/,i3(105)/ 1/ data cf(106)/-.9573266412991872E+06/ data i1(106)/ 1/,i2(106)/ 0/,i3(106)/ 7/ data cf(107)/0.3660006027798350E+07/ data i1(107)/ 1/,i2(107)/ 1/,i3(107)/ 6/ data cf(108)/-.3579182406816544E+07/ data i1(108)/ 1/,i2(108)/ 2/,i3(108)/ 5/ data cf(109)/0.4908434108074186E+06/ data i1(109)/ 1/,i2(109)/ 3/,i3(109)/ 4/ data cf(110)/0.8841257482936078E+06/ data i1(110)/ 1/,i2(110)/ 4/,i3(110)/ 3/ data cf(111)/-.3945427971954461E+06/ data i1(111)/ 1/,i2(111)/ 5/,i3(111)/ 2/ data cf(112)/-.3641577385141246E+06/ data i1(112)/ 1/,i2(112)/ 6/,i3(112)/ 1/ data cf(113)/0.2787908297523257E+06/ data i1(113)/ 1/,i2(113)/ 7/,i3(113)/ 0/ data cf(114)/-.4265748672606891E+06/ data i1(114)/ 2/,i2(114)/ 0/,i3(114)/ 6/ data cf(115)/-.1090531816860121E+07/ data i1(115)/ 2/,i2(115)/ 1/,i3(115)/ 5/ data cf(116)/0.3268588544650054E+07/ data i1(116)/ 2/,i2(116)/ 2/,i3(116)/ 4/ data cf(117)/-.2210235187784999E+07/ data i1(117)/ 2/,i2(117)/ 3/,i3(117)/ 3/ data cf(118)/-.6539848534896662E+06/ data i1(118)/ 2/,i2(118)/ 4/,i3(118)/ 2/ data cf(119)/0.1348075644418475E+07/ data i1(119)/ 2/,i2(119)/ 5/,i3(119)/ 1/ data cf(120)/-.4278198059377750E+06/ data i1(120)/ 2/,i2(120)/ 6/,i3(120)/ 0/ data cf(121)/0.1079607076438976E+07/ data i1(121)/ 3/,i2(121)/ 0/,i3(121)/ 5/ data cf(122)/-.2118255095290509E+07/ data i1(122)/ 3/,i2(122)/ 1/,i3(122)/ 4/ data cf(123)/0.7094243438289305E+06/ data i1(123)/ 3/,i2(123)/ 2/,i3(123)/ 3/ data cf(124)/0.1335049302898237E+07/ data i1(124)/ 3/,i2(124)/ 3/,i3(124)/ 2/ data cf(125)/-.1290360026217662E+07/ data i1(125)/ 3/,i2(125)/ 4/,i3(125)/ 1/ data cf(126)/0.1121839752548341E+06/ data i1(126)/ 3/,i2(126)/ 5/,i3(126)/ 0/ data cf(127)/-.1485373714724166E+06/ data i1(127)/ 4/,i2(127)/ 0/,i3(127)/ 4/ data cf(128)/0.1239338069622415E+06/ data i1(128)/ 4/,i2(128)/ 1/,i3(128)/ 3/ data cf(129)/0.2828839682310345E+06/ data i1(129)/ 4/,i2(129)/ 2/,i3(129)/ 2/ data cf(130)/0.6671708723270607E+06/ data i1(130)/ 4/,i2(130)/ 3/,i3(130)/ 1/ data cf(131)/0.4367879112322491E+05/ data i1(131)/ 4/,i2(131)/ 4/,i3(131)/ 0/ data cf(132)/-.3914877266825115E+05/ data i1(132)/ 5/,i2(132)/ 0/,i3(132)/ 3/ data cf(133)/0.3146696214940414E+06/ data i1(133)/ 5/,i2(133)/ 1/,i3(133)/ 2/ data cf(134)/-.2665905341176484E+07/ data i1(134)/ 5/,i2(134)/ 2/,i3(134)/ 1/ data cf(135)/0.5161357529993323E+06/ data i1(135)/ 5/,i2(135)/ 3/,i3(135)/ 0/ data cf(136)/-.9414447415030508E+06/ data i1(136)/ 6/,i2(136)/ 0/,i3(136)/ 2/ data cf(137)/0.3294477282977226E+07/ data i1(137)/ 6/,i2(137)/ 1/,i3(137)/ 1/ data cf(138)/-.8712783503435561E+06/ data i1(138)/ 6/,i2(138)/ 2/,i3(138)/ 0/ data cf(139)/-.9755615008296034E+06/ data i1(139)/ 7/,i2(139)/ 0/,i3(139)/ 1/ data cf(140)/-.1091424280758889E+07/ data i1(140)/ 7/,i2(140)/ 1/,i3(140)/ 0/ data cf(141)/0.2176725063317629E+06/ data i1(141)/ 0/,i2(141)/ 1/,i3(141)/ 8/ data cf(142)/0.3304380343655173E+05/ data i1(142)/ 0/,i2(142)/ 2/,i3(142)/ 7/ data cf(143)/-.5969937478955020E+06/ data i1(143)/ 0/,i2(143)/ 3/,i3(143)/ 6/ data cf(144)/0.4718972126363198E+05/ data i1(144)/ 0/,i2(144)/ 4/,i3(144)/ 5/ data cf(145)/0.3875686871973745E+06/ data i1(145)/ 0/,i2(145)/ 5/,i3(145)/ 4/ data cf(146)/-.4977319580586850E+05/ data i1(146)/ 0/,i2(146)/ 6/,i3(146)/ 3/ data cf(147)/-.3984450229005935E+06/ data i1(147)/ 0/,i2(147)/ 7/,i3(147)/ 2/ data cf(148)/0.5095579546026827E+05/ data i1(148)/ 0/,i2(148)/ 8/,i3(148)/ 1/ data cf(149)/0.6966047860760380E+06/ data i1(149)/ 1/,i2(149)/ 0/,i3(149)/ 8/ data cf(150)/-.3084480680366338E+07/ data i1(150)/ 1/,i2(150)/ 1/,i3(150)/ 7/ data cf(151)/0.3089309329754936E+07/ data i1(151)/ 1/,i2(151)/ 2/,i3(151)/ 6/ data cf(152)/0.1386203780239782E+06/ data i1(152)/ 1/,i2(152)/ 3/,i3(152)/ 5/ data cf(153)/-.1622601471072440E+07/ data i1(153)/ 1/,i2(153)/ 4/,i3(153)/ 4/ data cf(154)/-.1311575733843850E+06/ data i1(154)/ 1/,i2(154)/ 5/,i3(154)/ 3/ data cf(155)/0.6370564960016095E+06/ data i1(155)/ 1/,i2(155)/ 6/,i3(155)/ 2/ data cf(156)/0.3414853437538099E+06/ data i1(156)/ 1/,i2(156)/ 7/,i3(156)/ 1/ data cf(157)/-.3093492631186083E+06/ data i1(157)/ 1/,i2(157)/ 8/,i3(157)/ 0/ data cf(158)/0.5829867369850066E+06/ data i1(158)/ 2/,i2(158)/ 0/,i3(158)/ 7/ data cf(159)/0.1248911043536329E+07/ data i1(159)/ 2/,i2(159)/ 1/,i3(159)/ 6/ data cf(160)/-.4824546541109537E+07/ data i1(160)/ 2/,i2(160)/ 2/,i3(160)/ 5/ data cf(161)/0.3363628897894940E+07/ data i1(161)/ 2/,i2(161)/ 3/,i3(161)/ 4/ data cf(162)/0.2269765240352807E+07/ data i1(162)/ 2/,i2(162)/ 4/,i3(162)/ 3/ data cf(163)/-.1179494012401426E+07/ data i1(163)/ 2/,i2(163)/ 5/,i3(163)/ 2/ data cf(164)/-.1169900528555702E+07/ data i1(164)/ 2/,i2(164)/ 6/,i3(164)/ 1/ data cf(165)/0.5193782584430051E+06/ data i1(165)/ 2/,i2(165)/ 7/,i3(165)/ 0/ data cf(166)/-.1283419503714328E+07/ data i1(166)/ 3/,i2(166)/ 0/,i3(166)/ 6/ data cf(167)/0.1639982315854683E+07/ data i1(167)/ 3/,i2(167)/ 1/,i3(167)/ 5/ data cf(168)/0.2893033951260524E+07/ data i1(168)/ 3/,i2(168)/ 2/,i3(168)/ 4/ data cf(169)/-.7634246503301354E+07/ data i1(169)/ 3/,i2(169)/ 3/,i3(169)/ 3/ data cf(170)/0.2772486812555089E+07/ data i1(170)/ 3/,i2(170)/ 4/,i3(170)/ 2/ data cf(171)/0.6050136381143964E+06/ data i1(171)/ 3/,i2(171)/ 5/,i3(171)/ 1/ data cf(172)/-.2050897067703426E+06/ data i1(172)/ 3/,i2(172)/ 6/,i3(172)/ 0/ data cf(173)/0.7049918510800622E+06/ data i1(173)/ 4/,i2(173)/ 0/,i3(173)/ 5/ data cf(174)/-.1730636384705670E+07/ data i1(174)/ 4/,i2(174)/ 1/,i3(174)/ 4/ data cf(175)/0.2301814624724537E+07/ data i1(175)/ 4/,i2(175)/ 2/,i3(175)/ 3/ data cf(176)/0.6164692759512253E+06/ data i1(176)/ 4/,i2(176)/ 3/,i3(176)/ 2/ data cf(177)/-.5630780999709301E+06/ data i1(177)/ 4/,i2(177)/ 4/,i3(177)/ 1/ data cf(178)/0.2064125706719125E+06/ data i1(178)/ 4/,i2(178)/ 5/,i3(178)/ 0/ data cf(179)/-.3307907615632712E+06/ data i1(179)/ 5/,i2(179)/ 0/,i3(179)/ 4/ data cf(180)/-.5743100408019097E+05/ data i1(180)/ 5/,i2(180)/ 1/,i3(180)/ 3/ data cf(181)/-.9931879174122943E+06/ data i1(181)/ 5/,i2(181)/ 2/,i3(181)/ 2/ data cf(182)/-.4592394628473503E+06/ data i1(182)/ 5/,i2(182)/ 3/,i3(182)/ 1/ data cf(183)/-.2942406774670298E+06/ data i1(183)/ 5/,i2(183)/ 4/,i3(183)/ 0/ data cf(184)/0.5589659158174123E+06/ data i1(184)/ 6/,i2(184)/ 0/,i3(184)/ 3/ data cf(185)/-.4402827985189819E+06/ data i1(185)/ 6/,i2(185)/ 1/,i3(185)/ 2/ data cf(186)/0.2739813274823615E+07/ data i1(186)/ 6/,i2(186)/ 2/,i3(186)/ 1/ data cf(187)/-.2402218343404057E+06/ data i1(187)/ 6/,i2(187)/ 3/,i3(187)/ 0/ data cf(188)/0.7371687603188831E+06/ data i1(188)/ 7/,i2(188)/ 0/,i3(188)/ 2/ data cf(189)/-.3277718057510877E+07/ data i1(189)/ 7/,i2(189)/ 1/,i3(189)/ 1/ data cf(190)/0.7050265708633160E+06/ data i1(190)/ 7/,i2(190)/ 2/,i3(190)/ 0/ data cf(191)/0.8273792313227295E+06/ data i1(191)/ 8/,i2(191)/ 0/,i3(191)/ 1/ data cf(192)/0.9422822573410947E+06/ data i1(192)/ 8/,i2(192)/ 1/,i3(192)/ 0/ vex1=0.1126748420000000E+01 vex2=0.8944411640000000E+00 vex3=0.9030453890000000E+00 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