inter_curv_computation.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "com04_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	inter_curv_computation (icurv, nrtm_t, irect, c_max_local, curv_max, x)

Function/Subroutine Documentation

inter_curv_computation()

subroutine inter_curv_computation	(	integer, intent(in)	icurv,
		integer, intent(in)	nrtm_t,
		integer, dimension(4,nrtm_t), intent(in)	irect,
		intent(inout)	c_max_local,
		intent(inout)	curv_max,
		intent(in)	x )

Definition at line 28 of file inter_curv_computation.F.

!$COMMENT
!       INTER_CURV_COMPUTATION description :
!       
!       INTER_CURV_COMPUTATION :
!$ENDCOMMENT
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        INTEGER, INTENT(in) :: ICURV,NRTM_T
        INTEGER, DIMENSION(4,NRTM_T), INTENT(in) :: IRECT
        my_real, INTENT(inout) :: c_max_local
        my_real, DIMENSION(NRTM_T), INTENT(inout) :: curv_max
        my_real, DIMENSION(3,NUMNOD), INTENT(in) :: x
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
        INTEGER :: I
        my_real :: xxx,yyy,zzz 
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
C
        IF(icurv/=0)THEN
            c_max_local = zero
            DO i=1,nrtm_t
                xxx=max(x(1,irect(1,i)),x(1,irect(2,i)),
     .                  x(1,irect(3,i)),x(1,irect(4,i)))
     .               -min(x(1,irect(1,i)),x(1,irect(2,i)),
     .                  x(1,irect(3,i)),x(1,irect(4,i)))
                yyy=max(x(2,irect(1,i)),x(2,irect(2,i)),
     .                  x(2,irect(3,i)),x(2,irect(4,i)))
     .               -min(x(2,irect(1,i)),x(2,irect(2,i)),
     .                  x(2,irect(3,i)),x(2,irect(4,i)))
                zzz=max(x(3,irect(1,i)),x(3,irect(2,i)),
     .                  x(3,irect(3,i)),x(3,irect(4,i)))
     .               -min(x(3,irect(1,i)),x(3,irect(2,i)),
     .                  x(3,irect(3,i)),x(3,irect(4,i)))
                curv_max(i) = half * max(xxx,yyy,zzz)
                c_max_local = max(c_max_local,curv_max(i))
            ENDDO
        ELSE
            DO i=1,nrtm_t
                curv_max(i)=zero
            ENDDO
        ENDIF
C
        RETURN