lecr3d.f

Go to the documentation of this file.

!                       *****************
                        SUBROUTINE lecr3d
!                       *****************
!
     &(irec,at,z,u,v,w,npoin3,npoin2,nplan,nres,fformat,nva3,tab)
!
!***********************************************************************
! POSTEL3D VERSION 5.1   01/09/99   T. DENOT (LNH) 01 30 87 74 89
! FORTRAN90
!***********************************************************************
!
!     FONCTION  : LIT LES ENREGISTREMENTS 3D D'UN PAS DE TEMPS
!
!-----------------------------------------------------------------------
!                             ARGUMENTS
! .________________.____.______________________________________________.
! !      NOM       !MODE!                   ROLE                       !
! !________________!____!______________________________________________!
! !   IREC         !--> ! PAS DE TEMPS TRAITE                          !
! !   AT           !<-- ! TEMPS CORRESPONDANT AU PAS TRAITE            !
! !   U,V,W        !<-- ! COMPOSANTES 3D DE LA VITESSE                 !
! !   TA,TP        !<-- ! CONCENTRATIONS DES TRACEURS                  !
! !   NUX,NUY,NUZ  !<-- ! COEFFICIENTS DE VISCOSITE POUR LES VITESSES  !
! !   NAX,NAY,NAZ  !<-- ! COEFFICIENTS DE VISCOSITE POUR LES TR.ACTIFS !
! !   NPX,NPY,NPZ  !<-- ! COEFFICIENTS DE VISCOSITE POUR LES TR.PASSIFS!
! !   RI           !<-- ! NOMBRE DE RICHARDSON                         !
! !   AK,EP        !<-- ! VARIABLES DU MODELE K-EPSILON                !
! !   RHO          !<-- ! ECARTS RELATIFS DE DENSITE                   !
! !   H            !<-- ! HAUTEUR D'EAU                                !
! !   Z            !<-- ! COTES DES NOEUDS                             !
! !   ZSTAR        ! -->! COTES RELATIVES DES NOEUDS                   !
! !   NPOIN2       ! -->! NOMBRE DE POINTS DU MAILLAGE 2D              !
! !   NPOIN3       ! -->! NOMBRE DE POINTS DU MAILLAGE 3D              !
! !   NRES         ! -->! NUMERO DE CANAL DU FICHIER DE RESULTAT 3D    !
! !   NPLAN        ! -->! NOMBRE DE PLANS                              !
! !   NPLINT       ! -->! NUMERO DU PLAN DE CHANGEMENT DE TRANSFORMATION
! !   NTRAC        ! -->! NOMBRE DE TRACEURS ACTIFS                    !
! !   NTRPA        ! -->! NOMBRE DE TRACEURS PASSIFS                   !
! !   SORG3D       ! -->! INDICATEUR DES VARIABLES ENREGISTREES        !
! !   R4           ! -->! TABLEAU DE REELS SIMPLE PRECISION POUR LIT   !
! !________________!____!______________________________________________!
! MODE : -->(DONNEE NON MODIFIEE), <--(RESULTAT), <-->(DONNEE MODIFIEE)
!
! SOUS-PROGRAMME APPELE PAR : POSTEL3D
! SOUS-PROGRAMME APPELES : LIT
!
!history Y AUDOUIN (LNHE)
!+       25/05/2015
!+       V7P0
!+       Modification to comply with the hermes module
!
!**********************************************************************
!
      USE bief
      USE declarations_postel3d, ONLY: textlu
      USE interface_hermes
      USE interface_postel3d, ex_lecr3d => lecr3d
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN) :: IREC
      INTEGER, INTENT(IN) :: NPOIN3,NPOIN2,NPLAN,NRES
      CHARACTER(LEN=8), INTENT(IN) :: FFORMAT
      INTEGER, INTENT(INOUT) :: NVA3
      type(bief_obj) , INTENT(INOUT) :: tab
      DOUBLE PRECISION , INTENT(INOUT) :: AT
      DOUBLE PRECISION , INTENT(INOUT) :: U(npoin3)
      DOUBLE PRECISION , INTENT(INOUT) :: V(npoin3)
      DOUBLE PRECISION , INTENT(INOUT) :: W(npoin3)
      DOUBLE PRECISION , INTENT(INOUT) :: Z(npoin2,nplan)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER I
      INTEGER IERR
!
!***********************************************************************
!
!
! LECTURE DU TEMPS DU DEBUT DU CALCUL
!
      CALL get_data_time(fformat,nres,irec,at,ierr)
      CALL check_call(ierr,'LECR3D:GET_DATA_TIME')
!
! LECTURE DES VITESSES U,V ET W
!
      CALL get_data_value(fformat,nres,irec,textlu(1),z,npoin3,ierr)
      CALL check_call(ierr,'LECR3D:GET_DATA_VALUE:Z')
      CALL get_data_value(fformat,nres,irec,textlu(2),u,npoin3,ierr)
      CALL check_call(ierr,'LECR3D:GET_DATA_VALUE:U')
      CALL get_data_value(fformat,nres,irec,textlu(3),v,npoin3,ierr)
      CALL check_call(ierr,'LECR3D:GET_DATA_VALUE:V')
      CALL get_data_value(fformat,nres,irec,textlu(4),w,npoin3,ierr)
      CALL check_call(ierr,'LECR3D:GET_DATA_VALUE:W')
!
      IF (nva3.GT.4) THEN
        DO i=1,nva3-4
          CALL get_data_value(fformat,nres,irec,textlu(i+4),
     &                        tab%ADR(i)%P%R,npoin3,ierr)
          CALL check_call(ierr,'LECR3D:GET_DATA_VALUE:W')
        ENDDO
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END