actuzf.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac3d\actuzf.f
00002 !
00049                      SUBROUTINE ACTUZF
00050 !                    *****************
00051 !
00052      & ( IVIDE , EPAI , ZF , NPOIN2, NPFMAX , NPF )
00053 !
00054 !***********************************************************************
00055 ! TELEMAC3D   V6P1                                   21/08/2010
00056 !***********************************************************************
00057 !
00058 !
00059 !
00060 !
00061 !
00062 !
00063 !
00064 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00065 !| EPAI           |-->| THICKNESS OF MESH ELEMENTS DISCRETISING THE BED
00066 !| IVIDE          |-->| INDEX OF EMPTY SPACES AT MESH POINTS
00067 !| NPF            |-->| NUMBER OF POINTS OF THE BOTTOM ON ONE VERTICAL
00068 !| NPFMAX         |-->| MAXIMUM NUMBER OF HORIZONTAL PLANES THAT
00069 !|                |   | DISCRETISE MUDDY BOTTOM
00070 !| NPOIN2         |-->| NUMBER OF POINTS IN 2D
00071 !| ZF             |<->| ELEVATION OF BOTTOM
00072 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00073 !
00074       IMPLICIT NONE
00075 !
00076 !-----------------------------------------------------------------------
00077 !
00078       INTEGER, INTENT(IN)             :: NPOIN2, NPFMAX
00079       DOUBLE PRECISION, INTENT(IN)    :: IVIDE(NPFMAX,NPOIN2)
00080       DOUBLE PRECISION, INTENT(IN)    :: EPAI(NPFMAX-1,NPOIN2)
00081       DOUBLE PRECISION, INTENT(INOUT) :: ZF(NPOIN2)
00082       INTEGER, INTENT(IN)             :: NPF(NPOIN2)
00083 !
00084 !-----------------------------------------------------------------------
00085 !
00086       INTEGER IPOIN , IPF
00087       DOUBLE PRECISION  ECOUCH
00088 !
00089 !-----------------------------------------------------------------------
00090 !
00091       DO IPOIN=1,NPOIN2
00092 !
00093 !       -----COMPUTES THE BOTTOM ELEVATION-----
00094 !
00095         DO IPF=1,NPF(IPOIN)-1
00096           ECOUCH=(IVIDE(IPF,IPOIN)+IVIDE(IPF+1,IPOIN))/2.D0
00097           ZF(IPOIN)=ZF(IPOIN)+(1.D0+ECOUCH)*EPAI(IPF,IPOIN)
00098         END DO
00099 !
00100       END DO
00101 !
00102 !-----------------------------------------------------------------------
00103 !
00104       RETURN
00105       END SUBROUTINE ACTUZF

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