Esope source files

   1 : C J3REDU    SOURCE    CHAT      05/01/13    00:47:20     5004           
   2 :       SUBROUTINE J3REDU(WWORK)
   3 : C----------------------------------------------------
   4 : C     REDUCTION DU NOMBRE DE TROU DANS WWORK
   5 : C
   6 : C     PP 6/97 refaite 12/98
   7 : C     Pierre Pegon/JRC Ispra
   8 : C----------------------------------------------------
   9 :       IMPLICIT INTEGER(I-N)
  10 : -INC  CCOPTIO

  11 :       SEGMENT WWORK
  12 :         REAL*8 PORIG(3),VNORM(3),VI(3),VJ(3)
  13 :         INTEGER FWORK
  14 :         INTEGER TWORK(NTROU)
  15 :       ENDSEGMENT
  16 : C
  17 :       NTROU1=TWORK(/1)
  18 : C
  19 :       I1=0
  20 :       NTROU=NTROU1
  21 :       DO IE1=1,NTROU1
  22 :         IF(TWORK(IE1).EQ.0)THEN
  23 :           NTROU=NTROU-1
  24 :         ELSE
  25 :           I1=I1+1
  26 :           IDUM=TWORK(IE1)
  27 :           TWORK(IE1)=TWORK(I1)
  28 :           TWORK(I1)=IDUM
  29 :         ENDIF
  30 :       ENDDO
  31 :       IF(NTROU.LT.NTROU1)THEN
  32 :         SEGADJ,WWORK
  33 :       ENDIF
  34 : C
  35 :       RETURN
  36 :       END
  37 :  
  38 :