ecran1.F File Reference

#include "implicit_f.inc"
#include "scr11_c.inc"
#include "com04_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	ecran1 (detonator, x, vdet)

Function/Subroutine Documentation

ecran1()

subroutine ecran1	(	type(detonator_wave_shaper_struct_), target	detonator,
			x,
		intent(in)	vdet )

Definition at line 35 of file ecran1.F.

C-----------------------------------------------
      USE message_mod
      USE detonators_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      my_real x(3,numnod)
      my_real,INTENT(IN)::vdet
      TYPE(DETONATOR_WAVE_SHAPER_STRUCT_),TARGET :: DETONATOR
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "scr11_c.inc"
#include      "com04_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, II, J, JJ, JJJ,NPE
      my_real dtomin, ddmx, d2
      INTEGER,POINTER,DIMENSION(:) :: IECR, IORDR, IFLG
      my_real,POINTER,DIMENSION(:) :: dtime
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
      INTEGER,EXTERNAL :: IOMBR
C-----------------------------------------------
 
      npe=detonator%NUMNOD
      iecr => detonator%NODES(1:)
      iordr => detonator%IORDR(1:)
      iflg => detonator%FLAG(1:)
      dtime => detonator%TIME(1:)
C---------------------------------------
C     COMPUTING DETONATION TIMES OF THE SCREEN POINTS
C---------------------------------------
C     1) LIGHTENING FROM DETONATION ORIGIN
C         check if the target can be directly lightened from detonation origin.
C            yes if no crossing with the screen (IOMBR returns 0)
C             no otherwise (IOMBR return 1)
C---------------------------------------
      dtomin=ep20
      iordr(1)=0
 
      DO i=1,npe
        ddmx  =ep20
        iflg(i)=0
        ii=iecr(i)
        yl=x(2,ii)
        zl=x(3,ii)
        IF(iombr(detonator,x,iecr,ddmx,vdet) == 0)THEN   !Vdet is optional velocity from origin (default is material velocity VDTO)
          dtime(i)=dto !DT0 is lightening time from function IOMBR
          IF(dtime(i) <= dtomin)THEN
            iordr(1)=i
            dtomin=dtime(i)
          ENDIF
        ELSE
          dtime(i)=ep20
        ENDIF
      END DO !I=1,NPE
 
      IF(iordr(1) == 0)THEN
        CALL ancmsg(msgid=422,
     .              msgtype=msgerror,
     .              anmode=aninfo)
      ENDIF
 
C-----------------------------------------
C     2) LIGHTENING SCREEN POINTS FROM SCREEN POINTS THEMSELVES
C-----------------------------------------
      DO j=1,npe-1
        jj=iordr(j)
        jjj=iecr(jj)
        !define J as a new source. (J from N-1, N-2, ..., 1 ; first point is the nearest one from detonation origin)
        yd=x(2,jjj)
        zd=x(3,jjj)
        dto0=dtime(jj)
        iflg(jj)=1
        dtomin=ep20
        DO i=1,npe
          ! check if I can be lightened from J (I from 1,2,...,N ; first node is the farthest one from detonation origin)
          IF(iflg(i) == 1)cycle !I
          ddmx  =ep20
          ii=iecr(i)
          yl=x(2,ii)
          zl=x(3,ii)
          IF(iabs(i-jj) == 1)THEN
            d2 =(yd-yl)**2+(zd-zl)**2
            dto=dto0+sqrt(d2)/vdto
            dtime(i)= min(dtime(i),dto)
          ELSEIF(iombr(detonator,x,iecr,ddmx,vdto) == 0)THEN
            dtime(i)= min(dtime(i),dto)
          ENDIF
          IF(dtime(i) <= dtomin)THEN
            iordr(j+1)=i
            dtomin=dtime(i)
          ENDIF
        END DO !I=1,NPE
      END DO !J=1,NPE-1
 
      RETURN