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.  
  6. -INC PPARAM
  7. -INC CCOPTIO
  8. -INC SMLREEL
  9. -INC SMEVOLL
  10. call lirobj('EVOLUTIO',mevol1,1,iret)
  11. if(ierr.ne.0) return
  12. segini,mevoll=mevol1
  13. n=ievoll(/1)
  14. do 1 i=1,n
  15. kevol1 = ievoll(i)
  16. segini,kevoll=kevol1
  17. ievoll(i)=kevoll
  18. if(typx.ne.'LISTREEL'.or.typy.ne.'LISTREEL') go to 1
  19. mlreel=iprogx
  20. mlree1=iprogy
  21. segini,mlree2=mlreel
  22. segini,mlree3=mlree1
  23. iprogx=mlree2
  24. iprogy=mlree3
  25. b=mlree2.prog(1)
  26. S = 0.d0
  27. do 2 j=2,mlree2.prog(/1)
  28. a= mlree2.prog(j)
  29. if(B.gt.a) go to 1000
  30. sp=(mlree3.prog(j)+mlree3.prog(j-1))*(a-b)/2.d0
  31. mlree3.prog(j-1)=s
  32. s = s + sp
  33. b = a
  34. 2 continue
  35. mlree3.prog(mlree3.prog(/1))=s
  36. segdes mlree2,mlree3, kevoll
  37. 1 continue
  38. segdes mevoll
  39. call ecrobj('EVOLUTIO',mevoll)
  40. return
  41. 1000 call erreur(872)
  42. segdes mevoll,kevoll,mlree2,mlree3
  43. return
  44. end
  45.  
  46.  
  47.  

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