preqt2.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE preqt2
!                   *****************
!
!***********************************************************************
! TOMAWAC   V8P0                                   32/12/2018
!***********************************************************************
!
!brief    SOURCE TERM RELATED TO NON-LINEAR INTERACTIONS
!+                BETWEEN FREQUENCY TRIPLETS.
!+                DEVELOPED FROM THE BOUSSINESQ EQUATIONS.
!
!history  EDF/DER/LNH
!+        11/06/98
!+        V5P0
!+
!
!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)
!+        23/06/2011
!+        V6P1
!+   Translation of French names of the variables in argument
!
      USE declarations_tomawac
      USE interface_tomawac, ex_preqt2 => preqt2
      IMPLICIT NONE
!
!.....VARIABLES FROM MODULE TOMAWAC
!     """""""""""""""""""""""""""""
!| QINDI          |<--| CONFIGURATION INDEX
!| NBD            |-->| NUMBER OF TRIAD CONFIGURATIONS
!| BDISPB         |-->| LOWER DIRECTIONAL BOUND. OF SPB TRIAD MODEL
!| BDSSPB         |-->| UPPER DIRECTIONAL BOUND. OF SPB TRIAD MODEL
!| KSPB           |-->| COEFFICIENT K OF SPB TRIAD INTERACTION MODEL

!
!.....LOCAL VARIABLES
!     """""""""""""""""
      INTEGER           IPL
      DOUBLE PRECISION  AP2 , DTETA
!
      INTEGER           NBPL , NBPU, NB1
!
!     """"""""""""""""""""""""""""""""""""""""""""""""""""""""""
      IF(.NOT.deja_qt2) THEN
        ALLOCATE(qindi(ndire))
        deja_qt2=.true.
      ENDIF
!
      dteta = teta(2)-teta(1)
      IF(bdsspb.GE.bdispb) THEN
        ap2  = (bdispb-teta(1))/dteta
        nbpl = nint(ap2) + 1
        ap2  = (bdsspb-teta(1))/dteta
        nbpu = nint(ap2)
        nbd=nbpu-nbpl+1
        DO ipl=1,nbd
          qindi(ipl)=nbpl+ipl-1
        END DO
      ELSE
        ap2  = (bdsspb-teta(1))/dteta
        nbpu = nint(ap2)  + 1
        ap2  = (bdispb-teta(1))/dteta
        nbpl = nint(ap2)
        IF(nbpl.GT.ndire) THEN
          nbpl = 1
          qindi(1) = 1
          nbd  = nbpu - nbpl + 1
          DO ipl = 2,nbd
            qindi(ipl)=ipl
          END DO
        ELSE
          nb1 = ndire - nbpl + 1
          nbd = nb1 + nbpu
          DO ipl = 1,nb1
            qindi(ipl)=nbpl+ipl-1
          END DO
          DO ipl = 1,nbpu
            qindi(ipl+nb1)=ipl
          END DO
        ENDIF
      ENDIF
!
      RETURN
      END