Télécharger xpost.eso

Retour à la liste

Numérotation des lignes :

xpost
  1. C XPOST SOURCE CB215821 24/04/12 21:17:31 11897
  2. C
  3. SUBROUTINE XPOST
  4. C
  5. C Post traitement des calculs avec elements de type X-FEM necessitant
  6. C une attention particuliere
  7. C
  8. IMPLICIT REAL*8 (A-H,O-Z)
  9. C
  10. C
  11. C SEGMENTS INCLUDE
  12.  
  13. -INC PPARAM
  14. -INC CCOPTIO
  15. -INC SMMODEL
  16. -INC SMCHPOI
  17. -INC SMCHAML
  18. C
  19. C
  20. PARAMETER (NMOT0=2)
  21. CHARACTER*4 MOTCLE(NMOT0)
  22. DATA MOTCLE/'RECO','FISS'/
  23. * DATA MOTCLE/'DEP2DEP','GRAD ','CHANCHAM','SSDECOUP',
  24. * $ 'CHAMPAUX','GRADCHAM'/
  25.  
  26. CALL LIRMOT(MOTCLE,NMOT0,IMOT0,1)
  27.  
  28. c XFEM 'RECO'
  29. IF(IMOT0.EQ.1) THEN
  30. CALL LIROBJ('CHPOINT ',IPCHP1,1,IRETOU)
  31. CALL ACTOBJ('CHPOINT ',IPCHP1,1)
  32. if (iretou.ne.1)
  33. & write(6,*) 'il faut fournir un chPO de déplacements'
  34. CALL LIROBJ('MMODEL',IPMOD1,1,IRETOU)
  35. CALL ACTOBJ('MMODEL',IPMOD1,1)
  36. if (iretou.ne.1)
  37. & write(6,*) 'il faut fournir le modèle'
  38. IF (IERR.NE.0) RETURN
  39.  
  40. c lecture eventuelle d'un maillage sur lequel faire la recombinaison
  41. CALL LIROBJ('MAILLAGE',IPT2in,0,IRETOU)
  42. IF(IERR.NE.0) RETURN
  43. IF(IRETOU.NE.1) THEN
  44. c on travaille sur les noeuds du chpoint (approximation simple)
  45. CALL XPOST1(IPCHP1,IPMOD1,IPCHP2)
  46. ELSE
  47. c on recombine sur les noeuds du maillage fourni (interpolation)
  48. CALL ACTOBJ('MAILLAGE',IPT2in,1)
  49. IF (IERR.NE.0) RETURN
  50. CALL XPOST3(IPCHP1,IPMOD1,IPT2in,IPCHP2)
  51. RETURN
  52. ENDIF
  53. IF (IERR.NE.0) RETURN
  54.  
  55. C ECRITURE
  56. CALL ACTOBJ('CHPOINT ',IPCHP2,1)
  57. CALL ECROBJ('CHPOINT ',IPCHP2)
  58.  
  59. c XFEM 'FISS'
  60. ELSEIF(IMOT0.EQ.2) THEN
  61. CALL XPOST2
  62.  
  63. c ELSEIF(IMOT0.EQ.3) THEN
  64. c CALL XPOST3
  65. c
  66. c ELSEIF(IMOT0.EQ.4) THEN
  67. c CALL XPOST4
  68. c
  69. c ELSEIF(IMOT0.EQ.5) THEN
  70. c call LIRENT(IVAL,1,IRETOU)
  71. c if(IVAL.eq.1) call XPOS51
  72. c if(IVAL.eq.2) call XPOS52
  73. C
  74. c ELSEIF(IMOT0.EQ.6) THEN
  75. c CALL XPOST6
  76.  
  77. ELSE
  78. WRITE(*,*) '!!! OPTION NON IMPLANTÉE !!!'
  79. RETURN
  80.  
  81. ENDIF
  82.  
  83. END
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  

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