NEW Sharing area
User documentation
Admin documentation
Development resources
Download
Private area
$(variable_name)
@{formula_expression}
&CRISTAL 'darwin4cig.cristal'and
&ISOTOPES 27 U234 U235 U236 U238 Pu238 Pu239 Pu240 Pu241 Pu242 Np237 Am241 Am243 Rh103 Cs133 Nd143 Nd145 Gd155 Mo95 Tc99 Ru101 Ag109 Sm147 Sm149 Sm150 Sm151 Sm152 Eu153to build a suitable output file with given isotopes
DEBUT_APOLLO2
**************************************************************
* *
* Procedure saphib pour DARWIN *
* *
**************************************************************
*
DEBPROC SAPHYB_DARWIN ;
***************************
*
ARGUMENT IMPRES/'ENTIER' TSTR*'TABLE'
TOPT*'TABLE' TRES*'TABLE' ;
*
WRITE: 'LISTING' 'EXECUTION DE LA SORTIE SAPHIB POUR DARWIN' ;
*********************************************
* initialisation de la saphyb (1er passage)
*********************************************
COMBUSTIBLE = MILIEU: &EQUI ;
IMPS: TSTR;
*nom_calc = TSTR.'indic' ;
WRITE 'LISTING' nom_calc ;
TSTR.'EQHOMO' = EQUIGEOM: &EDIT 3 TSTR.'GEO'
&EVOL &GLOB COMBUSTIBLE
;
SI ( NON ( EXISTE TSTR.'SAPH' ) ) ;
TOPT.'N_CYCL' = 1 ;
TSTR.'SAPH' = SAPHYB: &EDIT 3 TSTR.'EQHOMO'
TOPT.'SAPHYB'.'ADRESSE'
'Saphyb'
&PARA 'Cycle No' 'CYCL' &VALE &ENTI
&PARA 'Puissance' 'PUIS' &VALE &FLOT
&PARA 'temps' 'DTIM' &VALE &FLOT
&PARA 'temp_comb'
'TCOM' &TEMP TSTR.'MILREF'.'FISSIL1_ 1'
&PARA 'temp_mode' 'TMOD' &TEMP
TSTR.'MILREF'.'STRUCT2 '
&PARA 'irrad_comb' 'BURN' &IRRA
&ISOT &FISS &PF
&MACR 'macro' &TOUT
&REAC &ABSO &FISS &CHI &NUFI &ENER &SNNN
&LOCA 'irrad_locale' 'BLOC' &IRRA
&LOCA 'puiss_locale' 'PLOC' &PUIS
&LOCA 'dens_noy_lrds' 'MASL' &MASL ;
FINSI ;
*
*-------------------------------------------------------------------
* Condensation des informations
*-------------------------------------------------------------------
*
IMPS: TSTR ;
FLXN FNORM = NORMALISER_FLUX 0 'PUIS_WG' TOPT.'PWPG'
TSTR.'FLX'
TSTR.'CALIN'
TSTR.'MAC'
TSTR.'GEO' ;
WRITE 'LISTING' 'tout va bien' ;
TSTR.'MAQ' = PONDFLUX: &EDIT 2
FLXN TSTR.'MAC'
TSTR.'MAIL' TSTR.'EQHOMO' ;
TSTR.'APOLIQ' = CONDENSE: &EDIT 2
TSTR.'APOLIB' FLXN
TSTR.'MAIL' TSTR.'EQHOMO'
TSTR.'GEO'
TSTR.'CALIN' ;
*
*--------------------------------------------------------------------
* Recuperation du temps
*--------------------------------------------------------------------
*
DERNIER = EXTRAIRE: &ENTI TSTR.'MILREF'.'FISSIL1_ 1'
&TYPM 'PMED' 1 'NTIME' ;
SI ( DERNIER EQ 1 ) ALORS ;
TEMPO = 0. ;
SINON ;
TEMPO = EXTRAIRE:
&FLOT TSTR.'MILREF'.'FISSIL1_ 1'
&TYPM 'PMED' 1 'PTIME'
&FIND 'NOMT' &ITEM 'NPT' 'TEMPS'
&INDI 'NTIME' DERNIER ;
FINSI ;
*
*-------------------------------------------------------------------
* Inscription des resultats dans la saphyb
*-------------------------------------------------------------------
*
TSTR.'SAPH' = SAPHYB: &EDIT 3
TSTR.'SAPH' TSTR.'EQHOMO'
TSTR.'MAQ' TSTR.'APOLIQ'
&CYCL TOPT.'N_CYCL'
&PUIS TOPT.'PWPG'
&DTIM TEMPO ;
*
WRITE: 'LISTING' 'POINT AJOUTE A LA SAPHYB POUR :' ' '
' PAS NUMERO :' DERNIER
' TEMPS :' TEMPO
' BURNUP :' TOPT.'BU'
' PUISSANCE :' TOPT.'PWPG' ;
*
DETRUIRE: TSTR.'MAQ' ;
DETRUIRE: TSTR.'APOLIQ' ;
DETRUIRE: TSTR.'EQHOMO' ;
*
WRITE: 'LISTING' 'COURONNE 1' ;
IMPS: TSTR.'MILREF'.'FISSIL1_ 1' ;
WRITE: 'LISTING' 'COURONNE 2' ;
IMPS: TSTR.'MILREF'.'FISSIL1_ 2' ;
WRITE: 'LISTING' 'COURONNE 3' ;
IMPS: TSTR.'MILREF'.'FISSIL1_ 3' ;
WRITE: 'LISTING' 'COURONNE 4' ;
IMPS: TSTR.'MILREF'.'FISSIL1_ 4' ;
*
RESPROC TRES TSTR TOPT ;
*****
FINPROC ;
***************************
DEBPROC GENERE_MILIEUX_ISA ;
*
* -----oO0 CRISTAL 0Oo-----
*
* Proc�dure de generation des milieux et de preparation des
* donnees d'autoprotection dans le cadre de CRISTAL
*
* AUTEUR INITIAL : S. LOUBIERE (02/1998) FMC 001 / FSC 005
* MODIFIEE PAR : S. LOUBIERE (04/1998) FMC 002 / FI
* MODIFIEE PAR : S. LOUBIERE (04/1998) AJOUT d'ISOTOPES
* AUTOPROTEGEABLES
* MODIFIEE PAR : S. MENGELLE (06/2002) AJOUT d'ISOTOPES SM151
* AUTOPROTEGEABLES en V6
* Option autop des PF du CBU
* MODIFICATION : S. MENGELLE (03/2003) Introduction de TRES
* MODIFICATION : S. MENGELLE (09/2003) Introduction de TOPT.'VERSION'
* TOPT.'AUTOP_PF' = NON par defaut
* ajout d'isotopes autoprotegeables
* 'ER162' 'ER168' 'ER170'
* 'EU151' 'EU153' 'SM147' 'SM149'
* 'SM152' 'AG109'
* MODIFICATION : S. MENGELLE (10/2003)
* Introduction de TOPT.'STMIL'.nmil.'NOMREG'
*
* VERIFICATRICE : I. GUIMIER
* VERSION : TOPT.'VERSION'.'GENERE_MILIEUX_S'
* _________________________________________________________________
ARGUMENT IMPR/'ENTIER' TSTR*'TABLE' TOPT*'TABLE' TRES*'TABLE' ;
SI ( NON ( EXISTE IMPR ) ) ;
IMPR = 0 ;
FINSI ;
SI ( IMPR GE 1 ) ALORS ;
WRITE: 'LISTING'
'Execution de la procedure GENERE_MILIEUX_S version'
TOPT.'VERSION'.'GENERE_MILIEUX_S' ;
FINSI ;
MESSAGE = 'Erreur dans la procedure GENERE_MILIEUX_S' ;
* _________________________________________________________________
*
* GESTION DES INDICES D'IMPRESSION
* _________________________________________________________________
*
SI ( NON ( EXISTE TOPT.'STIMP' ) ) ;
TOPT.'STIMP' = TABLE: ;
FINSI ;
SI ( NON ( EXISTE TOPT.'STIMP'.'MILIEU' ) ) ;
TOPT.'STIMP'.'MILIEU' = 0 ;
FINSI ;
SI ( NON ( EXISTE TOPT.'AUTOP_PF' ) ) ;
TOPT.'AUTOP_PF' = 'NON' ;
FINSI ;
* _________________________________________________________________
*
* VARIABLES A FOURNIR EN ENTREE
* _________________________________________________________________
*
SI ( NON ( EXISTE TOPT.'STMIL') ) ;
WRITE: 'LISTING'
'IL MANQUE LA TABLE TOPT.STMIL' ' ' ;
ARRET: MESSAGE ;
FINSI ;
SI ( NON ( EXISTE TSTR.'MILREF') ) ;
WRITE: 'LISTING'
'IL MANQUE LA TABLE TSTR.MILREF' ' ' ;
ARRET: MESSAGE ;
FINSI ;
* _________________________________________________________________
*
* BOUCLE SUR LES MILIEUX ET LES ISOTOPES LES COMPOSANT
* _________________________________________________________________
*
nbmil = DIMTABLE: TOPT.'STMIL' ;
imil = 1 ;
tm = INITABLE: ( TOUTABLE: TOPT.'STMIL' 'INDICE' ) ;
*
tempo = TABLE: ;
ireg = 0 ;
SI ( TOPT.'BIBLIO' EQ 'CEA93_V6' ) ;
SI ( TOPT.'AUTOP_PF' EQ 'NON' ) ;
tauto = INITABLE:
'U238' 'PU240' 'U235' 'PU239' 'PU241'
'PU242' 'U236' 'PU238' 'U233' 'U234'
'NP237' 'AM241' 'AM243' 'ER162'
'ER166' 'ER167' 'ER168' 'ER170' 'EU151'
'EU153' 'GD154' 'GD156' 'GD155' 'GD157'
'GD158' 'GD160' 'SM147' 'SM149' 'SM152'
'AG107' 'AG109' 'IN115' 'ZRNAT' 'HF174'
'HF176' 'HF177' 'HF178' 'HF179' 'HF180'
'CM243' 'DY160' 'DY161' 'DY162' 'DY163'
'DY164'
;
SINON ;
tauto = INITABLE:
'U238' 'PU240' 'U235' 'PU239' 'PU241'
'PU242' 'U236' 'PU238' 'U233' 'U234'
'NP237' 'AM241' 'AM243' 'ER162'
'ER166' 'ER167' 'ER168' 'ER170' 'EU151'
'EU153' 'GD154' 'GD155' 'GD156' 'GD157'
'GD158' 'GD160' 'SM147' 'SM149' 'SM152'
'AG107' 'AG109' 'IN115' 'ZRNAT' 'HF174'
'HF176' 'HF177' 'HF178' 'HF179' 'HF180'
'CM243' 'TC99' 'CS133' 'MO95' 'RH103'
'ND145' 'DY160' 'DY161' 'DY162' 'DY163'
'DY164' 'SM151'
;
FINSI ;
SINON ;
SI ( TOPT.'AUTOP_PF' EQ 'NON' ) ;
tauto = INITABLE:
'U238' 'PU240' 'U235' 'PU239' 'PU241'
'PU242' 'U236' 'PU238' 'U233' 'U234'
'NP237' 'AM241' 'AM243' 'ER166'
'ER167' 'EU151' 'EU153' 'GD154' 'GD156'
'GD155' 'GD157' 'GD158' 'GD160' 'SM147'
'SM149' 'SM152' 'AG107' 'AG109' 'IN115'
'ZRNAT' 'HF174' 'HF176' 'HF177' 'HF178'
'HF179' 'HF180' 'CM243' 'DY160' 'DY161'
'DY162' 'DY163' 'DY164'
;
SINON ;
tauto = INITABLE:
'U238' 'PU240' 'U235' 'PU239' 'PU241'
'PU242' 'U236' 'PU238' 'U233' 'U234'
'NP237' 'AM241' 'AM243' 'ER166'
'ER167' 'EU151' 'EU153' 'GD154' 'GD155'
'GD156' 'GD157' 'GD158' 'GD160' 'SM147'
'SM149' 'SM152' 'AG107' 'AG109' 'IN115'
'ZRNAT' 'HF174' 'HF176' 'HF177' 'HF178'
'HF179' 'HF180' 'CM243' 'TC99' 'CS133'
'MO95' 'RH103' 'ND145' 'DY160' 'DY161'
'DY162' 'DY163' 'DY164'
;
FINSI ;
FINSI ;
tauto2 = TABLE: ;
nauto = DIMTABLE: tauto ;
i = 1 ;
REPETER B nauto ;
j = tauto.i ;
tauto2.j = TABLE: ;
i = i + 1 ;
FIN B ;
*
REPETER BOUCLE nbmil ;
nmil = tm.imil ;
TOPT.'STMIL'.nmil = CRISTAL_ISO_NAT TOPT.'STMIL'.nmil ;
nbiso = DIMTABLE: TOPT.'STMIL'.nmil ;
iiso = 1 ;
ti = INITABLE: ( TOUTABLE: TOPT.'STMIL'.nmil 'INDICE' ) ;
txt_mil = ' ' ;
txt_deb = ' ' ;
REPETER SOUS_BOUCLE nbiso ;
niso = ti.iiso ;
SI ( ( niso NE 'TEMPERATURE' ) ET ( niso NE 'NOMREG' ) ) ;
SI ( niso EQ 'DEBYE' ) ;
txt_deb = CONCAT: ' &DEBYE TOPT.''STMIL''.''' nmil
'''.''DEBYE''' ;
SINON ;
ciso = FORMAT: TOPT.'STMIL'.nmil.niso 'E10.4' ;
txt_mil = CONCAT: txt_mil ' TSTR.''TISO''.''' niso ''' '
ciso ;
FINSI ;
FINSI ;
*
* preparation de l'autoprotection
*
SI ( EXISTE tauto2.niso ) ;
SI ( EXISTE tempo.niso ) ;
ii = ( DIMTABLE: tempo.niso.'MILIEUX' + 1 ) ;
SINON ;
tempo.niso = TABLE: ;
ii = 1 ;
tempo.niso.'MILIEUX' = TABLE: ;
tempo.niso.'REGIONS' = TABLE: ;
tempo.niso.'ISO' = niso ;
FINSI ;
SI ( NON ( EXISTE creg ) ) ;
ireg = ireg + 1 ;
creg = FORMAT: ireg 'I2.2' ;
SI ( NON ( EXISTE TOPT.'STMIL'.nmil.'NOMREG' ) ) ;
nomreg = CONCAT: 'BB' creg ;
SINON ;
nomreg = TOPT.'STMIL'.nmil.'NOMREG' ;
FINSI ;
FINSI ;
tempo.niso.'MILIEUX'.ii = CONCAT: 'TSTR.''MILREF''.'''
nmil '''' ;
tempo.niso.'REGIONS'.ii = nomreg ;
FINSI ;
*
iiso = iiso + 1 ;
FIN SOUS_BOUCLE ;
txt_evol = ' ' ;
SI ( nmil EQ 'FISSIL1_ 1' ) ;
txt_evol = ' &CHAIN TSTR.''CHEVO'' ' ;
FINSI ;
SI ( nmil EQ 'FISSIL1_ 2' ) ;
txt_evol = ' &CHAIN TSTR.''CHEVO'' ' ;
FINSI ;
SI ( nmil EQ 'FISSIL1_ 3' ) ;
txt_evol = ' &CHAIN TSTR.''CHEVO'' ' ;
FINSI ;
SI ( nmil EQ 'FISSIL1_ 4' ) ;
txt_evol = ' &CHAIN TSTR.''CHEVO'' ' ;
FINSI ;
TSTR.'MILREF'.nmil = MILIEU: &EDIT TOPT.'STIMP'.'MILIEU'
( TEXTE txt_mil )
&TEMP TOPT.'STMIL'.nmil.'TEMPERATURE'
( TEXTE txt_deb )
( TEXTE txt_evol )
;
SI ( IMPR GE 2 ) ;
IMPS: TSTR.'MILREF'.nmil ;
FINSI ;
imil = imil + 1 ;
SI ( EXISTE creg ) ;
DETRUIRE: creg ;
FINSI ;
FIN BOUCLE ;
*
* Mise dans l'ordre des isotopes a autoproteger
*
nauto = DIMTABLE: tauto ;
iauto = 1 ;
iordre = 1 ;
TOPT.'STAUTO' = TABLE: ;
TOPT.'STAUTO'.'NBITER' = 1 ;
TOPT.'STAUTO'.'TORDRE' = TABLE: ;
REPETER B1 nauto ;
cauto = tauto.iauto ;
SI ( EXISTE tempo.cauto ) ;
TOPT.'STAUTO'.'TORDRE'.iordre = tempo.cauto ;
iordre = iordre + 1 ;
FINSI ;
iauto = iauto + 1 ;
FIN B1 ;
*
* Cas des calculs homogenes
*
SI ( nbmil EQ 1 ) ;
TOPT.'STAUTO'.'HOMOG' = 'OUI' ;
SINON ;
TOPT.'STAUTO'.'HOMOG' = 'NON' ;
FINSI ;
*
RESPROC TRES TSTR TOPT ;
FINPROC ;
*
*=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+
* CIGALES version 3.1 en date du 06/04/2006
* Creation du Fichier le 04/05/2007 15:43:34
*=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+
* -=- INITIALISATION - CALCUL 1 -=-
TOPT = TABLE: ;
TRES = TABLE: ;
TSTR = TABLE: ;
TOPT.'CALCUL_CRISTAL' = 1 ;
REPPROC = OUVRIR: 22 'VARIABLE' 1024 10000
'ADRESSE' 'aprocristal' ;
CHARGE_APROCRISTAL = LIRE: REPPROC 'APROC' 'CHARGE_APROCRISTAL' ;
FERMER: REPPROC ;
EXECUTER CHARGE_APROCRISTAL ;
TOPT.'BIBLIO' = 'CEA93_V9' ;
TRES TSTR TOPT = INITIALISER_CRISTAL 1 TSTR TOPT TRES ;
* -=- OPTIONS -=-
*
TOPT.'STCRI'.'NGROUP_FINAL' = 172 ;
TOPT.'STCRI'.'ANISOTROPIE' = 'P5' ;
TSTR.'CHEVO' = CHAINE_STD_CEA93_V9 3 TSTR.'TISO' ;
IMPS: TSTR.'CHEVO' ;
*
***********************************************************************
* d�but de la recopie de la partie obtenue par CIGALES
* (penser � remplacer GENERE_MILIEUX_S par GENERE_MILIEUX_ISA)
*
************************************************************************
*==============================================================
* APOLLO PIJ CALCUL 1
*
ANISO = CONCAT: '&' TOPT.'STCRI'.'ANISOTROPIE' ;
*==============================================================
*
*Crayon2 CAS 1
TITRE: 'Crayon2 CAS 1 ' ;
CALCUL_AP2 = 1 ;
WRITE: TOPT.'RESU' 'UOX_4% '
' CAS 1' ;
*
* -=- Description des milieux -=-
'CELLUL1 NZ=3 C1=.4127 C2=.4744 C3=.7448' ;
nom_calc = 'CELLUL1' ;
TOPT.'STCRI'.'CALCUL_INITIAL' = nom_calc ;
TOPT.'STCRI'.'CALCULS_INITIAUX'.nom_calc = TABLE: ;
TSTR.nom_calc = TABLE: ;
*
*UOX_4%
nom_mil = 'FISSIL1_ 1' ;
TOPT.'STMIL'.nom_mil = TABLE: ;
TOPT.'STMIL'.nom_mil.'U234 ' = 7.8342E-06 ;
TOPT.'STMIL'.nom_mil.'U235 ' = 9.0971E-04 ;
TOPT.'STMIL'.nom_mil.'U236 ' = 2.9666E-06 ;
TOPT.'STMIL'.nom_mil.'U238 ' = 2.1547E-02 ;
TOPT.'STMIL'.nom_mil.'O16 ' = 4.4934E-02 ;
TOPT.'STMIL'.nom_mil.'TEMPERATURE' = 600. ;
*
TOPT.'STMIL'.'FISSIL1_ 2' = TOPT.'STMIL'.'FISSIL1_ 1' ;
TOPT.'STMIL'.'FISSIL1_ 3' = TOPT.'STMIL'.'FISSIL1_ 1' ;
TOPT.'STMIL'.'FISSIL1_ 4' = TOPT.'STMIL'.'FISSIL1_ 1' ;
* Gaine
nom_mil = 'STRUCT1 ' ;
TOPT.'STMIL'.nom_mil = TABLE: ;
TOPT.'STMIL'.nom_mil.'FENAT ' = 1.383E-04 ;
TOPT.'STMIL'.nom_mil.'CRNAT ' = 7.073E-05 ;
TOPT.'STMIL'.nom_mil.'O16 ' = 2.874E-04 ;
TOPT.'STMIL'.nom_mil.'ZRNAT ' = 3.956E-02 ;
TOPT.'STMIL'.nom_mil.'TEMPERATURE' = 335. ;
*
* Moderateur
nom_mil = 'STRUCT2 ' ;
TOPT.'STMIL'.nom_mil = TABLE: ;
TOPT.'STMIL'.nom_mil.'H2O ' = 2.430E-02 ;
TOPT.'STMIL'.nom_mil.'B10 ' = 3.692E-06 ;
TOPT.'STMIL'.nom_mil.'B11 ' = 1.477E-05 ;
TOPT.'STMIL'.nom_mil.'TEMPERATURE' = 306.3 ;
*
TRES TSTR TOPT = GENERE_MILIEUX_ISA 2 TSTR TOPT TRES ;
*
* -=- Creation de la geometrie -=-
*
TSTR.nom_calc.'GEO' = GEOM: &CYLI &MAIL
4
.291795 .369094 .402211 .41266
1 &EQD .475364
3 &EQD .750416
&MILI
TSTR.'MILREF'.'FISSIL1_ 1' 1
TSTR.'MILREF'.'FISSIL1_ 2' 2
TSTR.'MILREF'.'FISSIL1_ 3' 3
TSTR.'MILREF'.'FISSIL1_ 4' 4
TSTR.'MILREF'.'STRUCT1 ' 5
TSTR.'MILREF'.'STRUCT2 ' 6 &A 8
;
*
***********************************************************************
* fin de la recopie de la partie obtenue par CIGALES
*
* (penser � remplacer GENERE_MILIEUX_S par GENERE_MILIEUX_ISA)
*
************************************************************************
*
*
TOPT.'SAPHYB' = TABLE: ;
*////////////////////////////////
* Adresse SAPHYB
*///////////////////////////////
TOPT.'SAPHYB'.'ADRESSE' =
'saphyb' ;
*
*
*######################################################"
*
*
TSTR.'GEO' = TSTR.nom_calc.'GEO' ;
* -=- Creation de la bibliotheque interne -=-
*
TSTR.'APOLIB' = BIBINT: &EDIT 1 TSTR.'IDB' TSTR.'GEO'
&SFIN &TP
( TEXTE TOPT.'REPBIB' ) ;
*
* -=- Autoprotection -=-
*
TSTR.'GEOAU' = TSTR.'GEO' ;
TRES TSTR TOPT = AUTOPROTECTION_S 3 TSTR TOPT TRES ;
*
*
* -=- Flux a B2 critique -=-
*
TOPT.'STPIJ' = TABLE: ;
TOPT.'STPIJ'.'UP' = 'LINEAIRE' ;
TOPT.'TYPE_B2' = 'CRITIQUE' ;
TRES TSTR TOPT = CALFLUX_PIJ_S 3 TSTR TOPT TRES ;
*
*
* -=- Sortie Saphyb -=-
*
*Puissance de normalisation (en W/g)
TOPT.'PWPG' = 36.6 ;
*
TOPT.'BU' = 0. ;
TRES TSTR TOPT = SAPHYB_DARWIN 3 TSTR TOPT TRES ;
*
* -=- Evolution -=-
*
TOPT.'TYPE_EVOL' = 'BURNUP' ;
TOPT.'STEVOL' = TABLE: ;
* TOPT.'STEVOL'.'NBITER' = 500 ;
* TOPT.'STEVOL'.'EPS1' = 0.05 ;
* TOPT.'STEVOL'.'EPS2' = 0.2 ;
* TOPT.'EV_EPSI' = 0.01 ;
TOPT.'STEVOL'.'BURNUP' = INITABLE: 9.375 1 18.75 1 37.5 1
75. 1 112.5 1 150. 1 325. 1 500. 1
750. 1 1000. 1 4000. 6 8000. 8
10000. 8 12000. 8 14000. 8 16000. 4
18000. 4 20000. 5 30000. 5 40000. 5
50000. 5 60000. 5 70000. 2 ;
TOPT.'STEVOL'.1 = INITABLE: 'OUI' AUTOPROTECTION_S 4000. 8000.
12000. 24000. 36000. 42000. ;
TOPT.'STEVOL'.2 = INITABLE: 'OUI' SAPHYB_DARWIN 37.5 75. 112.5 150.
325. 500. 750. 1000. 1500. 2000.
2500. 3000. 4000. 5000. 6000. 7000.
8000.
10000. 12000. 14000. 16000. 18000.
20000.
22000. 24000. 26000. 28000. 30000.
32000.
34000. 36000. 38000. 40000. 42000.
44000.
46000. 48000. 50000. 52000.
54000. 56000. 58000. 60000.
65000. 70000. ;
TSTR.'EQGEV' = EQUIGEOM: &EDIT 3 TSTR.'GEO' &EVOL &INDI FISSI ;
TRES TSTR TOPT = EVOLUTION_PIJ_S 2 TSTR TOPT TRES ;
EDTIME: ;
ARRET: ;
FIN_APOLLO2
DEBUT_DARWIN2
************************************************************
* PEPIN2 SAPHYR : Saisie interactive *
************************************************************
CODE: PEPIN ;
************************************************************
* Execute PSAPHY
************************************************************
PEP_obj = PSAPHY:
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
* Nom complet du fichier Saphyb
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
'saphyb'
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
* Saisie du mot cl� &PF ou &NL ou &PFNL ou &NLPF
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
&PFNL
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
* D�finition des paliers utilisateur par la donn�e d'un
* historique de puissance (en utilisant la fonction inte
* rpollation
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
&INTERPOL &PALIER 1
0.0 SECONDE 1912.57 JOUR 36.6
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
* Saisie �ventuelle du mot cl� de la temp�rature du mod�
* rateur
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
&TMOD
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
* Saisie des noms utilis�s � la constitution de la SAPHY
* B pour d�finir les param�tres
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
&PARA &CYCL 'CYCL'
&PARA &TIME 'DTIM'
&PARA &BURN 'BURN'
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
* S�lection des milieux � faire �voluer
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
&MILIEU 1 COMBUSTIBLE &TCOM
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
* Saisie du nom du fichier CHAINE DE FILIATIONS
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
'/produits/sec/AUTRES_CODES/DARWIN//biblio/CHAINE_PFNL.JEF2.BBL'
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
* Saisir du nom du fichier des CONSTANTES NUCL�AIRES
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
'/produits/sec/AUTRES_CODES/DARWIN//biblio/DARWIN_PFNL.JEF2.BBL'
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
* Saisie �ventuelle des d�finitions des mod�les de calcu
* l des taux de r�action pour les isotopes trait�s par A
* PPOLO et ceux non trait�s par APPOLO
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-==-=-=-=-
&MODELE
&APOLLO 3
&BIBLIO 2
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
* Saisie �ventuelle des noms de fichiers de type GENDF c
* ontenant les �valuations des sections efficaces micros
* copiques multigroupe
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
&FIGENDF 1
'/produits/sec/AUTRES_CODES/DARWIN//biblio/jef2.neutron.PFNL.bcd.gendf.923.172'
&EDITION
&PALIER &TOUS
&PUISS &NUPAL &TOUS
&NUMIL &TOUS
&FIN
*
&CRISTAL
'darwin4cig.cristal'
;
*
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
* Soumission du module PEPIN2
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
PEPIN2: PEP_obj
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
* Les options de calcul
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
&NOPT 3
&CALCUL DEBUT
&NBTREF 1
&PALIMP ELEMENTAIRE
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
* Titre identificateur du calcul
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
'benchmark_JAEA'
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
* base de r�sultats des concentrations
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
'CONC.resbrut'
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
* Nom du fichier de reprise ou de sauvegarde de reprise
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
'CONC.REPRISE'
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
* crit�re d'�volution
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
&EVOLUTION
&METHODE NUMERIQUE
&EPSILON 0.1
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
* temps de refroidissement en fin d'irradiation
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
&REFROIDISSEMENT
90.0 JOUR
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
* Saisie �ventuelle des application � r�aliser
*=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
&APPLICATION 2
CONCENTRATION
MASSE
&FIN
*
&MULTIREF
2
0.0 SECONDE
90. JOUR
&BUTAB 1 70.0000E+03
&ISOTOPES 27
U234 U235 U236 U238
Pu238 Pu239 Pu240 Pu241 Pu242 Np237 Am241 Am243
Rh103 Cs133 Nd143 Nd145 Gd155 Mo95 Tc99 Ru101 Ag109
Sm147 Sm149 Sm150 Sm151 Sm152 Eu153
;
ARRET: ;
FIN_DARWIN2
<CODE name="Darwin" cplugin="file:./plugins/calc/Darwin.cplugin.jar" command="**Promethee_Install_Path**/scripts/darwin2" />
line in calculator.xml file inside ”<calculator>” tag,