dxyzsect.F File Reference

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

Go to the source code of this file.

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

Function/Subroutine Documentation

dxyzsect()

subroutine dxyzsect	(	integer, dimension(*)	nstrf,
			rwbuf,
		integer, dimension(*)	nprw,
			x,
			xmin,
			ymin,
			zmin,
			xmax,
			ymax,
			zmax,
		integer, dimension(nspmd+1,*)	fr_sec,
		integer, dimension(nspmd+2,*)	fr_wall,
		integer, dimension(*)	weight,
		integer, dimension(*)	itab )

Definition at line 38 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"
#include      "task_c.inc"
#include      "com01_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
      integer
     .   fr_sec(nspmd+1,*),fr_wall(nspmd+2,*),weight(*)
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-----------------------------------------------
      loc_proc=ispmd+1
      k1=1
      IF (nspmd == 1) THEN
        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
      ELSE
        IF (nsect > 0) CALL spmd_aget_sect(nstrf,x,xsec,weight,itab)
 
          IF (ispmd==0) THEN
 
            DO i=1,nsect
 
             xx1 = xsec(1,1,i)
             yy1 = xsec(1,2,i)
             zz1 = xsec(1,3,i)
 
             xx2 = xsec(2,1,i)
             yy2 = xsec(2,2,i)
             zz2 = xsec(2,3,i)
 
             xx3 = xsec(3,1,i)
             yy3 = xsec(3,2,i)
             zz3 = xsec(3,3,i)
C  
             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)
            ENDDO
          ENDIF
      ENDIF
C
      k=1
      DO n=1,nrwall
        n2=n +nrwall
        n3=n2+nrwall
        n4=n3+nrwall
        msr = nprw(n3)
        IF (nspmd == 1) THEN
          IF(msr==0)THEN
            xwl=rwbuf(4,n)
            ywl=rwbuf(5,n)
            zwl=rwbuf(6,n)
          ELSE
C verifier que ce noeud est sur proc0 !
            xwl=x(1,msr)
            ywl=x(2,msr)
            zwl=x(3,msr)
          ENDIF
        ELSE
          CALL spmd_agetmsr(fr_wall(1,n),x,msr,xwl,ywl,zwl,rwbuf(1,n))
        END IF
        IF (ispmd==0) THEN
          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)
        END IF
      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,fr_wall(1,n))
        ELSEIF(ityp==2)THEN
           CALL dxwalc(x,rwbuf(1,n),nprw(n3),xmin ,ymin ,
     .                 zmin  ,xmax ,ymax  , zmax,fr_wall(1,n))
        ELSEIF(ityp==3)THEN
           CALL dxwals(x,rwbuf(1,n),nprw(n3),fr_wall(1,n))
        ELSEIF(ityp==4)THEN
           CALL dxwalp(x,rwbuf(1,n),nprw(n3),fr_wall(1,n))
        ENDIF
        k=k+nprw(n)
        IF(nprw(n4)==-1)k=k+nint(rwbuf(8,n))
      ENDDO
C
      RETURN