lecr3d.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\postel3d\lecr3d.f
00002 !
00032                         SUBROUTINE LECR3D
00033 !                       *****************
00034 !
00035      &(AT,Z,U,V,W,NPOIN3,NPOIN2,NPLAN,NRES,BINRES,RB,NVA3,TAB,VARSUB)
00036 !
00037 !***********************************************************************
00038 ! POSTEL3D VERSION 5.1   01/09/99   T. DENOT (LNH) 01 30 87 74 89
00039 ! FORTRAN90
00040 !***********************************************************************
00041 !
00042 !     FONCTION  : LIT LES ENREGISTREMENTS 3D D'UN PAS DE TEMPS
00043 !
00044 !-----------------------------------------------------------------------
00045 !                             ARGUMENTS
00046 ! .________________.____.______________________________________________.
00047 ! !      NOM       !MODE!                   ROLE                       !
00048 ! !________________!____!______________________________________________!
00049 ! !   AT           !<-- ! TEMPS CORRESPONDANT AU PAS TRAITE            !
00050 ! !   U,V,W        !<-- ! COMPOSANTES 3D DE LA VITESSE                 !
00051 ! !   TA,TP        !<-- ! CONCENTRATIONS DES TRACEURS                  !
00052 ! !   NUX,NUY,NUZ  !<-- ! COEFFICIENTS DE VISCOSITE POUR LES VITESSES  !
00053 ! !   NAX,NAY,NAZ  !<-- ! COEFFICIENTS DE VISCOSITE POUR LES TR.ACTIFS !
00054 ! !   NPX,NPY,NPZ  !<-- ! COEFFICIENTS DE VISCOSITE POUR LES TR.PASSIFS!
00055 ! !   RI           !<-- ! NOMBRE DE RICHARDSON                         !
00056 ! !   AK,EP        !<-- ! VARIABLES DU MODELE K-EPSILON                !
00057 ! !   RHO          !<-- ! ECARTS RELATIFS DE DENSITE                   !
00058 ! !   H            !<-- ! HAUTEUR D'EAU                                !
00059 ! !   Z            !<-- ! COTES DES NOEUDS                             !
00060 ! !   ZSTAR        ! -->! COTES RELATIVES DES NOEUDS                   !
00061 ! !   NPOIN2       ! -->! NOMBRE DE POINTS DU MAILLAGE 2D              !
00062 ! !   NPOIN3       ! -->! NOMBRE DE POINTS DU MAILLAGE 3D              !
00063 ! !   NRES         ! -->! NUMERO DE CANAL DU FICHIER DE RESULTAT 3D    !
00064 ! !   NPLAN        ! -->! NOMBRE DE PLANS                              !
00065 ! !   NPLINT       ! -->! NUMERO DU PLAN DE CHANGEMENT DE TRANSFORMATION
00066 ! !   NTRAC        ! -->! NOMBRE DE TRACEURS ACTIFS                    !
00067 ! !   NTRPA        ! -->! NOMBRE DE TRACEURS PASSIFS                   !
00068 ! !   SORG3D       ! -->! INDICATEUR DES VARIABLES ENREGISTREES        !
00069 ! !   R4           ! -->! TABLEAU DE REELS SIMPLE PRECISION POUR LIT   !
00070 ! !________________!____!______________________________________________!
00071 ! MODE : -->(DONNEE NON MODIFIEE), <--(RESULTAT), <-->(DONNEE MODIFIEE)
00072 !
00073 ! SOUS-PROGRAMME APPELE PAR : POSTEL3D
00074 ! SOUS-PROGRAMME APPELES : LIT
00075 !
00076 !**********************************************************************
00077 !
00078       USE BIEF
00079 !
00080       IMPLICIT NONE
00081       INTEGER LNG,LU
00082       COMMON/INFO/LNG,LU
00083 !
00084       INTEGER NPOIN3,NPOIN2,NPLAN,NRES
00085       INTEGER NVA3
00086 !      INTEGER , INTENT(IN) :: NELEM3
00087       CHARACTER*3 BINRES
00088 !
00089       DOUBLE PRECISION , INTENT(INOUT) :: AT
00090       DOUBLE PRECISION , INTENT(INOUT) :: U(NPOIN3)
00091       DOUBLE PRECISION , INTENT(INOUT) :: V(NPOIN3)
00092       DOUBLE PRECISION , INTENT(INOUT) :: W(NPOIN3)
00093       DOUBLE PRECISION , INTENT(INOUT) :: Z(NPOIN2,NPLAN)
00094       LOGICAL , INTENT(IN) :: VARSUB
00095       TYPE (BIEF_OBJ) , INTENT(INOUT) :: TAB
00096 !
00097       INTEGER I,ISTAT
00098       DOUBLE PRECISION XB(2)
00099       INTEGER IB(2)
00100       CHARACTER(LEN=1) CB
00101       REAL RB(NPOIN3)
00102 !
00103 !***********************************************************************
00104 !
00105 !
00106 ! LECTURE DU TEMPS DU DEBUT DU CALCUL
00107 !
00108       CALL LIT(XB,RB,IB,CB,1,'R4',NRES,BINRES,ISTAT)
00109       AT=XB(1)
00110 !
00111 ! LECTURE DES VITESSES U,V ET W
00112 !
00113       CALL LIT(Z,RB,IB,CB,NPOIN3,'R4',NRES,BINRES,ISTAT)
00114       CALL LIT(U,RB,IB,CB,NPOIN3,'R4',NRES,BINRES,ISTAT)
00115       CALL LIT(V,RB,IB,CB,NPOIN3,'R4',NRES,BINRES,ISTAT)
00116       CALL LIT(W,RB,IB,CB,NPOIN3,'R4',NRES,BINRES,ISTAT)
00117 !
00118       IF (NVA3.GT.4) THEN
00119       DO I=1,NVA3-4
00120         CALL LIT(TAB%ADR(I)%P%R,RB,IB,CB,NPOIN3,'R4',NRES,BINRES,ISTAT)
00121       ENDDO
00122       ENDIF
00123 !
00124 !-----------------------------------------------------------------------
00125 !
00126       RETURN
00127       END

Generated on Fri Aug 31 2013 18:12:58 by S.E.Bourban (HRW) using doxygen 1.7.0