Télécharger geoqua.eso

Retour à la liste

Numérotation des lignes :

geoqua
  1. C GEOQUA SOURCE GOUNAND 21/06/02 21:16:14 11022
  2. SUBROUTINE GEOQUA(ITQUAF,
  3. $ JDIAMA,
  4. $ IMPR,IRET)
  5. IMPLICIT REAL*8 (A-H,O-Z)
  6. IMPLICIT INTEGER (I-N)
  7. C***********************************************************************
  8. C NOM : GEOQUA
  9. C DESCRIPTION :
  10. C
  11. * Calcul d'une propriété géométrique d'un QUAF régulier de côté 1 :
  12. * ici le diamètre du cercle circonscrit.
  13. * Cela sert pour le decentrement.
  14. C
  15. C LANGAGE : ESOPE
  16. C AUTEUR : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  17. C mél : gounand@semt2.smts.cea.fr
  18. C***********************************************************************
  19. C APPELES :
  20. C APPELES (E/S) :
  21. C APPELES (BLAS) :
  22. C APPELES (CALCUL) :
  23. C APPELE PAR :
  24. C***********************************************************************
  25. C SYNTAXE GIBIANE :
  26. C ENTREES :
  27. C ENTREES/SORTIES :
  28. C SORTIES :
  29. C CODE RETOUR (IRET) : = 0 si tout s'est bien passé
  30. C***********************************************************************
  31. C VERSION : v1, 04/10/2005, version initiale
  32. C HISTORIQUE : v1, 04/10/2005, création
  33. C HISTORIQUE :
  34. C HISTORIQUE :
  35. C***********************************************************************
  36. C Prière de PRENDRE LE TEMPS de compléter les commentaires
  37. C en cas de modification de ce sous-programme afin de faciliter
  38. C la maintenance !
  39. C***********************************************************************
  40.  
  41. -INC PPARAM
  42. -INC CCOPTIO
  43. -INC CCGEOME
  44. -INC TNLIN
  45. *-INC SMCHAEL
  46. INTEGER NBLIG,NBCOL,N2LIG,N2COL,NBPOI,NBELM
  47. POINTEUR JDIAMA.MCHEVA
  48. *
  49. CHARACTER*4 CQUAF
  50. LOGICAL LBID
  51. INTEGER IMPR,IRET
  52.  
  53. *
  54. * Executable statements
  55. *
  56. IF (IMPR.GT.2) WRITE(IOIMP,*) 'Entrée dans geoqua.eso'
  57. CQUAF=NOMS(ITQUAF)
  58. *
  59. IF (CQUAF.EQ.'SEG3') THEN
  60. XDIAMA=1.D0
  61. ELSEIF (CQUAF.EQ.'TRI7') THEN
  62. XDIAMA=2.D0/(SQRT(3.D0))
  63. ELSEIF (CQUAF.EQ.'QUA9') THEN
  64. XDIAMA=SQRT(2.D0)
  65. ELSEIF (CQUAF.EQ.'TE15') THEN
  66. XDIAMA=SQRT(3.D0/2.D0)
  67. ELSEIF (CQUAF.EQ.'PY19') THEN
  68. XDIAMA=SQRT(2.D0)
  69. ELSEIF (CQUAF.EQ.'PR21') THEN
  70. XDIAMA=SQRT(7.D0/3.D0)
  71. ELSEIF (CQUAF.EQ.'CU27') THEN
  72. XDIAMA=SQRT(3.D0)
  73. ELSE
  74. WRITE(IOIMP,*) 'Diametre interne de ',CQUAF,' non implemente'
  75. GOTO 9999
  76. ENDIF
  77. *
  78. NBLIG=1
  79. NBCOL=1
  80. N2LIG=1
  81. N2COL=1
  82. NBPOI=1
  83. NBELM=1
  84. SEGINI JDIAMA
  85. JDIAMA.WELCHE(1,1,1,1,1,1)=XDIAMA
  86. SEGDES JDIAMA
  87. *
  88. * Normal termination
  89. *
  90. IRET=0
  91. RETURN
  92. *
  93. * Format handling
  94. *
  95. 9999 CONTINUE
  96. IRET=1
  97. WRITE(IOIMP,*) 'An error was detected in subroutine geoqua'
  98. RETURN
  99. *
  100. * End of subroutine GEOQUA
  101. *
  102. END
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  

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