1 : C CHMTMD SOURCE CHAT 05/01/12 22:00:30 5004 2 : SUBROUTINE CHMTMD(IDSCHI,LGKMOD,ICOTY3,TMP,TMPNEW) 3 : C===================================================================== 4 : C 5 : C ISSU DE TRIO-EF (TRTMPD) 6 : C 7 : C PREMIERE APPROXIMATION D ULICH 8 : C 9 : C 10 : C====================================================================== 11 : IMPLICIT INTEGER(I-N) 12 : IMPLICIT REAL*8(A-H,O-Z) 13 : -INC SMLENTI 14 : POINTEUR ICOTY3.MLENTI 15 : SEGMENT IDSCHI 16 : REAL*8 GK(NYDIM),AA(NYDIM,NXDIM),FF(NZDIM,NPDIM) 17 : INTEGER IDX(NXDIM),IDY(NYDIM),IDZ(NZDIM),IDP(NPDIM),NN(6) 18 : INTEGER IDECY(NYDIM),IONZ(NXDIM) 19 : CHARACTER*32 NAME(NXDIM),NAMESP(NYDIM) 20 : ENDSEGMENT 21 : SEGMENT LGKMOD 22 : REAL*8 DELH0(NYDIM),DELCP0(NYDIM) 23 : ENDSEGMENT 24 : * 25 : C 26 : TMP=25.D0 27 : 28 : IF (ICOTY3.NE.0) NO3=ICOTY3.LECT(/1) 29 : IF (TMPNEW.EQ.TMP) RETURN 30 : NX=NN(1)+1 31 : NC=NN(1)+NN(2)+NN(3)+NN(4)+NN(5)+NN(6) 32 : T1=273.16D0+TMP 33 : T2=273.16D0+TMPNEW 34 : VI=219.30D0*(T2-T1)/T1/T2 35 : 36 : CBRUNO 37 : DO 200 I=NX,NN(1)+NN(2) 38 : GK(I)=GK(I)+VI*DELH0(I) 39 : 200 CONTINUE 40 : 41 : NI3=NN(1)+NN(2)+1 42 : NF3=NN(1)+NN(2)+NN(3) 43 : NI4=NN(1)+NN(2)+NN(3)+1 44 : DO 210 I=NI4,NC 45 : GK(I)=GK(I)+VI*DELH0(I) 46 : 210 CONTINUE 47 : 48 : DO 220 I=NI3,NF3 49 : IF (ICOTY3.NE.0) THEN 50 : DO J=1,NO3 51 : IF (IDY(I).EQ.ICOTY3.LECT(J)) GOTO 240 52 : ENDDO 53 : ENDIF 54 : GK(I)=GK(I)+VI*DELH0(I) 55 : 220 CONTINUE 56 : 57 : 240 CONTINUE 58 : * 59 : TMP=TMPNEW 60 : 61 : 62 : RETURN 63 : END 64 : C 65 : C 66 : 67 : 68 : 69 : 70 :
© Cast3M 2003 - All rights reserved.
Disclaimer