s8sav12.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "scr05_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	s8sav12 (offg, offg0, sav, x, xdp, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8, nel, jcvt)

Function/Subroutine Documentation

s8sav12()

subroutine s8sav12	(		offg,
			offg0,
		double precision, dimension(nel,21)	sav,
			x,
		double precision, dimension(3,*)	xdp,
		integer, dimension(*)	nc1,
		integer, dimension(*)	nc2,
		integer, dimension(*)	nc3,
		integer, dimension(*)	nc4,
		integer, dimension(*)	nc5,
		integer, dimension(*)	nc6,
		integer, dimension(*)	nc7,
		integer, dimension(*)	nc8,
		integer	nel,
		integer, intent(in)	jcvt )

Definition at line 31 of file s8sav12.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "scr05_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: JCVT
      INTEGER NEL
C     REAL
      INTEGER 
     .        NC1(*),NC2(*),NC3(*),NC4(*),NC5(*),NC6(*),NC7(*),NC8(*)
      my_real
     .   offg(*),offg0(*),x(3,*)
      DOUBLE PRECISION 
     .   XDP(3,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I ,IH    
      DOUBLE PRECISION 
     .   XD1(MVSIZ), XD2(MVSIZ), XD3(MVSIZ), XD4(MVSIZ),
     .   XD5(MVSIZ), XD6(MVSIZ), XD7(MVSIZ), XD8(MVSIZ),
     .   YD1(MVSIZ), YD2(MVSIZ), YD3(MVSIZ), YD4(MVSIZ),
     .   YD5(MVSIZ), YD6(MVSIZ), YD7(MVSIZ), YD8(MVSIZ),
     .   ZD1(MVSIZ), ZD2(MVSIZ), ZD3(MVSIZ), ZD4(MVSIZ),
     .   ZD5(MVSIZ), ZD6(MVSIZ), ZD7(MVSIZ), ZD8(MVSIZ), SAV(NEL,21)
C=======================================================================
C-----Always in global system and fixed then (not like Ismstr=11
C-----for JCVT>0 Fixed SAV is rotated w/ local system not like Ismstr=1
        ih=0
        DO i=1,nel
         IF(offg(i)/=offg0(i).AND.abs(offg(i)) > one ) ih=1
          END DO
          IF (ih==0) RETURN
C
      IF(iresp==1)THEN 
        DO i=1,nel
          xd1(i)=xdp(1,nc1(i))  
          yd1(i)=xdp(2,nc1(i))  
          zd1(i)=xdp(3,nc1(i))  
          xd2(i)=xdp(1,nc2(i))  
          yd2(i)=xdp(2,nc2(i))  
          zd2(i)=xdp(3,nc2(i))  
          xd3(i)=xdp(1,nc3(i))  
          yd3(i)=xdp(2,nc3(i))  
          zd3(i)=xdp(3,nc3(i))  
          xd4(i)=xdp(1,nc4(i))  
          yd4(i)=xdp(2,nc4(i))  
          zd4(i)=xdp(3,nc4(i))  
          xd5(i)=xdp(1,nc5(i))  
          yd5(i)=xdp(2,nc5(i))  
          zd5(i)=xdp(3,nc5(i))  
          xd6(i)=xdp(1,nc6(i))  
          yd6(i)=xdp(2,nc6(i))  
          zd6(i)=xdp(3,nc6(i))  
          xd7(i)=xdp(1,nc7(i))  
          yd7(i)=xdp(2,nc7(i))  
          zd7(i)=xdp(3,nc7(i))  
          xd8(i)=xdp(1,nc8(i))  
          yd8(i)=xdp(2,nc8(i))  
          zd8(i)=xdp(3,nc8(i))              
        ENDDO
      ELSE
        DO i=1,nel
          xd1(i)=x(1,nc1(i))  
          yd1(i)=x(2,nc1(i))  
          zd1(i)=x(3,nc1(i))  
          xd2(i)=x(1,nc2(i))  
          yd2(i)=x(2,nc2(i))  
          zd2(i)=x(3,nc2(i))  
          xd3(i)=x(1,nc3(i))  
          yd3(i)=x(2,nc3(i))  
          zd3(i)=x(3,nc3(i))  
          xd4(i)=x(1,nc4(i))  
          yd4(i)=x(2,nc4(i))  
          zd4(i)=x(3,nc4(i))  
          xd5(i)=x(1,nc5(i))  
          yd5(i)=x(2,nc5(i))  
          zd5(i)=x(3,nc5(i))  
          xd6(i)=x(1,nc6(i))  
          yd6(i)=x(2,nc6(i))  
          zd6(i)=x(3,nc6(i))  
          xd7(i)=x(1,nc7(i))  
          yd7(i)=x(2,nc7(i))  
          zd7(i)=x(3,nc7(i))  
          xd8(i)=x(1,nc8(i))  
          yd8(i)=x(2,nc8(i))  
          zd8(i)=x(3,nc8(i))  
        ENDDO      
      ENDIF
C   -----attention, format is different to Ismstr=10,11
        DO i=1,nel
          IF(offg(i)/=offg0(i).AND.abs(offg(i)) > one )THEN
           sav(i,1)=xd1(i)-xd8(i)
           sav(i,2)=yd1(i)-yd8(i)
           sav(i,3)=zd1(i)-zd8(i)
           sav(i,4)=xd2(i)-xd8(i)
           sav(i,5)=yd2(i)-yd8(i)
           sav(i,6)=zd2(i)-zd8(i)
           sav(i,7)=xd3(i)-xd8(i)
           sav(i,8)=yd3(i)-yd8(i)
           sav(i,9)=zd3(i)-zd8(i)
           sav(i,10)=xd4(i)-xd8(i)
           sav(i,11)=yd4(i)-yd8(i)
           sav(i,12)=zd4(i)-zd8(i)
           sav(i,13)=xd5(i)-xd8(i)
           sav(i,14)=yd5(i)-yd8(i)
           sav(i,15)=zd5(i)-zd8(i)
           sav(i,16)=xd6(i)-xd8(i)
           sav(i,17)=yd6(i)-yd8(i)
           sav(i,18)=zd6(i)-zd8(i)
           sav(i,19)=xd7(i)-xd8(i)
           sav(i,20)=yd7(i)-yd8(i)
           sav(i,21)=zd7(i)-zd8(i)
          ENDIF
        ENDDO
C-----------
      RETURN