twcclosest.f

Go to the documentation of this file.

!                   *********************
                    SUBROUTINE twcclosest
!                   *********************
!
!***********************************************************************
! ARTEMIS   V7P4                                     Nov 2017
!***********************************************************************
!
!brief    IDENTIFIES ARTEMIS BOUNDARY NODE CLOSEST TO REFERENCE POINT
!+               (ONLY REQUIRED FOR CALL TO ENTART, artemis.f line 489)
!
!history  N.DURAND (HRW)
!+        Nov 2017
!+        V7P4
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief
      USE declarations_artemis
!
      USE declarations_special
      IMPLICIT NONE
!
!-----------------------------------------------------------------------
!
      INTEGER                       :: IPTFR
      DOUBLE PRECISION              :: EPS, DIST
!
!=======================================================================
!
      eps = 3000.d0
      iptfr_ref = 0
!
      DO iptfr = 1 , nptfr
        dist = dsqrt((x(mesh%NBOR%I(iptfr))-x_sfref)**2 +
     &               (y(mesh%NBOR%I(iptfr))-y_sfref)**2 )
        IF(dist.LT.eps) THEN
          eps = dist
          iptfr_ref = iptfr
        ENDIF
      ENDDO
!
      IF(iptfr_ref.EQ.0) THEN
        WRITE(lu,201) x_sfref, y_sfref
!       CALL PLANTE(1)
!       STOP
      ENDIF
!
!-----------------------------------------------------------------------
!
!     PRINTOUT FORMATS:
!
201   FORMAT(/,1x,'CLOSEST : PLEASE REVIEW COORDINATES FOR THE',
     &       ' REFERENCE POINT (',1f9.2,';',1f9.2,')')
!
!-----------------------------------------------------------------------
!
      RETURN
      END