mt01aa.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\utils\bief\mt01aa.f
00002 !
00065                      SUBROUTINE MT01AA
00066 !                    *****************
00067 !
00068      &( A11 , A12 , A13 ,
00069      &        A22 , A23 ,
00070      &              A33 ,
00071      &  XMUL,SURFAC,NELEM,NELMAX)
00072 !
00073 !***********************************************************************
00074 ! BIEF   V6P1                                   21/08/2010
00075 !***********************************************************************
00076 !
00077 !
00078 !
00079 !
00080 !
00081 !
00082 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00083 !| A11            |<--| ELEMENTS OF MATRIX
00084 !| A12            |<--| ELEMENTS OF MATRIX
00085 !| A13            |<--| ELEMENTS OF MATRIX
00086 !| A22            |<--| ELEMENTS OF MATRIX
00087 !| A23            |<--| ELEMENTS OF MATRIX
00088 !| A33            |<--| ELEMENTS OF MATRIX
00089 !| NELEM          |-->| NUMBER OF ELEMENTS
00090 !| NELMAX         |-->| MAXIMUM NUMBER OF ELEMENTS
00091 !| SURFAC         |-->| AREA OF TRIANGLES
00092 !| XMUL           |-->| MULTIPLICATION FACTOR
00093 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00094 !
00095       USE BIEF, EX_MT01AA => MT01AA
00096 !
00097       IMPLICIT NONE
00098       INTEGER LNG,LU
00099       COMMON/INFO/LNG,LU
00100 !
00101 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00102 !
00103       INTEGER, INTENT(IN) :: NELEM,NELMAX
00104 !
00105       DOUBLE PRECISION, INTENT(INOUT) :: A11(*),A12(*),A13(*)
00106       DOUBLE PRECISION, INTENT(INOUT) ::        A22(*),A23(*)
00107       DOUBLE PRECISION, INTENT(INOUT) ::               A33(*)
00108 !
00109       DOUBLE PRECISION, INTENT(IN) :: XMUL
00110       DOUBLE PRECISION, INTENT(IN) :: SURFAC(NELMAX)
00111 !
00112 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00113 !
00114 !     DECLARATIONS SPECIFIC TO THIS SUBROUTINE
00115 !
00116       INTEGER IELEM
00117       DOUBLE PRECISION SUR12,DET
00118 !
00119 !=======================================================================
00120 !
00121       SUR12 = XMUL/12.D0
00122 !
00123 !-----------------------------------------------------------------------
00124 !
00125 !   LOOP ON THE ELEMENTS
00126 !
00127       DO IELEM = 1 , NELEM
00128 !
00129       DET = SURFAC(IELEM) * SUR12
00130 !
00131 !  ELEMENTS OFF THE DIAGONAL
00132 !
00133       A12(IELEM) = DET
00134       A13(IELEM) = DET
00135       A23(IELEM) = DET
00136 !
00137 !  DIAGONAL TERMS
00138 !
00139       A11(IELEM) = DET + DET
00140       A22(IELEM) = DET + DET
00141       A33(IELEM) = DET + DET
00142 !
00143 !   END OF THE LOOP ON THE ELEMENTS
00144 !
00145       ENDDO ! IELEM
00146 !
00147 !-----------------------------------------------------------------------
00148 !
00149       RETURN
00150       END

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