Télécharger clebl3.eso

Retour à la liste

Numérotation des lignes :

clebl3
  1. C CLEBL3 SOURCE PV 17/12/08 21:15:48 9660
  2. SUBROUTINE CLEBL3(iwrk52,iele,igau)
  3. C
  4. C
  5. C
  6. C Modele de Leblond et Devaux
  7. C descente au niveau des points de Gauss
  8. C 06/97
  9. C
  10. C SOURCE BRIOTTET
  11. C
  12. C
  13. C
  14. IMPLICIT INTEGER(I-N)
  15. IMPLICIT REAL*8(A-H,O-Z)
  16. C
  17.  
  18. -INC PPARAM
  19. -INC CCOPTIO
  20. -INC SMCHAML
  21. -INC CCHAMP
  22. -INC DECHE
  23. C
  24. REAL*8 VWRK1(7),VWRK2(7)
  25.  
  26. wrk52 = iwrk52
  27. * DTPS fourni dans caracteristiques
  28. DTPS = valma0(14)
  29. 10 continue
  30. C
  31. T0=ture0(1)
  32. C
  33. if ((tempf-temp0).gt.0) then
  34. TPOINT0=(turef(1) - ture0(1))/(tempf-temp0)
  35. else
  36. TPOINT0 = 0.
  37. endif
  38. C
  39. ZA0=valma0(9)
  40. C
  41. ZF0=valma0(10)
  42. C
  43. ZB0=valma0(11)
  44. C
  45. ZM0=valma0(12)
  46. C
  47. VMS0=valma0(13)
  48. C
  49. VWRK1(1)=T0
  50. VWRK1(2)=TPOINT0
  51. VWRK1(3)=ZA0
  52. VWRK1(4)=ZF0
  53. VWRK1(5)=ZB0
  54. VWRK1(6)=ZM0
  55. VWRK1(7)=VMS0
  56. C
  57. C CALCUL DU NOUVEAU CHAMP DE PHASES
  58. C
  59. CALL CLEBL4(iwrk52,VWRK1,DTPS,VWRK2)
  60. C
  61. C sortie
  62. xmatf(9)=VWRK2(3)
  63. C
  64. xmatf(10)=VWRK2(4)
  65. C
  66. xmatf(11)=VWRK2(5)
  67. C
  68. xmatf(12)=VWRK2(6)
  69. C
  70. xmatf(13)=VWRK2(7)
  71. C
  72. do ii = 1,8
  73. xmatf(ii) = xmat(ii)
  74. enddo
  75. C
  76. RETURN
  77. END
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  

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