smtrac.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac2d\smtrac.f
00002 !
00046                      SUBROUTINE SMTRAC
00047 !                    *****************
00048 !
00049      &(NPOIN,DIMT,AT,DT,SMTR,SMH,NREJET,ISCE,TSCE2,MAXSCE,MAXTRA,ITRAC)
00050 !
00051 !***********************************************************************
00052 ! TELEMAC2D   V6P1                                   21/08/2010
00053 !***********************************************************************
00054 !
00055 !
00056 !
00057 !
00058 !
00059 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00060 !| AT             |-->| TIME
00061 !| DIMT           |-->| DIMENSION OF THE TRACER
00062 !| DT             |-->| HYDRO TIME STEP
00063 !| ISCE           |-->| GLOBAL INDICES OF SOURCE POINTS
00064 !| ITRAC          |-->| TRCER INDEX
00065 !| MAXSCE         |-->| MAXIMUM NUMBER OF SOURCES
00066 !| MAXTRA         |-->| MAXIMUM NUMER OF TRACERS
00067 !| NPOIN          |-->| TOTAL NUMBER OF NODES IN THE MESH
00068 !| NREJET         |-->| NUMBER OF SOURCE/SINK
00069 !| SMH            |-->| SOURCE TERMS FOR CONTINUITY EQUATION
00070 !| SMTR           |-->| SOURCE TERMS FOR TRACER
00071 !| TSCE2          |-->| VALUES OF TRACERS AT SOURCES
00072 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00073 !
00074       IMPLICIT NONE
00075       INTEGER LNG,LU
00076       COMMON/INFO/LNG,LU
00077 !
00078 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00079 !
00080       INTEGER, INTENT(IN) :: NPOIN,NREJET,ISCE(*),DIMT,ITRAC
00081       INTEGER, INTENT(IN) :: MAXSCE,MAXTRA
00082       DOUBLE PRECISION, INTENT(IN)    :: AT,DT,SMH(NPOIN)
00083       DOUBLE PRECISION, INTENT(IN)    :: TSCE2(MAXSCE,MAXTRA)
00084       DOUBLE PRECISION, INTENT(INOUT) :: SMTR(DIMT)
00085 !
00086 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00087 !
00088       INTEGER I,IS
00089 !
00090 !-----------------------------------------------------------------------
00091 !
00092       IF(NREJET.NE.0) THEN
00093         DO I=1,NREJET
00094           IS =ISCE(I)
00095           IF(IS.GT.0) THEN
00096             SMTR(IS) = SMTR(IS) + DT*SMH(IS) * TSCE2(I,ITRAC)
00097           ENDIF
00098         ENDDO
00099       ENDIF
00100 !
00101 !-----------------------------------------------------------------------
00102 !
00103       RETURN
00104       END

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