reinit.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac2d\reinit.f
00002 !
00045                      SUBROUTINE REINIT
00046 !                    *****************
00047 !
00048      &(NS,NSEG,NPTFR,H,SMTR,HSTOK,HC,HCSTOK,FLUXT,FLUHBOR,DTT,NTRAC)
00049 !
00050 !***********************************************************************
00051 ! TELEMAC2D   V6P1                                   21/08/2010
00052 !***********************************************************************
00053 !
00054 !
00055 !
00056 !
00057 !
00058 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00059 !| DTT            |<--| TIME STEP FOR TRACER
00060 !| FLUHBOR        |<--| TRACER FLUX AT HE BOUNDARY (UPDATED)
00061 !| FLUXT          |<--| TRACER FLUX (UPDATED)
00062 !| H              |-->| WATER DEPTH
00063 !| HC             |-->| RECONSRUCTED H FOR ORDRE 2 (CORRECTED)
00064 !| HCSTOK         |<--| RECONSTRUCTED H FOR ORDRE 2(CORRECTED AND STOCKE)
00065 !| HSTOK          |<--| STOCKED WATER DEPTHS
00066 !| NPTFR          |-->| NUMBER OF BOUNDARY POINTS
00067 !| NS             |-->| TOTAL NUMER OF POINTS IN THE MESH
00068 !| NSEG           |-->| NUMBER OF EDGES IN THE MESH
00069 !| NTRAC          |-->| NUMBER OF TRACERS
00070 !| SMTR           |<--| TRACER SOURCE TERMS
00071 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00072 !
00073       USE BIEF
00074 !
00075       IMPLICIT NONE
00076 !
00077 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00078 !
00079       INTEGER, INTENT(IN)             :: NS,NSEG,NPTFR,NTRAC
00080       DOUBLE PRECISION, INTENT(INOUT) :: DTT
00081       DOUBLE PRECISION, INTENT(INOUT) :: HSTOK(*),HCSTOK(2,*)
00082       DOUBLE PRECISION, INTENT(IN)    :: H(*),HC(2,*)
00083       TYPE(BIEF_OBJ), INTENT(INOUT)   :: SMTR,FLUXT,FLUHBOR
00084 !
00085 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00086 !
00087       INTEGER IS,NSG,ITRAC
00088 !
00089       DTT = 0.D0
00090 !
00091       DO IS=1,NS
00092         HSTOK(IS)=H(IS)
00093       ENDDO
00094 !
00095       DO ITRAC=1,NTRAC
00096         DO IS=1,NS
00097           SMTR%ADR(ITRAC)%P%R(IS)=0.D0
00098         ENDDO
00099       ENDDO
00100 !
00101       DO NSG=1,NSEG
00102         HCSTOK(1,NSG) = HC(1,NSG)
00103         HCSTOK(2,NSG) = HC(2,NSG)
00104       ENDDO
00105 !
00106       DO ITRAC=1,NTRAC
00107         DO NSG=1,NSEG
00108           FLUXT%ADR(ITRAC)%P%R(NSG)=0.D0
00109         ENDDO
00110       ENDDO
00111 !
00112       DO ITRAC=1,NTRAC
00113         DO IS=1,NPTFR
00114           FLUHBOR%ADR(ITRAC)%P%R(IS)=0.D0
00115         ENDDO
00116       ENDDO
00117 !
00118 !-----------------------------------------------------------------------
00119 !
00120       RETURN
00121       END

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