Numérotation des lignes :

C PB302     SOURCE    MAGN      10/05/18    21:16:33     6675      SUBROUTINE PB302(XREF,X,Y,PG,FN,GR,FM,GM,ND,NP,MP,NPG,NOM2,ITYPI)      IMPLICIT INTEGER(I-N)      IMPLICIT REAL*8 (A-H,O-Z)C************************************************************************CC     CALCULE LES FONCTIONS DE FORME D'UN : TRI3CC     ^ etaC     |C   a |n3C     |\C     | \     a=sqrt(2)C     |  \C     |___\ _____>ksiC     0   aC     n1  n2CCC************************************************************************       REAL*8 XREF(ND,NP),X(NPG),Y(NPG)      DIMENSION FN(NP,NPG),GR(ND,NP,NPG),PG(NPG)      DIMENSION FM(MP,NPG),GM(ND,MP,NPG)      REAL*8 A,B,C,D,U(5),H(5)      CHARACTER*4 NOM2C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>       NP2=NP*ND       R2=SQRT(2.D0)       XREF(1,1)=0.D0      XREF(2,1)=0.D0      XREF(1,2)=R2      XREF(2,2)=0.D0      XREF(1,3)=0.D0      XREF(2,3)=R2       IF(ITYPI.EQ.2)THEN      X(1)=0.D0      Y(1)=0.D0      X(2)=R2      Y(2)=0.D0      X(3)=0.D0      Y(3)=R2      DO 2 L=1,3      PG(L)=1.D0/3.D0 2    CONTINUE      ELSE      CALL CALUHH(X,Y,PG,NPG)      ENDIF      DO 1 L=1,NPGC      FN(1,L)=-R2/2.D0*(X(L)+Y(L)-R2)      FN(2,L)=R2/2.D0*X(L)      FN(3,L)=R2/2.D0*Y(L)C      GR(1,1,L)=-R2/2.D0      GR(2,1,L)=-R2/2.D0      GR(1,2,L)=R2/2.D0      GR(2,2,L)=0.D0      GR(1,3,L)=0.D0      GR(2,3,L)=R2/2.D0C       IF(NOM2.EQ.'P1P1')THEN      FM(1,L)=FN(1,L)      FM(2,L)=FN(2,L)      FM(3,L)=FN(3,L)      GM(1,1,L)=-R2/2.D0      GM(2,1,L)=-R2/2.D0      GM(1,2,L)=R2/2.D0      GM(2,2,L)=0.D0      GM(1,3,L)=0.D0      GM(2,3,L)=R2/2.D0      ELSE      FM(1,L)=1.D0      GM(1,1,L)=0.D0      GM(2,1,L)=0.D0      ENDIF  1    CONTINUE C     WRITE(6,101)C     WRITE(6,1002)FMC     WRITE(6,1002)GMC     WRITE(6,1002)FNC     WRITE(6,1002)GRC     WRITE(6,101)      RETURN 1002 FORMAT(10(1X,1PD11.4)) 1001 FORMAT(20(1X,I5)) 101  FORMAT(1X,'... SUBPB302 ... FM,GM,FN,GR ',9(10H..........)/)C      END

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