dismoy.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE dismoy
!                   *****************
!
     &(npoin,nelem,x,y,ikle,k,lishho)
!
!***********************************************************************
! ARTEMIS   V7P3                                     Aug 2017
!***********************************************************************
!
!brief    CALCULATES THE NUMBER OF SMOOTHINGS ON THE WAVE
!+                HEIGHT (LISHHO), A PRIORI NECESSARY TO FILTER OUT
!+                THE PARASITIC OSCILLATIONS (REGULAR WAVES).
!+                ESTIMATED FROM THE AVERAGE DISTANCE BETWEEN
!+                NODES AND THE AVERAGE NUMBER OF NODES IN HALF A
!+                WAVELENGTH.
!
!history  J-M HERVOUET (LNH)
!+
!+
!+   LINKED TO BIEF 5.0
!
!history  D. AELBRECHT (LNH) ; P. THELLIER (LNH)
!+        04/06/1999
!+        V5P1
!+   THANK YOU PAUL
!
!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  N.DURAND (HRW)
!+        August 2017
!+        V7P3
!+   PI now defined in DECLARATIONS_ARTEMIS
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| IKLE           |---| CONNECTIVITY TABLE
!| K              |-->| WAVE NUMBER
!| LISHHO         |<--| NUMBER OF SMOOTHING FOR HHO
!| NELEM          |-->| NUMBER OF ELEMENTS
!| NPOIN          |-->| NUMBER OF POINT
!| X,Y            |-->| MESH COORDINATES
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE interface_artemis, ex_dismoy => dismoy
!
      USE declarations_special
      USE declarations_artemis, ONLY : pi
      IMPLICIT NONE
!
!
! VARIABLES FOR THE COMPUTATION OF AVERAGE DISTANCES
!
      INTEGER NG, ISOM, IELEM
      DOUBLE PRECISION SOM, D1, D2, DMOY, SOMD
!
!
      INTEGER NPOIN,NELEM,LISHHO
      INTEGER IKLE(nelem,*)
!
      DOUBLE PRECISION X(npoin),Y(npoin),K(npoin)
!
!-----------------------------------------------------------------------
!
!     COMPUTES THE AVERAGE DISTANCE BETWEEN A NODE AND ITS NEIGHBOURS
!     -----------------------------------------------------------
!
      somd = 0.d0
      DO ng=1,npoin
        isom = 0
        som = 0.d0
!
        DO ielem=1,nelem
!
          IF (ikle(ielem,1).EQ.ng) THEN
!         --                       ---
!
!           PT 1 IS COMMON, COMPUTES THE DISTANCE TO PTS 2 AND 3
            d1 = sqrt( (x(ng)-x(ikle(ielem,2))) **2.d0 +
     &                  (y(ng)-y(ikle(ielem,2))) **2.d0 )
            d2 = sqrt( (x(ng)-x(ikle(ielem,3))) **2.d0 +
     &                  (y(ng)-y(ikle(ielem,3))) **2.d0 )
            som = som + d1 + d2
            isom = isom + 2
!
          ELSEIF (ikle(ielem,2).EQ.ng) THEN
!         ------                       ----
!
!           PT 2 IS COMMON, COMPUTES THE DISTANCE TO PTS 1 AND 3
            d1 = sqrt( (x(ng)-x(ikle(ielem,1))) **2.d0 +
     &                  (y(ng)-y(ikle(ielem,1))) **2.d0 )
            d2 = sqrt( (x(ng)-x(ikle(ielem,3))) **2.d0 +
     &                  (y(ng)-y(ikle(ielem,3))) **2.d0 )
            som = som + d1 + d2
            isom = isom + 2
!
          ELSEIF (ikle(ielem,3).EQ.ng) THEN
!         ------                       ----
!
!           PT 3 IS COMMON, COMPUTES THE DISTANCE TO PTS 1 AND 2
            d1 = sqrt( (x(ng)-x(ikle(ielem,1))) **2.d0 +
     &                  (y(ng)-y(ikle(ielem,1))) **2.d0 )
            d2 = sqrt( (x(ng)-x(ikle(ielem,2))) **2.d0 +
     &                  (y(ng)-y(ikle(ielem,2))) **2.d0 )
            som = som + d1 + d2
            isom = isom + 2
!
          ENDIF
!         -----
!
        ENDDO ! IELEM
!
        dmoy = som / float(isom)
        somd = somd + (pi / (k(ng)*dmoy))
!
      ENDDO ! NG
!
!     ESTIMATES THE NUMBER OF SMOOTHINGS FROM THE AVERAGE DISTANCE
!     -------------------------------------------------------------------
!
      lishho = int((somd/float(npoin))) * 10
!
      RETURN
      END