mt01cc.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\utils\bief\mt01cc.f
00002 !
00049                      SUBROUTINE MT01CC
00050 !                    *****************
00051 !
00052      &( A11 , A12 , A13 , A14 , A15 , A16 ,
00053      &        A22 , A23 , A24 , A25 , A26 ,
00054      &              A33 , A34 , A35 , A36 ,
00055      &                    A44 , A45 , A46 ,
00056      &                          A55 , A56 ,
00057      &                                A66 ,
00058      &  XMUL,SURFAC,NELEM,NELMAX )
00059 !
00060 !***********************************************************************
00061 ! BIEF   V6P1                                   21/08/2010
00062 !***********************************************************************
00063 !
00064 !
00065 !
00066 !
00067 !
00068 !
00069 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00070 !| A11            |<--| ELEMENTS OF MATRIX
00071 !| A12            |<--| ELEMENTS OF MATRIX
00072 !| A13            |<--| ELEMENTS OF MATRIX
00073 !| A14            |<--| ELEMENTS OF MATRIX
00074 !| A15            |<--| ELEMENTS OF MATRIX
00075 !| A16            |<--| ELEMENTS OF MATRIX
00076 !| A22            |<--| ELEMENTS OF MATRIX
00077 !| A23            |<--| ELEMENTS OF MATRIX
00078 !| A24            |<--| ELEMENTS OF MATRIX
00079 !| A25            |<--| ELEMENTS OF MATRIX
00080 !| A26            |<--| ELEMENTS OF MATRIX
00081 !| A33            |<--| ELEMENTS OF MATRIX
00082 !| A34            |<--| ELEMENTS OF MATRIX
00083 !| A35            |<--| ELEMENTS OF MATRIX
00084 !| A36            |<--| ELEMENTS OF MATRIX
00085 !| A44            |<--| ELEMENTS OF MATRIX
00086 !| A45            |<--| ELEMENTS OF MATRIX
00087 !| A46            |<--| ELEMENTS OF MATRIX
00088 !| A55            |<--| ELEMENTS OF MATRIX
00089 !| A56            |<--| ELEMENTS OF MATRIX
00090 !| A66            |<--| ELEMENTS OF MATRIX
00091 !| NELEM          |-->| NUMBER OF ELEMENTS
00092 !| NELMAX         |-->| MAXIMUM NUMBER OF ELEMENTS
00093 !| SURFAC         |-->| AREA OF TRIANGLES
00094 !| XMUL           |-->| MULTIPLICATION FACTOR
00095 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00096 !
00097       USE BIEF !, EX_MT01CC => MT01CC
00098 !
00099       IMPLICIT NONE
00100       INTEGER LNG,LU
00101       COMMON/INFO/LNG,LU
00102 !
00103 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00104 !
00105       INTEGER, INTENT(IN) :: NELEM,NELMAX
00106 !
00107       DOUBLE PRECISION, INTENT(INOUT) :: A11(*),A12(*),A13(*),A14(*)
00108       DOUBLE PRECISION, INTENT(INOUT) :: A15(*),A16(*)
00109       DOUBLE PRECISION, INTENT(INOUT) ::        A22(*),A23(*),A24(*)
00110       DOUBLE PRECISION, INTENT(INOUT) :: A25(*),A26(*)
00111       DOUBLE PRECISION, INTENT(INOUT) :: A33(*),A34(*),A35(*),A36(*)
00112       DOUBLE PRECISION, INTENT(INOUT) :: A44(*),A45(*),A46(*)
00113       DOUBLE PRECISION, INTENT(INOUT) :: A55(*),A56(*),A66(*)
00114 !
00115       DOUBLE PRECISION, INTENT(IN) :: XMUL
00116 !
00117       DOUBLE PRECISION, INTENT(IN) :: SURFAC(NELMAX)
00118 !
00119 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00120 !
00121 !     DECLARATIONS SPECIFIC TO THIS SUBROUTINE
00122 !
00123       INTEGER IELEM
00124       DOUBLE PRECISION XSUR180
00125 !
00126 !=======================================================================
00127 !
00128       XSUR180 = XMUL / 180.D0
00129 !
00130       DO IELEM = 1 , NELEM
00131 !
00132 !  EXTRADIAGONAL TERMS
00133 !
00134         A12(IELEM) = - SURFAC(IELEM) * XSUR180
00135         A13(IELEM) =          A12(IELEM)
00136         A14(IELEM) =   0.D0
00137         A15(IELEM) =   4.D0 * A12(IELEM)
00138         A16(IELEM) =   0.D0
00139         A23(IELEM) =          A12(IELEM)
00140         A24(IELEM) =   0.D0
00141         A25(IELEM) =   0.D0
00142         A26(IELEM) =          A15(IELEM)
00143         A34(IELEM) =          A15(IELEM)
00144         A35(IELEM) =   0.D0
00145         A36(IELEM) =   0.D0
00146         A45(IELEM) = - 4.D0 * A15(IELEM)
00147         A46(IELEM) =          A45(IELEM)
00148         A56(IELEM) =          A45(IELEM)
00149 !
00150 !  DIAGONAL TERMS
00151 !
00152         A11(IELEM) = - 6.D0 * A12(IELEM)
00153         A22(IELEM) =          A11(IELEM)
00154         A33(IELEM) =          A11(IELEM)
00155         A44(IELEM) = - 8.D0 * A15(IELEM)
00156         A55(IELEM) =          A44(IELEM)
00157         A66(IELEM) =          A44(IELEM)
00158 !
00159       ENDDO
00160 !
00161 !-----------------------------------------------------------------------
00162 !
00163       RETURN
00164       END

Generated on Fri Aug 31 2013 18:12:58 by S.E.Bourban (HRW) using doxygen 1.7.0