Télécharger crech1.eso

Retour à la liste

Numérotation des lignes :

  1. C CRECH1 SOURCE CB215821 16/09/30 21:15:05 9105
  2. * preconditionnement des meleme cree par crechp
  3. subroutine crech1(meleme,idet)
  4. -INC CCPRECO
  5. -INC SMELEME
  6. if=nbemel
  7. do 10 i=1,nbemel
  8. if (premel(i).eq.0) goto 20
  9. ipt1=premel(i)
  10. if (ipt1.eq.meleme) return
  11. if (nbso(i).ne.num(/2)) goto 10
  12. segact ipt1
  13. if (ipt1.num(/2).ne.num(/2)) goto 10
  14. do j=1,num(/2)
  15. if (num(1,j).ne.ipt1.num(1,j)) goto 10
  16. if (icolor(j).ne.ipt1.icolor(j)) goto 10
  17. enddo
  18. * maillage identiques
  19. * write (6,*) ' crech1 ',meleme,' remplace par ',ipt1,
  20. * > 'position ',i
  21. *** if (idet.eq.1) segsup meleme
  22. meleme=ipt1
  23. segact meleme
  24. if = i
  25. goto 20
  26. return
  27. 10 continue
  28. 20 continue
  29. if (if.lt.nbemel/3) return
  30. * on rajoute le maillage courant en tête
  31. call ooogll(1)
  32. do j=if,2,-1
  33. premel(j)=premel(j-1)
  34. nbso(j)=nbso(j-1)
  35. enddo
  36. premel(1)=meleme
  37. nbso(1)=num(/2)
  38. call ooogll(0)
  39. return
  40. end
  41.  
  42.  
  43.  
  44.  

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