Télécharger genrd1.eso

Retour à la liste

Numérotation des lignes :

genrd1
  1. C GENRD1 SOURCE CHAT 05/01/13 00:17:37 5004
  2. SUBROUTINE GENRD1(FAC1,FAC2,NUM2,NPTO,TOL2,IRET)
  3. C--------------------------------------------------------------------
  4. C
  5. C ON ORDONNE FAC2 SUIVANT FAC1 EN FAISANT SUIVRE LES PERMUTATIONS
  6. C SUR NUM2
  7. C
  8. C PP 9/97
  9. C Pierre Pegon/JRC Ispra
  10. C--------------------------------------------------------------------
  11. IMPLICIT INTEGER(I-N)
  12. IMPLICIT REAL*8 (A-H,O-Z)
  13. C
  14.  
  15. -INC PPARAM
  16. -INC CCOPTIO
  17. C
  18. DIMENSION FAC1(3,NPTO),FAC2(3,NPTO),NUM2(NPTO)
  19. C
  20. DO IE1=1,NPTO-1
  21. DO IE2=IE1,NPTO
  22. DIST=0.D0
  23. DO IE3=1,3
  24. DIST=DIST+(FAC1(IE3,IE1)-FAC2(IE3,IE2))**2
  25. ENDDO
  26. IF(DIST.LT.TOL2)GOTO 1
  27. ENDDO
  28. IRET=IRET+1
  29. WRITE(IOIMP,*)'GENRD1: Impossible de ce retrouver la!'
  30. GOTO 2
  31. 1 CONTINUE
  32. IF(IE2.NE.IE1)THEN
  33. DO IE3=1,3
  34. PDUM=FAC2(IE3,IE2)
  35. FAC2(IE3,IE2)=FAC2(IE3,IE1)
  36. FAC2(IE3,IE1)=PDUM
  37. ENDDO
  38. NDUM=NUM2(IE2)
  39. NUM2(IE2)=NUM2(IE1)
  40. NUM2(IE1)=NDUM
  41. ENDIF
  42. ENDDO
  43. C
  44. 2 CONTINUE
  45. RETURN
  46. END
  47.  
  48.  
  49.  

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