C PRIMI0    SOURCE    OF166741  25/02/20    21:17:26     12165          
      subroutine primi0
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)

-INC PPARAM
-INC CCOPTIO
-INC SMLREEL
-INC SMEVOLL
      call lirobj('EVOLUTIO',mevol1,1,iret)
      if(ierr.ne.0) return
      segini,mevoll=mevol1
      n=ievoll(/1)
      do 1 i=1,n
      kevol1 = ievoll(i)
      segini,kevoll=kevol1
      ievoll(i)=kevoll
      if(typx.ne.'LISTREEL'.or.typy.ne.'LISTREEL')  go to 1
      mlreel=iprogx
      mlree1=iprogy
      segini,mlree2=mlreel
      segini,mlree3=mlree1
      iprogx=mlree2
      iprogy=mlree3
      b=mlree2.prog(1)
      S = 0.d0
      do 2 j=2,mlree2.prog(/1)
      a= mlree2.prog(j)
      if(B.gt.a) go to 1000
      sp=(mlree3.prog(j)+mlree3.prog(j-1))*(a-b)/2.d0
      mlree3.prog(j-1)=s
      s = s + sp
      b = a
    2 continue
      mlree3.prog(mlree3.prog(/1))=s
      segdes mlree2,mlree3, kevoll
    1 continue
      segdes mevoll
      call ecrobj('EVOLUTIO',mevoll)
      return
 1000 call erreur(872)
      segdes mevoll,kevoll,mlree2,mlree3
      return
      end