1 : C CRECH1B SOURCE PV 14/01/10 21:15:04 7903 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 : segact mcoord 10 : ** write (6,*) ' premel avant ',(premel(i),i=1,nbemel) 11 : do 30 i=1,nbemel 12 : if (premel(i).eq.0) goto 40 13 : meleme=premel(i) 14 : segact meleme 15 : 30 continue 16 : 40 continue 16 : call ooogll(1) 17 : nbpts=xcoor(/1)/(idim+1) 18 : do 10 i=1,nbemel 19 : if (premel(i).eq.0) goto 20 20 : meleme=premel(i) 21 : do j=1,num(/2) 22 : if (num(1,j).gt.nbpts) premel(i)=0 23 : enddo 24 : 10 continue 25 : 20 continue 26 : il=0 27 : do i=1,nbemel 28 : if (premel(i).ne.0) then 29 : il=il+1 30 : premel(il)=premel(i) 31 : endif 32 : enddo 33 : do i=il+1,nbemel 34 : premel(i)=0 35 : enddo 35 : call ooogll(0) 36 : segdes mcoord 37 : ** write (6,*) ' premel apres ',(premel(i),i=1,nbemel) 38 : return 39 : end 40 : 41 : 42 : 43 : 44 :
© Cast3M 2003 - All rights reserved.
Disclaimer