Télécharger la plaquette Cast3M


fr en

Télécharger bb302.eso

Retour à la liste

Numérotation des lignes : 

bb302

  1. C BB302     SOURCE    MAGN      10/05/31    21:15:00     6679
  2.       SUBROUTINE BB302(XREF,XCOPG,XPOPG,
  3.      $     FFPGV,DFFPGV,
  4.      $     FFPGP,DFFPGP,
  5.      $     X,Y,PG,FN,GR,FM,GM,ND,NP,MP,NPG)
  6.       IMPLICIT INTEGER(I-N)
  7.       IMPLICIT REAL*8 (A-H,O-Z)
  8. C************************************************************************
  9. C
  10. C     CALCULE LES FONCTIONS DE FORME D'UN : TRI4   P1 + bulle P1
  11. C                                              ou  P1 + bulle P0
  12. C     ^ eta
  13. C     |
  14. C   a |n3
  15. C     |\
  16. C     | \     a=sqrt(2)
  17. C     |n4\
  18. C     |___\ _____>ksi
  19. C     0   a
  20. C     n1  n2
  21. C
  22. C
  23. C************************************************************************
  24. * Entrées
  25.       INTEGER ND,NP,MP,NPG
  26. *
  27.       REAL*8 XREF(ND,NP)
  28.       REAL*8 XCOPG(ND,NPG)
  29.       REAL*8 XPOPG(NPG)
  30.       REAL*8 FFPGV(NP,NPG)
  31.       REAL*8 DFFPGV(NP,ND,NPG)
  32.       REAL*8 FFPGP(MP,NPG)
  33.       REAL*8 DFFPGP(MP,ND,NPG)
  34. * Sorties
  35.       REAL*8 X(NPG),Y(NPG)
  36.       REAL*8 PG(NPG)
  37.       REAL*8 FN(NP,NPG),GR(ND,NP,NPG)
  38.       REAL*8 FM(MP,NPG),GM(ND,MP,NPG)
  39. * Travail
  40.       INTEGER IPG,IP,ID
  41.       REAL*8 UNSVOL,USVUSD
  42. C***
  43.       A=SQRT(2.D0)
  44.       XREF(1,1)=0.D0
  45.       XREF(2,1)=0.D0
  46.       XREF(1,2)=A
  47.       XREF(2,2)=0.D0
  48.       XREF(1,3)=0.D0
  49.       XREF(2,3)=A
  50.       XREF(1,4)=A/3.D0
  51.       XREF(2,4)=A/3.D0
  52.  
  53. * magn veut que les volumes de ses éléments fasse 1
  54.       UNSVOL=2.D0
  55.       USVUSD=UNSVOL**(1.D0/DBLE(ND))
  56. * Recopie des points de Gauss
  57.       DO 1 IPG=1,NPG
  58.          X(IPG)=XCOPG(1,IPG)*USVUSD
  59.          Y(IPG)=XCOPG(2,IPG)*USVUSD
  60.          PG(IPG)=XPOPG(IPG)*UNSVOL
  61.  1    CONTINUE
  62. * Recopie des fns. de forme vitesse
  63.       DO 3 IPG=1,NPG
  64.          DO 32 IP=1,NP
  65.             FN(IP,IPG)=FFPGV(IP,IPG)
  66.             DO 322 ID=1,ND
  67.                GR(ID,IP,IPG)=DFFPGV(IP,ID,IPG)
  68.  322        CONTINUE
  69.  32      CONTINUE
  70.  3    CONTINUE
  71. * Recopie des fns. de forme pression
  72.       DO 5 IPG=1,NPG
  73.          DO 52 IP=1,MP
  74.             FM(IP,IPG)=FFPGP(IP,IPG)
  75.             DO 522 ID=1,ND
  76.                GM(ID,IP,IPG)=DFFPGP(IP,ID,IPG)
  77.  522        CONTINUE
  78.  52      CONTINUE
  79.  5    CONTINUE
  80. C     WRITE(6,101)
  81. C     WRITE(6,1002)FM
  82. C     WRITE(6,1002)GM
  83. C     WRITE(6,1002)FN
  84. C     WRITE(6,1002)GR
  85. C     WRITE(6,101)
  86.       RETURN
  87.  1002 FORMAT(10(1X,1PD11.4))
  88.  1001 FORMAT(20(1X,I5))
  89.  101  FORMAT(1X,'... SUBPB302 ... FM,GM,FN,GR ',9(10H..........)/)
  90. C
  91.       END
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  

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