Télécharger part6i.eso

Retour à la liste

Numérotation des lignes :

  1. C PART6I SOURCE CB215821 16/04/11 12:39:59 8899
  2. SUBROUTINE PART6I(ithr)
  3. C Cette subroutine est l'interface de PART6 et peut etre appelee en
  4. C parallele pour le partitionnement ARLEQUIN d'un MMODEL
  5.  
  6. C Creation 18/09/2015
  7. C Createur CB215821
  8. C Historique des Corrections apportees :
  9. C -
  10. C -
  11. C -
  12.  
  13. IMPLICIT INTEGER(I-N)
  14. IMPLICIT REAL*8 (A-H,O-Z)
  15.  
  16. -INC CCOPTIO
  17. -INC SMELEME
  18.  
  19. SEGMENT SSREF
  20. C NBZON : Nombre de ZONES demandees
  21. C NBSM : Nombre de IMODEL du MMODEL
  22. C ISREF : Tableau indiquant si la sous-zone contient ou non des elements
  23. C IMELE : Tableau avec les references des MELEME SIMPLES de chaque ZONE
  24. C IMODE : Tableau avec les references des IMODEL ELEMENTAIRES de chaque ZONE
  25. C NBELE : Tableau indiquant le nombre d'elements reellement present dans le MELEME SIMPLE en COURS
  26. C MTA : Tableau avec les references des MMODEL de chaque ZONE
  27. C IPOS : Tableau indiquant le nombre de IMODEL contenus dans les MMODEL de chaque ZONE
  28. INTEGER ISREF(NBZON,NBSM)
  29. INTEGER IMELE(NBZON,NBSM)
  30. INTEGER IMODE(NBZON,NBSM)
  31. INTEGER NBELE(NBZON)
  32. INTEGER MTA (NBZON)
  33. INTEGER IPOS (NBZON)
  34. ENDSEGMENT
  35.  
  36. C Declaration du COMMON pour le travail en parallele
  37. COMMON/part6c/NBTHR,SSREF,JA,IPT2,NBNN,NBELEM,ISOUS,NBZONE
  38.  
  39. DO 40 J = ithr,NBELEM,NBTHR
  40. C Boucle sur les elements du MELEME a partitionner par pas de SOUS-ZONES
  41. C Determination de la ZONE d'appartenance
  42. IZ=MOD(J+JA,NBZONE)
  43. IF (IZ .EQ. 0) THEN
  44. IZ = NBZONE
  45. ENDIF
  46.  
  47. C Determination du numero d'element dans la ZONE d'appartenance
  48. JF = ((J - 1) / NBZONE) + 1
  49.  
  50. C Chargement du bon MELEME SIMPLE dans la bonne ZONE
  51. IPT5 = IMELE(IZ,ISOUS)
  52. IPT5.ICOLOR(JF)=IPT2.ICOLOR(J)
  53.  
  54. C Boucle sur la connectivite des elements
  55. DO 50 II=1,IPT2.NUM(/1)
  56. IPT5.NUM(II,JF)=IPT2.NUM(II,J)
  57. 50 CONTINUE
  58.  
  59. C Enregistrement du nombre maximum d'element de la Sous Ref pour toutes les Zones
  60. IF (J .GT. (NBELEM - NBZONE)) THEN
  61. NBELE(IZ) = JF
  62. ENDIF
  63. 40 CONTINUE
  64.  
  65. RETURN
  66. END
  67.  
  68.  
  69.  
  70.  

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