1 : C KVERI2 SOURCE CHAT 05/01/13 01:11:28 5004 2 : SUBROUTINE KVERI2 (SHC2D,SKBUF2,SKRESO) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8 (A-H,O-Z) 5 : C 6 : C VERIFICATION DU BILAN AUTOUR D'UNE DIRECTION DONNEE 7 : C 8 : C----------------------------------------------------------------------- 9 : SEGMENT SKRESO 10 : INTEGER KFC,NRES,KES,KIMP 11 : ENDSEGMENT 12 : C KFC : NOMBRE DE FACES H.C 13 : C NRES: RESOLUTION 14 : C KES : DIM ESPACE 15 : C KIMP: IMPRESSION 16 : C----------------------------------------------------------------------- 17 : C----------------------------------------------------------------------- 18 : SEGMENT SHC2D 19 : INTEGER IR(NR),KA(NFC),IM(NFC,NFC) 20 : INTEGER KRO(NFC,NES),KSI(NFC,NES) 21 : REAL*8 V(NES,NR),G(NR) 22 : ENDSEGMENT 23 : 24 : C DESCRIPTION DU H.C DE PROJECTION 25 : C -------------------------------- 26 : C V : DIRECTION UNITAIRE DES CELLULES 27 : C G : FACTEUR DE FORME ASSOCIE 28 : C IR: CORRESPONDANCE 29 : C KRO , KSI : POUR LE CHANGEMENT DE REPERE 30 : C IM : REFERENCE 31 : C NR : RESOLUTION 32 : C NFC : NOMBRE DE FACES 33 : C----------------------------------------------------------------------- 34 : SEGMENT SKBUF2 35 : INTEGER NUMF(NFC,NOC,NR),NTYP(NFC,NR) 36 : REAL*8 ZB(NFC,NR),PSC(NFC,NR) 37 : ENDSEGMENT 38 : C 39 : C BUFFER ASSOCIE AU H.C 40 : C --------------------- 41 : C NUMF : INDICE DE LE DERNIERE FACE RENCONTREE 42 : C NTYP : TYPES ASSOCIES 43 : C ZB : PROFONDEUR 44 : C PSC : PRODUIT SCALAIRE (NORMALE.DIRECTION CELLULE) 45 : C----------------------------------------------------------------------- 46 : 47 : F = 0 48 : NC = 0 49 : DO 500 NF = 1,KFC 50 : DO 501 I = 1,NRES 51 : IF (PSC(NF,I).GT.-1.) THEN 52 : NC = NC + 1 53 : F = F + PSC(NF,I)*G(I) 54 : ENDIF 55 : 501 CONTINUE 56 : 500 CONTINUE 57 : ERR = DABS(1.-F) 58 : IF (ERR.GT.0.01.OR.KIMP.GE.4) THEN 59 : WRITE(6,1000) NC,F 60 : 1000 FORMAT(1X,'VERIF: NCEL ',I6,' F ',F10.5) 61 : ENDIF 62 : RETURN 63 : END 64 :
© Cast3M 2003 - All rights reserved.
Disclaimer