complete_ebe_fluint.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac3d\complete_ebe_fluint.f
00002 !
00026                      SUBROUTINE COMPLETE_EBE_FLUINT
00027 !                    ******************************
00028 !
00029      &(EBE_FLUINT,NELEM2,NPLAN)
00030 !
00031 !***********************************************************************
00032 ! TELEMAC3D
00033 !***********************************************************************
00034 !
00035 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00036 !| EBE_FLUINT     |<->| INTERIOR FLUXES TO BE COMPLETED
00037 !| NELEM2         |-->| NUMBER OF 2D ELEMENTS
00038 !| NPLAN          |-->| NUMBER OF PLANES ON THE VERTICAL
00039 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00040 !
00041       USE BIEF
00042 !
00043       IMPLICIT NONE
00044       INTEGER LNG,LU
00045       COMMON/INFO/LNG,LU
00046 !
00047 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00048 !
00049       INTEGER, INTENT(IN)             :: NELEM2,NPLAN
00050       DOUBLE PRECISION, INTENT(INOUT) :: EBE_FLUINT(NELEM2,NPLAN-1,6)
00051 !
00052 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00053 !
00054       INTEGER ILAYER,NLAYER,IELEM2
00055 !
00056 !-----------------------------------------------------------------------
00057 !
00058       NLAYER=NPLAN-1
00059 !
00060 !-----------------------------------------------------------------------
00061 !
00062 !     UPPER LAYER: INITIALISING WITH COEFFICIENTS OF LOWER PLANE
00063 !
00064       DO IELEM2=1,NELEM2
00065         EBE_FLUINT(IELEM2,NLAYER,4)=-EBE_FLUINT(IELEM2,NLAYER,1)
00066         EBE_FLUINT(IELEM2,NLAYER,5)=-EBE_FLUINT(IELEM2,NLAYER,2)
00067         EBE_FLUINT(IELEM2,NLAYER,6)=-EBE_FLUINT(IELEM2,NLAYER,3)
00068       ENDDO
00069 !
00070 !     OTHER LAYERS
00071 !
00072       IF(NLAYER.GE.2) THEN
00073         DO ILAYER=1,NLAYER-1
00074           DO IELEM2=1,NELEM2
00075             EBE_FLUINT(IELEM2,NLAYER,4)=EBE_FLUINT(IELEM2,NLAYER,4)
00076      &                                 -EBE_FLUINT(IELEM2,ILAYER,1)
00077      &                                 -EBE_FLUINT(IELEM2,ILAYER,4)
00078             EBE_FLUINT(IELEM2,NLAYER,5)=EBE_FLUINT(IELEM2,NLAYER,5)
00079      &                                 -EBE_FLUINT(IELEM2,ILAYER,2)
00080      &                                 -EBE_FLUINT(IELEM2,ILAYER,5)
00081             EBE_FLUINT(IELEM2,NLAYER,6)=EBE_FLUINT(IELEM2,NLAYER,6)
00082      &                                 -EBE_FLUINT(IELEM2,ILAYER,3)
00083      &                                 -EBE_FLUINT(IELEM2,ILAYER,6)
00084           ENDDO
00085         ENDDO
00086       ENDIF
00087 !
00088 !-----------------------------------------------------------------------
00089 !
00090       RETURN
00091       END

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