read_spectra_coords.f

Go to the documentation of this file.

!                   ******************************
                    SUBROUTINE read_spectra_coords
!                   ******************************
!
     &     (fid,np,xp,yp)
!
!***********************************************************************
!     TOMAWAC   V7P3                                   21/02/2017
!***********************************************************************
!
!     brief    READS A LIST OF COORDINATES FROM AN EXTERNAL TEXT FILE
!     +        AND ADDS THEM TO THE EXISTING COORDINATES
!
!     history  A. JOLY (EDF - LNHE)
!     +        21/02/2017
!     +        V7P3
!     +   CREATED

!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| FID            |-->| ID OF THE TEXT FILE TO BE READ
!| NP             |<->| NUMPER OF POINTS READ
!| XP,YP          |<->| COORDINATES OF THE SPECTRUM POINTS
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_special
!
      IMPLICIT NONE
!
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN)                         :: FID
      INTEGER, INTENT(INOUT)                      :: NP
      DOUBLE PRECISION,ALLOCATABLE, INTENT(INOUT) :: XP(:),YP(:)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER I
      INTEGER I_GLO
      DOUBLE PRECISION NE,ZP
!
!-----------------------------------------------------------------------
!
! CHECK TO SEE IF THE PRINTOUT POINTS ARE ONLY DEFINED ONE WAY
!
      IF(np.GT.0)THEN
        WRITE(lu,*) ''
        WRITE(lu,*) '**************************************'
        WRITE(lu,*) '   ERROR : SPECTRA PRINTOUT POINTS'
        WRITE(lu,*) '   THE PRINTOUT POINTS ARE DEFINED'
        WRITE(lu,*) '   USING BOTH KEYWORDS AND A FILE'
        WRITE(lu,*) '**************************************'
        CALL plante(1)
        stop
      ENDIF
!
! READ THE NUMBER OF POINTS
      READ(fid,*) np,ne
      IF(ALLOCATED(xp))DEALLOCATE(xp)
      ALLOCATE(xp(np))
      IF(ALLOCATED(yp))DEALLOCATE(yp)
      ALLOCATE(yp(np))
! READ ALL THE POINTS
      DO i=1,np
        READ(fid,*) i_glo,xp(i),yp(i),zp
      ENDDO
!
      RETURN
      END