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 : -INC CCOPTIO 37 : -INC SMCHAML 38 : * 39 : IRET=0 40 : MELVA1=IELVA1 41 : MELVA2=IELVA2 42 : SEGACT MELVA1 43 : SEGACT MELVA2 44 : * 45 : NBP1=MELVA1.VELCHE(/1) 46 : NBP2=MELVA2.VELCHE(/1) 47 : NEL1=MELVA1.VELCHE(/2) 48 : NEL2=MELVA2.VELCHE(/2) 49 : N1PTEL=MAX(NBP1,NBP2) 50 : N1EL =MAX(NEL1,NEL2) 51 : N2PTEL=0 52 : N2EL =0 53 : SEGINI MELVAL 54 : * 55 : DO 5 IB=1,N1EL 56 : IBMN1=MIN(IB ,MELVA1.VELCHE(/2)) 57 : IBMN2=MIN(IB ,MELVA2.VELCHE(/2)) 58 : * 59 : DO 6 IGAU=1,N1PTEL 60 : IGMN1=MIN(IGAU,MELVA1.VELCHE(/1)) 61 : IGMN2=MIN(IGAU,MELVA2.VELCHE(/1)) 62 : * 63 : * 64 : VELCHE(IGAU,IB)=MELVA1.VELCHE(IGMN1,IBMN1)- 65 : & MELVA2.VELCHE(IGMN2,IBMN2) 66 : MOYS=VELCHE(IGAU,IB)/N1PTEL+MOYS 67 : 6 CONTINUE 68 : 5 CONTINUE 69 : IMELVA=MELVAL 70 : SEGDES MELVAL 71 : SEGDES MELVA1 72 : SEGDES MELVA2 73 : RETURN 74 : END 75 : 76 : 77 :
© Cast3M 2003 - All rights reserved.
Disclaimer