1 subroutine transit(compo,dt0)
3 dimension q(24),temp(24)
34 c initialisation puissance a t=0
36 call cpeen(compo,CP_TEMPS,ti,1,"topo",20,topo,info)
37 CALL cpeRE(compo,CP_TEMPS, ti, 1, 'puissance', 24,
39 IF( info.NE. CPOK )GO TO 9000
42 c do while( ti.lT.100. )
46 c lecture de la temperature combustible entre ti et tf
48 CALL cplRE(compo,CP_TEMPS, ti,tf, npas, 'temperature', 24,
50 IF( info.NE. CPOK )GO TO 9000
57 c calcul de la puissance degagee en fonction de la position
58 c des barres et de la temperature
67 q(npt)=100.*(1.-0.0001*(temp(npt)-1000.))
71 write(6,*)"NEUTRO:","temps=",tf
73 c ecriture de la puissance a tf
75 call cpeen(compo,CP_TEMPS,ti,1,"topo",20,topo,info)
76 CALL cpeRE(compo,CP_TEMPS, tf, npas+1, 'puissance', 24,
78 IF( info.NE. CPOK )GO TO 9000
83 CALL cpfin(compo,CP_ARRET, info)
86 subroutine perma(compo)
88 dimension q(24),temp(24)
112 c initialisation puissance a iter=0
114 CALL cpeRE(compo,CP_ITERATION, ti, iter, 'puissi', 24,
116 IF( info.NE. CPOK )GO TO 9000
118 do while( iconv .eq. 0)
120 c lecture de la temperature combustible a iter
122 CALL cplRE(compo,CP_ITERATION, ti,tf, iter, 'tempi', 24,
124 IF( info.NE. CPOK )GO TO 9000
132 c calcul de la puissance degagee en fonction de la position
133 c des barres et de la temperature
142 q(npt)=100.*(1.-0.0001*(temp(npt)-1000.))
147 c ecriture de la puissance a iter+1
150 CALL cpeRE(compo,CP_ITERATION, ti, iter, 'puissi', 24,
152 IF( info.NE. CPOK )GO TO 9000
154 c lecture du flag de convergence iconv
156 CALL cplEN(compo,CP_ITERATION,ti, tf, iter , 'iconv', 1,
157 & nval, iconv , info)
158 write(6,*)"info:",info
159 write(6,*)"NEUTRO:",iter,iconv
161 IF( info.NE. CPOK )GO TO 9000
163 if(iconv.eq.1)go to 9000
168 CALL cpfin(compo,CP_ARRET, info)