1 : C ELFBL4 SOURCE CHAT 05/01/12 23:31:23 5004 2 : SUBROUTINE ELFBL4(N,DC,DNORM,P,VPOCHA) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8(A-H,O-Z) 5 : C 6 : C ===================================================================== 7 : C SOUS-PROGRAMME FORTRAN APPELE PAR ELFRES 8 : C CALCUL DES COMPOSANTES DNCN DANS LE REPERE GENERAL 9 : C N*2 : NOMBRE D'ELEMENTS (2 PTS PAR ELEMENT) 10 : C DNORM : VECTEUR DE NORMALISATION 11 : C P : MATRICES DE CHANGEMENT DE REPERE 12 : C DC : VECTEUR DNCN LOCAL 13 : C VPOCHA: VECTEUR DNCN GLOBAL 14 : C CREATION : 4/11/86 15 : C PROGRAMMEUR : GUILBAUD 16 : C ===================================================================== 17 : C 18 : DIMENSION DC(1),DNORM(1),P(1),DN(3),VPOCHA(N,1) 19 : IP=0 20 : NSTR=N/2 21 : NN=0 22 : DO 60 NS=1,NSTR 23 : MP=9*(NS-1) 24 : DO 50 NP=1,2 25 : NN=NN+1 26 : DO 40 KI=1,4 27 : I3=3*(KI-1) 28 : I3P=I3+IP 29 : DO 10 J=1,3 30 : DN(J)=DNORM(I3P+J)*DC(I3P+J) 31 : 10 CONTINUE 32 : DO 30 I=1,3 33 : AVI=0.D0 34 : MPI=MP+3*(I-1) 35 : DO 20 J=1,3 36 : AVI=AVI+P(MPI+J)*DN(J) 37 : 20 CONTINUE 38 : VPOCHA(NN,I3+I)=AVI 39 : 30 CONTINUE 40 : 40 CONTINUE 41 : IP=IP+12 42 : 50 CONTINUE 43 : 60 CONTINUE 44 : RETURN 45 : END 46 :
© Cast3M 2003 - All rights reserved.
Disclaimer