Télécharger la plaquette Cast3M


fr en

Télécharger mpmag1.eso

Retour à la liste

Numérotation des lignes : 

mpmag1

  1. C MPMAG1    SOURCE    CB215821  23/11/02    21:15:07     11779          
  2.       SUBROUTINE MPMAG1(IMAIL,XEPAIS,ILREE)
  3.  
  4. * ---------------------------------------------------------------------
  5. *
  6. *                               ROUTINE MPMAGV
  7. *
  8. * ---------------------------------------------------------------------
  9. * Auteur : Nikola JERANCE
  10. *
  11. * Historique :
  12. *
  13. *       20/09/2022 - Premiere ecriture (N. Jerance)
  14. *
  15. *
  16. * Descriptif :
  17. *
  18. * Calcul de la matrice  qui exprime le potentiel vecteur magnetique
  19. * à partir de la densite de courant J. En entree on a un maillage 2D
  20. * ayant une epaisseur donnee et en sortie on a une liste de reels
  21. * (qui est en fait une matrice M).
  22. * Pour calculer le potentiel vecteur (en 3D) par la suite, il faut multiplier M
  23. *   par la densite de courant J, en appelant la procedure JPMAGV.
  24. *  Ce calcul correspond, par exemple, à la geometrie d'une galette supraconductrice
  25. * La densite de courant n'a que 2 composantes.
  26. *
  27. *
  28. * Arguments :
  29. *
  30. *       (E)     IMAIL = Pointeur sur un MAILLAGE
  31. *       (E)     XEPAIS = Epaisseur
  32. **      (S)     ILREE1 = Pointeur sur la liste de reels (resultat)
  33. *
  34. * ---------------------------------------------------------------------
  35.  
  36. * ----------------------------------------------------------------------
  37. *
  38. *                         0 - DECLARATIONS ET IMPORTS
  39. *
  40. * ----------------------------------------------------------------------
  41.  
  42.       IMPLICIT INTEGER(I-N)
  43.       IMPLICIT REAL*8(A-H,O-Z)
  44.  
  45. -INC PPARAM
  46. -INC CCOPTIO
  47. -INC CCGEOME
  48. -INC CCREEL
  49. -INC SMCOORD
  50. -INC SMELEME
  51. -INC SMLREEL
  52.  
  53.  
  54.       INTEGER NBELEM,NBNN
  55.       POINTEUR IMAIL1.MELEME
  56.       POINTEUR ILREE1.MLREEL
  57.  
  58.       IMAIL1 = IMAIL
  59.  
  60. *     Verification que le maillage est un maillage simple, c-à-d avec un
  61. *     seul type d'elements (pour l'instant)
  62.       IF (IMAIL1.LISOUS(/1) .NE. 0) THEN
  63.          CALL ERREUR(426)
  64.          RETURN
  65.       ENDIF
  66.  
  67.       NBNN   = IMAIL1.NUM(/1)
  68.       NBELEM = IMAIL1.NUM(/2)
  69.  
  70.       CALL MATPMV(IMAIL1,ILREE1,NBELEM,NBNN,XEPAIS)
  71.  
  72.       ILREE = ILREE1
  73. * On calcule les valeurs du champ par element
  74.  
  75.       RETURN
  76.       END
  77.  
  78.  

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