Télécharger pb201.eso

Retour à la liste

Numérotation des lignes :

pb201
  1. C PB201 SOURCE MAGN 10/05/18 21:16:30 6675
  2. SUBROUTINE PB201(XREF,X,PG,FN,GR,FM,GM,ND,NP,MP,NG,NPG,NOM2,ITYPI)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. C************************************************************************
  6. C
  7. C CALCULE LES FONCTIONS DE FORME D'UN : SEG2
  8. C
  9. C
  10. C************************************************************************
  11.  
  12. REAL*8 X(NPG),XREF(ND,NP)
  13. REAL*8 FN(NP,NPG),GR(ND,NP,NPG),PG(NPG)
  14. REAL*8 FM(MP,NPG),GM(ND,MP,NPG)
  15. REAL*8 A,B,C,D,U(5),H(5)
  16. CHARACTER*4 NOM2
  17. C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  18.  
  19. XREF(1,1)=0.D0
  20. XREF(1,2)=1.D0
  21.  
  22. NP2=NP*ND
  23.  
  24. CALL CALUHG(U,H,NG)
  25. A=0.D0
  26. B=1.D0
  27. IF(ITYPI.EQ.2)THEN
  28. X(1)=0.D0
  29. X(2)=1.D0
  30. PG(1)=0.5D0
  31. PG(2)=0.5D0
  32. ELSE
  33. CALL CALG1(A,B,NG,H,U,X,PG)
  34. ENDIF
  35. DO 1 L=1,NPG
  36. C
  37. FN(1,L)=1.D0-X(L)
  38. FN(2,L)=X(L)
  39. GR(1,1,L)=-1.D0
  40. GR(1,2,L)=1.D0
  41. 1 CONTINUE
  42.  
  43. IF(MP.EQ.1)THEN
  44. DO 2 L=1,NPG
  45. FM(1,L)=1.D0
  46. GM(1,1,L)=0.D0
  47. 2 CONTINUE
  48.  
  49. ENDIF
  50. IF(NOM2.EQ.'P1P1')THEN
  51. DO 3 L=1,NPG
  52. FM(1,L)=FN(1,L)
  53. FM(2,L)=FN(2,L)
  54. GM(1,1,L)=GR(1,1,L)
  55. GM(2,1,L)=GR(2,1,L)
  56. GM(1,2,L)=GR(1,2,L)
  57. GM(1,2,L)=GR(1,2,L)
  58. 3 CONTINUE
  59. ENDIF
  60.  
  61. C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  62. C WRITE(6,101)
  63. C WRITE(6,1002)FN
  64. C WRITE(6,1002)GR
  65. C WRITE(6,101)
  66. RETURN
  67. 1002 FORMAT(10(1X,1PD11.4))
  68. 1001 FORMAT(20(1X,I5))
  69. 101 FORMAT(1X,'... SUB PB201 ... FN,GR,FX,GX ',9(10H..........)/)
  70. C
  71. END
  72.  
  73.  
  74.  
  75.  
  76.  
  77.  
  78.  
  79.  

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