dricv.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac3d\dricv.f
00002 !
00050                      SUBROUTINE DRICV
00051 !                    ****************
00052 !
00053      & (FRI,FRT,RI,NPOIN3)
00054 !
00055 !***********************************************************************
00056 ! TELEMAC3D   V6P1                                   21/08/2010
00057 !***********************************************************************
00058 !
00059 !
00060 !reference  JOURNAL OF MARINE RESEARCH VOLUME 1  1948.
00061 !
00062 !
00063 !
00064 !
00065 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00066 !| FRI            |<->| DAMPING FUNCTION FOR VELOCITIES
00067 !| FRT            |<->| DAMPING FUNCTION FOR TRACERS
00068 !| NPOIN3         |-->| NUMBER OF 3D POINTS
00069 !| RI             |-->| RICHARDSON NUMBER
00070 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00071 !
00072       IMPLICIT NONE
00073       INTEGER LNG,LU
00074       COMMON/INFO/LNG,LU
00075 !
00076 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00077 !
00078       INTEGER, INTENT(IN)             :: NPOIN3
00079       DOUBLE PRECISION, INTENT(INOUT) :: FRI(NPOIN3),FRT(NPOIN3)
00080       DOUBLE PRECISION, INTENT(IN)    :: RI(NPOIN3)
00081 !
00082 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00083 !
00084       INTEGER I
00085       DOUBLE PRECISION A,B,EPS
00086 !
00087 !-----------------------------------------------------------------------
00088 !
00089 ! DAMPING FUNCTION FOR VELOCITIES: (1+A*RI)**B
00090 !
00091       A=10.D0
00092       B=-0.5D0
00093       EPS=1.D-8
00094 !
00095       DO I=1,NPOIN3
00096         IF(RI(I).GT.EPS) THEN
00097           FRI(I)=(1.D0+A*RI(I))**B
00098         ELSE
00099           FRI(I)=1.D0
00100         ENDIF
00101       ENDDO
00102 !
00103 !-----------------------------------------------------------------------
00104 !
00105 ! DAMPING FUNCTION FOR TRACERS: (1+A*RI)**B
00106 !
00107       A=3.33D0
00108       B=-1.5D0
00109       EPS=1.D-8
00110 !
00111       DO I=1,NPOIN3
00112         IF(RI(I).GT.EPS) THEN
00113           FRT(I)=(1.D0+A*RI(I))**B
00114         ELSE
00115           FRT(I)=1.D0
00116         ENDIF
00117       ENDDO
00118 !
00119 !-----------------------------------------------------------------------
00120 !
00121       RETURN
00122       END

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