ifabtom.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\tomawac\ifabtom.f
00002 !
00039                      SUBROUTINE IFABTOM
00040 !                    ******************
00041 !
00042      &(IFABOR,NELEM2,NETAGE)
00043 !
00044 !***********************************************************************
00045 ! TELEMAC3D   V6P3                                   21/08/2010
00046 !***********************************************************************
00047 !
00048 !
00049 !
00050 !
00051 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00052 !| IFABOR         |<->| CORRESPONDENCE BOUNDARY FACE - 2D ELEMENT 2D
00053 !| NELEM2         |-->| NUMBER OF ELEMENTS IN 2D
00054 !| NETAGE         |-->| NUMBER OF PLANES - 1
00055 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00056 !
00057       USE BIEF
00058 !
00059       IMPLICIT NONE
00060       INTEGER LNG,LU
00061       COMMON/INFO/LNG,LU
00062 !
00063 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00064 !
00065       INTEGER, INTENT(IN)          :: NELEM2,NETAGE
00066       INTEGER, INTENT(INOUT)       :: IFABOR(NELEM2,5,NETAGE)
00067 !
00068 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00069 !
00070       INTEGER IELEM2,IETAGE,LOC
00071 !
00072 !=======================================================================
00073 !  IFABOR ARRAYS
00074 !=======================================================================
00075 !
00076 !     SEE SCHAR41 FOR THE MEANING OF LOC
00077 !     1 IS : DO NOT RECOMPUTE VELOCITIES WHEN CROSSING A MESH PLANE
00078 !     2 IS : RECOMPUTE VELOCITIES WHEN CROSSING A MESH PLANE
00079 !
00080       LOC = 1
00081 !
00082 !  IFABOR ON THE VERTICAL FACES OF THE ELEMENTS OF THE FIRST LAYER
00083 !
00084       DO IELEM2 = 1,NELEM2
00085 !       ALREADY DONE BY 2D
00086 !       IFABOR(IELEM2,1,IETAGE) = IFABOR(IELEM2,1,1)
00087 !       IFABOR(IELEM2,2,IETAGE) = IFABOR(IELEM2,2,1)
00088 !       IFABOR(IELEM2,3,IETAGE) = IFABOR(IELEM2,3,1)
00089         IFABOR(IELEM2,4,1) = LOC
00090         IFABOR(IELEM2,5,1) = LOC
00091       ENDDO
00092       IF(NETAGE.GE.2) THEN
00093         DO IETAGE = 2,NETAGE
00094           DO IELEM2 = 1,NELEM2
00095             IFABOR(IELEM2,1,IETAGE) = IFABOR(IELEM2,1,1)
00096             IFABOR(IELEM2,2,IETAGE) = IFABOR(IELEM2,2,1)
00097             IFABOR(IELEM2,3,IETAGE) = IFABOR(IELEM2,3,1)
00098             IFABOR(IELEM2,4,IETAGE) = LOC
00099             IFABOR(IELEM2,5,IETAGE) = LOC
00100           ENDDO
00101         ENDDO
00102       ENDIF
00103 !
00104 !-----------------------------------------------------------------------
00105 !
00106 !     IFABOR ON THE BOTTOM AND SURFACE
00107 !
00108       DO IELEM2 = 1,NELEM2
00109 !       -1 : SOLID
00110 !        0 : LIQUID
00111 !        1 : PERIODICITY (TREATED AS INTERNAL BOUNDARY)
00112         IFABOR(IELEM2,4,1)      = LOC
00113         IFABOR(IELEM2,5,NETAGE) = LOC
00114 !
00115       ENDDO
00116 !
00117 !-----------------------------------------------------------------------
00118 !
00119       RETURN
00120       END

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