Télécharger lekdiv.eso

Retour à la liste

Numérotation des lignes :

  1. C LEKDIV SOURCE BP208322 16/11/18 21:18:40 9177
  2. SUBROUTINE LEKDIV(MTABD,NU,MTABX,KINC,MCHELM,KPOIND)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. C***********************************************************************
  6. C
  7. C Ce sous programme lit le NUième coefficient d'un opérateur, stocké
  8. C à l'indice ARGnu de la table KIZX de pointeur MTABX associée à
  9. C l'opérateur en cours de traitement.
  10. C Le coefficient est de type ENTIER, FLOTTANT, POINT, CHPOINT ou MOT.
  11. C Dans le cas MOT, on récupère la donnée se trouvant à l'indice MOT
  12. C de la table INCO.
  13. C Ce Sp lit lit le coefficient et rend un MCHAML quoi qu'il arrive
  14. C Le MCHAML en retour est jetable il est calcule aux pts d'integrations
  15. C----------------------------------------------------------------------
  16. C HISTORIQUE : 20/10/01 : Création
  17. C
  18. C HISTORIQUE :
  19. C
  20. C
  21. C---------------------------
  22. C Paramètres Entrée/Sortie :
  23. C---------------------------
  24. C
  25. C E/ MTABD : Objet model de la zone
  26. C E/ NU : Rang du coefficient à aller chercher (ENTIER)
  27. C E/ MTABX : Pointeur de la table contenant le coefficient
  28. C E/ KINC : Pointeur de la table INCO dans laquelle on va chercher
  29. C les valeurs des indices lorsque ceux-ci sont des MOTS.
  30. C /S MCHELM : Chamelem pts d'intégration pour le COEF
  31. C E/ KPOIND : ENTIER type du support GÉométrique DUAL du shéma
  32. C d'intégration => pour la création du Chamelem
  33. C----------------------------------------------------------------------
  34. C KPOIN = 0->SOMMET 1-> FACE 2-> CENTRE 3-> CENTREP0 4-> CENTREP1 5-> MSOMMET
  35. C----------------------------------------------------------------------
  36. C************************************************************************
  37.  
  38. -INC SMCHAML
  39. -INC SMCHPOI
  40. -INC SMELEME
  41. -INC SMCOORD
  42. -INC CCOPTIO
  43. -INC CCGEOME
  44. CHARACTER*8 TYPE,NOM,NARG
  45. LOGICAL LOGI
  46. DIMENSION XPOI(3)
  47. C*****************************************************************************
  48. CLEKDIV
  49. C
  50. C- Récupération du TYPE du NUième argument de l'opérateur
  51. C
  52. c write(6,*)' DEBUT LEKDIV NU=',NU
  53.  
  54. IF (NU.LE.0 .OR. NU.GE.100) THEN
  55. INTERR(1) = NU
  56. INTERR(2) = 1
  57. INTERR(3) = 99
  58. CALL ERREUR(789)
  59. RETURN
  60. ENDIF
  61. IF (NU.LE.9) THEN
  62. WRITE(NARG,FMT='(A3,I1)')'ARG',Nu
  63. ELSE
  64. WRITE(NARG,FMT='(A3,I2)')'ARG',Nu
  65. ENDIF
  66. TYPE = ' '
  67. CALL LENCHA(NARG,LC)
  68. CALL ACCTAB(MTABX,'MOT', 0,0.D0,NARG(1:LC),LOGI, 0,
  69. & TYPE ,IVAL,COEF, NOM,LOGI,MCHPOI)
  70.  
  71. c write(6,*)' TYPE ,MCHPOI =',TYPE ,MCHPOI
  72. IF(TYPE.EQ.'MOT')THEN
  73. TYPE=' '
  74. CALL ACMO(KINC,NOM,TYPE,MCHPOI)
  75. ENDIF
  76.  
  77.  
  78. IF(TYPE.NE.'CHPOINT ')THEN
  79. C On ne veut pas d'objet de type %m1:8
  80. MOTERR(1: 8) = TYPE
  81. CALL ERREUR(39)
  82. RETURN
  83. ENDIF
  84.  
  85. CALL MELDIV(MTABD,MCHPOI,MCHELM,KPOIND)
  86.  
  87. c write(6,*)' FIN LEKDIV '
  88. RETURN
  89. 1001 FORMAT(20(1X,I5))
  90. 1002 FORMAT(10(1X,1PE11.4))
  91. END
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  

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