gettri.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\utils\bief\gettri.f
00002 !
00066                      SUBROUTINE GETTRI
00067 !                    *****************
00068 !
00069      &(M,MDIFF,TETA,MESH3D,NPLAN,NPOIN2,NSEG2D,IELM3,NELEM2)
00070 !
00071 !***********************************************************************
00072 ! BIEF   V6P2                                   21/08/2010
00073 !***********************************************************************
00074 !
00075 !
00076 !
00077 !
00078 !
00079 !
00080 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00081 !| IELM3          |<--| TYPE OF ELEMENT
00082 !| M              |<--| TRIDIAGONAL MATRIX
00083 !| MDIFF          |-->| ORIGINAL DIFFUSION MATRIX
00084 !| MESH3D         |-->| 3D MESH STRUCTURE
00085 !| NELEM2         |-->| NUMBER OF TRIANGLES OF ORIGINAL 2D MESH
00086 !| NPLAN          |-->| NUMBER OF PLANES
00087 !| NPOIN2         |-->| NUMBER OF POINTS OF 2D MESH
00088 !| NSEG2D         |-->| NUMBER OF SEGMENTS IN 2D
00089 !| TETA           |-->| COEFFICIENT USED IN THE RESULT
00090 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00091 !
00092       USE BIEF, EX_GETTRI => GETTRI
00093 !
00094       IMPLICIT NONE
00095       INTEGER LNG,LU
00096       COMMON/INFO/LNG,LU
00097 !
00098 !-----------------------------------------------------------------------
00099 !
00100       INTEGER, INTENT(IN) :: NPLAN,NPOIN2,NSEG2D,IELM3,NELEM2
00101 !
00102       DOUBLE PRECISION, INTENT(IN)    :: TETA
00103       DOUBLE PRECISION, INTENT(INOUT) :: M(NPOIN2*NPLAN,*)
00104 !
00105       TYPE(BIEF_OBJ),  INTENT(INOUT)  :: MDIFF
00106       TYPE(BIEF_MESH), INTENT(IN)     :: MESH3D
00107 !
00108 !-----------------------------------------------------------------------
00109 !
00110       IF(MDIFF%STO.EQ.1) THEN
00111 !
00112         CALL GETTRIEBE(M,MDIFF%D%R,MDIFF%X%R,TETA,
00113      &                 MESH3D%IKLE%I,MESH3D%NPOIN,MESH3D%NELEM,
00114      &                 MESH3D%NELMAX,MESH3D,IELM3,NELEM2,NPLAN,
00115      &                 MESH3D%KNOLG%I)
00116 !
00117       ELSEIF(MDIFF%STO.EQ.3) THEN
00118 !
00119         CALL GETTRISEG(M,MDIFF%D%R,MDIFF%X%R,TETA,
00120      &                 MESH3D%NPOIN,MESH3D,
00121      &                 MESH3D%NSEG,NSEG2D,NPLAN,NPOIN2,IELM3)
00122 !
00123       ELSE
00124 !
00125         WRITE(LU,*) 'UNKNOWN STORAGE FOR MDIFF IN GETTRI'
00126         CALL PLANTE(1)
00127         STOP
00128 !
00129       ENDIF
00130 !
00131 !-----------------------------------------------------------------------
00132 !
00133       RETURN
00134       END

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