Esope source files

   1 : C CRECH1    SOURCE    PV        13/04/17    21:15:05     7767           
   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 :         do j=1,num(/2)
  14 :          if (num(1,j).ne.ipt1.num(1,j)) goto 10
  15 :         enddo
  16 : *  maillage identiques
  17 : *       write (6,*) ' crech1 ',meleme,' remplace par ',ipt1,
  18 : *    >   'position ',i
  19 :         if (idet.eq.1) segsup meleme
  20 :         meleme=ipt1
  21 :         segact meleme
  22 :         if = i
  23 :         goto 20
  24 :         return
  25 :   10  continue
  26 :   20  continue
  27 :       if (if.lt.nbemel/3) return
  28 : *  on rajoute le maillage courant en tête
  28 :       call ooogll(1)
  29 :       do j=if,2,-1
  30 :         premel(j)=premel(j-1)
  31 :         nbso(j)=nbso(j-1)
  32 :       enddo
  33 :       premel(1)=meleme
  34 :       nbso(1)=num(/2)
  34 :       call ooogll(0)
  35 :       return
  36 :       end
  37 :  
  38 :  
  39 :  
  40 :