Numérotation des lignes :

C FFPR15    SOURCE    CHAT      05/01/12    23:59:06     5004      SUBROUTINE FFPR15(XA)      IMPLICIT INTEGER(I-N)      IMPLICIT REAL*8 (A-H,O-Z)C************************************************************************CCCC************************************************************************      DIMENSION FN(15),COOR(3,3),XA(3,*)C      A=SQRT(2.D0)      AS2=A*0.5D0C?    CA2=2.D0/A/A      CA2=1.D0C?    CA4=4.D0/A/A      CA4=2.D0      A4=4.D0/A       COOR(1,1)=AS2      COOR(2,1)=0.D0      COOR(3,1)=0.5D0       COOR(1,2)=AS2      COOR(2,2)=AS2      COOR(3,2)=0.5D0       COOR(1,3)=0.D0      COOR(2,3)=AS2      COOR(3,3)=0.5D0       DO 1 I=1,3      X=COOR(1,I)      Y=COOR(2,I)      Z=COOR(3,I)       FN(1)= CA2*(X+Y-A)*(1.D0-Z)*(X+Y+A*Z-AS2)      FN(2)=-CA4*X*(X+Y-A)*(1.D0-Z)      FN(3)=-CA2*X*(1.D0-Z)*(-X+A*Z+AS2)      FN(4)= CA4*X*Y*(1.D0-Z)      FN(5)=-CA2*Y*(1.D0-Z)*(-Y+A*Z+AS2)      FN(6)=-CA4*Y*(X+Y-A)*(1.D0-Z)      FN(7)=-A4*Z*(1.D0-Z)*(X+Y-A)      FN(8)= A4*X*Z*(1.D0-Z)      FN(9) = A4*Y*Z*(1.D0-Z)      FN(10)= CA2*(X+Y-A)*Z*(X+Y+A*(1.D0-Z)-AS2)      FN(11)=-CA4*X*(X+Y-A)*Z      FN(12)=-CA2*X*Z*(-X+A*(1.D0-Z)+AS2)      FN(13)= CA4*X*Y*Z      FN(14)=-CA2*Y*Z*(-Y+A*(1.D0-Z)+AS2)      FN(15)=-CA4*Y*Z*(X+Y-A)       DO 2 N=1,3      XA(N,I+15)=     &         FN(1)*XA(N,1)+FN(2)*XA(N,2)+FN(3)*XA(N,3)+FN(4)*XA(N,4)     &+        FN(5)*XA(N,5)+FN(6)*XA(N,6)+FN(7)*XA(N,7)+FN(8)*XA(N,8)     &+  FN(9)*XA(N,9)+FN(10)*XA(N,10)+FN(11)*XA(N,11)+FN(12)*XA(N,12)     &+FN(13)*XA(N,13)+FN(14)*XA(N,14)+FN(15)*XA(N,15) 2    CONTINUE 1    CONTINUE    1002  format(10(1x,1pe11.4))      RETURN      END

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