Télécharger chmchl.eso

Retour à la liste

Numérotation des lignes :

  1. C CHMCHL SOURCE CHAT 05/01/12 21:58:41 5004
  2. SUBROUTINE CHMCHL(IDSCHI,SP2,IZTR,PE)
  3. C
  4. C Issu de TRIOEF (TRCHEL)
  5. C
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8 (A-H,O-Z)
  8. SEGMENT IDSCHI
  9. REAL*8 GK(NYDIM),AA(NYDIM,NXDIM),FF(NZDIM,NPDIM)
  10. INTEGER IDX(NXDIM),IDY(NYDIM),IDZ(NZDIM),IDP(NPDIM),NN(6)
  11. INTEGER IDECY(NYDIM),IONZ(NXDIM)
  12. CHARACTER*32 NAME(NXDIM),NAMESP(NYDIM)
  13. ENDSEGMENT
  14. SEGMENT SP2
  15. REAL*8 GX(NXDIM),XX(NXDIM),GS(NZDIM),SS(NZDIM)
  16. REAL*8 TOT(NXDIM),TOTAQ(NXDIM),TOTFIX(NXDIM),GKS(NZDIM)
  17. REAL*8 YY(NXDIM),ZZ(NXDIM,NXDIM),CC(NYDIM),GC(NYDIM)
  18. ENDSEGMENT
  19. SEGMENT IZTR
  20. REAL*8 A0(NXDIM)
  21. ENDSEGMENT
  22. NXDIM=IDX(/1)
  23. NYDIM=IDY(/1)
  24. NZDIM=IDZ(/1)
  25. NPDIM=IDP(/1)
  26. CALL CHIADY(IDY,NYDIM,99,IDY1)
  27. NN6=NN(1)+NN(2)+NN(3)+NN(4)+NN(5)
  28. N5F=NN6-1
  29.  
  30. C SEGINI IZTR
  31. C0=CC(IDY1)
  32. GC0=GC(IDY1)
  33. IDYC0=IDECY(IDY1)
  34. DO KJ=1,NXDIM
  35. A0(KJ) = AA(IDY1,KJ)
  36. ENDDO
  37.  
  38. DO I=IDY1,N5F
  39. IDY(I)=IDY(I+1)
  40. CC(I)=CC(I+1)
  41. GC(I)=GC(I+1)
  42. GK(I)=GK(I+1)
  43. IDECY(I)=IDECY(I+1)
  44. DO KJ=1,NXDIM
  45. AA(I,KJ)=AA(I+1,KJ)
  46. ENDDO
  47. ENDDO
  48. IDY(NN6)= 99
  49. CC(NN6) = 10.D0**(-PE)
  50. GC(NN6) = -PE
  51. GK(NN6) =0.D0
  52. IDECY(NN6)=IDYC0
  53. DO KJ=1,NXDIM
  54. AA(NN6,KJ)=A0(KJ)
  55. ENDDO
  56. NN(3)=NN(3)-1
  57. NN(6)=NN(6)+1
  58. C SEGSUP IZTR
  59.  
  60. RETURN
  61. END
  62.  
  63.  
  64.  
  65.  
  66.  

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