Télécharger la plaquette Cast3M


fr en

Télécharger chmchl.eso

Retour à la liste

Numérotation des lignes : 

chmchl

  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