Télécharger la plaquette Cast3M


fr en

Télécharger jpmm1.eso

Retour à la liste

Numérotation des lignes : 

jpmm1

  1. C JPMM1     SOURCE    CB215821  24/04/12    21:16:28     11897          
  2.       SUBROUTINE JPMM1(ILREE,ICHAMJ,ILREE2,IMOD5,MINTE1,NBELEM)
  3.  
  4. * ---------------------------------------------------------------------
  5. * 
  6. *                               ROUTINE JPMAG1
  7. * 
  8. * ---------------------------------------------------------------------
  9. * Auteur : Nikola JERANCE
  10. * 
  11. * Historique : 
  12. * 
  13. *       20/09/2022 - Première écriture (N. Jerance)
  14. *     
  15. * 
  16. * Descriptif :
  17. * 
  18. *  Le potentiel magnétique vecteur est calculé à partir de J(x,y)
  19. *  et la matrice obtenue par MPMAGV. Le potentiel est en 3D.
  20. * 
  21. *  
  22. *  
  23. * Arguments :
  24. * 
  25. *       (E)     ILREE = Pointeur sur la liste de réels d'entrée
  26. *       (E)     ICHAMJ = J(x,y) en tant que MCHAML
  27. *       (S)     ILREE2 = Pointeur sur la liste de réels de sortie
  28. *       (E)     IMOD5 = pointeur sur le modèle
  29. *       (E)     MINTE1 = pointeur sur MINTE
  30. *       (E)     NBELEM = nb d'éléments
  31. * 
  32. * ---------------------------------------------------------------------
  33.  
  34. * ----------------------------------------------------------------------
  35. * 
  36. *                         0 - DECLARATIONS ET IMPORTS
  37. * 
  38. * ----------------------------------------------------------------------
  39.  
  40.       IMPLICIT INTEGER(I-N)
  41.       IMPLICIT REAL*8(A-H,O-Z)
  42.  
  43. -INC PPARAM
  44. -INC CCOPTIO
  45. -INC CCGEOME
  46. -INC CCREEL
  47. -INC SMCOORD
  48. -INC SMELEME
  49. -INC SMLREEL
  50. -INC SMCHAML
  51. -INC SMMODEL
  52. -INC SMINTE
  53.  
  54.  
  55.  
  56.       INTEGER NBEL1,NBEL2
  57.       POINTEUR IA.MCHELM
  58.       POINTEUR ILREE1.MLREEL
  59.       POINTEUR ILREE2.MLREEL
  60.       POINTEUR IJ.MCHELM
  61.       POINTEUR IJX.MELVAL
  62.       POINTEUR IJY.MELVAL
  63.       POINTEUR IMOD1.IMODEL
  64.       POINTEUR IJ1.MCHAML
  65.  
  66. C      WRITE(IOIMP,*) '    ***  CALCUL DE A  *** '
  67.  
  68.       IJ = ICHAMJ
  69.       SEGACT,IJ
  70.  
  71.       NB1 = IJ.ICHAML(/1)
  72.  
  73.       IJ1 = IJ.ICHAML(1)
  74.       SEGACT,IJ1
  75.  
  76.       IMOD1 = IMOD5
  77.  
  78.  
  79.  
  80.       ILREE1 = ILREE
  81.       SEGACT,ILREE1
  82.       JG = NBELEM*NBELEM
  83.       SEGINI,ILREE2
  84.  
  85.       NCOMPO = IJ1.IELVAL(/1)
  86.  
  87. *      DO ICOMPO = 1,NCOMPO
  88.         ICOMPO = 1
  89.          IJX = IJ1.IELVAL(ICOMPO)
  90.          SEGACT,IJX
  91.  
  92.           NBEL1 = IJX.VELCHE(/2)
  93.           NBEL2 = IJX.VELCHE(/1)
  94.          IF (NBEL1.EQ.1) THEN
  95.            V = 0.0
  96.              DO I1 = 1,NBEL2
  97.               VX = IJX.VELCHE(I1,1)
  98.               V = V + VX/NBEL2
  99.              END DO
  100.  
  101.           K=0
  102.           DO I = 1, NBELEM
  103.  
  104.  
  105.            DO J1=1,NBELEM
  106.              K = K + 1
  107.                ILREE2.PROG(K) = ILREE1.PROG(K)*V
  108.            END DO
  109.  
  110.           END DO
  111.  
  112.  
  113.          ELSE
  114.           K=0
  115.           DO I = 1, NBELEM
  116.            V = 0.0
  117.              DO I1 = 1,NBEL2
  118.               VX = IJX.VELCHE(I1,I)
  119.               V = V + VX/NBEL2
  120.              END DO
  121.            DO J=1,NBELEM
  122.              K = K + 1
  123.              ILREE2.PROG(K) =  ILREE1.PROG(K)*V
  124.            END DO
  125.  
  126.          END DO   
  127.  
  128.         END IF
  129.  
  130.  
  131.  
  132.       ICHAMA = IA
  133.  
  134.       SEGDES,ILREE1
  135.       SEGDES,ILREE2
  136.       SEGDES,IJ
  137.       SEGDES,IJX
  138.  
  139.  
  140.  
  141.       RETURN
  142.       END 
  143.  
  144.  

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