Télécharger la plaquette Cast3M


fr en

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