limite.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\tomawac\limite.f
00002 !
00046                         SUBROUTINE LIMITE
00047 !                       *****************
00048      &( F     , DEPTH , FREQ  , NPOIN2, NPLAN , NF    )
00049 !
00050 !***********************************************************************
00051 ! TOMAWAC   V7P0                                 30/07/2014
00052 !***********************************************************************
00053 !
00054 !
00055 !
00056 !
00057 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~!
00058 !| DEPTH          |-->| WATER DEPTH
00059 !| F              |-->| DIRECTIONAL SPECTRUM
00060 !| FREQ           |-->| DISCRETIZED FREQUENCIES
00061 !| NF             |-->| NUMBER OF FREQUENCIES
00062 !| NPLAN          |-->| NUMBER OF DIRECTIONS
00063 !| NPOIN2         |-->| NUMBER OF POINTS IN 2D MESH
00064 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~!
00065 !
00066 !
00067       USE DECLARATIONS_TOMAWAC, ONLY : DEUPI,GRAVIT
00068 !
00069       IMPLICIT NONE
00070 !
00071       INTEGER LNG,LU
00072       COMMON/INFO/ LNG,LU
00073 !
00074 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00075 !
00076       INTEGER, INTENT(IN)             :: NF,NPLAN,NPOIN2
00077       DOUBLE PRECISION, INTENT(IN)    :: FREQ(NF),DEPTH(NPOIN2)
00078       DOUBLE PRECISION, INTENT(INOUT) :: F(NPOIN2,NPLAN,NF)
00079 !
00080 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00081 !
00082       INTEGER          IP    , JP    , JF
00083       DOUBLE PRECISION EDEF  , COEF  , DTETAR, EMAX  , REDUC
00084 !
00085 !-----------------------------------------------------------------------
00086 !
00087       DTETAR=DEUPI/DBLE(NPLAN)
00088       COEF=0.0081D0*GRAVIT**2/DEUPI**4
00089 !
00090       DO IP=1,NPOIN2
00091         DO JF=1,NF
00092           EDEF=0.0D0
00093           DO JP=1,NPLAN
00094             EDEF=EDEF+F(IP,JP,JF)
00095           ENDDO
00096           EDEF=EDEF*DTETAR
00097           EMAX=COEF/FREQ(JF)**5
00098           IF (EDEF.GT.EMAX) THEN
00099             REDUC=EMAX/EDEF
00100             DO JP=1,NPLAN
00101               F(IP,JP,JF)=F(IP,JP,JF)*REDUC
00102             ENDDO
00103           ENDIF
00104         ENDDO
00105       ENDDO
00106 !
00107       RETURN
00108       END

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