source: trunk/fsource/spsubs/sgrmat.for @ 211

Last change on this file since 211 was 211, checked in by vondreele, 13 years ago

unfix includes

File size: 1.4 KB
Line 
1      SUBROUTINE SGRMAT(IOP,RT,N,A,B,C,D,E,F,G,H,O)
2
3!Purpose:      S.R. to create a 3,3 matrix from 9 scalers
4
5      INCLUDE       '../INCLDS/COPYRIGT.FOR' 
6
7C       THIS PROGRAM WAS DEVELOPED FOR
8C                    THE DIVISION OF CHEMISTRY
9C                               OF
10C               THE NATIONAL RESEARCH COUNCIL OF CANADA
11C                               BY
12C       ALLEN C. LARSON, 14 CERRADO LOOP, SANTA FE, NM 87505, USA
13
14!Calling sequence parameters:
15
16      INTEGER*4     IOP                 !Matrix generator count
17      REAL*4        RT(5,4,25)          !Matrix to be generated
18      INTEGER*4     N                   !Number of the matrix to be generated
19      REAL*4        A,B,C,D,E,F,G,H,O   !Matrix terms
20
21!Local varaibles:
22
23!Code:
24
25      RT(1,1,N) = A
26      RT(1,2,N) = B
27      RT(1,3,N) = C
28      RT(1,4,N) = 0.0
29      RT(2,1,N) = D
30      RT(2,2,N) = E
31      RT(2,3,N) = F
32      RT(2,4,N) = 0.0
33      RT(3,1,N) = G
34      RT(3,2,N) = H
35      RT(3,3,N) = O
36      RT(3,4,N) = 0.0
37      RT(4,1,N) = 0.0
38      RT(4,2,N) = 0.0
39      RT(4,3,N) = 0.0
40      RT(4,4,N) = 1.0
41      RT(5,1,N) = 81*(2*RT(1,1,N)+3*RT(1,2,N)+4*RT(1,3,N))
42     1  +9*(2*RT(2,1,N)+3*RT(2,2,N)+4*RT(2,3,N))
43     1  +2*RT(3,1,N)+3*RT(3,2,N)+4*RT(3,3,N)
44      RT(5,2,N) = 0.0                                          !Clear the translation info
45      RT(5,3,N) = IOP
46      RT(5,4,N) = 20.
47
48      RETURN
49      END
Note: See TracBrowser for help on using the repository browser.