Télécharger la plaquette Cast3M


fr en

Télécharger pilep2.eso

Retour à la liste

Numérotation des lignes : 

pilep2

  1. C PILEP2    SOURCE    PV        09/03/12    21:30:10     6325
  2.        SUBROUTINE PILEP2(IELVA1,IELVA2,XX,IRET,CRIT)
  3. *
  4.       IMPLICIT INTEGER(I-N)
  5.       IMPLICIT REAL*8(A-H,O-Z)
  6. *
  7.  
  8. -INC PPARAM
  9. -INC CCOPTIO
  10. -INC SMCHAML
  11. -INC SMLREEL
  12. -INC SMCOORD
  13. *
  14.       IRET=0
  15.       MELVA1=IELVA1
  16.       MELVA2=IELVA2
  17.       SEGACT MELVA1
  18.       SEGACT MELVA2
  19. *
  20.       NBP1=MELVA1.VELCHE(/1)
  21.       IF (NBP1.EQ.0) THEN
  22.         IRET = 197
  23.       ELSE
  24.          NBP2=MELVA2.VELCHE(/1)
  25.          NEL1=MELVA1.VELCHE(/2)
  26.          NEL2=MELVA2.VELCHE(/2)
  27.          N1PTEL=MAX(NBP1,NBP2)
  28.          N1EL  =MAX(NEL1,NEL2)
  29.          N2PTEL=0
  30.          N2EL  =0
  31.          DO 5 IGAU=1,N1PTEL
  32.             IGMN1=MIN(IGAU,MELVA1.VELCHE(/1))
  33.             IGMN2=MIN(IGAU,MELVA2.VELCHE(/1))
  34.             XLA = 1.E30
  35.             DO 5 IB=1,N1EL
  36. *
  37.                IBMN1=MIN(IB  ,MELVA1.VELCHE(/2))
  38.                IBMN2=MIN(IB  ,MELVA2.VELCHE(/2))
  39. *
  40.         IF ( MELVA2.VELCHE(IGMN2,IBMN2) . NE . 0.)
  41.      *  XLA =(SIGN(CRIT,MELVA2.VELCHE(IGMN2,IBMN2))
  42.      *          - MELVA1.VELCHE(IGMN1,IBMN1))/
  43.      *                        MELVA2.VELCHE(IGMN2,IBMN2)
  44.          IF( XLA . LE . XX ) XX = XLA
  45. *
  46.   5      CONTINUE
  47.       ENDIF
  48.       SEGDES MELVA1
  49.       SEGDES MELVA2
  50. *
  51.       RETURN
  52.       END
  53.  
  54.  
  55.  
  56.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales