Numérotation des lignes :

C KREMPA    SOURCE    CB215821  16/04/21    21:17:36     8920      SUBROUTINE KREMPA (K2,KG,NG,NF,NA,NAL,C,U2,     -         SHC3D,SKCEL,SPROJA,SKBUFF)      IMPLICIT INTEGER(I-N)      IMPLICIT REAL*8 (A-H,O-Z)      DIMENSION U2(1),KG(2)C----------------------------------------------------------------------C Calcul des facteurs de forme en 3DC Sp appele par la famille KREMxGCC     CELLULES-ARETESC     ---------------CC     DETERMINATION DES CELLULES SITUEES SUR LA FACE 'NF' APPARTENANTC     A LA PROJECTION DE L'ARETE 'NA' DE L'ELEMENT 'K2'CC     NA  : INDICE DE L'ARETEC     U2  : NORMALE DE L'ELEMENT K2CC      NAL : INDICE DU NUMERO DE FACE DANS L'ARETE 1,2 OU 3C            SI 3  ON CALCULE NFCC      SI NFA = 3 : SEUL CAS OU ON CALCULE NFC---------------------------------------------------------------------- -INC TFFOR3D      SEGMENT SKCEL      INTEGER  KBCEL(NR,NR),IINT(2,NSTAC),IS(NSTAC),JS(NSTAC)      REAL*8   RMAX      ENDSEGMENTC        IF(NAL.NE.3) THEN                IF(NF.EQ.IFA(NAL,NA)) THEN                  LF = 1               ELSE                  LF = 2               ENDIF        ELSE           IF (NFA(NA).EQ.3) THEN               LF = 3               NF = IFA(3,NA)           ELSE               IF (IFA(3,NA).EQ.NF) THEN                   LF = 3               ELSE                   LF = 4               ENDIF           ENDIF        ENDIFC        KG(1) = KG(1) + IG(LF,1,NA)        KG(2) = KG(2) + IG(LF,2,NA)        NG = NG + 1C        NOC = NUMF(/2)        DO 403 KC = 1,NCEL(LF,NA)            I = ICEL(LF,1,KC,NA)            J = ICEL(LF,2,KC,NA)            KBCEL(I,J) = 1C          IF (PSC(NF,I,J).GT.-1.) THEN            B = 0.            DO 40   IES = 1,KES            B = B + U2(IES)*KSI(NF,IES)*V(KRO(NF,IES),I,J) 40         CONTINUECC           CELLULE SOMMET OU ARETEC            IF (ABS(B).GT.0.0001) THEN               Z = -  C / B               DFF = Z - ZB(NF,I,J)               DIFF = ABS(DFF)               NTY = NTYP(NF,I,J)               IF (DIFF.LT.1E-4.AND.NTY.LT.NOC) THEN                 DO 100 KT=1,NTY                 K = NUMF(NF,KT,I,J)                 IF (K.EQ.K2) GOTO 101 100             CONTINUE                  NTY = NTY + 1                  NUMF(NF,NTY,I,J) = K2                  NTYP(NF,I,J) = NTY 101             CONTINUE               ELSE                   IF (DFF.LT.-1E-3.AND.Z.GT.1E-4)  THEN                       ZB(NF,I,J) = Z                       NUMF(NF,1,I,J) = K2                       NTYP(NF,I,J) = 1                    ENDIF               ENDIF            ENDIF         ENDIF 403  CONTINUE       RETURN      END

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