hpropa.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac2d\hpropa.f
00002 !
00053                      SUBROUTINE HPROPA
00054 !                    *****************
00055 !
00056      &(HPROP ,HN,H,PROLIN,HAULIN,TETA,NSOUSI)
00057 !
00058 !***********************************************************************
00059 ! TELEMAC2D   V6P1                                   21/08/2010
00060 !***********************************************************************
00061 !
00062 !         OF CONTINUITY EQUATION.
00063 !
00064 !
00065 !
00066 !
00067 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00068 !| H              |<--| WATER DEPTH AT TIME T(N+1)
00069 !| HAULIN         |-->| MEAN DEPTH FOR LINEARISATION
00070 !| HN             |<--| WATER DEPTH AT TIME T(N)
00071 !| HPROP          |<--| WATER DEPTH FOR PROPAGATION
00072 !| NSOUSI         |-->| NUMBER OF SUB-ITERATIONS
00073 !| PROLIN         |-->| KEY-WORD "LINEARISED PROPAGATON"
00074 !| TETA           |-->| IMPLICITATION COEFFICIENT ON DEPTH
00075 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00076 !
00077       USE BIEF
00078 !
00079       IMPLICIT NONE
00080       INTEGER LNG,LU
00081       COMMON/INFO/LNG,LU
00082 !
00083 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00084 !
00085       INTEGER, INTENT(IN)           :: NSOUSI
00086       LOGICAL, INTENT(IN)           :: PROLIN
00087       DOUBLE PRECISION, INTENT(IN)  :: TETA,HAULIN
00088       TYPE(BIEF_OBJ), INTENT(IN)    :: HN,H
00089       TYPE(BIEF_OBJ), INTENT(INOUT) :: HPROP
00090 !
00091 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00092 !
00093       IF(PROLIN) THEN
00094         CALL OS( 'X=C     ' , X=HPROP , C=HAULIN    )
00095       ELSEIF(NSOUSI.EQ.1) THEN
00096         CALL OS( 'X=Y     ' , X=HPROP , Y=HN )
00097       ELSE
00098         CALL OS( 'X=CY    ' , X=HPROP , Y=HN , C=1.D0-TETA )
00099         CALL OS( 'X=X+CY  ' , X=HPROP , Y=H  , C= TETA )
00100       ENDIF
00101 !
00102 !-----------------------------------------------------------------------
00103 !
00104 !     CLIPS HPROP
00105 !
00106       IF(.NOT.PROLIN) THEN
00107         CALL OS('X=+(Y,C)',X=HPROP,Y=HPROP,C=0.D0)
00108       ENDIF
00109 !
00110 !-----------------------------------------------------------------------
00111 !
00112       RETURN
00113       END

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