Télécharger primi0.eso

Retour à la liste

Numérotation des lignes :

  1. C PRIMI0 SOURCE CHAT 05/01/13 02:27:02 5004
  2. subroutine primi0
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. -INC CCOPTIO
  6. -INC SMLREEL
  7. -INC SMEVOLL
  8. call lirobj('EVOLUTIO',mevol1,1,iret)
  9. if(ierr.ne.0) return
  10. segini,mevoll=mevol1
  11. n=ievoll(/1)
  12. do 1 i=1,n
  13. kevol1 = ievoll(i)
  14. segini,kevoll=kevol1
  15. ievoll(i)=kevoll
  16. if(typx.ne.'LISTREEL'.or.typy.ne.'LISTREEL') go to 1
  17. mlreel=iprogx
  18. mlree1=iprogy
  19. segini,mlree2=mlreel
  20. segini,mlree3=mlree1
  21. iprogx=mlree2
  22. iprogy=mlree3
  23. b=mlree2.prog(1)
  24. S = 0.d0
  25. do 2 j=2,mlree2.prog(/1)
  26. a= mlree2.prog(j)
  27. if(B.gt.a) go to 1000
  28. sp=(mlree3.prog(j)+mlree3.prog(j-1))*(a-b)/2.d0
  29. mlree3.prog(j-1)=s
  30. s = s + sp
  31. b = a
  32. 2 continue
  33. mlree3.prog(mlree3.prog(/1))=s
  34. segdes mlree2,mlree3, kevoll
  35. 1 continue
  36. segdes mevoll
  37. call ecrobj('EVOLUTIO',mevoll)
  38. return
  39. 1000 call erreur(872)
  40. segdes mevoll,kevoll,mlree2,mlree3
  41. return
  42. end
  43.  
  44.  
  45.  

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