Télécharger kveri2.eso

Retour à la liste

Numérotation des lignes :

  1. C KVERI2 SOURCE CB215821 16/04/21 21:17:38 8920
  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 = ABS(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.  
  65.  
  66.  
  67.  

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