cqueue.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE cqueue
!                   *****************
!
     &( jfre  , jbis  , coef1 )
!
!***********************************************************************
! TOMAWAC   V6P1                                   14/06/2011
!***********************************************************************
!
!brief    ADJUSTS FREQUENCY INDICES AND COMPUTES TAIL
!
!note     THE SPECTRUM IS ASSUMED TO BE 0 FOR FREQUENCIES LOWER THAN
!+          THE FIRST DISCRETISED FREQUENCY.
!note   BEYOND THE LAST DISCRETISED FREQUENCY THE SPECTRUM IS
!+          ASSUMED TO DECREASE FOLLOWING A FREQ**(-TAILF) LAW.
!
!history  M. BENOIT
!+        26/06/96
!+        V1P2
!+
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        13/07/2010
!+        V6P0
!+   Translation of French comments within the FORTRAN sources into
!+   English comments
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        21/08/2010
!+        V6P0
!+   Creation of DOXYGEN tags for automated documentation and
!+   cross-referencing of the FORTRAN sources
!
!history  G.MATTAROLO (EDF - LNHE)
!+        14/06/2011
!+        V6P1
!+   Translation of French names of the variables in argument
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| COEF1          |---| MULTIPLYING COEF. F(JFRE)=COEF1*F(JBIS)
!| JBIS           |---| ADJUSTED INDEX IN THE INTERVAL [1;NF]
!| JFRE           |-->| FREQUENCY INDEX
!| NF             |-->| NUMBER OF FREQUENCIES
!| RAISF          |-->| FREQUENTIAL RATIO
!| TAILF          |-->| SPECTRUM QUEUE FACTOR
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE interface_tomawac, ex_cqueue => cqueue
      USE declarations_tomawac, ONLY : nf , raisf , tailf
      IMPLICIT NONE
!
!.....VARIABLES IN ARGUMENT
!     """""""""""""""""""""
      INTEGER, INTENT(IN)             :: JFRE
      DOUBLE PRECISION, INTENT(INOUT) :: COEF1
      INTEGER, INTENT(INOUT)          :: JBIS
!
!
      IF (jfre.GT.nf) THEN
        jbis = nf
        coef1= 1.d0/raisf**(dble(jfre-nf)*tailf)
      ELSEIF (jfre.LT.1) THEN
        jbis = 1
        coef1= 0.d0
      ELSE
        jbis = jfre
        coef1= 1.d0
      ENDIF
!
      RETURN
      END