fdiss3d.f

Go to the documentation of this file.

!#######################################################################
                        SUBROUTINE fdiss3d
     &  (fdx, fdy, npoin2, xk, ndire, fs,nf)
!  SURFACE STRESS DUE TO DEPTH INDUCED WAVE BREAKING
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| FS             |-->| VARIANCE DENSITY DIRECTIONAL SPECTRUM
!| FDX            |<--| SURFACE STRESS ALONG X
!| FDY            |<--| SURFACE STRESS ALONG Y 
!| NF             |-->| NUMBER OF FREQUENCIES
!| NDIRE          |-->| NUMBER OF DIRECTIONS
!| NPOIN2         |-->| NUMBER OF POINTS IN 2D MESH
!| XK             |-->| DISCRETIZED WAVE NUMBER
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
      USE bief
      USE declarations_tomawac, ONLY: freq, dfreq, sintet, costet, 
     &                                betabr, deupi 
      USE interface_tomawac, ex_fdiss3d => fdiss3d
    
      IMPLICIT NONE
      
!.....VARIABLES IN ARGUMENT
!     """"""""""""""""""""
      INTEGER, INTENT(IN) :: NPOIN2, NDIRE,NF
      DOUBLE PRECISION, INTENT(IN) :: FS(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN) :: XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: FDX(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: FDY(npoin2)  

!     """""""""""""""""
!.....LOCAL VARIABLES
!     """""""""""""""""
      INTEGER  JP    , JF    , IP
      DOUBLE PRECISION SIGMA
      DOUBLE PRECISION DTETAR, AUX1
!
      dtetar=deupi/dble(ndire)
      DO ip=1,npoin2
        fdx(ip) = 0.d0
        fdy(ip) = 0.d0
      ENDDO
      DO jf=1,nf
        sigma=deupi*freq(jf)
        aux1=dfreq(jf)*dtetar        
        DO jp=1,ndire
          DO ip=1,npoin2
            fdx(ip)=fdx(ip)+((xk(ip,jf)/sigma)*sintet(jp)
     & *betabr(ip)*fs(ip,jp,jf))*aux1
            fdy(ip)=fdy(ip)+((xk(ip,jf)/sigma)*costet(jp)
     & *betabr(ip)*fs(ip,jp,jf))*aux1
          ENDDO
        ENDDO
      ENDDO

      RETURN
      END