user_corrxy.f

Go to the documentation of this file.

!                   **********************
                    SUBROUTINE user_corrxy
!                   **********************
!
     & (x,y,npoin)
!
!***********************************************************************
! BIEF
!***********************************************************************
!
!brief    MODIFIES THE COORDINATES OF THE POINTS IN THE MESH.
!
!warning  DO NOT PERFORM ROTATIONS AS IT WILL CHANGE
!+            THE NUMBERING OF THE LIQUID BOUNDARIES
!
!history  EMILE RAZAFINDRAKOTO (LNHE)
!+        17/10/05
!+        V5P6
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| NPOIN          |-->| NUMBER OF POINTS IN THE MESH
!| X              |<->| ABSCISSAE OF POINTS IN THE MESH
!| Y              |<->| ORDINATES OF POINTS IN THE MESH
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_user_corrxy => user_corrxy
!
!     OTHER DATA ARE AVAILABLE WITH THE DECLARATIONS OF EACH PROGRAM
!
!     USE DECLARATIONS_TELEMAC2D
!     USE DECLARATIONS_TELEMAC3D
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN) :: NPOIN
      DOUBLE PRECISION, INTENT(INOUT) :: X(npoin),Y(npoin)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER I
      DOUBLE PRECISION DEGTORAD
      INTRINSIC acos
!
!-----------------------------------------------------------------------
!
!     THIS SUBROUTINE MUST BE MODIFIED ACCORDING TO
!     THE CALLING PROGRAM AND THE NEEDED MODIFICATION
!     BY ADDING USE DECLARATIONS_"NAME OF CALLING CODE"
!     ALL THE DATA STRUCTURE OF THIS CODE IS AVAILABLE
!
!-----------------------------------------------------------------------
!
!      EXAMPLE 1: MULTIPLIES BY A CONSTANT (SCALES THE MESH)
!                 AND CHANGES THE ORIGIN
!
!
      IF(.false.) THEN
        DO i = 1 , npoin
          x(i) = 3.d0 * x(i) + 100.d0
          y(i) = 5.d0 * y(i) - 50.d0
        ENDDO
      ENDIF
!
!
!-----------------------------------------------------------------------
!
!      EXAMPLE 2: CHANGING LATITUDE-LONGITUDE IN DEGREES TO RADIANS
!
!
      IF(.false.) THEN
        degtorad=acos(-1.d0)/180.d0
        DO i = 1 , npoin
          x(i) = x(i) * degtorad
          y(i) = y(i) * degtorad
        ENDDO
      ENDIF
!
!-----------------------------------------------------------------------
!
!     THIS SHOULD BE CHANGED IF MODIFICATIONS ARE DONE
!
      WRITE(lu,*)'USER_CORRXY (BIEF):NO MODIFICATION OF COORDINATES'
      WRITE(lu,*)
!
!-----------------------------------------------------------------------
!
      RETURN
      END