Télécharger cneq.eso

Retour à la liste

Numérotation des lignes :

  1. C CNEQ SOURCE CB215821 19/08/01 21:15:24 10279
  2. SUBROUTINE CNEQ
  3. *_______________________________________________________________________
  4. *
  5. * OPERATEUR DE CONVERSION DE CHAMPS VOLUMIQUES EN CHAMPS NODAUX
  6. *
  7. * ZOZO = CNEQ MOD CHP1 (CAR1) ;
  8. * CHAM
  9. *
  10. * MOD objet de type MMODEL
  11. * CHP1 objet de type CHPOINT
  12. * CHAM objet de type MCHAML
  13. * CAR1 CHAMP PAR ELEMENT DE CARACTERISTIQUES
  14. * GEOMETRIQUES ET MATERIELLES, TYPE MCHAML (facultatif)
  15. * ZOZO CHPOINT donnant les VALEURS nodales
  16. *
  17. * Francois CAFFIN 1992
  18. *_______________________________________________________________________
  19. *
  20. IMPLICIT INTEGER(I-N)
  21. IMPLICIT REAL*8(A-H,O-Z)
  22. -INC CCOPTIO
  23. -INC SMCHAML
  24. -INC SMCHPOI
  25. -INC SMCOORD
  26. *
  27. IPCHE1=0
  28. IPCHE2=0
  29. IPCHPO=0
  30. IPCHP2=0
  31. *
  32. * LECTURE DU MODELE
  33. *
  34. CALL LIROBJ('MMODEL ',IPMODL,1,IRT1)
  35. CALL ACTOBJ('MMODEL ',IPMODL,1)
  36. IF (IERR.NE.0) RETURN
  37. *
  38. * LECTURE D'UN CHPOINT OU D'UN MCHAML
  39. *
  40. CALL LIROBJ('CHPOINT ',IPCHPO,0,IRT1)
  41. IF (IERR.NE.0) RETURN
  42. IF (IRT1.EQ.1) THEN
  43. CALL ACTOBJ('CHPOINT ',IPCHPO,1)
  44. ELSE
  45. CALL LIROBJ('MCHAML ',IPIN,1,IRT1)
  46. CALL ACTOBJ('MCHAML ',IPIN,1)
  47. IF (IERR.NE.0) RETURN
  48.  
  49. CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  50. IF(IR .NE. 1) CALL ERREUR(KER)
  51. IF(IERR .NE. 0) RETURN
  52. ENDIF
  53. *
  54. * LECTURE EVENTUELLE D'UN CHAMP DE CARACTERISTIQUES
  55. *
  56. CALL LIROBJ('MCHAML ',IPIN,0,IRT1)
  57. IF (IERR .NE. 0) RETURN
  58. IPCHE2=0
  59. IF (IRT1 .EQ. 1) THEN
  60. CALL ACTOBJ('MCHAML ',IPIN,1)
  61. CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
  62. IF(IR .NE. 1) CALL ERREUR(KER)
  63. IF(IERR .NE. 0) RETURN
  64. ENDIF
  65. *
  66. CALL CNEQP(IPMODL,IPCHE1,IPCHPO,IPCHE2,IPCHP2,IRET)
  67. *
  68. IF (IRET.EQ.1) THEN
  69. *
  70. * on attribue une nature discrete au champ resultant
  71. MCHPOI = IPCHP2
  72. SEGACT MCHPOI*MOD
  73. JATTRI(1) = 2
  74.  
  75. CALL ACTOBJ('CHPOINT ',IPCHP2,1)
  76. CALL ECROBJ('CHPOINT ',IPCHP2)
  77. ENDIF
  78.  
  79. END
  80.  
  81.  
  82.  

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