diff options
Diffstat (limited to 'Donjon/data/rep900EnrichCOMPO.x2m')
| -rw-r--r-- | Donjon/data/rep900EnrichCOMPO.x2m | 230 |
1 files changed, 230 insertions, 0 deletions
diff --git a/Donjon/data/rep900EnrichCOMPO.x2m b/Donjon/data/rep900EnrichCOMPO.x2m new file mode 100644 index 0000000..57343bc --- /dev/null +++ b/Donjon/data/rep900EnrichCOMPO.x2m @@ -0,0 +1,230 @@ +* +* Input file to compute diffusion flux on a homogenize assembly +* in an infinite space +* Homogenization can be homogeneous or heterogeneous +* +* Author : R. Chambon (05/2014) +* modified : R. Chambon (09/2015) +* +SEQ_ASCII XCOMPOU :: FILE '_CompoAss_1BU' ; +* +SEQ_ASCII AMF Amap AFLU ATRK :: FILE 'AMF' 'Amap' 'AFLU' 'ATRK' ; +SEQ_ASCII ACPOnew :: FILE 'ACPOnew' ; +LINKED_LIST Cpo DONNEES GeoInf GeoRes Matex Flux Track Micro + MacroF Macro2 System Fmap Macro ; +* +MODULE GEO: USPLIT: TRIVAT: FLUD: TRIVAA: NCR: DELETE: GREP: NAP: UTL: ABORT: END: ; +* +*REAL pbore_0 := 600. ; ! natural boron concentration in ppm + +INTEGER maxstep istep npin npinHt ifx ; +INTEGER MaxR := 16 ; +REAL Cote Lame mxel mxhl mxhr mxer Size ; +REAL burnup ppmBore crod ; +INTEGER iBoreMax := 1 ; !5 +INTEGER iTfMax := 1 ; +INTEGER iTcaMax := 1 ; +INTEGER iDcaMax := 1 ; +INTEGER iBarrMax := 2 ; +* +EVALUATE Cote := 1.262082 ; +EVALUATE Lame := .07792594 ; +* +STRING TypeCompo := "Heter0" ; +STRING Dir ; +* +IF TypeCompo "Heter0" = THEN + EVALUATE npinHt := 1 ; +* EVALUATE Dir := "EDI2A " ; + EVALUATE Dir := "Ass " ; +ELSEIF TypeCompo "Heter1" = THEN + EVALUATE npinHt := 1 ; + EVALUATE Dir := "EDI2C " ; +ENDIF ; +IF TypeCompo "Heter0" = THEN + EVALUATE ifx := 0 ; +ELSE + EVALUATE ifx := npinHt ; +ENDIF ; +* +REAL keff2 ; +* +EVALUATE npin := 17 ; +EVALUATE mxer := npin I_TO_R Cote * Lame 2. * + 2. / ; +EVALUATE mxhr := npin npinHt 2 * - I_TO_R Cote * 2. / ; +EVALUATE mxhl := mxhr -1. * ; +EVALUATE mxel := mxer -1. * ; +EVALUATE Size := npin I_TO_R Cote * Lame 2. * + ; +ECHO mxel mxhl mxhr mxer ; +* +Cpo := XCOMPOU ; + +IF ifx 0 = THEN + GeoInf := GEO: :: CAR2D 1 1 + EDIT 1 + X- REFL X+ REFL + Y- REFL Y+ REFL + MIX 1 + MESHX <<mxel>> <<mxer>> + MESHY <<mxel>> <<mxer>> + SPLITX 1 + SPLITY 1 + ; + Track := TRIVAT: GeoInf :: + MAXR <<MaxR>> MCFD 1 ; +ELSE + GeoInf := GEO: :: CAR2D 3 3 + EDIT 1 + X- REFL X+ REFL + Y- REFL Y+ REFL + MIX 3 2 3 + 2 1 2 + 3 2 3 + MESHX <<mxel>> <<mxhl>> <<mxhr>> <<mxer>> + MESHY <<mxel>> <<mxhl>> <<mxhr>> <<mxer>> + SPLITX 1 1 1 + SPLITY 1 1 1 + ; + Track := TRIVAT: GeoInf :: + MAXR <<MaxR>> DUAL 3 3 ; +ENDIF ; +ATRK := Track ; + + EVALUATE maxstep := 1 ; +* EVALUATE maxstep := 6 ; + DONNEES := UTL: :: + CREA + BURN <<maxstep>> = +0.0 +! 0.0 9.37498 18.7500 37.4999 74.9999 150.000 + ; + +DONNEES := UTL: DONNEES :: +* EDIT 10 + CREA + CBOR <<iBoreMax>> = +* 0.0 600.0 1200.0 1800.0 2400.0 +* 0.0 600.0 + 0.0 + CREA + TCOM <<iTfMax>> = + 5.58150024E+02 + CREA + FLUE <<iTcaMax>> = + 0.0 + CREA + DCOM <<iDcaMax>> = + 7.01012015E-01 + CREA + BARR <<iBarrMax>> = + 1.00000000E+00 2.00000000E+00 + ; +INTEGER iBoreStep := 0 ; + +ECHO "$$$ Crod loop" ; +INTEGER iCrodStep := 0 ; +WHILE iCrodStep iBoreMax < DO + EVALUATE iCrodStep := iCrodStep 1 + ; + ECHO "iCrodStep = " iCrodStep " of " iBarrMax ; + GREP: DONNEES :: GETVAL 'BARR' <<iCrodStep>> >>crod<< ; + +ECHO "$$$ Boron loop" ; +EVALUATE iBoreStep := 0 ; +WHILE iBoreStep iBoreMax < DO + EVALUATE iBoreStep := iBoreStep 1 + ; + ECHO "iBoreStep = " iBoreStep " of " iBoreMax ; + GREP: DONNEES :: GETVAL 'CBOR' <<iBoreStep>> >>ppmBore<< ; + +ECHO "$$$ Burnup loop" ; + +EVALUATE istep := 0 ; +WHILE istep maxstep < DO + EVALUATE istep := istep 1 + ; + GREP: DONNEES :: GETVAL 'BURN' <<istep>> >>burnup<< ; + +IF ifx 0 = THEN + Micro := NCR: Cpo :: + EDIT 0 + NMIX 1 + COMPO Cpo <<Dir>> + MIX 1 FROM 1 SET 'BURN' <<burnup>> + SET 'CBOR' <<ppmBore>> + SET 'TCOM' 5.58150024E+02 + SET 'DMOD' 7.01012015E-01 + SET 'FLUE' 0.0 + SET 'BARR' <<crod>> + ENDMIX + ; +ELSE + Micro := NCR: Cpo :: + EDIT 0 + NMIX 3 + COMPO Cpo <<Dir>> + MIX 1 FROM 1 SET 'BURN' <<burnup>> + SET 'CBOR' <<ppmBore>> + SET 'TCOM' 5.58150024E+02 + SET 'DMOD' 7.01012015E-01 + SET 'FLUE' 0.0 + SET 'BARR' <<crod>> + ENDMIX + MIX 2 FROM 2 SET 'BURN' <<burnup>> + SET 'CBOR' <<ppmBore>> + SET 'TCOM' 5.58150024E+02 + SET 'DMOD' 7.01012015E-01 + SET 'FLUE' 0.0 + SET 'BARR' <<crod>> + ENDMIX + MIX 3 FROM 3 SET 'BURN' <<burnup>> + SET 'CBOR' <<ppmBore>> + SET 'TCOM' 5.58150024E+02 + SET 'DMOD' 7.01012015E-01 + SET 'FLUE' 0.0 + SET 'BARR' <<crod>> + ENDMIX + ; +ENDIF ; +MacroF := Micro :: STEP UP MACROLIB ; + +*AMF := MacroF ; + +Macro2 := MacroF ; + +System := TRIVAA: Macro2 Track ; + +Flux := FLUD: System Track :: + EDIT 0 + ADI 5 +; +*AFLU := Flux ; +GREP: Flux :: + GETVAL 'K-EFFECTIVE ' 1 >>keff2<< ; +ECHO "ppmBore=" ppmBore "K-effective 2 = " keff2 ; + +Cpo := NAP: Cpo Track Flux :: + EDIT 10 + PROJECTION + STEP <<Dir>> + IFX <<ifx>> + SET 'BURN' <<burnup>> + SET 'CBOR' <<ppmBore>> + SET 'TCOM' 5.58150024E+02 + SET 'DMOD' 7.01012015E-01 + SET 'FLUE' 0.0 + SET 'BARR' <<crod>> +; + +Flux MacroF Macro2 Micro System + := DELETE: +Flux MacroF Macro2 Micro System ; + +ENDWHILE ; ! reference burnup loop + +ENDWHILE ; ! boron loop + +ENDWHILE ; ! crod loop + +ACPOnew := Cpo ; + +ECHO "test rep900EnrichCOMPO.x2m completed" ; +END: ; +QUIT . |