Télécharger h8crsp.eso

Retour à la liste

Numérotation des lignes :

  1. C H8CRSP SOURCE CHAT 06/03/29 21:22:55 5360
  2. C
  3. SUBROUTINE H8CRSP( NB1,NB2,NB3, ITRNOE,
  4. > NBNMAX,NBE,NBN,iarr)
  5. C **********************************************************************
  6. C OBJET H8CRSP : CREATION DE LA TOPOLOGIE D'UNE GRILLE 3D (CAS SIMPLE)
  7. C
  8. C EN ENTREE :
  9. C NB1,NB2,NB3 : NOMBRE D'ELEMENTS DE LA GRILLE SUR LA LIGNE,
  10. C LA COLONNE, LA HAUTEUR
  11. C ITRNOE : TABLEAU DES ELEMENTS (A REMPLIR)
  12. C SA TAILLE DOIT ETRE DE (NB1*NB2*NB3*NBNMAX)
  13. C NBNMAX : NOMBRE MAXIMUM DE NOEUDS PAR ELEMENTS (8 OU +)
  14. C
  15. C EN SORTIE :
  16. C ITRNOE : DEFINITION DES ELEMENTS (HEXA)
  17. C NBE, NBN : NOMBRE D'ELEMENTS ET DE NOEUDS
  18. C iarr : CODE D'ERREUR
  19. C **********************************************************************
  20. IMPLICIT INTEGER(I-N)
  21. INTEGER NB1,NB2,NB3
  22. INTEGER ITRNOE(*),NBNMAX,NBE,NBN,iarr
  23. C
  24. INTEGER I,J,K,IP,IE
  25. C
  26. iarr = -1
  27. NBE = 0
  28. NBN = 0
  29. IF( NBNMAX.LT.4 )GOTO 999
  30. IF((NB1.LE.0).OR.(NB2.LE.0))GOTO 999
  31. iarr = 0
  32. NBE = NB1*NB2*NB3
  33. NBN = (NB1+1)*(NB2+1)*(NB3+1)
  34. DO 20 K=1,NB3
  35. DO 10 J=1,NB2
  36. DO 5 I=1,NB1
  37. IE = I + (J-1)*NB1+(K-1)*NB1*NB2
  38. IP = I + (J-1)*(NB1+1)+(K-1)*(NB1+1)*(NB2+1)
  39. C
  40. ITRNOE((IE-1)*NBNMAX+1)= IP
  41. ITRNOE((IE-1)*NBNMAX+2)= IP+1
  42. ITRNOE((IE-1)*NBNMAX+3)= IP+(NB1+1)+1
  43. ITRNOE((IE-1)*NBNMAX+4)= IP+(NB1+1)
  44. ITRNOE((IE-1)*NBNMAX+5)= IP+(NB1+1)*(NB2+1)
  45. ITRNOE((IE-1)*NBNMAX+6)= IP+(NB1+1)*(NB2+1)+ 1
  46. ITRNOE((IE-1)*NBNMAX+7)= IP+(NB1+1)*(NB2+1)+(NB1+1)+1
  47. ITRNOE((IE-1)*NBNMAX+8)= IP+(NB1+1)*(NB2+1)+(NB1+1)
  48. 5 CONTINUE
  49. 10 CONTINUE
  50. 20 CONTINUE
  51. 999 END
  52.  
  53.  
  54.  
  55.  

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