read_bin_2d.f

Go to the documentation of this file.

!                    *************************
                      SUBROUTINE read_bin_2d
!                    *************************
!
     &( q ,varname,at,nfic,fformat,npoin,found,offset)
!
!***********************************************************************
! BIEF   V7P2                                       01/09/2016
!***********************************************************************
!
!brief    READS AND INTERPOLATES VALUES FROM A BINARY
!+        FILE ON THE WHOLE 2D MESH.
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| AT             |-->| TIME IN SECONDS
!| NFIC           |-->| LOGICAL UNIT OF FILE
!| Q              |<--| ARRAY WHERE THE VARIABLE IS STORED
!| FOUND          |<--| IF FALSE: VARIABLE NOT FOUND
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_read_bin_2d => read_bin_2d
      USE interface_hermes
!
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      CHARACTER(LEN=16), INTENT(IN)      :: VARNAME
      INTEGER         , INTENT(IN)       :: NFIC,NPOIN
      DOUBLE PRECISION, INTENT(IN)       :: AT
      DOUBLE PRECISION, INTENT(INOUT)    :: Q(npoin)
      CHARACTER(LEN=8), INTENT(IN)       :: FFORMAT
      LOGICAL         , INTENT(OUT)      :: FOUND
      DOUBLE PRECISION, INTENT(IN)       :: OFFSET
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER :: J,NPOINATM,IERR
!
!-----------------------------------------------------------------------
!
!     INITIALISE THE ARRAY TO BE FILLED
      DO j=1,npoin
        q(j) = 0.d0
      ENDDO
!
!-----------------------------------------------------------------------
!
!     READ THE FILE
!     AND STORE VARIABLE VALUE AT ALL 2D POINTS
!
!-----------------------------------------------------------------------
!
!     GET THE NUMBER OF POINTS IN THE FILE
!
      CALL get_mesh_npoin(fformat,nfic,triangle_elt_type,
     &                    npoinatm,ierr)
      CALL check_call(ierr,"READ_BIN_FRILQ: GET_MESH_NPOIN")
      IF(npoinatm.NE.npoin) THEN
        WRITE(lu,*) "READ_BIN_2D: ",
     &              "INCORRECT NUMBER OF POINTS IN THE FILE"
        CALL plante(1)
        stop
      ENDIF
!
!-----------------------------------------------------------------------
!
!     STORE THE VARIABLE VALUES IN Q
!
      CALL find_variable(fformat,nfic,varname,q,npoin,ierr,
     &                 time=at,eps_time=1.d-6,offset=offset)
      IF(ierr.EQ.hermes_var_unknown_err) THEN
        found = .false.
      ELSE
        found = .true.
      ENDIF
      CALL check_call(ierr,"READ_BIN_2D: FIND_VARIABLE")

      RETURN
      END