c
c Aguado/Paniagua H3 surface
c
c From: "Dr. A. Aguado" <ALFREDO@vm1.sdi.uam.es>
c Subject: H3 PES
c Status: R
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C POTENTIAL ENERGY SURFACE FOR THE GROUND STATE OF H3
C See: A. Aguado and M. Paniagua. JCP 96, 1265 (1992)
C
C This main is only a test for the PES.
C
C R1,R2,ANG = INTERNUCLEAR DISTANCES FOR HH, HH and ANGLE HHH,
C These are converted to
C R1,R2,R3 = INTERNUCLEAR DISTANCES FOR HH,HH,HH RESPECTIVELY
C
c CALL H3aguado7 (R1,R2,R3,E7, D7 )
c CALL H3aguado11(R1,R2,R3,E11,D11)
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C H3 SURFACE EVALUATION ROUTINES from Aguado
C ==============================
C
C
SUBROUTINE H3aguado11(X,Y,Z,ENER,DER)
IMPLICIT REAL * 8 (A-H,O-Z)
DIMENSION DER(3)
CALL H2LIU3(X,E1,D1)
CALL H2LIU3(Y,E2,D2)
CALL H2LIU3(Z,E3,D3)
CALL V311(X,Y,Z,ENER,DER)
ENER = ENER + E1+E2+E3
DER(1) = DER(1)+D1
DER(2) = DER(2)+D2
DER(3) = DER(3)+D3
RETURN
END
************************************************************************
* *
* SUBROUTINE AJUSTA *
* *
* ESTA SUBRUTINA CALCULA LAS ENERGIAS DE UNA SUPERFICIE BIDIMENSIONAL *
* CON UN RMS DE 0.07303625 KCAL/MOL *
* CON UN EMAX DE 0.24606775 KCAL/MOL
* *
************************************************************************
SUBROUTINE V311(X,Y,Z,ENER,DER)
IMPLICIT REAL * 8 (A-H,O-Z)
DIMENSION I(330), J(330), K(330), CF(330)
DIMENSION XXX(-1:10), YYY(-1:10), ZZZ(-1:10)
DIMENSION DER(3)
DATA CF( 1)/0.14677087656214488780899781394333D+02/
DATA I( 1)/ 1/,J( 1)/ 1/,K( 1)/ 0/
DATA CF( 2)/0.14677087656214488780899781394333D+02/
DATA I( 2)/ 1/,J( 2)/ 0/,K( 2)/ 1/
DATA CF( 3)/0.14677087656214488780899781394333D+02/
DATA I( 3)/ 0/,J( 3)/ 1/,K( 3)/ 1/
DATA CF( 4)/-.48326064166157666607775390716017D+03/
DATA I( 4)/ 1/,J( 4)/ 1/,K( 4)/ 1/
DATA CF( 5)/-.16244877244465470633845548163692D+03/
DATA I( 5)/ 2/,J( 5)/ 1/,K( 5)/ 0/
DATA CF( 6)/-.16244877244465470633845548163692D+03/
DATA I( 6)/ 2/,J( 6)/ 0/,K( 6)/ 1/
DATA CF( 7)/-.16244877244465470633845548163692D+03/
DATA I( 7)/ 1/,J( 7)/ 2/,K( 7)/ 0/
DATA CF( 8)/-.16244877244465470633845548163692D+03/
DATA I( 8)/ 1/,J( 8)/ 0/,K( 8)/ 2/
DATA CF( 9)/-.16244877244465470633845548163692D+03/
DATA I( 9)/ 0/,J( 9)/ 2/,K( 9)/ 1/
DATA CF( 10)/-.16244877244465470633845548163692D+03/
DATA I( 10)/ 0/,J( 10)/ 1/,K( 10)/ 2/
DATA CF( 11)/0.63437717785836300593518383365217D+04/
DATA I( 11)/ 2/,J( 11)/ 1/,K( 11)/ 1/
DATA CF( 12)/0.63437717785836300593518383365217D+04/
DATA I( 12)/ 1/,J( 12)/ 2/,K( 12)/ 1/
DATA CF( 13)/0.63437717785836300593518383365217D+04/
DATA I( 13)/ 1/,J( 13)/ 1/,K( 13)/ 2/
DATA CF( 14)/0.73221190661019997088698509021166D+04/
DATA I( 14)/ 2/,J( 14)/ 2/,K( 14)/ 0/
DATA CF( 15)/0.73221190661019997088698509021166D+04/
DATA I( 15)/ 2/,J( 15)/ 0/,K( 15)/ 2/
DATA CF( 16)/0.73221190661019997088698509021166D+04/
DATA I( 16)/ 0/,J( 16)/ 2/,K( 16)/ 2/
DATA CF( 17)/-.10064602936430782171849387178899D+04/
DATA I( 17)/ 3/,J( 17)/ 1/,K( 17)/ 0/
DATA CF( 18)/-.10064602936430782171849387178899D+04/
DATA I( 18)/ 3/,J( 18)/ 0/,K( 18)/ 1/
DATA CF( 19)/-.10064602936430782171849387178899D+04/
DATA I( 19)/ 1/,J( 19)/ 3/,K( 19)/ 0/
DATA CF( 20)/-.10064602936430782171849387178899D+04/
DATA I( 20)/ 1/,J( 20)/ 0/,K( 20)/ 3/
DATA CF( 21)/-.10064602936430782171849387178899D+04/
DATA I( 21)/ 0/,J( 21)/ 3/,K( 21)/ 1/
DATA CF( 22)/-.10064602936430782171849387178899D+04/
DATA I( 22)/ 0/,J( 22)/ 1/,K( 22)/ 3/
DATA CF( 23)/-.10765866433345379136925737824945D+06/
DATA I( 23)/ 2/,J( 23)/ 2/,K( 23)/ 1/
DATA CF( 24)/-.10765866433345379136925737824945D+06/
DATA I( 24)/ 2/,J( 24)/ 1/,K( 24)/ 2/
DATA CF( 25)/-.10765866433345379136925737824945D+06/
DATA I( 25)/ 1/,J( 25)/ 2/,K( 25)/ 2/
DATA CF( 26)/-.97680403659657906145924550804058D+04/
DATA I( 26)/ 3/,J( 26)/ 1/,K( 26)/ 1/
DATA CF( 27)/-.97680403659657906145924550804058D+04/
DATA I( 27)/ 1/,J( 27)/ 3/,K( 27)/ 1/
DATA CF( 28)/-.97680403659657906145924550804058D+04/
DATA I( 28)/ 1/,J( 28)/ 1/,K( 28)/ 3/
DATA CF( 29)/-.55610159332873447244746232721176D+05/
DATA I( 29)/ 3/,J( 29)/ 2/,K( 29)/ 0/
DATA CF( 30)/-.55610159332873447244746232721176D+05/
DATA I( 30)/ 3/,J( 30)/ 0/,K( 30)/ 2/
DATA CF( 31)/-.55610159332873447244746232721176D+05/
DATA I( 31)/ 2/,J( 31)/ 3/,K( 31)/ 0/
DATA CF( 32)/-.55610159332873447244746232721176D+05/
DATA I( 32)/ 2/,J( 32)/ 0/,K( 32)/ 3/
DATA CF( 33)/-.55610159332873447244746232721176D+05/
DATA I( 33)/ 0/,J( 33)/ 3/,K( 33)/ 2/
DATA CF( 34)/-.55610159332873447244746232721176D+05/
DATA I( 34)/ 0/,J( 34)/ 2/,K( 34)/ 3/
DATA CF( 35)/0.28044719015568109673381955463300D+05/
DATA I( 35)/ 4/,J( 35)/ 1/,K( 35)/ 0/
DATA CF( 36)/0.28044719015568109673381955463300D+05/
DATA I( 36)/ 4/,J( 36)/ 0/,K( 36)/ 1/
DATA CF( 37)/0.28044719015568109673381955463300D+05/
DATA I( 37)/ 1/,J( 37)/ 4/,K( 37)/ 0/
DATA CF( 38)/0.28044719015568109673381955463300D+05/
DATA I( 38)/ 1/,J( 38)/ 0/,K( 38)/ 4/
DATA CF( 39)/0.28044719015568109673381955463300D+05/
DATA I( 39)/ 0/,J( 39)/ 4/,K( 39)/ 1/
DATA CF( 40)/0.28044719015568109673381955463300D+05/
DATA I( 40)/ 0/,J( 40)/ 1/,K( 40)/ 4/
DATA CF( 41)/0.11083828565178803149309834176145D+07/
DATA I( 41)/ 2/,J( 41)/ 2/,K( 41)/ 2/
DATA CF( 42)/0.59297165758948566629941986071239D+06/
DATA I( 42)/ 3/,J( 42)/ 2/,K( 42)/ 1/
DATA CF( 43)/0.59297165758948566629941986071239D+06/
DATA I( 43)/ 3/,J( 43)/ 1/,K( 43)/ 2/
DATA CF( 44)/0.59297165758948566629941986071239D+06/
DATA I( 44)/ 2/,J( 44)/ 3/,K( 44)/ 1/
DATA CF( 45)/0.59297165758948566629941986071239D+06/
DATA I( 45)/ 2/,J( 45)/ 1/,K( 45)/ 3/
DATA CF( 46)/0.59297165758948566629941986071239D+06/
DATA I( 46)/ 1/,J( 46)/ 3/,K( 46)/ 2/
DATA CF( 47)/0.59297165758948566629941986071239D+06/
DATA I( 47)/ 1/,J( 47)/ 2/,K( 47)/ 3/
DATA CF( 48)/0.65828190802810608011809540589227D+06/
DATA I( 48)/ 3/,J( 48)/ 3/,K( 48)/ 0/
DATA CF( 49)/0.65828190802810608011809540589227D+06/
DATA I( 49)/ 3/,J( 49)/ 0/,K( 49)/ 3/
DATA CF( 50)/0.65828190802810608011809540589227D+06/
DATA I( 50)/ 0/,J( 50)/ 3/,K( 50)/ 3/
DATA CF( 51)/-.22545143708228407284186245202466D+06/
DATA I( 51)/ 4/,J( 51)/ 1/,K( 51)/ 1/
DATA CF( 52)/-.22545143708228407284186245202466D+06/
DATA I( 52)/ 1/,J( 52)/ 4/,K( 52)/ 1/
DATA CF( 53)/-.22545143708228407284186245202466D+06/
DATA I( 53)/ 1/,J( 53)/ 1/,K( 53)/ 4/
DATA CF( 54)/0.12460954469096832559595394013825D+06/
DATA I( 54)/ 4/,J( 54)/ 2/,K( 54)/ 0/
DATA CF( 55)/0.12460954469096832559595394013825D+06/
DATA I( 55)/ 4/,J( 55)/ 0/,K( 55)/ 2/
DATA CF( 56)/0.12460954469096832559595394013825D+06/
DATA I( 56)/ 2/,J( 56)/ 4/,K( 56)/ 0/
DATA CF( 57)/0.12460954469096832559595394013825D+06/
DATA I( 57)/ 2/,J( 57)/ 0/,K( 57)/ 4/
DATA CF( 58)/0.12460954469096832559595394013825D+06/
DATA I( 58)/ 0/,J( 58)/ 4/,K( 58)/ 2/
DATA CF( 59)/0.12460954469096832559595394013825D+06/
DATA I( 59)/ 0/,J( 59)/ 2/,K( 59)/ 4/
DATA CF( 60)/-.27311873046773076105553203577757D+06/
DATA I( 60)/ 5/,J( 60)/ 1/,K( 60)/ 0/
DATA CF( 61)/-.27311873046773076105553203577757D+06/
DATA I( 61)/ 5/,J( 61)/ 0/,K( 61)/ 1/
DATA CF( 62)/-.27311873046773076105553203577757D+06/
DATA I( 62)/ 1/,J( 62)/ 5/,K( 62)/ 0/
DATA CF( 63)/-.27311873046773076105553203577757D+06/
DATA I( 63)/ 1/,J( 63)/ 0/,K( 63)/ 5/
DATA CF( 64)/-.27311873046773076105553203577757D+06/
DATA I( 64)/ 0/,J( 64)/ 5/,K( 64)/ 1/
DATA CF( 65)/-.27311873046773076105553203577757D+06/
DATA I( 65)/ 0/,J( 65)/ 1/,K( 65)/ 5/
DATA CF( 66)/-.56485441456054714255055820418740D+07/
DATA I( 66)/ 3/,J( 66)/ 2/,K( 66)/ 2/
DATA CF( 67)/-.56485441456054714255055820418740D+07/
DATA I( 67)/ 2/,J( 67)/ 3/,K( 67)/ 2/
DATA CF( 68)/-.56485441456054714255055820418740D+07/
DATA I( 68)/ 2/,J( 68)/ 2/,K( 68)/ 3/
DATA CF( 69)/-.36187475096593855297382628456987D+07/
DATA I( 69)/ 3/,J( 69)/ 3/,K( 69)/ 1/
DATA CF( 70)/-.36187475096593855297382628456987D+07/
DATA I( 70)/ 3/,J( 70)/ 1/,K( 70)/ 3/
DATA CF( 71)/-.36187475096593855297382628456987D+07/
DATA I( 71)/ 1/,J( 71)/ 3/,K( 71)/ 3/
DATA CF( 72)/-.40181535619630801734866576836888D+06/
DATA I( 72)/ 4/,J( 72)/ 2/,K( 72)/ 1/
DATA CF( 73)/-.40181535619630801734866576836888D+06/
DATA I( 73)/ 4/,J( 73)/ 1/,K( 73)/ 2/
DATA CF( 74)/-.40181535619630801734866576836888D+06/
DATA I( 74)/ 2/,J( 74)/ 4/,K( 74)/ 1/
DATA CF( 75)/-.40181535619630801734866576836888D+06/
DATA I( 75)/ 2/,J( 75)/ 1/,K( 75)/ 4/
DATA CF( 76)/-.40181535619630801734866576836888D+06/
DATA I( 76)/ 1/,J( 76)/ 4/,K( 76)/ 2/
DATA CF( 77)/-.40181535619630801734866576836888D+06/
DATA I( 77)/ 1/,J( 77)/ 2/,K( 77)/ 4/
DATA CF( 78)/-.26135417638587459421695613435314D+07/
DATA I( 78)/ 4/,J( 78)/ 3/,K( 78)/ 0/
DATA CF( 79)/-.26135417638587459421695613435314D+07/
DATA I( 79)/ 4/,J( 79)/ 0/,K( 79)/ 3/
DATA CF( 80)/-.26135417638587459421695613435314D+07/
DATA I( 80)/ 3/,J( 80)/ 4/,K( 80)/ 0/
DATA CF( 81)/-.26135417638587459421695613435314D+07/
DATA I( 81)/ 3/,J( 81)/ 0/,K( 81)/ 4/
DATA CF( 82)/-.26135417638587459421695613435314D+07/
DATA I( 82)/ 0/,J( 82)/ 4/,K( 82)/ 3/
DATA CF( 83)/-.26135417638587459421695613435314D+07/
DATA I( 83)/ 0/,J( 83)/ 3/,K( 83)/ 4/
DATA CF( 84)/0.37516486872187164322777082881040D+06/
DATA I( 84)/ 5/,J( 84)/ 1/,K( 84)/ 1/
DATA CF( 85)/0.37516486872187164322777082881040D+06/
DATA I( 85)/ 1/,J( 85)/ 5/,K( 85)/ 1/
DATA CF( 86)/0.37516486872187164322777082881040D+06/
DATA I( 86)/ 1/,J( 86)/ 1/,K( 86)/ 5/
DATA CF( 87)/-.14825443019185591607764716682947D+06/
DATA I( 87)/ 5/,J( 87)/ 2/,K( 87)/ 0/
DATA CF( 88)/-.14825443019185591607764716682947D+06/
DATA I( 88)/ 5/,J( 88)/ 0/,K( 88)/ 2/
DATA CF( 89)/-.14825443019185591607764716682947D+06/
DATA I( 89)/ 2/,J( 89)/ 5/,K( 89)/ 0/
DATA CF( 90)/-.14825443019185591607764716682947D+06/
DATA I( 90)/ 2/,J( 90)/ 0/,K( 90)/ 5/
DATA CF( 91)/-.14825443019185591607764716682947D+06/
DATA I( 91)/ 0/,J( 91)/ 5/,K( 91)/ 2/
DATA CF( 92)/-.14825443019185591607764716682947D+06/
DATA I( 92)/ 0/,J( 92)/ 2/,K( 92)/ 5/
DATA CF( 93)/0.19898213723338149128222361551379D+07/
DATA I( 93)/ 6/,J( 93)/ 1/,K( 93)/ 0/
DATA CF( 94)/0.19898213723338149128222361551379D+07/
DATA I( 94)/ 6/,J( 94)/ 0/,K( 94)/ 1/
DATA CF( 95)/0.19898213723338149128222361551379D+07/
DATA I( 95)/ 1/,J( 95)/ 6/,K( 95)/ 0/
DATA CF( 96)/0.19898213723338149128222361551379D+07/
DATA I( 96)/ 1/,J( 96)/ 0/,K( 96)/ 6/
DATA CF( 97)/0.19898213723338149128222361551379D+07/
DATA I( 97)/ 0/,J( 97)/ 6/,K( 97)/ 1/
DATA CF( 98)/0.19898213723338149128222361551379D+07/
DATA I( 98)/ 0/,J( 98)/ 1/,K( 98)/ 6/
DATA CF( 99)/0.18739308703982654194817220130164D+08/
DATA I( 99)/ 3/,J( 99)/ 3/,K( 99)/ 2/
DATA CF(100)/0.18739308703982654194817220130164D+08/
DATA I(100)/ 3/,J(100)/ 2/,K(100)/ 3/
DATA CF(101)/0.18739308703982654194817220130164D+08/
DATA I(101)/ 2/,J(101)/ 3/,K(101)/ 3/
DATA CF(102)/0.20132189123118299532126527673379D+08/
DATA I(102)/ 4/,J(102)/ 2/,K(102)/ 2/
DATA CF(103)/0.20132189123118299532126527673379D+08/
DATA I(103)/ 2/,J(103)/ 4/,K(103)/ 2/
DATA CF(104)/0.20132189123118299532126527673379D+08/
DATA I(104)/ 2/,J(104)/ 2/,K(104)/ 4/
DATA CF(105)/0.11155480577023401095782070039961D+08/
DATA I(105)/ 4/,J(105)/ 3/,K(105)/ 1/
DATA CF(106)/0.11155480577023401095782070039961D+08/
DATA I(106)/ 4/,J(106)/ 1/,K(106)/ 3/
DATA CF(107)/0.11155480577023401095782070039961D+08/
DATA I(107)/ 3/,J(107)/ 4/,K(107)/ 1/
DATA CF(108)/0.11155480577023401095782070039961D+08/
DATA I(108)/ 3/,J(108)/ 1/,K(108)/ 4/
DATA CF(109)/0.11155480577023401095782070039961D+08/
DATA I(109)/ 1/,J(109)/ 4/,K(109)/ 3/
DATA CF(110)/0.11155480577023401095782070039961D+08/
DATA I(110)/ 1/,J(110)/ 3/,K(110)/ 4/
DATA CF(111)/0.16624368438223255207626155212816D+08/
DATA I(111)/ 4/,J(111)/ 4/,K(111)/ 0/
DATA CF(112)/0.16624368438223255207626155212816D+08/
DATA I(112)/ 4/,J(112)/ 0/,K(112)/ 4/
DATA CF(113)/0.16624368438223255207626155212816D+08/
DATA I(113)/ 0/,J(113)/ 4/,K(113)/ 4/
DATA CF(114)/-.14977679449302366501725995228152D+08/
DATA I(114)/ 5/,J(114)/ 2/,K(114)/ 1/
DATA CF(115)/-.14977679449302366501725995228152D+08/
DATA I(115)/ 5/,J(115)/ 1/,K(115)/ 2/
DATA CF(116)/-.14977679449302366501725995228152D+08/
DATA I(116)/ 2/,J(116)/ 5/,K(116)/ 1/
DATA CF(117)/-.14977679449302366501725995228152D+08/
DATA I(117)/ 2/,J(117)/ 1/,K(117)/ 5/
DATA CF(118)/-.14977679449302366501725995228152D+08/
DATA I(118)/ 1/,J(118)/ 5/,K(118)/ 2/
DATA CF(119)/-.14977679449302366501725995228152D+08/
DATA I(119)/ 1/,J(119)/ 2/,K(119)/ 5/
DATA CF(120)/0.19943965720153235321866127593855D+07/
DATA I(120)/ 5/,J(120)/ 3/,K(120)/ 0/
DATA CF(121)/0.19943965720153235321866127593855D+07/
DATA I(121)/ 5/,J(121)/ 0/,K(121)/ 3/
DATA CF(122)/0.19943965720153235321866127593855D+07/
DATA I(122)/ 3/,J(122)/ 5/,K(122)/ 0/
DATA CF(123)/0.19943965720153235321866127593855D+07/
DATA I(123)/ 3/,J(123)/ 0/,K(123)/ 5/
DATA CF(124)/0.19943965720153235321866127593855D+07/
DATA I(124)/ 0/,J(124)/ 5/,K(124)/ 3/
DATA CF(125)/0.19943965720153235321866127593855D+07/
DATA I(125)/ 0/,J(125)/ 3/,K(125)/ 5/
DATA CF(126)/0.13626363297493910760960565233653D+08/
DATA I(126)/ 6/,J(126)/ 1/,K(126)/ 1/
DATA CF(127)/0.13626363297493910760960565233653D+08/
DATA I(127)/ 1/,J(127)/ 6/,K(127)/ 1/
DATA CF(128)/0.13626363297493910760960565233653D+08/
DATA I(128)/ 1/,J(128)/ 1/,K(128)/ 6/
DATA CF(129)/0.26591250219558287007352416247288D+07/
DATA I(129)/ 6/,J(129)/ 2/,K(129)/ 0/
DATA CF(130)/0.26591250219558287007352416247288D+07/
DATA I(130)/ 6/,J(130)/ 0/,K(130)/ 2/
DATA CF(131)/0.26591250219558287007352416247288D+07/
DATA I(131)/ 2/,J(131)/ 6/,K(131)/ 0/
DATA CF(132)/0.26591250219558287007352416247288D+07/
DATA I(132)/ 2/,J(132)/ 0/,K(132)/ 6/
DATA CF(133)/0.26591250219558287007352416247288D+07/
DATA I(133)/ 0/,J(133)/ 6/,K(133)/ 2/
DATA CF(134)/0.26591250219558287007352416247288D+07/
DATA I(134)/ 0/,J(134)/ 2/,K(134)/ 6/
DATA CF(135)/-.10603362125543017112298076225028D+08/
DATA I(135)/ 7/,J(135)/ 1/,K(135)/ 0/
DATA CF(136)/-.10603362125543017112298076225028D+08/
DATA I(136)/ 7/,J(136)/ 0/,K(136)/ 1/
DATA CF(137)/-.10603362125543017112298076225028D+08/
DATA I(137)/ 1/,J(137)/ 7/,K(137)/ 0/
DATA CF(138)/-.10603362125543017112298076225028D+08/
DATA I(138)/ 1/,J(138)/ 0/,K(138)/ 7/
DATA CF(139)/-.10603362125543017112298076225028D+08/
DATA I(139)/ 0/,J(139)/ 7/,K(139)/ 1/
DATA CF(140)/-.10603362125543017112298076225028D+08/
DATA I(140)/ 0/,J(140)/ 1/,K(140)/ 7/
DATA CF(141)/-.46586804047421649167456709691511D+07/
DATA I(141)/ 3/,J(141)/ 3/,K(141)/ 3/
DATA CF(142)/-.50072086671134631451525240221034D+08/
DATA I(142)/ 4/,J(142)/ 3/,K(142)/ 2/
DATA CF(143)/-.50072086671134631451525240221034D+08/
DATA I(143)/ 4/,J(143)/ 2/,K(143)/ 3/
DATA CF(144)/-.50072086671134631451525240221034D+08/
DATA I(144)/ 3/,J(144)/ 4/,K(144)/ 2/
DATA CF(145)/-.50072086671134631451525240221034D+08/
DATA I(145)/ 3/,J(145)/ 2/,K(145)/ 4/
DATA CF(146)/-.50072086671134631451525240221034D+08/
DATA I(146)/ 2/,J(146)/ 4/,K(146)/ 3/
DATA CF(147)/-.50072086671134631451525240221034D+08/
DATA I(147)/ 2/,J(147)/ 3/,K(147)/ 4/
DATA CF(148)/0.63067299559141288404858221452481D+08/
DATA I(148)/ 4/,J(148)/ 4/,K(148)/ 1/
DATA CF(149)/0.63067299559141288404858221452481D+08/
DATA I(149)/ 4/,J(149)/ 1/,K(149)/ 4/
DATA CF(150)/0.63067299559141288404858221452481D+08/
DATA I(150)/ 1/,J(150)/ 4/,K(150)/ 4/
DATA CF(151)/-.46183485511361133054452833472168D+08/
DATA I(151)/ 5/,J(151)/ 2/,K(151)/ 2/
DATA CF(152)/-.46183485511361133054452833472168D+08/
DATA I(152)/ 2/,J(152)/ 5/,K(152)/ 2/
DATA CF(153)/-.46183485511361133054452833472168D+08/
DATA I(153)/ 2/,J(153)/ 2/,K(153)/ 5/
DATA CF(154)/-.72895989122727108071150503318796D+08/
DATA I(154)/ 5/,J(154)/ 3/,K(154)/ 1/
DATA CF(155)/-.72895989122727108071150503318796D+08/
DATA I(155)/ 5/,J(155)/ 1/,K(155)/ 3/
DATA CF(156)/-.72895989122727108071150503318796D+08/
DATA I(156)/ 3/,J(156)/ 5/,K(156)/ 1/
DATA CF(157)/-.72895989122727108071150503318796D+08/
DATA I(157)/ 3/,J(157)/ 1/,K(157)/ 5/
DATA CF(158)/-.72895989122727108071150503318796D+08/
DATA I(158)/ 1/,J(158)/ 5/,K(158)/ 3/
DATA CF(159)/-.72895989122727108071150503318796D+08/
DATA I(159)/ 1/,J(159)/ 3/,K(159)/ 5/
DATA CF(160)/-.33502851126217904274051721498598D+08/
DATA I(160)/ 5/,J(160)/ 4/,K(160)/ 0/
DATA CF(161)/-.33502851126217904274051721498598D+08/
DATA I(161)/ 5/,J(161)/ 0/,K(161)/ 4/
DATA CF(162)/-.33502851126217904274051721498598D+08/
DATA I(162)/ 4/,J(162)/ 5/,K(162)/ 0/
DATA CF(163)/-.33502851126217904274051721498598D+08/
DATA I(163)/ 4/,J(163)/ 0/,K(163)/ 5/
DATA CF(164)/-.33502851126217904274051721498598D+08/
DATA I(164)/ 0/,J(164)/ 5/,K(164)/ 4/
DATA CF(165)/-.33502851126217904274051721498598D+08/
DATA I(165)/ 0/,J(165)/ 4/,K(165)/ 5/
DATA CF(166)/0.10072719455528817262310606069143D+09/
DATA I(166)/ 6/,J(166)/ 2/,K(166)/ 1/
DATA CF(167)/0.10072719455528817262310606069143D+09/
DATA I(167)/ 6/,J(167)/ 1/,K(167)/ 2/
DATA CF(168)/0.10072719455528817262310606069143D+09/
DATA I(168)/ 2/,J(168)/ 6/,K(168)/ 1/
DATA CF(169)/0.10072719455528817262310606069143D+09/
DATA I(169)/ 2/,J(169)/ 1/,K(169)/ 6/
DATA CF(170)/0.10072719455528817262310606069143D+09/
DATA I(170)/ 1/,J(170)/ 6/,K(170)/ 2/
DATA CF(171)/0.10072719455528817262310606069143D+09/
DATA I(171)/ 1/,J(171)/ 2/,K(171)/ 6/
DATA CF(172)/0.56681657887384393627269408358825D+07/
DATA I(172)/ 6/,J(172)/ 3/,K(172)/ 0/
DATA CF(173)/0.56681657887384393627269408358825D+07/
DATA I(173)/ 6/,J(173)/ 0/,K(173)/ 3/
DATA CF(174)/0.56681657887384393627269408358825D+07/
DATA I(174)/ 3/,J(174)/ 6/,K(174)/ 0/
DATA CF(175)/0.56681657887384393627269408358825D+07/
DATA I(175)/ 3/,J(175)/ 0/,K(175)/ 6/
DATA CF(176)/0.56681657887384393627269408358825D+07/
DATA I(176)/ 0/,J(176)/ 6/,K(176)/ 3/
DATA CF(177)/0.56681657887384393627269408358825D+07/
DATA I(177)/ 0/,J(177)/ 3/,K(177)/ 6/
DATA CF(178)/-.95925361906341497286942014904531D+08/
DATA I(178)/ 7/,J(178)/ 1/,K(178)/ 1/
DATA CF(179)/-.95925361906341497286942014904531D+08/
DATA I(179)/ 1/,J(179)/ 7/,K(179)/ 1/
DATA CF(180)/-.95925361906341497286942014904531D+08/
DATA I(180)/ 1/,J(180)/ 1/,K(180)/ 7/
DATA CF(181)/-.13424883050104978878474408465758D+08/
DATA I(181)/ 7/,J(181)/ 2/,K(181)/ 0/
DATA CF(182)/-.13424883050104978878474408465758D+08/
DATA I(182)/ 7/,J(182)/ 0/,K(182)/ 2/
DATA CF(183)/-.13424883050104978878474408465758D+08/
DATA I(183)/ 2/,J(183)/ 7/,K(183)/ 0/
DATA CF(184)/-.13424883050104978878474408465758D+08/
DATA I(184)/ 2/,J(184)/ 0/,K(184)/ 7/
DATA CF(185)/-.13424883050104978878474408465758D+08/
DATA I(185)/ 0/,J(185)/ 7/,K(185)/ 2/
DATA CF(186)/-.13424883050104978878474408465758D+08/
DATA I(186)/ 0/,J(186)/ 2/,K(186)/ 7/
DATA CF(187)/0.35080556353049139044777189322872D+08/
DATA I(187)/ 8/,J(187)/ 1/,K(187)/ 0/
DATA CF(188)/0.35080556353049139044777189322872D+08/
DATA I(188)/ 8/,J(188)/ 0/,K(188)/ 1/
DATA CF(189)/0.35080556353049139044777189322872D+08/
DATA I(189)/ 1/,J(189)/ 8/,K(189)/ 0/
DATA CF(190)/0.35080556353049139044777189322872D+08/
DATA I(190)/ 1/,J(190)/ 0/,K(190)/ 8/
DATA CF(191)/0.35080556353049139044777189322872D+08/
DATA I(191)/ 0/,J(191)/ 8/,K(191)/ 1/
DATA CF(192)/0.35080556353049139044777189322872D+08/
DATA I(192)/ 0/,J(192)/ 1/,K(192)/ 8/
DATA CF(193)/0.21313165250776255116174382605637D+08/
DATA I(193)/ 4/,J(193)/ 3/,K(193)/ 3/
DATA CF(194)/0.21313165250776255116174382605637D+08/
DATA I(194)/ 3/,J(194)/ 4/,K(194)/ 3/
DATA CF(195)/0.21313165250776255116174382605637D+08/
DATA I(195)/ 3/,J(195)/ 3/,K(195)/ 4/
DATA CF(196)/0.33715830673581714383277368128765D+09/
DATA I(196)/ 4/,J(196)/ 4/,K(196)/ 2/
DATA CF(197)/0.33715830673581714383277368128765D+09/
DATA I(197)/ 4/,J(197)/ 2/,K(197)/ 4/
DATA CF(198)/0.33715830673581714383277368128765D+09/
DATA I(198)/ 2/,J(198)/ 4/,K(198)/ 4/
DATA CF(199)/-.13047887760177294841027911534426D+09/
DATA I(199)/ 5/,J(199)/ 3/,K(199)/ 2/
DATA CF(200)/-.13047887760177294841027911534426D+09/
DATA I(200)/ 5/,J(200)/ 2/,K(200)/ 3/
DATA CF(201)/-.13047887760177294841027911534426D+09/
DATA I(201)/ 3/,J(201)/ 5/,K(201)/ 2/
DATA CF(202)/-.13047887760177294841027911534426D+09/
DATA I(202)/ 3/,J(202)/ 2/,K(202)/ 5/
DATA CF(203)/-.13047887760177294841027911534426D+09/
DATA I(203)/ 2/,J(203)/ 5/,K(203)/ 3/
DATA CF(204)/-.13047887760177294841027911534426D+09/
DATA I(204)/ 2/,J(204)/ 3/,K(204)/ 5/
DATA CF(205)/-.98725001335457966924651673228532D+08/
DATA I(205)/ 5/,J(205)/ 4/,K(205)/ 1/
DATA CF(206)/-.98725001335457966924651673228532D+08/
DATA I(206)/ 5/,J(206)/ 1/,K(206)/ 4/
DATA CF(207)/-.98725001335457966924651673228532D+08/
DATA I(207)/ 4/,J(207)/ 5/,K(207)/ 1/
DATA CF(208)/-.98725001335457966924651673228532D+08/
DATA I(208)/ 4/,J(208)/ 1/,K(208)/ 5/
DATA CF(209)/-.98725001335457966924651673228532D+08/
DATA I(209)/ 1/,J(209)/ 5/,K(209)/ 4/
DATA CF(210)/-.98725001335457966924651673228532D+08/
DATA I(210)/ 1/,J(210)/ 4/,K(210)/ 5/
DATA CF(211)/0.12313020517064100209904542848769D+09/
DATA I(211)/ 5/,J(211)/ 5/,K(211)/ 0/
DATA CF(212)/0.12313020517064100209904542848769D+09/
DATA I(212)/ 5/,J(212)/ 0/,K(212)/ 5/
DATA CF(213)/0.12313020517064100209904542848769D+09/
DATA I(213)/ 0/,J(213)/ 5/,K(213)/ 5/
DATA CF(214)/0.21360958570731033908361223452011D+09/
DATA I(214)/ 6/,J(214)/ 2/,K(214)/ 2/
DATA CF(215)/0.21360958570731033908361223452011D+09/
DATA I(215)/ 2/,J(215)/ 6/,K(215)/ 2/
DATA CF(216)/0.21360958570731033908361223452011D+09/
DATA I(216)/ 2/,J(216)/ 2/,K(216)/ 6/
DATA CF(217)/0.18402806647355499508765877624796D+09/
DATA I(217)/ 6/,J(217)/ 3/,K(217)/ 1/
DATA CF(218)/0.18402806647355499508765877624796D+09/
DATA I(218)/ 6/,J(218)/ 1/,K(218)/ 3/
DATA CF(219)/0.18402806647355499508765877624796D+09/
DATA I(219)/ 3/,J(219)/ 6/,K(219)/ 1/
DATA CF(220)/0.18402806647355499508765877624796D+09/
DATA I(220)/ 3/,J(220)/ 1/,K(220)/ 6/
DATA CF(221)/0.18402806647355499508765877624796D+09/
DATA I(221)/ 1/,J(221)/ 6/,K(221)/ 3/
DATA CF(222)/0.18402806647355499508765877624796D+09/
DATA I(222)/ 1/,J(222)/ 3/,K(222)/ 6/
DATA CF(223)/-.18754201886738159269921471865807D+08/
DATA I(223)/ 6/,J(223)/ 4/,K(223)/ 0/
DATA CF(224)/-.18754201886738159269921471865807D+08/
DATA I(224)/ 6/,J(224)/ 0/,K(224)/ 4/
DATA CF(225)/-.18754201886738159269921471865807D+08/
DATA I(225)/ 4/,J(225)/ 6/,K(225)/ 0/
DATA CF(226)/-.18754201886738159269921471865807D+08/
DATA I(226)/ 4/,J(226)/ 0/,K(226)/ 6/
DATA CF(227)/-.18754201886738159269921471865807D+08/
DATA I(227)/ 0/,J(227)/ 6/,K(227)/ 4/
DATA CF(228)/-.18754201886738159269921471865807D+08/
DATA I(228)/ 0/,J(228)/ 4/,K(228)/ 6/
DATA CF(229)/-.26983444083488089431413304065887D+09/
DATA I(229)/ 7/,J(229)/ 2/,K(229)/ 1/
DATA CF(230)/-.26983444083488089431413304065887D+09/
DATA I(230)/ 7/,J(230)/ 1/,K(230)/ 2/
DATA CF(231)/-.26983444083488089431413304065887D+09/
DATA I(231)/ 2/,J(231)/ 7/,K(231)/ 1/
DATA CF(232)/-.26983444083488089431413304065887D+09/
DATA I(232)/ 2/,J(232)/ 1/,K(232)/ 7/
DATA CF(233)/-.26983444083488089431413304065887D+09/
DATA I(233)/ 1/,J(233)/ 7/,K(233)/ 2/
DATA CF(234)/-.26983444083488089431413304065887D+09/
DATA I(234)/ 1/,J(234)/ 2/,K(234)/ 7/
DATA CF(235)/0.10003922186620932158069727344562D+08/
DATA I(235)/ 7/,J(235)/ 3/,K(235)/ 0/
DATA CF(236)/0.10003922186620932158069727344562D+08/
DATA I(236)/ 7/,J(236)/ 0/,K(236)/ 3/
DATA CF(237)/0.10003922186620932158069727344562D+08/
DATA I(237)/ 3/,J(237)/ 7/,K(237)/ 0/
DATA CF(238)/0.10003922186620932158069727344562D+08/
DATA I(238)/ 3/,J(238)/ 0/,K(238)/ 7/
DATA CF(239)/0.10003922186620932158069727344562D+08/
DATA I(239)/ 0/,J(239)/ 7/,K(239)/ 3/
DATA CF(240)/0.10003922186620932158069727344562D+08/
DATA I(240)/ 0/,J(240)/ 3/,K(240)/ 7/
DATA CF(241)/0.24850791774710503554465144114161D+09/
DATA I(241)/ 8/,J(241)/ 1/,K(241)/ 1/
DATA CF(242)/0.24850791774710503554465144114161D+09/
DATA I(242)/ 1/,J(242)/ 8/,K(242)/ 1/
DATA CF(243)/0.24850791774710503554465144114161D+09/
DATA I(243)/ 1/,J(243)/ 1/,K(243)/ 8/
DATA CF(244)/0.21706723533830221321729370428596D+08/
DATA I(244)/ 8/,J(244)/ 2/,K(244)/ 0/
DATA CF(245)/0.21706723533830221321729370428596D+08/
DATA I(245)/ 8/,J(245)/ 0/,K(245)/ 2/
DATA CF(246)/0.21706723533830221321729370428596D+08/
DATA I(246)/ 2/,J(246)/ 8/,K(246)/ 0/
DATA CF(247)/0.21706723533830221321729370428596D+08/
DATA I(247)/ 2/,J(247)/ 0/,K(247)/ 8/
DATA CF(248)/0.21706723533830221321729370428596D+08/
DATA I(248)/ 0/,J(248)/ 8/,K(248)/ 2/
DATA CF(249)/0.21706723533830221321729370428596D+08/
DATA I(249)/ 0/,J(249)/ 2/,K(249)/ 8/
DATA CF(250)/-.62285853147928862367989342510437D+08/
DATA I(250)/ 9/,J(250)/ 1/,K(250)/ 0/
DATA CF(251)/-.62285853147928862367989342510437D+08/
DATA I(251)/ 9/,J(251)/ 0/,K(251)/ 1/
DATA CF(252)/-.62285853147928862367989342510437D+08/
DATA I(252)/ 1/,J(252)/ 9/,K(252)/ 0/
DATA CF(253)/-.62285853147928862367989342510437D+08/
DATA I(253)/ 1/,J(253)/ 0/,K(253)/ 9/
DATA CF(254)/-.62285853147928862367989342510437D+08/
DATA I(254)/ 0/,J(254)/ 9/,K(254)/ 1/
DATA CF(255)/-.62285853147928862367989342510437D+08/
DATA I(255)/ 0/,J(255)/ 1/,K(255)/ 9/
DATA CF(256)/0.71914062149624329879935815273289D+09/
DATA I(256)/ 4/,J(256)/ 4/,K(256)/ 3/
DATA CF(257)/0.71914062149624329879935815273289D+09/
DATA I(257)/ 4/,J(257)/ 3/,K(257)/ 4/
DATA CF(258)/0.71914062149624329879935815273289D+09/
DATA I(258)/ 3/,J(258)/ 4/,K(258)/ 4/
DATA CF(259)/-.12258942322137567155438157196262D+10/
DATA I(259)/ 5/,J(259)/ 3/,K(259)/ 3/
DATA CF(260)/-.12258942322137567155438157196262D+10/
DATA I(260)/ 3/,J(260)/ 5/,K(260)/ 3/
DATA CF(261)/-.12258942322137567155438157196262D+10/
DATA I(261)/ 3/,J(261)/ 3/,K(261)/ 5/
DATA CF(262)/-.21798005283565026935713712646780D+09/
DATA I(262)/ 5/,J(262)/ 4/,K(262)/ 2/
DATA CF(263)/-.21798005283565026935713712646780D+09/
DATA I(263)/ 5/,J(263)/ 2/,K(263)/ 4/
DATA CF(264)/-.21798005283565026935713712646780D+09/
DATA I(264)/ 4/,J(264)/ 5/,K(264)/ 2/
DATA CF(265)/-.21798005283565026935713712646780D+09/
DATA I(265)/ 4/,J(265)/ 2/,K(265)/ 5/
DATA CF(266)/-.21798005283565026935713712646780D+09/
DATA I(266)/ 2/,J(266)/ 5/,K(266)/ 4/
DATA CF(267)/-.21798005283565026935713712646780D+09/
DATA I(267)/ 2/,J(267)/ 4/,K(267)/ 5/
DATA CF(268)/0.50624308448637753756469430350211D+09/
DATA I(268)/ 5/,J(268)/ 5/,K(268)/ 1/
DATA CF(269)/0.50624308448637753756469430350211D+09/
DATA I(269)/ 5/,J(269)/ 1/,K(269)/ 5/
DATA CF(270)/0.50624308448637753756469430350211D+09/
DATA I(270)/ 1/,J(270)/ 5/,K(270)/ 5/
DATA CF(271)/0.63673139780850884861127040932408D+09/
DATA I(271)/ 6/,J(271)/ 3/,K(271)/ 2/
DATA CF(272)/0.63673139780850884861127040932408D+09/
DATA I(272)/ 6/,J(272)/ 2/,K(272)/ 3/
DATA CF(273)/0.63673139780850884861127040932408D+09/
DATA I(273)/ 3/,J(273)/ 6/,K(273)/ 2/
DATA CF(274)/0.63673139780850884861127040932408D+09/
DATA I(274)/ 3/,J(274)/ 2/,K(274)/ 6/
DATA CF(275)/0.63673139780850884861127040932408D+09/
DATA I(275)/ 2/,J(275)/ 6/,K(275)/ 3/
DATA CF(276)/0.63673139780850884861127040932408D+09/
DATA I(276)/ 2/,J(276)/ 3/,K(276)/ 6/
DATA CF(277)/-.30448335666417391836922042075543D+09/
DATA I(277)/ 6/,J(277)/ 4/,K(277)/ 1/
DATA CF(278)/-.30448335666417391836922042075543D+09/
DATA I(278)/ 6/,J(278)/ 1/,K(278)/ 4/
DATA CF(279)/-.30448335666417391836922042075543D+09/
DATA I(279)/ 4/,J(279)/ 6/,K(279)/ 1/
DATA CF(280)/-.30448335666417391836922042075543D+09/
DATA I(280)/ 4/,J(280)/ 1/,K(280)/ 6/
DATA CF(281)/-.30448335666417391836922042075543D+09/
DATA I(281)/ 1/,J(281)/ 6/,K(281)/ 4/
DATA CF(282)/-.30448335666417391836922042075543D+09/
DATA I(282)/ 1/,J(282)/ 4/,K(282)/ 6/
DATA CF(283)/-.75560911535432422981760748915256D+08/
DATA I(283)/ 6/,J(283)/ 5/,K(283)/ 0/
DATA CF(284)/-.75560911535432422981760748915256D+08/
DATA I(284)/ 6/,J(284)/ 0/,K(284)/ 5/
DATA CF(285)/-.75560911535432422981760748915256D+08/
DATA I(285)/ 5/,J(285)/ 6/,K(285)/ 0/
DATA CF(286)/-.75560911535432422981760748915256D+08/
DATA I(286)/ 5/,J(286)/ 0/,K(286)/ 6/
DATA CF(287)/-.75560911535432422981760748915256D+08/
DATA I(287)/ 0/,J(287)/ 6/,K(287)/ 5/
DATA CF(288)/-.75560911535432422981760748915256D+08/
DATA I(288)/ 0/,J(288)/ 5/,K(288)/ 6/
DATA CF(289)/-.51966065093426703154268264157805D+09/
DATA I(289)/ 7/,J(289)/ 2/,K(289)/ 2/
DATA CF(290)/-.51966065093426703154268264157805D+09/
DATA I(290)/ 2/,J(290)/ 7/,K(290)/ 2/
DATA CF(291)/-.51966065093426703154268264157805D+09/
DATA I(291)/ 2/,J(291)/ 2/,K(291)/ 7/
DATA CF(292)/-.70629011340840606134339365502511D+08/
DATA I(292)/ 7/,J(292)/ 3/,K(292)/ 1/
DATA CF(293)/-.70629011340840606134339365502511D+08/
DATA I(293)/ 7/,J(293)/ 1/,K(293)/ 3/
DATA CF(294)/-.70629011340840606134339365502511D+08/
DATA I(294)/ 3/,J(294)/ 7/,K(294)/ 1/
DATA CF(295)/-.70629011340840606134339365502511D+08/
DATA I(295)/ 3/,J(295)/ 1/,K(295)/ 7/
DATA CF(296)/-.70629011340840606134339365502511D+08/
DATA I(296)/ 1/,J(296)/ 7/,K(296)/ 3/
DATA CF(297)/-.70629011340840606134339365502511D+08/
DATA I(297)/ 1/,J(297)/ 3/,K(297)/ 7/
DATA CF(298)/0.73883357514484293938215349416806D+08/
DATA I(298)/ 7/,J(298)/ 4/,K(298)/ 0/
DATA CF(299)/0.73883357514484293938215349416806D+08/
DATA I(299)/ 7/,J(299)/ 0/,K(299)/ 4/
DATA CF(300)/0.73883357514484293938215349416806D+08/
DATA I(300)/ 4/,J(300)/ 7/,K(300)/ 0/
DATA CF(301)/0.73883357514484293938215349416806D+08/
DATA I(301)/ 4/,J(301)/ 0/,K(301)/ 7/
DATA CF(302)/0.73883357514484293938215349416806D+08/
DATA I(302)/ 0/,J(302)/ 7/,K(302)/ 4/
DATA CF(303)/0.73883357514484293938215349416806D+08/
DATA I(303)/ 0/,J(303)/ 4/,K(303)/ 7/
DATA CF(304)/0.26574545795392667263357823250913D+09/
DATA I(304)/ 8/,J(304)/ 2/,K(304)/ 1/
DATA CF(305)/0.26574545795392667263357823250913D+09/
DATA I(305)/ 8/,J(305)/ 1/,K(305)/ 2/
DATA CF(306)/0.26574545795392667263357823250913D+09/
DATA I(306)/ 2/,J(306)/ 8/,K(306)/ 1/
DATA CF(307)/0.26574545795392667263357823250913D+09/
DATA I(307)/ 2/,J(307)/ 1/,K(307)/ 8/
DATA CF(308)/0.26574545795392667263357823250913D+09/
DATA I(308)/ 1/,J(308)/ 8/,K(308)/ 2/
DATA CF(309)/0.26574545795392667263357823250913D+09/
DATA I(309)/ 1/,J(309)/ 2/,K(309)/ 8/
DATA CF(310)/-.41282706219359481969277501766754D+08/
DATA I(310)/ 8/,J(310)/ 3/,K(310)/ 0/
DATA CF(311)/-.41282706219359481969277501766754D+08/
DATA I(311)/ 8/,J(311)/ 0/,K(311)/ 3/
DATA CF(312)/-.41282706219359481969277501766754D+08/
DATA I(312)/ 3/,J(312)/ 8/,K(312)/ 0/
DATA CF(313)/-.41282706219359481969277501766754D+08/
DATA I(313)/ 3/,J(313)/ 0/,K(313)/ 8/
DATA CF(314)/-.41282706219359481969277501766754D+08/
DATA I(314)/ 0/,J(314)/ 8/,K(314)/ 3/
DATA CF(315)/-.41282706219359481969277501766754D+08/
DATA I(315)/ 0/,J(315)/ 3/,K(315)/ 8/
DATA CF(316)/-.23007097631534204110452149063500D+09/
DATA I(316)/ 9/,J(316)/ 1/,K(316)/ 1/
DATA CF(317)/-.23007097631534204110452149063500D+09/
DATA I(317)/ 1/,J(317)/ 9/,K(317)/ 1/
DATA CF(318)/-.23007097631534204110452149063500D+09/
DATA I(318)/ 1/,J(318)/ 1/,K(318)/ 9/
DATA CF(319)/-.85393983340225713105855823357489D+07/
DATA I(319)/ 9/,J(319)/ 2/,K(319)/ 0/
DATA CF(320)/-.85393983340225713105855823357489D+07/
DATA I(320)/ 9/,J(320)/ 0/,K(320)/ 2/
DATA CF(321)/-.85393983340225713105855823357489D+07/
DATA I(321)/ 2/,J(321)/ 9/,K(321)/ 0/
DATA CF(322)/-.85393983340225713105855823357489D+07/
DATA I(322)/ 2/,J(322)/ 0/,K(322)/ 9/
DATA CF(323)/-.85393983340225713105855823357489D+07/
DATA I(323)/ 0/,J(323)/ 9/,K(323)/ 2/
DATA CF(324)/-.85393983340225713105855823357489D+07/
DATA I(324)/ 0/,J(324)/ 2/,K(324)/ 9/
DATA CF(325)/0.45163512742840575520571476720490D+08/
DATA I(325)/10/,J(325)/ 1/,K(325)/ 0/
DATA CF(326)/0.45163512742840575520571476720490D+08/
DATA I(326)/10/,J(326)/ 0/,K(326)/ 1/
DATA CF(327)/0.45163512742840575520571476720490D+08/
DATA I(327)/ 1/,J(327)/10/,K(327)/ 0/
DATA CF(328)/0.45163512742840575520571476720490D+08/
DATA I(328)/ 1/,J(328)/ 0/,K(328)/10/
DATA CF(329)/0.45163512742840575520571476720490D+08/
DATA I(329)/ 0/,J(329)/10/,K(329)/ 1/
DATA CF(330)/0.45163512742840575520571476720490D+08/
DATA I(330)/ 0/,J(330)/ 1/,K(330)/10/
DATA XXX(-1)/ 0.D0/,YYY(-1)/ 0.D0/,ZZZ(-1)/ 0.D0/
DATA XXX(0)/ 1.D0/,YYY(0)/ 1.D0/,ZZZ(0)/ 1.D0/
VEX = 0.11660900D+01
X1 = DEXP(-VEX*X)
Y1 = DEXP(-VEX*Y)
Z1 = DEXP(-VEX*Z)
XX = X*X1
YY = Y*Y1
ZZ = Z*Z1
DO 1 L=1,10
XXX(L) = XX*XXX(L-1)
YYY(L) = YY*YYY(L-1)
ZZZ(L) = ZZ*ZZZ(L-1)
1 CONTINUE
***** CALCULO DEL POTENCIAL ****
ENER = 0.D0
DER1 = 0.D0
DER2 = 0.D0
DER3 = 0.D0
DO 2 L=1,330
ENER=ENER+CF(L)*XXX(I(L))*YYY(J(L))*ZZZ(K(L))
DER1=DER1+CF(L)*I(L)*XXX(I(L)-1)*YYY(J(L))*ZZZ(K(L))
DER2=DER2+CF(L)*J(L)*XXX(I(L))*YYY(J(L)-1)*ZZZ(K(L))
DER3=DER3+CF(L)*K(L)*XXX(I(L))*YYY(J(L))*ZZZ(K(L)-1)
2 CONTINUE
DER(1) = DER1 * (1.D0-VEX*X)*X1
DER(2) = DER2 * (1.D0-VEX*Y)*Y1
DER(3) = DER3 * (1.D0-VEX*Z)*Z1
RETURN
END
************************************************************************
* *
* SUBROUTINE AJUSTA *
* *
* ESTA SUBRUTINA CALCULA LAS ENERGIAS DE UNA CURVA DIATOMICA *
* CON UN RMS DE .0027 KCAL/MOL *
* Y UN EMAX DE .0056 KCAL/MOL *
* *
************************************************************************
SUBROUTINE H2LIU(X,ENER)
IMPLICIT REAL*8 (A-H,O-Z)
DIMENSION CF( 0: 6)
DATA CF( 0)/ .105458829431888D+01/
DATA CF( 1)/-.246365037221090D-01/
DATA CF( 2)/-.759893316614319D+01/
DATA CF( 3)/ .419034406717732D+02/
DATA CF( 4)/-.132708934232196D+03/
DATA CF( 5)/ .228198527326218D+03/
DATA CF( 6)/-.166018452889034D+03/
ENER = 0.000000000000000D+00
ENER=ENER+CF( 0)*DEXP(- .230433887268222D+01*X)/X
RHO =X*DEXP(- .109344220007699D+01*X)
DO 1 I= 1, 6
ENER=ENER+CF(I)*RHO**I
1 CONTINUE
RETURN
END
************************************************************************
* *
* SUBROUTINE AJUSTA *
* *
* ESTA SUBRUTINA CALCULA LAS ENERGIAS DE UNA CURVA DIATOMICA *
* CON UN RMS DE .0027 KCAL/MOL *
* Y UN EMAX DE .0056 KCAL/MOL *
* *
************************************************************************
SUBROUTINE H2LIU3(X,ENER,DER)
IMPLICIT REAL*8 (A-H,O-Z)
* DIMENSION CF( 0: 6)
PARAMETER (CF0= .105458829431888D+01)
PARAMETER (CF1=-.246365037221090D-01)
PARAMETER (CF2=-.759893316614319D+01)
PARAMETER (CF3= .419034406717732D+02)
PARAMETER (CF4=-.132708934232196D+03)
PARAMETER (CF5= .228198527326218D+03)
PARAMETER (CF6=-.166018452889034D+03)
PARAMETER (ALFA=.230433887268222D+01)
PARAMETER (BETA=.109344220007699D+01)
AUX = DEXP(- ALFA*X)
BUX = 1.D0/X
CUX = DEXP(- BETA*X)
***** CALCULO DE LA ENERGIA *****
* ENER = -.100000000000000D+01
ENER = 0.D0
ENER=ENER+CF0*AUX*BUX
RHO =X*CUX
ENER=ENER+CF1*RHO
ENER=ENER+CF2*RHO**2
ENER=ENER+CF3*RHO**3
ENER=ENER+CF4*RHO**4
ENER=ENER+CF5*RHO**5
ENER=ENER+CF6*RHO**6
***** CALCULO DE LA DERIVADA *****
DER =CF1 +2.D0*CF2*RHO
DER =DER +3.D0*CF3*RHO**2
DER =DER +4.D0*CF4*RHO**3
DER =DER +5.D0*CF5*RHO**4
DER =DER +6.D0*CF6*RHO**5
DER = DER*(1.D0-BETA*X)*CUX
DER = DER - CF0*BUX*(ALFA+BUX)*AUX
RETURN
END
SUBROUTINE H3aguado7(X,Y,Z,ENER,DER)
IMPLICIT REAL * 8 (A-H,O-Z)
DIMENSION DER(3)
CALL H2LIU3(X,E1,D1)
CALL H2LIU3(Y,E2,D2)
CALL H2LIU3(Z,E3,D3)
CALL V37(X,Y,Z,ENER,DER)
ENER = ENER + E1+E2+E3
DER(1) = DER(1)+D1
DER(2) = DER(2)+D2
DER(3) = DER(3)+D3
RETURN
END
************************************************************************
* *
* SUBROUTINE AJUSTA *
* *
* ESTA SUBRUTINA CALCULA LAS ENERGIAS DE UNA SUPERFICIE BIDIMENSIONAL *
* CON UN RMS DE 1.21071460 KCAL/MOL *
* CON UN EMAX DE 7.94354219 KCAL/MOL
* *
************************************************************************
SUBROUTINE V37(X,Y,Z,ENER,DER)
IMPLICIT REAL * 8 (A-H,O-Z)
DIMENSION I( 98), J( 98), K( 98), CF( 98)
DIMENSION XXX(-1:6), YYY(-1:6), ZZZ(-1:6)
DIMENSION DER(3)
DATA CF( 1)/ .9939012129422977D+01/
DATA I( 1)/ 1/,J( 1)/ 1/,K( 1)/ 0/
DATA CF( 2)/ .9939012129422977D+01/
DATA I( 2)/ 1/,J( 2)/ 0/,K( 2)/ 1/
DATA CF( 3)/ .9939012129422977D+01/
DATA I( 3)/ 0/,J( 3)/ 1/,K( 3)/ 1/
DATA CF( 4)/-.1761366669128080D+02/
DATA I( 4)/ 1/,J( 4)/ 1/,K( 4)/ 1/
DATA CF( 5)/-.4017676962559256D+02/
DATA I( 5)/ 2/,J( 5)/ 1/,K( 5)/ 0/
DATA CF( 6)/-.4017676962559256D+02/
DATA I( 6)/ 2/,J( 6)/ 0/,K( 6)/ 1/
DATA CF( 7)/-.4017676962559256D+02/
DATA I( 7)/ 1/,J( 7)/ 2/,K( 7)/ 0/
DATA CF( 8)/-.4017676962559256D+02/
DATA I( 8)/ 1/,J( 8)/ 0/,K( 8)/ 2/
DATA CF( 9)/-.4017676962559256D+02/
DATA I( 9)/ 0/,J( 9)/ 2/,K( 9)/ 1/
DATA CF( 10)/-.4017676962559256D+02/
DATA I( 10)/ 0/,J( 10)/ 1/,K( 10)/ 2/
DATA CF( 11)/ .2925204854657702D+03/
DATA I( 11)/ 2/,J( 11)/ 1/,K( 11)/ 1/
DATA CF( 12)/ .2925204854657702D+03/
DATA I( 12)/ 1/,J( 12)/ 2/,K( 12)/ 1/
DATA CF( 13)/ .2925204854657702D+03/
DATA I( 13)/ 1/,J( 13)/ 1/,K( 13)/ 2/
DATA CF( 14)/ .1033015975755931D+04/
DATA I( 14)/ 2/,J( 14)/ 2/,K( 14)/ 0/
DATA CF( 15)/ .1033015975755931D+04/
DATA I( 15)/ 2/,J( 15)/ 0/,K( 15)/ 2/
DATA CF( 16)/ .1033015975755931D+04/
DATA I( 16)/ 0/,J( 16)/ 2/,K( 16)/ 2/
DATA CF( 17)/-.3828824120557198D+03/
DATA I( 17)/ 3/,J( 17)/ 1/,K( 17)/ 0/
DATA CF( 18)/-.3828824120557198D+03/
DATA I( 18)/ 3/,J( 18)/ 0/,K( 18)/ 1/
DATA CF( 19)/-.3828824120557198D+03/
DATA I( 19)/ 1/,J( 19)/ 3/,K( 19)/ 0/
DATA CF( 20)/-.3828824120557198D+03/
DATA I( 20)/ 1/,J( 20)/ 0/,K( 20)/ 3/
DATA CF( 21)/-.3828824120557198D+03/
DATA I( 21)/ 0/,J( 21)/ 3/,K( 21)/ 1/
DATA CF( 22)/-.3828824120557198D+03/
DATA I( 22)/ 0/,J( 22)/ 1/,K( 22)/ 3/
DATA CF( 23)/ .4067566006634343D+03/
DATA I( 23)/ 2/,J( 23)/ 2/,K( 23)/ 1/
DATA CF( 24)/ .4067566006634343D+03/
DATA I( 24)/ 2/,J( 24)/ 1/,K( 24)/ 2/
DATA CF( 25)/ .4067566006634343D+03/
DATA I( 25)/ 1/,J( 25)/ 2/,K( 25)/ 2/
DATA CF( 26)/-.3367514367668495D+04/
DATA I( 26)/ 3/,J( 26)/ 1/,K( 26)/ 1/
DATA CF( 27)/-.3367514367668495D+04/
DATA I( 27)/ 1/,J( 27)/ 3/,K( 27)/ 1/
DATA CF( 28)/-.3367514367668495D+04/
DATA I( 28)/ 1/,J( 28)/ 1/,K( 28)/ 3/
DATA CF( 29)/-.2015701268716309D+04/
DATA I( 29)/ 3/,J( 29)/ 2/,K( 29)/ 0/
DATA CF( 30)/-.2015701268716309D+04/
DATA I( 30)/ 3/,J( 30)/ 0/,K( 30)/ 2/
DATA CF( 31)/-.2015701268716309D+04/
DATA I( 31)/ 2/,J( 31)/ 3/,K( 31)/ 0/
DATA CF( 32)/-.2015701268716309D+04/
DATA I( 32)/ 2/,J( 32)/ 0/,K( 32)/ 3/
DATA CF( 33)/-.2015701268716309D+04/
DATA I( 33)/ 0/,J( 33)/ 3/,K( 33)/ 2/
DATA CF( 34)/-.2015701268716309D+04/
DATA I( 34)/ 0/,J( 34)/ 2/,K( 34)/ 3/
DATA CF( 35)/ .2267829805462643D+04/
DATA I( 35)/ 4/,J( 35)/ 1/,K( 35)/ 0/
DATA CF( 36)/ .2267829805462643D+04/
DATA I( 36)/ 4/,J( 36)/ 0/,K( 36)/ 1/
DATA CF( 37)/ .2267829805462643D+04/
DATA I( 37)/ 1/,J( 37)/ 4/,K( 37)/ 0/
DATA CF( 38)/ .2267829805462643D+04/
DATA I( 38)/ 1/,J( 38)/ 0/,K( 38)/ 4/
DATA CF( 39)/ .2267829805462643D+04/
DATA I( 39)/ 0/,J( 39)/ 4/,K( 39)/ 1/
DATA CF( 40)/ .2267829805462643D+04/
DATA I( 40)/ 0/,J( 40)/ 1/,K( 40)/ 4/
DATA CF( 41)/-.3169171521097543D+05/
DATA I( 41)/ 2/,J( 41)/ 2/,K( 41)/ 2/
DATA CF( 42)/ .1358009279941601D+05/
DATA I( 42)/ 3/,J( 42)/ 2/,K( 42)/ 1/
DATA CF( 43)/ .1358009279941601D+05/
DATA I( 43)/ 3/,J( 43)/ 1/,K( 43)/ 2/
DATA CF( 44)/ .1358009279941601D+05/
DATA I( 44)/ 2/,J( 44)/ 3/,K( 44)/ 1/
DATA CF( 45)/ .1358009279941601D+05/
DATA I( 45)/ 2/,J( 45)/ 1/,K( 45)/ 3/
DATA CF( 46)/ .1358009279941601D+05/
DATA I( 46)/ 1/,J( 46)/ 3/,K( 46)/ 2/
DATA CF( 47)/ .1358009279941601D+05/
DATA I( 47)/ 1/,J( 47)/ 2/,K( 47)/ 3/
DATA CF( 48)/ .1329216151296855D+05/
DATA I( 48)/ 3/,J( 48)/ 3/,K( 48)/ 0/
DATA CF( 49)/ .1329216151296855D+05/
DATA I( 49)/ 3/,J( 49)/ 0/,K( 49)/ 3/
DATA CF( 50)/ .1329216151296855D+05/
DATA I( 50)/ 0/,J( 50)/ 3/,K( 50)/ 3/
DATA CF( 51)/-.8021803662639745D+04/
DATA I( 51)/ 4/,J( 51)/ 1/,K( 51)/ 1/
DATA CF( 52)/-.8021803662639745D+04/
DATA I( 52)/ 1/,J( 52)/ 4/,K( 52)/ 1/
DATA CF( 53)/-.8021803662639745D+04/
DATA I( 53)/ 1/,J( 53)/ 1/,K( 53)/ 4/
DATA CF( 54)/-.6968357580082175D+04/
DATA I( 54)/ 4/,J( 54)/ 2/,K( 54)/ 0/
DATA CF( 55)/-.6968357580082175D+04/
DATA I( 55)/ 4/,J( 55)/ 0/,K( 55)/ 2/
DATA CF( 56)/-.6968357580082175D+04/
DATA I( 56)/ 2/,J( 56)/ 4/,K( 56)/ 0/
DATA CF( 57)/-.6968357580082175D+04/
DATA I( 57)/ 2/,J( 57)/ 0/,K( 57)/ 4/
DATA CF( 58)/-.6968357580082175D+04/
DATA I( 58)/ 0/,J( 58)/ 4/,K( 58)/ 2/
DATA CF( 59)/-.6968357580082175D+04/
DATA I( 59)/ 0/,J( 59)/ 2/,K( 59)/ 4/
DATA CF( 60)/-.9993104153667776D+03/
DATA I( 60)/ 5/,J( 60)/ 1/,K( 60)/ 0/
DATA CF( 61)/-.9993104153667776D+03/
DATA I( 61)/ 5/,J( 61)/ 0/,K( 61)/ 1/
DATA CF( 62)/-.9993104153667776D+03/
DATA I( 62)/ 1/,J( 62)/ 5/,K( 62)/ 0/
DATA CF( 63)/-.9993104153667776D+03/
DATA I( 63)/ 1/,J( 63)/ 0/,K( 63)/ 5/
DATA CF( 64)/-.9993104153667776D+03/
DATA I( 64)/ 0/,J( 64)/ 5/,K( 64)/ 1/
DATA CF( 65)/-.9993104153667776D+03/
DATA I( 65)/ 0/,J( 65)/ 1/,K( 65)/ 5/
DATA CF( 66)/ .5583449663686610D+04/
DATA I( 66)/ 3/,J( 66)/ 2/,K( 66)/ 2/
DATA CF( 67)/ .5583449663686610D+04/
DATA I( 67)/ 2/,J( 67)/ 3/,K( 67)/ 2/
DATA CF( 68)/ .5583449663686610D+04/
DATA I( 68)/ 2/,J( 68)/ 2/,K( 68)/ 3/
DATA CF( 69)/ .9258544110748456D+05/
DATA I( 69)/ 3/,J( 69)/ 3/,K( 69)/ 1/
DATA CF( 70)/ .9258544110748456D+05/
DATA I( 70)/ 3/,J( 70)/ 1/,K( 70)/ 3/
DATA CF( 71)/ .9258544110748456D+05/
DATA I( 71)/ 1/,J( 71)/ 3/,K( 71)/ 3/
DATA CF( 72)/-.8202176930392115D+05/
DATA I( 72)/ 4/,J( 72)/ 2/,K( 72)/ 1/
DATA CF( 73)/-.8202176930392115D+05/
DATA I( 73)/ 4/,J( 73)/ 1/,K( 73)/ 2/
DATA CF( 74)/-.8202176930392115D+05/
DATA I( 74)/ 2/,J( 74)/ 4/,K( 74)/ 1/
DATA CF( 75)/-.8202176930392115D+05/
DATA I( 75)/ 2/,J( 75)/ 1/,K( 75)/ 4/
DATA CF( 76)/-.8202176930392115D+05/
DATA I( 76)/ 1/,J( 76)/ 4/,K( 76)/ 2/
DATA CF( 77)/-.8202176930392115D+05/
DATA I( 77)/ 1/,J( 77)/ 2/,K( 77)/ 4/
DATA CF( 78)/-.9568392519274544D+04/
DATA I( 78)/ 4/,J( 78)/ 3/,K( 78)/ 0/
DATA CF( 79)/-.9568392519274544D+04/
DATA I( 79)/ 4/,J( 79)/ 0/,K( 79)/ 3/
DATA CF( 80)/-.9568392519274544D+04/
DATA I( 80)/ 3/,J( 80)/ 4/,K( 80)/ 0/
DATA CF( 81)/-.9568392519274544D+04/
DATA I( 81)/ 3/,J( 81)/ 0/,K( 81)/ 4/
DATA CF( 82)/-.9568392519274544D+04/
DATA I( 82)/ 0/,J( 82)/ 4/,K( 82)/ 3/
DATA CF( 83)/-.9568392519274544D+04/
DATA I( 83)/ 0/,J( 83)/ 3/,K( 83)/ 4/
DATA CF( 84)/ .6490807101725352D+05/
DATA I( 84)/ 5/,J( 84)/ 1/,K( 84)/ 1/
DATA CF( 85)/ .6490807101725352D+05/
DATA I( 85)/ 1/,J( 85)/ 5/,K( 85)/ 1/
DATA CF( 86)/ .6490807101725352D+05/
DATA I( 86)/ 1/,J( 86)/ 1/,K( 86)/ 5/
DATA CF( 87)/ .1384326017766477D+05/
DATA I( 87)/ 5/,J( 87)/ 2/,K( 87)/ 0/
DATA CF( 88)/ .1384326017766477D+05/
DATA I( 88)/ 5/,J( 88)/ 0/,K( 88)/ 2/
DATA CF( 89)/ .1384326017766477D+05/
DATA I( 89)/ 2/,J( 89)/ 5/,K( 89)/ 0/
DATA CF( 90)/ .1384326017766477D+05/
DATA I( 90)/ 2/,J( 90)/ 0/,K( 90)/ 5/
DATA CF( 91)/ .1384326017766477D+05/
DATA I( 91)/ 0/,J( 91)/ 5/,K( 91)/ 2/
DATA CF( 92)/ .1384326017766477D+05/
DATA I( 92)/ 0/,J( 92)/ 2/,K( 92)/ 5/
DATA CF( 93)/-.6677360564525097D+04/
DATA I( 93)/ 6/,J( 93)/ 1/,K( 93)/ 0/
DATA CF( 94)/-.6677360564525097D+04/
DATA I( 94)/ 6/,J( 94)/ 0/,K( 94)/ 1/
DATA CF( 95)/-.6677360564525097D+04/
DATA I( 95)/ 1/,J( 95)/ 6/,K( 95)/ 0/
DATA CF( 96)/-.6677360564525097D+04/
DATA I( 96)/ 1/,J( 96)/ 0/,K( 96)/ 6/
DATA CF( 97)/-.6677360564525097D+04/
DATA I( 97)/ 0/,J( 97)/ 6/,K( 97)/ 1/
DATA CF( 98)/-.6677360564525097D+04/
DATA I( 98)/ 0/,J( 98)/ 1/,K( 98)/ 6/
DATA XXX(-1)/ 0.D0/,YYY(-1)/ 0.D0/,ZZZ(-1)/ 0.D0/
DATA XXX(0)/ 1.D0/,YYY(0)/ 1.D0/,ZZZ(0)/ 1.D0/
VEX = .1194369889579498D+01
X1 = DEXP(-VEX*X)
Y1 = DEXP(-VEX*Y)
Z1 = DEXP(-VEX*Z)
XX = X*X1
YY = Y*Y1
ZZ = Z*Z1
DO 1 L=1,6
XXX(L) = XX*XXX(L-1)
YYY(L) = YY*YYY(L-1)
ZZZ(L) = ZZ*ZZZ(L-1)
1 CONTINUE
***** CALCULO DEL POTENCIAL ****
ENER = 0.D0
DER1 = 0.D0
DER2 = 0.D0
DER3 = 0.D0
DO 2 L=1,98
ENER=ENER+CF(L)*XXX(I(L))*YYY(J(L))*ZZZ(K(L))
DER1=DER1+CF(L)*I(L)*XXX(I(L)-1)*YYY(J(L))*ZZZ(K(L))
DER2=DER2+CF(L)*J(L)*XXX(I(L))*YYY(J(L)-1)*ZZZ(K(L))
DER3=DER3+CF(L)*K(L)*XXX(I(L))*YYY(J(L))*ZZZ(K(L)-1)
2 CONTINUE
DER(1) = DER1 * (1.D0-VEX*X)*X1
DER(2) = DER2 * (1.D0-VEX*Y)*Y1
DER(3) = DER3 * (1.D0-VEX*Z)*Z1
RETURN
END