genjlf
C GENJLF SOURCE CHAT 05/01/13 00:16:57 5004 C-------------------------------------------------------------------- C C C PP 9/97 C Pierre Pegon/JRC Ispra C-------------------------------------------------------------------- IMPLICIT INTEGER(I-N) IMPLICIT REAL*8 (A-H,O-Z) C INTEGER P,SZ,FLA ,SZONE DIMENSION P(NPTO,1),SZ(1),NELM(1),FLA(1) DIMENSION NUM(1) DIMENSION N14TO4(4,6),N16TO4(4,3),N25TO4(4) DIMENSION N16TO3(3,2),N23TO3(3,4),N25TO3(3,4) DATA N14TO4/1,2,6,5, 2,3,7,6, 3,4,8,7, 4,1,5,8, 1,2,3,4, 5,6,7,8/ DATA N16TO4/1,2,5,4, 2,3,6,5, 3,1,4,6/ DATA N25TO4/1,2,3,4/ DATA N16TO3/1,2,3, 4,5,6/ DATA N23TO3/1,2,4, 2,3,4, 3,1,4, 1,2,3/ DATA N25TO3/1,2,5, 2,3,5, 3,4,5, 4,1,5/ C RETURN C ENTRY GF1424(P,SZ,NELM,FLA, IFA,NPTO, NUM,SZONE,NUMEL) DO IE1=1,6 IFA=IFA+1 SZ (IFA)=SZONE NELM(IFA)=NUMEL FLA (IFA)=0 DO IE2=1,4 P(IE2,IFA)=NUM(N14TO4(IE2,IE1)) FLA(IFA)=FLA(IFA)+P(IE2,IFA) ENDDO ENDDO RETURN C ENTRY GF1624(P,SZ,NELM,FLA, IFA,NPTO, NUM,SZONE,NUMEL) DO IE1=1,3 IFA=IFA+1 SZ (IFA)=SZONE NELM(IFA)=NUMEL FLA (IFA)=0 DO IE2=1,4 P(IE2,IFA)=NUM(N16TO4(IE2,IE1)) FLA(IFA)=FLA(IFA)+P(IE2,IFA) ENDDO ENDDO RETURN C ENTRY GF2524(P,SZ,NELM,FLA, IFA,NPTO, NUM,SZONE,NUMEL) IFA=IFA+1 SZ (IFA)=SZONE NELM(IFA)=NUMEL FLA (IFA)=0 DO IE2=1,4 P(IE2,IFA)=NUM(N25TO4(IE2)) FLA(IFA)=FLA(IFA)+P(IE2,IFA) ENDDO RETURN C ENTRY GF1623(P,SZ,NELM,FLA, IFA,NPTO, NUM,SZONE,NUMEL) DO IE1=1,2 IFA=IFA+1 SZ (IFA)=SZONE NELM(IFA)=NUMEL FLA (IFA)=0 DO IE2=1,3 P(IE2,IFA)=NUM(N16TO3(IE2,IE1)) FLA(IFA)=FLA(IFA)+P(IE2,IFA) ENDDO ENDDO RETURN C ENTRY GF2323(P,SZ,NELM,FLA, IFA,NPTO, NUM,SZONE,NUMEL) DO IE1=1,4 IFA=IFA+1 SZ (IFA)=SZONE NELM(IFA)=NUMEL FLA (IFA)=0 DO IE2=1,3 P(IE2,IFA)=NUM(N23TO3(IE2,IE1)) FLA(IFA)=FLA(IFA)+P(IE2,IFA) ENDDO ENDDO RETURN C ENTRY GF2523(P,SZ,NELM,FLA, IFA,NPTO, NUM,SZONE,NUMEL) DO IE1=1,4 IFA=IFA+1 SZ (IFA)=SZONE NELM(IFA)=NUMEL FLA (IFA)=0 DO IE2=1,3 P(IE2,IFA)=NUM(N25TO3(IE2,IE1)) FLA(IFA)=FLA(IFA)+P(IE2,IFA) ENDDO ENDDO RETURN C END
© Cast3M 2003 - Tous droits réservés.
Mentions légales