dxyzsect.F File Reference

#include "implicit_f.inc"
#include "param_c.inc"
#include "com04_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	dxyzsect (nstrf, rwbuf, nprw, x, xmin, ymin, zmin, xmax, ymax, zmax, itab)

Function/Subroutine Documentation

dxyzsect()

subroutine dxyzsect	(	integer, dimension(*)	nstrf,
			rwbuf,
		integer, dimension(*)	nprw,
			x,
			xmin,
			ymin,
			zmin,
			xmax,
			ymax,
			zmax,
		integer, dimension(*)	itab )

Definition at line 33 of file dxyzsect.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "param_c.inc"
#include      "com04_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NSTRF(*),NPRW(*),ITAB(*)
      my_real
     .  rwbuf(nrwlp,*),x(3,*),xmin ,ymin ,zmin  ,xmax ,ymax, zmax
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER J, I, K, K0, K1, N, NSEG, N1, N2, N3, N4,MSR, ITYP
      my_real
     .   xx1, yy1, zz1, xx2, yy2, zz2, xx3, yy3, zz3,
     .   xx4, yy4, zz4, d13, xxc, yyc, zzc, al4,xwl,ywl,zwl,
     .   pmain,loc_proc
 
      my_real
     .   xsec(3,3,nsect)
      REAL R4,SBUF(3*NSECT)
CC-----------------------------------------------
      k1=1
C
       k1 = 33
       DO i=1,nsect
        n1 = nstrf(k1+1)
        n2 = nstrf(k1+2)
        n3 = nstrf(k1+3)
        xx1=x(1,n1)
        yy1=x(2,n1)
        zz1=x(3,n1)
        xx2=x(1,n2)
        yy2=x(2,n2)
        zz2=x(3,n2)
        xx3=x(1,n3)
        yy3=x(2,n3)
        zz3=x(3,n3)
        xx4=xx2-xx1
        yy4=yy2-yy1
        zz4=zz2-zz1
        al4=sqrt(xx4**2+yy4**2+zz4**2)
        xx4=xx4/max(al4,em20)
        yy4=yy4/max(al4,em20)
        zz4=zz4/max(al4,em20)
C
        d13=(xx3-xx1)*xx4+(yy3-yy1)*yy4+(zz3-zz1)*zz4
        xxc=xx1+d13*xx4
        yyc=yy1+d13*yy4
        zzc=zz1+d13*zz4
C
        r4 = xxc
        CALL write_r_c(r4,1)
        r4 = yyc
        CALL write_r_c(r4,1)
        r4 = zzc
        CALL write_r_c(r4,1)
C
 
         k1= nstrf(k1+22)+2
       ENDDO
C
        k=1
        DO n=1,nrwall
         n2=n +nrwall
         n3=n2+nrwall
         n4=n3+nrwall
         msr = nprw(n3)
         IF(msr==0)THEN
          xwl=rwbuf(4,n)
          ywl=rwbuf(5,n)
          zwl=rwbuf(6,n)
         ELSE
          xwl=x(1,msr)
          ywl=x(2,msr)
          zwl=x(3,msr)
         ENDIF
         ityp= nprw(n4)
         IF(ityp==4)THEN
           xwl = xwl + half*(rwbuf(7,n)+rwbuf(10,n))
           ywl = ywl + half*(rwbuf(8,n)+rwbuf(11,n))
           zwl = zwl + half*(rwbuf(9,n)+rwbuf(12,n)) 
         ENDIF
         k=k+nprw(n)
         IF(nprw(n4)==-1)k=k+nint(rwbuf(8,n))
         r4 = xwl
         CALL write_r_c(r4,1)
         r4 = ywl
         CALL write_r_c(r4,1)
         r4 = zwl
         CALL write_r_c(r4,1)
        ENDDO
C
      k=1
      DO n=1,nrwall
         n2=n +nrwall
         n3=n2+nrwall
         n4=n3+nrwall
         ityp= nprw(n4)
 
         IF(iabs(ityp)==1)THEN
            CALL dxwall(x,rwbuf(1,n),nprw(n3),xmin ,ymin ,
     .                  zmin  ,xmax ,ymax  , zmax)
         ELSEIF(ityp==2)THEN
            CALL dxwalc(x,rwbuf(1,n),nprw(n3),xmin ,ymin ,
     .                  zmin  ,xmax ,ymax  , zmax)
         ELSEIF(ityp==3)THEN
            CALL dxwals(x,rwbuf(1,n),nprw(n3))
         ELSEIF(ityp==4)THEN
            CALL dxwalp(x,rwbuf(1,n),nprw(n3))
         ENDIF
         k=k+nprw(n)
         IF(nprw(n4)==-1)k=k+nint(rwbuf(8,n))
      ENDDO
C
      RETURN