Sources Esope | Cast3M

morayo
C MORAYO    SOURCE    MB234859  26/11/04    21:15:08     12691          
      SUBROUTINE MORAYO(IPT1,IPT2,IPT3)
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
 
-INC PPARAM
-INC CCOPTIO
-INC SMCOORD
-INC SMELEME
-INC SMLENTI
-INC CCGEOME
C
      LOGICAL bRACC,bLIAI
C
      CALL ACTOBJ('MAILLAGE',IPT1,1)
C
      NBMAI1=IPT1.LISOUS(/1)
      NIPT1 =NBMAI1
      IF (NBMAI1.EQ.0) THEN
        NIPT1=1
      ELSE
        NBNN=0
        NBELEM=0
        NBREF=0
        NBSOUS=NBMAI1
        SEGINI,IPT2,IPT3
      ENDIF
C
      IPT4=IPT1
      DO 1 IMAI=1,NIPT1
C
        IF (NBMAI1.NE.0) IPT4=IPT1.LISOUS(IMAI)
C
        NBN4=IPT4.NUM(/1)
        NBNN=NBN4/2
        NBELEM=IPT4.NUM(/2)
        NBREF=0
        NBSOUS=0
        SEGINI,IPT5,IPT6
        IPT5.ITYPEL=NBNN
        IPT6.ITYPEL=NBNN
        bRACC=((IPT4.ITYPEL.EQ.12).OR. (IPT4.ITYPEL.EQ.13))
        bLIAI=((IPT4.ITYPEL.GE.18).AND.(IPT4.ITYPEL.LE.21))
        IF (.NOT.(bRACC.OR.bLIAI)) THEN
          CALL ERREUR(16)
          RETURN
        ENDIF
        DO IELT=1,NBELEM
          DO INOE=1,NBNN
            INO2=NBNN+INOE
            IF (bRACC) INO2=NBN4-INOE+1
            IPT5.NUM(INOE,IELT)=IPT4.NUM(INOE,IELT)
            IPT6.NUM(INOE,IELT)=IPT4.NUM(INO2,IELT)
          ENDDO
        ENDDO
C
        IF (NBMAI1.EQ.0) THEN
          IPT2=IPT5
          IPT3=IPT6
        ELSE
          IPT2.LISOUS(IMAI)=IPT5
          IPT3.LISOUS(IMAI)=IPT6
        ENDIF
C
 1    CONTINUE
C
      END