Télécharger mamu53.eso

Retour à la liste

Numérotation des lignes :

mamu53
  1. C MAMU53 SOURCE PV 20/02/22 21:15:19 10531
  2. subroutine mamu53(long,val11,lpl1,val1,lpl,pt)
  3. implicit real*8 (a-h,o-z)
  4. real*8 val11(*),val1(*)
  5. real*8 pt(36)
  6. id12=lpl1
  7. id13=lpl1+ lpl1+1
  8. id2=lpl
  9. id3=lpl+ lpl+1
  10. id4=lpl+ lpl+1 +lpl+2
  11. id5=lpl+ lpl+1 +lpl+2 +lpl+3
  12. pt1=pt(1)
  13. pt2=pt(2)
  14. pt3=pt(3)
  15. pt7=pt(4)
  16. pt8=pt(5)
  17. pt9=pt(6)
  18. pt13=pt(7)
  19. pt14=pt(8)
  20. pt15=pt(9)
  21. pt19=pt(10)
  22. pt20=pt(11)
  23. pt21=pt(12)
  24. pt25=pt(13)
  25. pt26=pt(14)
  26. pt27=pt(15)
  27. do 110 i=1,long
  28. xval1=val1(i)
  29. xval2=val1(i+id2)
  30. xval3=val1(i+id3)
  31. xval4=val1(i+id4)
  32. xval5=val1(i+id5)
  33. xval11=val11(i)
  34. xval12=val11(i+id12)
  35. xval13=val11(i+id13)
  36. pt1=pt1+xval1*xval11
  37. pt2=pt2+xval1*xval12
  38. pt3=pt3+xval1*xval13
  39. pt7=pt7+xval2*xval11
  40. pt8=pt8+xval2*xval12
  41. pt9=pt9+xval2*xval13
  42. pt13=pt13+xval3*xval11
  43. pt14=pt14+xval3*xval12
  44. pt15=pt15+xval3*xval13
  45. pt19=pt19+xval4*xval11
  46. pt20=pt20+xval4*xval12
  47. pt21=pt21+xval4*xval13
  48. pt25=pt25+xval5*xval11
  49. pt26=pt26+xval5*xval12
  50. pt27=pt27+xval5*xval13
  51. 110 continue
  52. pt(1)=pt1
  53. pt(2)=pt2
  54. pt(3)=pt3
  55. pt(4)=pt7
  56. pt(5)=pt8
  57. pt(6)=pt9
  58. pt(7)=pt13
  59. pt(8)=pt14
  60. pt(9)=pt15
  61. pt(10)=pt19
  62. pt(11)=pt20
  63. pt(12)=pt21
  64. pt(13)=pt25
  65. pt(14)=pt26
  66. pt(15)=pt27
  67.  
  68. end
  69.  
  70.  
  71.  
  72.  
  73.  
  74.  
  75.  
  76.  

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