Télécharger deco.eso

Retour à la liste

Numérotation des lignes :

  1. C DECO SOURCE CB215821 16/12/05 21:15:16 9237
  2. SUBROUTINE DECO
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. *--------------------------------------------------------------------*
  6. * *
  7. * Op?rateur DEnsite de COurant *
  8. * ____________________________ *
  9. * *
  10. * Calcul de la densite de courant de Foucault *
  11. * *
  12. * Syntaxe n?1: DNS1 = DECO MOD1 DEP1 ( CAR1 ) ; *
  13. * *
  14. * Entr{es: *
  15. * *
  16. * MOD1 objet MMODEL *
  17. * DEP1 CHPOINT (FONCTION DE COURANT) *
  18. * CAR1 MCHAML de sous-type CARACTERISTIQUES *
  19. * *
  20. * Sortie: *
  21. * *
  22. * DNS1 MCHAML de sous-type GRADIENT *
  23. * *
  24. * *
  25. * Auteurs, date de cr{ation: *
  26. * *
  27. * *
  28. * Y. STEPHAN le 22/09/97 (copie de GRAD) *
  29. *--------------------------------------------------------------------*
  30. -INC CCOPTIO
  31. -INC CCHAMP
  32. *
  33. * NINFOS est le nombre d'informations contenues dans INFELE
  34. * ( VOIR LE S-P ELQUOI )
  35. *
  36. * PARAMETER ( NINFOS=15 )
  37. *
  38. * MC 19/01/98 : Fonction indisponible en d?fo planes g?n?ralis?es
  39. IF(IFOUR.EQ.-3) THEN
  40. CALL ERREUR(710)
  41. RETURN
  42. ENDIF
  43. *
  44. IPCHL1=0
  45. IPCHE1=0
  46. IRET1 =0
  47. IRET3 =0
  48. ICONV =0
  49. *
  50. *
  51. * LECTURE D'UN MMODEL
  52. *
  53. CALL LIROBJ('MMODEL',IPMODL,1,IRET)
  54. IF (IERR.NE.0) RETURN
  55. *
  56. IF (IRET.NE.0) THEN
  57. *
  58. * 1-}RE SYNTAXE
  59. * _____________
  60. *
  61. *
  62. * LECTURE D'UN CHPOINT
  63. *
  64. CALL LIROBJ('CHPOINT',IPCHP1,0,IRET1)
  65. IF(IERR.NE.0) RETURN
  66. IF (IRET1.NE.0) THEN
  67. *
  68. * CHPOINT ---> CHAMELEM AUX NOEUDS
  69. *
  70. CALL CHAME1(0,IPMODL,IPCHP1,' ',IPCHE2,1)
  71. IF (IERR.NE.0) RETURN
  72. ENDIF
  73. *
  74. * LECTURE D'UN MCHAML (CARACTERISTIQUE)
  75. *
  76. CALL LIROBJ('MCHAML',IPIN,0,IRET3)
  77. IF (IERR.NE.0) RETURN
  78. IPCHE1 = 0
  79. IF (IRET3 .EQ. 1) THEN
  80. CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  81. IF(IR .NE. 1) CALL ERREUR(KER)
  82. IF(IERR .NE. 0) RETURN
  83. ENDIF
  84. *
  85. ENDIF
  86. *
  87. * GESTION DES OP{RANDES
  88. * IPCHL1 = 0 (MCHAML calcule)
  89. *
  90. CALL DECO1(IPMODL,IPCHE2,IPCHE1,IPCHL1,IRET)
  91. *
  92. IF (IRET.EQ.1) THEN
  93. IF (ICONV.EQ.0) THEN
  94. IF (IRET1.NE.0) CALL DTCHAM(IPCHE2)
  95. CALL ECROBJ('MCHAML',IPCHL1)
  96. ELSE
  97. IF (IPCHE1.NE.0) CALL DTCHAM(IPCHE1)
  98. CALL DTCHAM(IPCHL1)
  99. CALL DTMODL(IPMODL)
  100. CALL ECROBJ('CHAMELEM',IPCHE2)
  101. ENDIF
  102. ENDIF
  103. RETURN
  104. END
  105.  
  106.  
  107.  

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