aleso2.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	aleso2 (ale_connectivity, pm, ixq, igeo)

Function/Subroutine Documentation

aleso2()

subroutine aleso2	(	type(t_ale_connectivity), intent(inout)	ale_connectivity,
		dimension(npropm,nummat), intent(in)	pm,
		integer, dimension(7,numelq), intent(in)	ixq,
		integer, dimension(npropgi,numgeo), intent(in)	igeo )

Definition at line 30 of file aleso2.F.

C-----------------------------------------------
C   D e s c r i p t i o n
C-----------------------------------------------
C This subroutine is displaying in Starter listing file
C  Connectivities element/element (ALE)
C  specific case of quads (4-node-solid-elems in 2D) : buffer is IXQ contrary to IXTG for triangles
C Precondition : IPRI >= 3 (/IOFLAG Starter option)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE ale_connectivity_mod
      use element_mod , only : nixq
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      "com04_c.inc"
#include      "units_c.inc"
#include      "param_c.inc"
#include      "scr03_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER,INTENT(IN) :: IXQ(7,NUMELQ),IGEO(NPROPGI,NUMGEO)
      my_real,INTENT(IN) :: pm(npropm,nummat)
      TYPE(t_ale_connectivity), INTENT(INOUT) :: ALE_CONNECTIVITY
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER IVPR(4), M, IAL, I, IAD1, LGTH
      INTEGER JALE_FROM_MAT, JALE_FROM_PROP, JTHE_FROM_MAT
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
      IF(ipri <= 3)THEN
        WRITE(iout,1200)
        DO m=1,numelq
          jthe_from_mat = nint(pm(71,iabs(ixq(1,m)))) !/HEAT/MAT
          jale_from_mat = nint(pm(72,iabs(ixq(1,m)))) !/EULER/MAT or /ALE/MAT
          jale_from_prop = igeo(62,iabs(ixq(6,m)))    !/PROP/TYPE14 (IALE_FLAG)
          ial = jthe_from_mat + jale_from_mat + jale_from_prop
          IF(ial == 0)cycle
          iad1 = ale_connectivity%ee_connect%iad_connect(m)
          lgth = ale_connectivity%ee_connect%iad_connect(m+1)-ale_connectivity%ee_connect%iad_connect(m)        
          DO i=1,lgth
            ivpr(i)=ale_connectivity%ee_connect%connected(iad1 + i - 1)
            IF(ivpr(i) /= 0)ivpr(i)=ixq(7,ivpr(i))
          ENDDO
          WRITE(iout,'(5I10)')ixq(7,m),(ivpr(i),i=1,4)
        ENDDO
      ENDIF
C-----------------------------------------------
      RETURN
C-----------------------------------------------    
 1200 FORMAT(//
     & 6x,'A.L.E.  ELEMENT CONNECTIVITIES                ',/
     & 6x,'------------------------------                ',/
     & 6x,'ELEM      ELE1      ELE2      ELE3      ELE4  ')