Télécharger cneq.eso

Retour à la liste

Numérotation des lignes :

cneq
  1. C CNEQ SOURCE CB215821 23/01/25 21:15:08 11573
  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.  
  23. -INC PPARAM
  24. -INC CCOPTIO
  25. -INC SMCHAML
  26. -INC SMCHPOI
  27. -INC SMCOORD
  28. *
  29. IPCHE1=0
  30. IPCHE2=0
  31. IPCHPO=0
  32. IPCHP2=0
  33. *
  34. * LECTURE DU MODELE
  35. *
  36. CALL LIROBJ('MMODEL ',IPMODL,1,IRT1)
  37. CALL ACTOBJ('MMODEL ',IPMODL,1)
  38. IF (IERR.NE.0) RETURN
  39. *
  40. * LECTURE D'UN CHPOINT OU D'UN MCHAML
  41. *
  42. CALL LIROBJ('CHPOINT ',IPCHPO,0,IRT1)
  43. IF (IERR.NE.0) RETURN
  44. IF (IRT1.EQ.1) THEN
  45. CALL ACTOBJ('CHPOINT ',IPCHPO,1)
  46. ELSE
  47. CALL LIROBJ('MCHAML ',IPIN,1,IRT1)
  48. CALL ACTOBJ('MCHAML ',IPIN,1)
  49. IF (IERR.NE.0) RETURN
  50.  
  51. CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  52. IF(IR .NE. 1) CALL ERREUR(KER)
  53. IF(IERR .NE. 0) RETURN
  54. ENDIF
  55. *
  56. * LECTURE EVENTUELLE D'UN CHAMP DE CARACTERISTIQUES
  57. *
  58. CALL LIROBJ('MCHAML ',IPIN,0,IRT1)
  59. IF (IERR .NE. 0) RETURN
  60. IPCHE2=0
  61. IF (IRT1 .EQ. 1) THEN
  62. CALL ACTOBJ('MCHAML ',IPIN,1)
  63. CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
  64. IF(IR .NE. 1) CALL ERREUR(KER)
  65. IF(IERR .NE. 0) RETURN
  66. ENDIF
  67. *
  68. SEGACT,MCOORD
  69. CALL CNEQP(IPMODL,IPCHE1,IPCHPO,IPCHE2,IPCHP2,IRET)
  70. SEGDES,MCOORD
  71. *
  72. IF (IRET.EQ.1) THEN
  73. *
  74. * on attribue une nature discrete au champ resultant
  75. MCHPOI = IPCHP2
  76. SEGACT MCHPOI*MOD
  77. JATTRI(1) = 2
  78.  
  79. CALL ACTOBJ('CHPOINT ',IPCHP2,1)
  80. CALL ECROBJ('CHPOINT ',IPCHP2)
  81. ENDIF
  82.  
  83. END
  84.  
  85.  
  86.  
  87.  
  88.  

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