Télécharger la plaquette Cast3M


fr en

Télécharger primi0.eso

Retour à la liste

Numérotation des lignes : 

primi0

  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