dxyzsect.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
!||====================================================================
!||    dxyzsect    ../starter/source/output/anim/dxyzsect.F
!||--- called by ------------------------------------------------------
!||    genani1     ../starter/source/output/anim/genani1.F
!||--- calls      -----------------------------------------------------
!||    dxwalc      ../starter/source/output/anim/dxwalc.F
!||    dxwall      ../starter/source/output/anim/dxwall.F
!||    dxwalp      ../starter/source/output/anim/dxwalp.F
!||    dxwals      ../starter/source/output/anim/dxwals.F
!||====================================================================
      SUBROUTINE dxyzsect (NSTRF,RWBUF,NPRW  ,X    , XMIN,
     2                     YMIN ,ZMIN ,XMAX ,YMAX  , ZMAX,
     3                     ITAB)
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
      END