Télécharger la plaquette Cast3M


fr en

Télécharger sensi3.eso

Retour à la liste

Numérotation des lignes : 

sensi3

  1. C SENSI3    SOURCE    PV        09/03/12    21:34:16     6325
  2. *
  3. *      BOUCLE SUR LES POINTS DE GAUSS DE L'ELEMENT DE LA SOUS ZONE
  4. *      EXTRAIT DE ADCHVE.ESO
  5. *      LELEU DIDIER HAZE FREDERIC
  6. *      03/03/93
  7. *
  8. *
  9.        SUBROUTINE SENSI3(IELVA1,IELVA2,IMELVA,IRET,MOYS)
  10. *
  11. *
  12. *
  13. *
  14. *
  15. *   ENTREE :
  16. *   ________
  17. *
  18. *       IELVA1       POINTEUR SUR LE MELVAL DU 1IER CHAMELEM
  19. *       IELVA2       POINTEUR SUR LE MELVA DU 2IEME CHAMELEM
  20. *
  21. *
  22. *   SORTIES :
  23. *   ---------
  24. *
  25. *       IMELVA       POINTEUR SUR LE MELVAL RESULTAT
  26. *       IRET         = 0 SI OK / =1 SINON
  27. *                    =104 POUR DECLENCHER LE MESSAGE D'ERREUR 104
  28. *                    =197 POUR DECLENCHER LE MESSAGE D'ERREUR 197
  29. *       MOYS         VALEUR DE LA SENSIBILITE (MOYENNE)
  30. *
  31. *
  32.       IMPLICIT INTEGER(I-N)
  33.       IMPLICIT REAL*8(A-H,O-Z)
  34.       REAL*8 MOYS
  35. *
  36.  
  37. -INC PPARAM
  38. -INC CCOPTIO
  39. -INC SMCHAML
  40. *
  41.       IRET=0
  42.       MELVA1=IELVA1
  43.       MELVA2=IELVA2
  44.       SEGACT MELVA1
  45.       SEGACT MELVA2
  46. *
  47.       NBP1=MELVA1.VELCHE(/1)
  48.       NBP2=MELVA2.VELCHE(/1)
  49.       NEL1=MELVA1.VELCHE(/2)
  50.       NEL2=MELVA2.VELCHE(/2)
  51.       N1PTEL=MAX(NBP1,NBP2)
  52.       N1EL  =MAX(NEL1,NEL2)
  53.       N2PTEL=0
  54.       N2EL  =0
  55.       SEGINI MELVAL
  56. *
  57.          DO 5 IB=1,N1EL
  58.             IBMN1=MIN(IB  ,MELVA1.VELCHE(/2))
  59.             IBMN2=MIN(IB  ,MELVA2.VELCHE(/2))
  60. *
  61.             DO 6 IGAU=1,N1PTEL
  62.                IGMN1=MIN(IGAU,MELVA1.VELCHE(/1))
  63.                IGMN2=MIN(IGAU,MELVA2.VELCHE(/1))
  64. *
  65. *
  66.                VELCHE(IGAU,IB)=MELVA1.VELCHE(IGMN1,IBMN1)-
  67.      &                         MELVA2.VELCHE(IGMN2,IBMN2)
  68.                MOYS=VELCHE(IGAU,IB)/N1PTEL+MOYS
  69.   6      CONTINUE
  70.   5      CONTINUE
  71.          IMELVA=MELVAL
  72.       SEGDES MELVAL
  73.       SEGDES MELVA1
  74.       SEGDES MELVA2
  75.       RETURN
  76.       END
  77.  
  78.  
  79.  
  80.  

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