Télécharger crech1b.eso

Retour à la liste

Numérotation des lignes :

  1. C CRECH1B SOURCE PV 19/07/12 21:15:05 10257
  2. * nettoyage du preconditionnement si le maillage est plus grand que les coordonnees
  3. subroutine crech1b
  4. -INC CCOPTIO
  5. -INC CCPRECO
  6. -INC SMELEME
  7. -INC SMCOORD
  8.  
  9. ith=oothrd
  10. segact mcoord
  11. ** write (6,*) ' premel avant ',(premel(i,ith),i=1,nbemel)
  12. do 30 i=1,nbemel
  13. if (premel(i,ith).eq.0) goto 40
  14. meleme=premel(i,ith)
  15. segact meleme
  16. 30 continue
  17. 40 continue
  18. nbpts=xcoor(/1)/(idim+1)
  19. do 10 i=1,nbemel
  20. if (premel(i,ith).eq.0) goto 20
  21. meleme=premel(i,ith)
  22. do j=1,num(/2)
  23. if (num(1,j).gt.nbpts) premel(i,ith)=0
  24. enddo
  25. 10 continue
  26. 20 continue
  27. il=0
  28. do i=1,nbemel
  29. if (premel(i,ith).ne.0) then
  30. il=il+1
  31. premel(il,ith)=premel(i,ith)
  32. nbso(il,ith)=nbso(i,ith)
  33. hash1(il,ith)=hash1(i,ith)
  34. hash2(il,ith)=hash2(i,ith)
  35. endif
  36. enddo
  37. do i=il+1,nbemel
  38. premel(i,ith)=0
  39. enddo
  40. segdes mcoord
  41. ** write (6,*) ' premel apres ',(premel(i,ith),i=1,nbemel)
  42. return
  43. end
  44.  
  45.  
  46.  
  47.  
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  

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