newstr.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac2d\newstr.f
00002 !
00066                      SUBROUTINE NEWSTR
00067 !                    *****************
00068 !
00069      &(SETSTR,SETSTR2,DESC,RO,RSTART,NPARAM,ESTIME,KFROT)
00070 !
00071 !***********************************************************************
00072 ! TELEMAC2D   V6P1                                   21/08/2010
00073 !***********************************************************************
00074 !
00075 !
00076 !
00077 !
00078 !
00079 !
00080 !
00081 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00082 !| DESC           |-->| DIRECTION OF DESCENT
00083 !| ESTIME         |---| ???????? NOT USED
00084 !| KFROT          |-->| LAW OF FRICTION
00085 !| NPARAM         |-->| NUMBER OF PARAMETERS TO ESTIMATE
00086 !| RO             |-->| SETSTR=SETSTR2+RO*DESC
00087 !| RSTART         |-->| LOGICAL, RESTART COMPUTATION BECAUSE OUT OF LIMITS
00088 !| SETSTR         |-->| SET OF FRICTION COEFFICIENTS IN A BIEF_OBJ STRUCTURE
00089 !| SETSTR2        |-->| OLD SET
00090 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00091 !
00092       USE BIEF
00093       IMPLICIT NONE
00094 !
00095       INTEGER LNG,LU
00096       COMMON/INFO/LNG,LU
00097 !
00098 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00099 !
00100       DOUBLE PRECISION , INTENT(IN)    :: RO
00101       TYPE (BIEF_OBJ)  , INTENT(IN)    :: DESC
00102       TYPE (BIEF_OBJ)  , INTENT(IN)    :: SETSTR2
00103       TYPE (BIEF_OBJ)  , INTENT(INOUT) :: SETSTR
00104       LOGICAL          , INTENT(INOUT) :: RSTART
00105       INTEGER          , INTENT(IN)    :: NPARAM,KFROT
00106       CHARACTER(LEN=72), INTENT(IN)    :: ESTIME
00107 !
00108 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00109 !
00110       INTEGER I
00111 !
00112 !---------------------------------------------------------------------
00113 !
00114       CALL OV( 'X=Y+CZ  ',SETSTR%R,SETSTR2%R,DESC%R,RO,NPARAM)
00115 !
00116 !     TESTS LIMITS
00117 !     LIMITS (1,100)
00118 !
00119       RSTART=.FALSE.
00120       DO I=1,NPARAM
00121         IF (SETSTR%R(I).LT.1.D0) THEN
00122           SETSTR%R(I)=1.D0
00123           RSTART=.TRUE.
00124         ENDIF
00125         IF (SETSTR%R(I).GT.100.D0) THEN
00126           SETSTR%R(I)=100.D0
00127           RSTART=.TRUE.
00128         ENDIF
00129       ENDDO
00130 !
00131 !---------------------------------------------------------------------
00132 !
00133       RETURN
00134       END

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