DATABASES/madrigal Files · source/madf/geolib/gmetbv.f

Fix bug in plotting

jespinoza - - Load All Authors

File last commit:

r0:b84e1135c2c4


                r21:781d2d915c68

Download file

             gmetbv.f
        
                    37 lines
            
             | 907 B
            
                | text/x-fortran
            
             |
                FortranFixedLexer
            
             / source / madf / geolib / gmetbv.f
          
                    History
                
                 |
                  Annotation
                 | Raw
                 |Copy content
                 |Copy permalink

      C     $Id: gmetbv.f 5093 2015-06-10 13:48:23Z brideout $

      C

            SUBROUTINE GMETBV(E1,E2,E3,G)

      C

      C     jmh - 5/06  ans fortran 66

      C

      C     GMETBV calculates the metric tensor G of a coordinate system given

      C     base vecors e1,e2,e3. See Margenau and Murphy, 5-86.

      C

      C       Input:

      C         E1,E2,E3 - base vectors

      C

      C      Output:

      C         G - metric tensor.

      C     $Id: gmetbv.f 5093 2015-06-10 13:48:23Z brideout $

      C

      C     .. Array Arguments ..

            DOUBLE PRECISION E1(3),E2(3),E3(3),G(3,3)

      C     ..

      C     .. Local Scalars ..

      C     ..

      C     .. External Functions ..

            DOUBLE PRECISION SPROD

            EXTERNAL SPROD

      C     ..

            G(1,1) = SPROD(E1,E1)

            G(2,1) = SPROD(E2,E1)

            G(3,1) = SPROD(E3,E1)

            G(1,2) = SPROD(E1,E2)

            G(2,2) = SPROD(E2,E2)

            G(3,2) = SPROD(E3,E2)

            G(1,3) = SPROD(E1,E3)

            G(2,3) = SPROD(E2,E3)

            G(3,3) = SPROD(E3,E3)

            RETURN

      C

            END

	Site-wide shortcuts
/	Use quick search box
g h	Goto home page
g g	Goto my private gists page
g G	Goto my public gists page
g 0-9	Goto bookmarked items from 0-9
n r	New repository page
n g	New gist page

	Repositories
g s	Goto summary page
g c	Goto changelog page
g f	Goto files page
g F	Goto files page with file search activated
g p	Goto pull requests page
g o	Goto repository settings
g O	Goto repository access permissions settings
t s	Toggle sidebar on some pages

				C $Id: gmetbv.f 5093 2015-06-10 13:48:23Z brideout $
				C
				SUBROUTINE GMETBV(E1,E2,E3,G)
				C
				C jmh - 5/06 ans fortran 66
				C
				C GMETBV calculates the metric tensor G of a coordinate system given
				C base vecors e1,e2,e3. See Margenau and Murphy, 5-86.
				C
				C Input:
				C E1,E2,E3 - base vectors
				C
				C Output:
				C G - metric tensor.
				C $Id: gmetbv.f 5093 2015-06-10 13:48:23Z brideout $
				C
				C .. Array Arguments ..
				DOUBLE PRECISION E1(3),E2(3),E3(3),G(3,3)
				C ..
				C .. Local Scalars ..
				C ..
				C .. External Functions ..
				DOUBLE PRECISION SPROD
				EXTERNAL SPROD
				C ..
				G(1,1) = SPROD(E1,E1)
				G(2,1) = SPROD(E2,E1)
				G(3,1) = SPROD(E3,E1)
				G(1,2) = SPROD(E1,E2)
				G(2,2) = SPROD(E2,E2)
				G(3,2) = SPROD(E3,E2)
				G(1,3) = SPROD(E1,E3)
				G(2,3) = SPROD(E2,E3)
				G(3,3) = SPROD(E3,E3)
				RETURN
				C
				END