detcord0.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	detcord0 (detonator_cord, alt, x, vdet2, iopt)

Function/Subroutine Documentation

detcord0()

subroutine detcord0	(	type(detonator_cord_struct_)	detonator_cord,
			alt,
			x,
			vdet2,
		integer	iopt )

Definition at line 30 of file detcord0.F.

C-----------------------------------------------
C   M o d u l e s
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-----------------------------------------------
      INTEGER :: IOPT
      my_real :: alt,vdet2,x(3,*)
      TYPE(DETONATOR_CORD_STRUCT_)::DETONATOR_CORD
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: J,I,NOD_ID,NNOD
      my_real :: xlp1,xlp2,ylp1,ylp2,zlp1,zlp2,d,xl0,yl0,zl0
C-----------------------------------------------
C   P r e - C o n d i t i o n
C-----------------------------------------------
      IF(iopt == 0)    RETURN
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
 
 
       nnod = detonator_cord%NUMNOD
 
       IF(iopt == 1)THEN
        detonator_cord%TDET_PATH(1) = alt
        DO j=1,nnod-1
          !first point
          nod_id   = detonator_cord%NODES(j)
          xlp1     = x(1,nod_id)
          ylp1     = x(2,nod_id)
          zlp1     = x(3,nod_id)
          !second point
          nod_id   = detonator_cord%NODES(j+1)
          xlp2     = x(1,nod_id)
          ylp2     = x(2,nod_id)
          zlp2     = x(3,nod_id)
          !distance P1-P2
          xl0      = (xlp1-xlp2)                                       
          yl0      = (ylp1-ylp2)                                       
          zl0      = (zlp1-zlp2)        
          d        = xl0**2 + yl0**2 + zl0**2
          d        = sqrt(d)
          detonator_cord%TDET_PATH(j+1) = detonator_cord%TDET_PATH(j) + d/vdet2
        ENDDO
       ELSE
          detonator_cord%TDET_PATH(1:nnod) = alt
       ENDIF
 
 
C-----------------------------------------------
      RETURN