ecran2.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	ecran2 (detonator_wave_shaper, x, vdet)

Function/Subroutine Documentation

ecran2()

subroutine ecran2	(	type(detonator_wave_shaper_struct_), target	detonator_wave_shaper,
			x,
		intent(in)	vdet )

Definition at line 33 of file ecran2.F.

C-----------------------------------------------
      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_WAVE_SHAPER
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,POINTER,DIMENSION(:) :: IECR, IORDR, IFLG
      my_real,POINTER,DIMENSION(:) ::  dtime
      INTEGER J, I, II, NPE
      my_real ddmx, dtos
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_wave_shaper%NUMNOD
      iecr => detonator_wave_shaper%NODES(1:npe)
      iordr => detonator_wave_shaper%IORDR(1:npe)
      dtime => detonator_wave_shaper%TIME(1:npe)
 
C---------------------------------------
C     COMPUTING DETONATION:TIME FOR ONE ELEMENT (CENTROID)
C---------------------------------------
C     1) DETONATION ORIGIN: USER DETONATION POINT
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---------------------------------------
      ddmx=ep20
      IF(iombr(detonator_wave_shaper,x,iecr,ddmx,vdet) == 0)RETURN   !Vdet is optional velocity from origin (default is material velocity VDTO)
 
C-----------------------------------------
C     2) DETONATION ORIGIN:EACH POINT OF THE SCREEN PATH
C-----------------------------------------
      dtos=ep20
 
      npe=detonator_wave_shaper%NUMNOD
      DO j=1,npe
        i=iordr(j)
        ii=iecr(i)
        yd=x(2,ii)
        zd=x(3,ii)
        dto0=dtime(i)
        IF(iombr(detonator_wave_shaper,x,iecr,ddmx,vdto) == 0)THEN
          dtos= min(dtos,dto) !keep minimum value (first arrival time)
        ENDIF
      END DO !J=1,NPE
 
      dto=dtos
 
      RETURN