Télécharger chmtmd.eso

Retour à la liste

Numérotation des lignes :

  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.  
  71.  

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