calcg.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac3d\calcg.f
00002 !
00053                      SUBROUTINE CALCG
00054 !                    ****************
00055 !
00056      & (TRAV2,TRAV3,U,V,UETCAR,NPOIN2,NPLAN)
00057 !
00058 !***********************************************************************
00059 ! TELEMAC3D   V6P1                                   21/08/2010
00060 !***********************************************************************
00061 !
00062 !
00063 !
00064 !
00065 !
00066 !
00067 !
00068 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00069 !| NPLAN          |-->| NUMBER OF PLANES IN THE 3D MESH OF PRISMS
00070 !| NPOIN2         |-->| NOMBER OF POINTS IN THE 2D MESH
00071 !| TRAV2          |-->| WORK ARRAYS
00072 !| TRAV3          |<->| WORK ARRAYS
00073 !| U              |-->| COMPONENT OF VELOCITY
00074 !| UETCAR         |-->| USTAR**2
00075 !| V              |-->| COMPONENT OF VELOCITY
00076 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00077 !
00078       IMPLICIT NONE
00079 !
00080 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00081 !
00082       INTEGER         , INTENT(IN)    :: NPOIN2,NPLAN
00083       DOUBLE PRECISION, INTENT(IN)    :: U(NPOIN2),V(NPOIN2)
00084       DOUBLE PRECISION, INTENT(IN)    :: UETCAR(NPOIN2)
00085       DOUBLE PRECISION, INTENT(IN)    :: TRAV2(NPOIN2*NPLAN)
00086       DOUBLE PRECISION, INTENT(INOUT) :: TRAV3(NPOIN2*NPLAN)
00087 !
00088 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00089 !
00090       INTEGER IPOIN2,IPLAN
00091 !
00092 !=======================================================================
00093 !
00094 !     COMPUTES AUBORF IN TRAV1
00095 !
00096       DO IPOIN2 = 1,NPOIN2
00097 !
00098 !        ACCORDING TO RICHARD SOULSBY (HRW):
00099 !        G=SQRT(EPSILON/NU) AND EPSILON=(TAU/RHO)*(DU/DZ)
00100 !        THEN G COULD BE ALSO SQRT(NUZ/NU)*DU/DZ
00101 !
00102 !
00103 !                              USTAR**2 * DU/DZ
00104 !        COMPUTES  G  =  SQRT( --------------------- ) : TRAV3
00105 !                                     NU
00106 !
00107         DO IPLAN = 1,NPLAN
00108           TRAV3(IPOIN2+(IPLAN-1)*NPOIN2) =
00109      &        SQRT(UETCAR(IPOIN2)*1.D06*TRAV2(IPOIN2+(IPLAN-1)*NPOIN2))
00110 !
00111         ENDDO
00112 !
00113       ENDDO
00114 !
00115 !=======================================================================
00116 !
00117       RETURN
00118       END

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