diff options
Diffstat (limited to 'Dragon/data/proc_LZC')
| -rw-r--r-- | Dragon/data/proc_LZC/DimCell.c2m | 29 | ||||
| -rw-r--r-- | Dragon/data/proc_LZC/DimFuel.c2m | 49 | ||||
| -rw-r--r-- | Dragon/data/proc_LZC/DimLZC.c2m | 52 | ||||
| -rwxr-xr-x | Dragon/data/proc_LZC/GeoCell.c2m | 306 | ||||
| -rwxr-xr-x | Dragon/data/proc_LZC/GeoLZC.c2m | 1020 | ||||
| -rwxr-xr-x | Dragon/data/proc_LZC/GeoLZCM.c2m | 1020 | ||||
| -rw-r--r-- | Dragon/data/proc_LZC/MicLZC.c2m | 98 | ||||
| -rw-r--r-- | Dragon/data/proc_LZC/MicWUC.c2m | 189 | ||||
| -rw-r--r-- | Dragon/data/proc_LZC/MicWULZC.c2m | 395 |
9 files changed, 3158 insertions, 0 deletions
diff --git a/Dragon/data/proc_LZC/DimCell.c2m b/Dragon/data/proc_LZC/DimCell.c2m new file mode 100644 index 0000000..e738481 --- /dev/null +++ b/Dragon/data/proc_LZC/DimCell.c2m @@ -0,0 +1,29 @@ +*DECK DimCell.c2m +*---- +* Name : DimCell.c2m +* Type : DRAGON procedure +* Use : Dimensions for Gentilly-2 cell +* Author : G. Marleau +* Date : 2006/03/01 +* Reference : Gentilly-2 . +* +* Procedure called as: +* +* DimCell :: >>rcool<< >>rpt<< >>rgap<< >>rct<< +* >>LPitch<< >>BLength<< ; +* +* Output data: +* rcool : outer coolant radius +* rpt : outer radius of pressure tube +* rgap : outer radius of gap +* rct : outer radius of calandria tube +* LPitch : Lattice pitch +* BLength : Bundle length +*---- +* Definition of the input and output data structures and +* of the input and output parameters for this procedure +*---- +REAL rcool rpt rgap rct LPitch BLength := + 5.168878 5.60320 6.44780 6.587482 28.575 49.53 ; +:: <<rcool>> <<rpt>> <<rgap>> <<rct>> <<LPitch>> <<BLength>> ; +QUIT . diff --git a/Dragon/data/proc_LZC/DimFuel.c2m b/Dragon/data/proc_LZC/DimFuel.c2m new file mode 100644 index 0000000..d79c7a3 --- /dev/null +++ b/Dragon/data/proc_LZC/DimFuel.c2m @@ -0,0 +1,49 @@ +*DECK DimFuel.c2m +*---- +* Name : DimFuel.c2m +* Type : DRAGON procedure +* Use : Fuel dimensions for Gentilly-2 cell +* Author : G. Marleau +* Date : 2006/03/01 +* Reference : Gentilly-2 . +* +* Description of the procedure: +* DimFuel :: <<Ring>> +* >>npinp<< >>rpinp<< >>apinp<< >>rradf<< >>rrads<< ; +* Input data: +* Ring : Ring number +* Output data: +* npinp : number of pins in this ring +* rpinp : ring radius +* apinp : first pin angle +* rradf : outer fuel radius +* rrads : outer sheat radius +*---- +* Definition of the input and output data structures and +* of the input and output parameters for this procedure +* +MODULE END: ; +INTEGER Ring ; + :: >>Ring<< ; +INTEGER npinp1 npinp2 npinp3 npinp4 := + 1 6 12 18 ; +REAL rradf1 rrads1 rpinp1 apinp1 := + 0.6122 0.6540 0.0 0.0 ; +REAL rradf2 rrads2 rpinp2 apinp2 := + rradf1 rrads1 1.4885 0.0 ; +REAL rradf3 rrads3 rpinp3 apinp3 := + rradf1 rrads1 2.8755 $Pi_R npinp3 I_TO_R / ; +REAL rradf4 rrads4 rpinp4 apinp4 := + rradf1 rrads1 4.3305 0.0 ; +IF Ring 4 = THEN + :: <<npinp4>> <<rpinp4>> <<apinp4>> <<rradf4>> <<rrads4>> ; +ELSEIF Ring 3 = THEN + :: <<npinp3>> <<rpinp3>> <<apinp3>> <<rradf3>> <<rrads3>> ; +ELSEIF Ring 2 = THEN + :: <<npinp2>> <<rpinp2>> <<apinp2>> <<rradf2>> <<rrads2>> ; +ELSEIF Ring 1 = THEN + :: <<npinp1>> <<rpinp1>> <<apinp1>> <<rradf1>> <<rrads1>> ; +ELSE + ECHO "WARNING: Invalid Ring number in DimFuel : " Ring ; +ENDIF ; +QUIT . diff --git a/Dragon/data/proc_LZC/DimLZC.c2m b/Dragon/data/proc_LZC/DimLZC.c2m new file mode 100644 index 0000000..af18d3d --- /dev/null +++ b/Dragon/data/proc_LZC/DimLZC.c2m @@ -0,0 +1,52 @@ +*DECK DimLZC.c2m +*---- +* Name : DimLZC.c2m +* Type : DRAGON procedure +* Use : Dimensions for Gentilly-2 Liquid Zone Controller +* (cluster model) +* Author : G. Marleau +* Date : 2006/05/05 + +* Description of the procedure: +* DimLZC :: +* >>RinnGT<< >>RoutGT<< +* >>RinnBa<< >>RoutBa<< +* >>RinnFe<< >>RoutFe<< +* >>RinnBu<< >>RoutBu<< +* >>RinnSc<< >>RoutSc<< +* >>Rcent<< ; +* Output data: +* RinnGT : inner radius of guide tube +* RoutGT : outer radius of guide tube +* RinnBa : inner Balance tube radius +* RoutBa : outer Balance tube radius +* RinnFe : inner Feeder tube radius +* RoutFe : outer Feeder tube radius +* RinnBu : inner Bubbler tube radius +* RoutBu : outer Bubbler tube radius +* RinnSc : inner Scavanger tube radius +* RoutSc : outer Scavanger tube radius +* Rcent : annulus radius where the LZC are located +*---- +* Definition of the input and output data structures and +* of the input and output parameters for this procedure +*---- +REAL RinnGT RoutGT := + 11.430 2.0 / 11.684 2.0 / ; +REAL RinnBa RoutBa := + 1.814 2.0 / 1.966 2.0 / ; +REAL RinnFe RoutFe := + 2.972 2.0 / 3.124 2.0 / ; +REAL RinnBu RoutBu := + 0.483 2.0 / 0.635 2.0 / ; +REAL RinnSc RoutSc := + 1.608 2.0 / 1.760 2.0 / ; +REAL Rcent := + RinnGT 2.0 + 2.0 / ; +:: <<RinnGT>> <<RoutGT>> + <<RinnBa>> <<RoutBa>> + <<RinnFe>> <<RoutFe>> + <<RinnBu>> <<RoutBu>> + <<RinnSc>> <<RoutSc>> + <<Rcent>> ; +QUIT . diff --git a/Dragon/data/proc_LZC/GeoCell.c2m b/Dragon/data/proc_LZC/GeoCell.c2m new file mode 100755 index 0000000..97f14f9 --- /dev/null +++ b/Dragon/data/proc_LZC/GeoCell.c2m @@ -0,0 +1,306 @@ +*DECK GeoLZC.c2m +*---------- +* Name : GeoLZC.c2m +* Mesh : DRAGON procedure +* Use : LZC geometry +* for G2 (cluster model) +* Author : G. Marleau +* Date : 2006/05/05 +* +* Procedure description: +* Geometry := GeoLZC :: +* <<PrtLev>> <<Type>> <<Mesh>> +* <<NbAngles>> <<Density>> ; +* +* Procedure parameters: +* Geometry : XSM_FILE containing the LZC geometry +* for G2 +* PrtLev : Print Level +* 0 -> Minimum print +* 1 -> Standard print +* 10 -> Debug +* Type : Geometry type +* 'XY' -> 2-D projection of geometry in XY plane +* 'XZ' -> 2-D projection of geometry in XZ plane +* '3D' -> 3-D geometry (default value if Type +* is not known) +* Mesh : Geometry Mesh +* Coarse -> Corase mesh 3-D geometry +* Fine -> Fine mesh 3-D geometry (default) +* NbAngles : Suggested angular quadrature for tracking +* Density : Suggested line density for tracking +* +*---------- +* Define procedure parameters and options and read +* the procedure input data +*---- +PARAMETER Geometry :: EDIT 0 ::: XSM_FILE Geometry ; ; +STRING Type Mesh ; +INTEGER PrtLev NbAngles ; +REAL Density ; +:: >>PrtLev<< >>Type<< >>Mesh<< ; +*---- +* Local procedures and modules used in this procedure +*---- +PROCEDURE DimFuel DimCell DimLZC ; +MODULE GEO: END: ; +*---- +* Get overall cell dimensions using procedure DimCell +*---- +REAL rcool4 rpt rgap rct LPitch BLength ; +DimCell :: + >>rcool4<< >>rpt<< >>rgap<< >>rct<< >>LPitch<< >>BLength<< ; +IF PrtLev 9 > THEN + ECHO "Overall cell dimensions" ; + ECHO rcool4 rpt rgap rct LPitch BLength ; +ENDIF ; +*---- +* Get ring by ring description of fuel using procedure DimFuel +*---- +INTEGER Ring ; +INTEGER npinp1 npinp2 npinp3 npinp4 ; +REAL rradf1 rrads1 rpinp1 apinp1 ; +REAL rradf2 rrads2 rpinp2 apinp2 ; +REAL rradf3 rrads3 rpinp3 apinp3 ; +REAL rradf4 rrads4 rpinp4 apinp4 ; +EVALUATE Ring := 1 ; +DimFuel :: <<Ring>> + >>npinp1<< >>rpinp1<< >>apinp1<< >>rradf1<< >>rrads1<< ; +EVALUATE Ring := 2 ; +DimFuel :: <<Ring>> + >>npinp2<< >>rpinp2<< >>apinp2<< >>rradf2<< >>rrads2<< ; +EVALUATE Ring := 3 ; +DimFuel :: <<Ring>> + >>npinp3<< >>rpinp3<< >>apinp3<< >>rradf3<< >>rrads3<< ; +EVALUATE Ring := 4 ; +DimFuel :: <<Ring>> + >>npinp4<< >>rpinp4<< >>apinp4<< >>rradf4<< >>rrads4<< ; +IF PrtLev 9 > THEN + ECHO "Fuel rings 1 to 4" ; + ECHO npinp1 npinp2 npinp3 npinp4 ; + ECHO rradf1 rrads1 rpinp1 apinp1 ; + ECHO rradf2 rrads2 rpinp2 apinp2 ; + ECHO rradf3 rrads3 rpinp3 apinp3 ; + ECHO rradf4 rrads4 rpinp4 apinp4 ; +ENDIF ; +*---- +* Get LZC using procedure DimLZC +*---- +REAL RinnGT RoutGT + RinnBa RoutBa + RinnFe RoutFe + RinnBu RoutBu + RinnSc RoutSc + Rcent ; +DimLZC :: + >>RinnGT<< >>RoutGT<< + >>RinnBa<< >>RoutBa<< + >>RinnFe<< >>RoutFe<< + >>RinnBu<< >>RoutBu<< + >>RinnSc<< >>RoutSc<< + >>Rcent<< ; +IF PrtLev 9 > THEN + ECHO "LZC dimensions" ; + ECHO "Guide tube " RinnGT RoutGT ; + ECHO "Feeder tubes " RinnBa RoutBa RinnFe RoutFe ; + ECHO "Scavenger tubes " RinnBu RoutBu RinnSc RoutSc ; + ECHO "LZC location " Rcent ; +ENDIF ; +*---- +* Evaluate derived dimensions for supercell description +*---- +REAL HLPitch := LPitch 2.0 / ; +REAL HBLength := BLength 2.0 / ; +REAL LimCoCell := HLPitch rct - RoutGT + 2.0 / ; +REAL IntFCell := LPitch LimCoCell - ; +REAL IntYCell := HLPitch LimCoCell - ; +REAL OffsetXL OffsetXR := LimCoCell 2.0 / LimCoCell -2.0 / ; +* Left and right fuel cell X position +REAL XLFmin XLFint XLFcent XLFmax + XRFmin XRFint XRFcent XRFmax := + LPitch -1.0 * IntFCell -1.0 * HLPitch -1.0 * LimCoCell -1.0 * + LimCoCell IntFCell HLPitch LPitch ; +* Cell Y position +REAL YFmin YFintb YFcent YFintt YFmax := + XLFcent IntYCell -1.0 * 0.0 IntYCell XRFcent ; +* Z position +REAL ZFmin ZFintf ZFcent ZFintb ZFmax := + HBLength -1.0 * XLFmax 0.0 XRFmin HBLength ; +IF PrtLev 9 > THEN + ECHO "Derived dimensions for supercell" ; + ECHO HLPitch HBLength LimCoCell IntFCell OffsetXL OffsetXR ; + ECHO XLFmin XLFint XLFcent XLFmax ; + ECHO XRFmin XRFint XRFcent XRFmax ; + ECHO YFmin YFintb YFcent YFintt YFmax ; + ECHO ZFmin ZFintf ZFcent ZFintb ZFmax ; +ENDIF ; +*---- +* Derived dimensions for fuel +*---- +REAL eps1 := 0.1 ; +REAL rcool1 rcool2 rcool3 := + rrads1 eps1 + rpinp2 rrads2 - + 2.0 / + rpinp2 rrads2 + rpinp3 rrads3 - + 2.0 / + rpinp3 rrads3 + rpinp4 rrads4 - + 2.0 / ; +REAL rmod := HLPitch LimCoCell - eps1 - ; +IF PrtLev 9 > THEN + ECHO "Derived dimensions for fuel" ; + ECHO eps1 rcool1 rcool2 rcool3 rmod ; +ENDIF ; +*---- +* Define mixtures +*---- +INTEGER mcL1 mpt mgap mct mmod := 1 2 3 4 5 ; +INTEGER mfL1 mfL2 mfL3 mfL4 := 6 7 8 9 ; +INTEGER mfR1 mfR2 mfR3 mfR4 := 10 11 12 13 ; +INTEGER ms := 14 ; +INTEGER mfill mgt := 15 16 ; +INTEGER ms11 ms12 ms13 ms14 := + 17 18 19 20 ; +INTEGER ms21 ms22 ms23 ms24 := + ms11 4 + ms12 4 + ms13 4 + ms14 4 + ; +INTEGER ms31 ms32 ms33 ms34 := + ms11 8 + ms12 8 + ms13 8 + ms14 8 + ; +INTEGER mf11 mf12 mf13 mf14 := + ms11 12 + ms12 12 + ms13 12 + ms14 12 + ; +INTEGER mf21 mf22 mf23 mf24 := + ms11 16 + ms12 16 + ms13 16 + ms14 16 + ; +REAL PIO2 PIO6 := $Pi_R 2.0 / $Pi_R 6.0 / ; +REAL ApinZs1 ApinZs2 ApinZs3 + ApinZf1 ApinZf2 := + 0.0 PIO6 4.0 * PIO6 8.0 * + PIO6 2.0 * PIO6 6.0 * ; +REAL ApinXs1 ApinXs2 ApinXs3 + ApinXf1 ApinXf2 := + PIO2 ApinZs1 - PIO2 ApinZs2 - PIO2 ApinZs3 - + PIO2 ApinZf1 - PIO2 ApinZf2 - ; +ECHO ApinZs1 ApinZs2 ApinZs3 ApinZf1 ApinZf2 ; +IF PrtLev 9 > THEN + ECHO "Mixtures for coolant, PT, gap, CT and moderator" ; + ECHO mcL1 mpt mgap mct mmod ; + ECHO "Mixtures for left fuel in each ring and for sheathing" ; + ECHO mfL1 mfL2 mfL3 mfL4 ms ; + ECHO "Mixtures for right fuel in each ring and for sheathing" ; + ECHO mfR1 mfR2 mfR3 mfR4 ms ; + ECHO "Mixtures for scavenger tube 1" ; + ECHO ms11 ms12 ms13 ms14 ; + ECHO "Mixtures for scavenger tube 2" ; + ECHO ms21 ms22 ms23 ms24 ; + ECHO "Mixtures for scavenger tube 3" ; + ECHO ms31 ms32 ms33 ms34 ; + ECHO "Mixtures for feeder tube 1" ; + ECHO mf11 mf12 mf13 mf14 ; + ECHO "Mixtures for feeder tube 12" ; + ECHO mf11 mf22 mf23 mf24 ; +ENDIF ; +*---- +* Geometry +*---- + IF Mesh "Coarse" = THEN +*---- +* Coarse mesh geometry +*---- + Geometry := GEO: :: CAR2D 1 1 + X- REFL X+ SYME Y- REFL Y+ SYME + CELL FuelG +*---- +* Fuel on the left +*---- + ::: FuelG := GEO: CARCEL 7 1 1 + MESHX <<YFmin>> <<YFmax>> + MESHY <<YFmin>> <<YFmax>> + RADIUS 0.00000 <<rcool1>> <<rcool2>> + <<rcool3>> <<rcool4>> <<rpt>> + <<rgap>> <<rct>> + MIX <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> + MIX <<mfL1>> <<ms>> ; + ::: ROD2 := GEO: TUBE 2 + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> + MIX <<mfL2>> <<ms>> ; + ::: ROD3 := GEO: TUBE 2 + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> + MIX <<mfL3>> <<ms>> ; + ::: ROD4 := GEO: TUBE 2 + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> + MIX <<mfL4>> <<ms>> ; + ; + ; + EVALUATE NbAngles Density := 15 30.0 ; + ELSE +*---- +* Fine mesh geometry (default) +*---- + Geometry := GEO: :: CAR2D 1 1 + X- REFL X+ REFL Y- REFL Y+ REFL + CELL FuelG +*---- +* Fuel on the left +*---- + ::: FuelG := GEO: CARCEL 8 4 4 + MESHX <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + RADIUS 0.00000 <<rcool1>> <<rcool2>> + <<rcool3>> <<rcool4>> <<rpt>> + <<rgap>> <<rct>> <<rmod>> + MIX <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> SPLITR -3 1 + MIX <<mfL1>> <<ms>> ; + ::: ROD2 := GEO: TUBE 2 + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> SPLITR -3 1 + MIX <<mfL2>> <<ms>> ; + ::: ROD3 := GEO: TUBE 2 + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> SPLITR -3 1 + MIX <<mfL3>> <<ms>> ; + ::: ROD4 := GEO: TUBE 2 + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> SPLITR -3 1 + MIX <<mfL4>> <<ms>> ; + ; + ; + EVALUATE NbAngles Density := 15 30.0 ; + ENDIF ; +:: <<NbAngles>> <<Density>> ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/proc_LZC/GeoLZC.c2m b/Dragon/data/proc_LZC/GeoLZC.c2m new file mode 100755 index 0000000..45b89fe --- /dev/null +++ b/Dragon/data/proc_LZC/GeoLZC.c2m @@ -0,0 +1,1020 @@ +*DECK GeoLZC.c2m +*---------- +* Name : GeoLZC.c2m +* Mesh : DRAGON procedure +* Use : LZC geometry +* for G2 (cluster model) +* Author : G. Marleau +* Date : 2006/05/05 +* +* Procedure description: +* Geometry := GeoLZC :: +* <<PrtLev>> <<Type>> <<Mesh>> +* <<NbAngles>> <<Density>> ; +* +* Procedure parameters: +* Geometry : XSM_FILE containing the LZC geometry +* for G2 +* PrtLev : Print Level +* 0 -> Minimum print +* 1 -> Standard print +* 10 -> Debug +* Type : Geometry type +* 'XY' -> 2-D projection of geometry in XY plane +* 'XZ' -> 2-D projection of geometry in XZ plane +* '3D' -> 3-D geometry (default value if Type +* is not known) +* Mesh : Geometry Mesh +* Coarse -> Corase mesh 3-D geometry +* Fine -> Fine mesh 3-D geometry (default) +* NbAngles : Suggested angular quadrature for tracking +* Density : Suggested line density for tracking +* +*---------- +* Define procedure parameters and options and read +* the procedure input data +*---- +PARAMETER Geometry :: EDIT 0 ::: XSM_FILE Geometry ; ; +STRING Type Mesh ; +INTEGER PrtLev NbAngles ; +REAL Density ; +:: >>PrtLev<< >>Type<< >>Mesh<< ; +*---- +* Local procedures and modules used in this procedure +*---- +PROCEDURE DimFuel DimCell DimLZC ; +MODULE GEO: END: ; +*---- +* Get overall cell dimensions using procedure DimCell +*---- +REAL rcool4 rpt rgap rct LPitch BLength ; +DimCell :: + >>rcool4<< >>rpt<< >>rgap<< >>rct<< >>LPitch<< >>BLength<< ; +IF PrtLev 9 > THEN + ECHO "Overall cell dimensions" ; + ECHO rcool4 rpt rgap rct LPitch BLength ; +ENDIF ; +*---- +* Get ring by ring description of fuel using procedure DimFuel +*---- +INTEGER Ring ; +INTEGER npinp1 npinp2 npinp3 npinp4 ; +REAL rradf1 rrads1 rpinp1 apinp1 ; +REAL rradf2 rrads2 rpinp2 apinp2 ; +REAL rradf3 rrads3 rpinp3 apinp3 ; +REAL rradf4 rrads4 rpinp4 apinp4 ; +EVALUATE Ring := 1 ; +DimFuel :: <<Ring>> + >>npinp1<< >>rpinp1<< >>apinp1<< >>rradf1<< >>rrads1<< ; +EVALUATE Ring := 2 ; +DimFuel :: <<Ring>> + >>npinp2<< >>rpinp2<< >>apinp2<< >>rradf2<< >>rrads2<< ; +EVALUATE Ring := 3 ; +DimFuel :: <<Ring>> + >>npinp3<< >>rpinp3<< >>apinp3<< >>rradf3<< >>rrads3<< ; +EVALUATE Ring := 4 ; +DimFuel :: <<Ring>> + >>npinp4<< >>rpinp4<< >>apinp4<< >>rradf4<< >>rrads4<< ; +IF PrtLev 9 > THEN + ECHO "Fuel rings 1 to 4" ; + ECHO npinp1 npinp2 npinp3 npinp4 ; + ECHO rradf1 rrads1 rpinp1 apinp1 ; + ECHO rradf2 rrads2 rpinp2 apinp2 ; + ECHO rradf3 rrads3 rpinp3 apinp3 ; + ECHO rradf4 rrads4 rpinp4 apinp4 ; +ENDIF ; +*---- +* Get LZC using procedure DimLZC +*---- +REAL RinnGT RoutGT + RinnBa RoutBa + RinnFe RoutFe + RinnBu RoutBu + RinnSc RoutSc + Rcent ; +DimLZC :: + >>RinnGT<< >>RoutGT<< + >>RinnBa<< >>RoutBa<< + >>RinnFe<< >>RoutFe<< + >>RinnBu<< >>RoutBu<< + >>RinnSc<< >>RoutSc<< + >>Rcent<< ; +IF PrtLev 9 > THEN + ECHO "LZC dimensions" ; + ECHO "Guide tube " RinnGT RoutGT ; + ECHO "Feeder tubes " RinnBa RoutBa RinnFe RoutFe ; + ECHO "Scavenger tubes " RinnBu RoutBu RinnSc RoutSc ; + ECHO "LZC location " Rcent ; +ENDIF ; +*---- +* Evaluate derived dimensions for supercell description +*---- +REAL HLPitch := LPitch 2.0 / ; +REAL HBLength := BLength 2.0 / ; +REAL LimCoCell := HLPitch rct - RoutGT + 2.0 / ; +REAL IntFCell := LPitch LimCoCell - ; +REAL IntYCell := HLPitch LimCoCell - ; +REAL OffsetXL OffsetXR := LimCoCell 2.0 / LimCoCell -2.0 / ; +* Left and right fuel cell X position +REAL XLFmin XLFint XLFcent XLFmax + XRFmin XRFint XRFcent XRFmax := + LPitch -1.0 * IntFCell -1.0 * HLPitch -1.0 * LimCoCell -1.0 * + LimCoCell IntFCell HLPitch LPitch ; +* Cell Y position +REAL YFmin YFintb YFcent YFintt YFmax := + XLFcent IntYCell -1.0 * 0.0 IntYCell XRFcent ; +* Z position +REAL ZFmin ZFintf ZFcent ZFintb ZFmax := + HBLength -1.0 * XLFmax 0.0 XRFmin HBLength ; +IF PrtLev 9 > THEN + ECHO "Derived dimensions for supercell" ; + ECHO HLPitch HBLength LimCoCell IntFCell OffsetXL OffsetXR ; + ECHO XLFmin XLFint XLFcent XLFmax ; + ECHO XRFmin XRFint XRFcent XRFmax ; + ECHO YFmin YFintb YFcent YFintt YFmax ; + ECHO ZFmin ZFintf ZFcent ZFintb ZFmax ; +ENDIF ; +*---- +* Derived dimensions for fuel +*---- +REAL eps1 := 0.1 ; +REAL rcool1 rcool2 rcool3 := + rrads1 eps1 + rpinp2 rrads2 - + 2.0 / + rpinp2 rrads2 + rpinp3 rrads3 - + 2.0 / + rpinp3 rrads3 + rpinp4 rrads4 - + 2.0 / ; +REAL rmod := HLPitch LimCoCell - eps1 - ; +IF PrtLev 9 > THEN + ECHO "Derived dimensions for fuel" ; + ECHO eps1 rcool1 rcool2 rcool3 rmod ; +ENDIF ; +*---- +* Define mixtures +*---- +INTEGER mcL1 mpt mgap mct mmod := 1 2 3 4 5 ; +INTEGER mfL1 mfL2 mfL3 mfL4 := 6 7 8 9 ; +INTEGER mfR1 mfR2 mfR3 mfR4 := 10 11 12 13 ; +INTEGER ms := 14 ; +INTEGER mfill mgt := 15 16 ; +INTEGER ms11 ms12 ms13 ms14 := + 17 18 19 20 ; +INTEGER ms21 ms22 ms23 ms24 := + ms11 4 + ms12 4 + ms13 4 + ms14 4 + ; +INTEGER ms31 ms32 ms33 ms34 := + ms11 8 + ms12 8 + ms13 8 + ms14 8 + ; +INTEGER mf11 mf12 mf13 mf14 := + ms11 12 + ms12 12 + ms13 12 + ms14 12 + ; +INTEGER mf21 mf22 mf23 mf24 := + ms11 16 + ms12 16 + ms13 16 + ms14 16 + ; +REAL PIO2 PIO6 := $Pi_R 2.0 / $Pi_R 6.0 / ; +REAL ApinZs1 ApinZs2 ApinZs3 + ApinZf1 ApinZf2 := + 0.0 PIO6 4.0 * PIO6 8.0 * + PIO6 2.0 * PIO6 6.0 * ; +REAL ApinXs1 ApinXs2 ApinXs3 + ApinXf1 ApinXf2 := + PIO2 ApinZs1 - PIO2 ApinZs2 - PIO2 ApinZs3 - + PIO2 ApinZf1 - PIO2 ApinZf2 - ; +ECHO ApinZs1 ApinZs2 ApinZs3 ApinZf1 ApinZf2 ; +IF PrtLev 9 > THEN + ECHO "Mixtures for coolant, PT, gap, CT and moderator" ; + ECHO mcL1 mpt mgap mct mmod ; + ECHO "Mixtures for left fuel in each ring and for sheathing" ; + ECHO mfL1 mfL2 mfL3 mfL4 ms ; + ECHO "Mixtures for right fuel in each ring and for sheathing" ; + ECHO mfR1 mfR2 mfR3 mfR4 ms ; + ECHO "Mixtures for scavenger tube 1" ; + ECHO ms11 ms12 ms13 ms14 ; + ECHO "Mixtures for scavenger tube 2" ; + ECHO ms21 ms22 ms23 ms24 ; + ECHO "Mixtures for scavenger tube 3" ; + ECHO ms31 ms32 ms33 ms34 ; + ECHO "Mixtures for feeder tube 1" ; + ECHO mf11 mf12 mf13 mf14 ; + ECHO "Mixtures for feeder tube 12" ; + ECHO mf11 mf22 mf23 mf24 ; +ENDIF ; +*---- +* Geometry +*---- +IF Type "XY" = THEN + IF Mesh "Coarse" = THEN +*---- +* Coarse mesh geometry +*---- + Geometry := GEO: :: CAR2D 3 1 + X- REFL X+ REFL Y- REFL Y+ SYME + CELL FuelL LZC FuelR +*---- +* Fuel on the left +*---- + ::: FuelL := GEO: CARCEL 7 1 1 + MESHX <<XLFmin>> <<XLFmax>> + MESHY <<YFmin>> <<YFmax>> + RADIUS 0.00000 <<rcool1>> <<rcool2>> + <<rcool3>> <<rcool4>> <<rpt>> + <<rgap>> <<rct>> + OFFCENTER <<OffsetXL>> 0.0 0.0 + MIX <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> + MIX <<mfL1>> <<ms>> ; + ::: ROD2 := GEO: TUBE 2 + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> + MIX <<mfL2>> <<ms>> ; + ::: ROD3 := GEO: TUBE 2 + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> + MIX <<mfL3>> <<ms>> ; + ::: ROD4 := GEO: TUBE 2 + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> + MIX <<mfL4>> <<ms>> ; + ; +*---- +* LZC +*---- + ::: LZC := GEO: CAR2D 1 1 + MESHX <<XLFmax>> <<XRFmin>> + MESHY <<YFmin>> <<YFmax>> + MIX <<mmod>> + ; +*---- +* Fuel on the right +*---- + ::: FuelR := GEO: CARCEL 7 1 1 + MESHX <<XRFmin>> <<XRFmax>> + MESHY <<YFmin>> <<YFmax>> + RADIUS 0.00000 <<rcool1>> <<rcool2>> + <<rcool3>> <<rcool4>> <<rpt>> + <<rgap>> <<rct>> + OFFCENTER <<OffsetXR>> 0.0 0.0 + MIX <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> + MIX <<mfR1>> <<ms>> ; + ::: ROD2 := GEO: TUBE 2 + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> + MIX <<mfR2>> <<ms>> ; + ::: ROD3 := GEO: TUBE 2 + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> + MIX <<mfR3>> <<ms>> ; + ::: ROD4 := GEO: TUBE 2 + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> + MIX <<mfR4>> <<ms>> ; + ; + ; + EVALUATE NbAngles Density := 15 30.0 ; + ELSE +*---- +* Fine mesh geometry (default) +*---- + Geometry := GEO: :: CAR2D 3 1 + X- REFL X+ REFL Y- REFL Y+ SYME + CELL FuelL LZC FuelR +*---- +* Fuel on the left +*---- + ::: FuelL := GEO: CARCEL 8 3 4 + MESHX <<XLFmin>> <<XLFint>> <<XLFcent>> <<XLFmax>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + RADIUS 0.00000 <<rcool1>> <<rcool2>> + <<rcool3>> <<rcool4>> <<rpt>> + <<rgap>> <<rct>> <<rmod>> + OFFCENTER <<OffsetXL>> 0.0 0.0 + MIX <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> SPLITR -3 1 + MIX <<mfL1>> <<ms>> ; + ::: ROD2 := GEO: TUBE 2 + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> SPLITR -3 1 + MIX <<mfL2>> <<ms>> ; + ::: ROD3 := GEO: TUBE 2 + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> SPLITR -3 1 + MIX <<mfL3>> <<ms>> ; + ::: ROD4 := GEO: TUBE 2 + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> SPLITR -3 1 + MIX <<mfL4>> <<ms>> ; + ; +*---- +* LZC +*---- + ::: LZC := GEO: CAR2D 2 4 + MESHX <<XLFmax>> 0.0 <<XRFmin>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + MIX <<mmod>> <<mmod>> + <<mmod>> <<mmod>> + <<mmod>> <<mmod>> + <<mmod>> <<mmod>> + ; +*---- +* Fuel on the right +*---- + ::: FuelR := GEO: CARCEL 8 3 4 + MESHX <<XRFmin>> <<XRFcent>> <<XRFint>> <<XRFmax>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + RADIUS 0.00000 <<rcool1>> <<rcool2>> + <<rcool3>> <<rcool4>> <<rpt>> + <<rgap>> <<rct>> <<rmod>> + OFFCENTER <<OffsetXR>> 0.0 0.0 + MIX <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> SPLITR -3 1 + MIX <<mfR1>> <<ms>> ; + ::: ROD2 := GEO: TUBE 2 + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> SPLITR -3 1 + MIX <<mfR2>> <<ms>> ; + ::: ROD3 := GEO: TUBE 2 + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> SPLITR -3 1 + MIX <<mfR3>> <<ms>> ; + ::: ROD4 := GEO: TUBE 2 + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> SPLITR -3 1 + MIX <<mfR4>> <<ms>> ; + ; + ; + EVALUATE NbAngles Density := 15 30.0 ; + ENDIF ; +ELSEIF Type "XZ" = THEN + IF Mesh "Coarse" = THEN +*---- +* Coarse mesh geometry +*---- + Geometry := GEO: :: CAR2D 3 1 + X- REFL X+ REFL Y- REFL Y+ REFL + CELL FuelL LZC FuelR +*---- +* Fuel on the left +*---- + ::: FuelL := GEO: CAR2D 1 1 + MESHX <<XLFmin>> <<XLFmax>> + MESHY <<ZFmin>> <<ZFmax>> + MIX <<mmod>> + ; +*---- +* LZC +*---- + ::: LZC := GEO: CARCEL 2 1 1 + MESHX <<XLFmax>> <<XRFmin>> + MESHY <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<RinnGT>> <<RoutGT>> + MIX <<mfill>> <<mgt>> <<mmod>> + CLUSTER S1 S2 S3 F1 F2 + ::: S1 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXs1>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms11>> <<ms12>> <<ms13>> <<ms14>> ; + ::: S2 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXs2>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms21>> <<ms22>> <<ms23>> <<ms24>> ; + ::: S3 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXs3>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms31>> <<ms32>> <<ms33>> <<ms34>> ; + ::: F1 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXf1>> + RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>> + MIX <<mf11>> <<mf12>> <<mf13>> <<mf14>> ; + ::: F2 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXf2>> + RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>> + MIX <<mf21>> <<mf22>> <<mf23>> <<mf24>> ; + ; +*---- +* Fuel on the right +*---- + ::: FuelR := GEO: CAR2D 1 1 + MESHX <<XRFmin>> <<XRFmax>> + MESHY <<ZFmin>> <<ZFmax>> + MIX <<mmod>> + ; + ; + EVALUATE NbAngles Density := 15 30.0 ; + ELSE +*---- +* Fine mesh geometry (default) +*---- + Geometry := GEO: :: CAR2D 3 1 + X- REFL X+ REFL Y- REFL Y+ REFL + CELL FuelL LZC FuelR +*---- +* Fuel on the left +*---- + ::: FuelL := GEO: CAR2D 3 4 + MESHX <<XLFmin>> <<XLFint>> <<XLFcent>> <<XLFmax>> + MESHY <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + MIX <<mmod>> <<mmod>> <<mmod>> + <<mmod>> <<mmod>> <<mmod>> + <<mmod>> <<mmod>> <<mmod>> + <<mmod>> <<mmod>> <<mmod>> + ; +*---- +* LZC +*---- + ::: LZC := GEO: CARCEL 2 2 4 + MESHX <<XLFmax>> 0.0 <<XRFmin>> + MESHY <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<RinnGT>> <<RoutGT>> + MIX <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + CLUSTER S1 S2 S3 F1 F2 + ::: S1 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXs1>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms11>> <<ms12>> <<ms13>> <<ms14>> ; + ::: S2 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXs2>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms21>> <<ms22>> <<ms23>> <<ms24>> ; + ::: S3 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXs3>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms31>> <<ms32>> <<ms33>> <<ms34>> ; + ::: F1 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXf1>> + RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>> + MIX <<mf11>> <<mf12>> <<mf13>> <<mf14>> ; + ::: F2 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXf2>> + RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>> + MIX <<mf21>> <<mf22>> <<mf23>> <<mf24>> ; + ; +*---- +* Fuel on the right +*---- + ::: FuelR := GEO: CAR2D 3 4 + MESHX <<XRFmin>> <<XRFcent>> <<XRFint>> <<XRFmax>> + MESHY <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + MIX <<mmod>> <<mmod>> <<mmod>> + <<mmod>> <<mmod>> <<mmod>> + <<mmod>> <<mmod>> <<mmod>> + <<mmod>> <<mmod>> <<mmod>> + ; + ; + EVALUATE NbAngles Density := 15 30.0 ; + ENDIF ; +ELSE + IF Mesh "Coarse" = THEN +*---- +* Coarse mesh geometry +*---- + Geometry := GEO: :: CAR3D 3 1 1 + X- REFL X+ REFL Y- REFL Y+ SYME Z- REFL Z+ REFL + CELL FuelL LZC FuelR +*---- +* Fuel on the left +*---- + ::: FuelL := GEO: CARCELZ 7 1 1 1 + MESHX <<XLFmin>> <<XLFmax>> + MESHY <<YFmin>> <<YFmax>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rcool1>> <<rcool2>> + <<rcool3>> <<rcool4>> <<rpt>> + <<rgap>> <<rct>> + OFFCENTER <<OffsetXL>> 0.0 0.0 + MIX <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBEZ 2 1 + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> + MIX <<mfL1>> <<ms>> ; + ::: ROD2 := GEO: TUBEZ 2 1 + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> + MIX <<mfL2>> <<ms>> ; + ::: ROD3 := GEO: TUBEZ 2 1 + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> + MIX <<mfL3>> <<ms>> ; + ::: ROD4 := GEO: TUBEZ 2 1 + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> + MIX <<mfL4>> <<ms>> ; + ; +*---- +* LZC +*---- + ::: LZC := GEO: CARCELY 2 1 1 1 + MESHX <<XLFmax>> <<XRFmin>> + MESHY <<YFmin>> <<YFmax>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<RinnGT>> <<RoutGT>> + MIX <<mfill>> <<mgt>> <<mmod>> + CLUSTER S1 S2 S3 F1 F2 + ::: S1 := GEO: TUBEY 4 1 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZs1>> + MESHY <<YFmin>> <<YFmax>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms11>> <<ms12>> <<ms13>> <<ms14>> ; + ::: S2 := GEO: TUBEY 4 1 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZs2>> + MESHY <<YFmin>> <<YFmax>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms21>> <<ms22>> <<ms23>> <<ms24>> ; + ::: S3 := GEO: TUBEY 4 1 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZs3>> + MESHY <<YFmin>> <<YFmax>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms31>> <<ms32>> <<ms33>> <<ms34>> ; + ::: F1 := GEO: TUBEY 4 1 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZf1>> + MESHY <<YFmin>> <<YFmax>> + RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>> + MIX <<mf11>> <<mf12>> <<mf13>> <<mf14>> ; + ::: F2 := GEO: TUBEY 4 1 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZf2>> + MESHY <<YFmin>> <<YFmax>> + RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>> + MIX <<mf21>> <<mf22>> <<mf23>> <<mf24>> ; + ; +*---- +* Fuel on the right +*---- + ::: FuelR := GEO: CARCELZ 7 1 1 1 + MESHX <<XRFmin>> <<XRFmax>> + MESHY <<YFmin>> <<YFmax>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rcool1>> <<rcool2>> + <<rcool3>> <<rcool4>> <<rpt>> + <<rgap>> <<rct>> + OFFCENTER <<OffsetXR>> 0.0 0.0 + MIX <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBEZ 2 1 + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> + MIX <<mfR1>> <<ms>> ; + ::: ROD2 := GEO: TUBEZ 2 1 + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> + MIX <<mfR2>> <<ms>> ; + ::: ROD3 := GEO: TUBEZ 2 1 + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> + MIX <<mfR3>> <<ms>> ; + ::: ROD4 := GEO: TUBEZ 2 1 + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> + MIX <<mfR4>> <<ms>> ; + ; + ; + EVALUATE NbAngles Density := 8 100.0 ; + ELSE +*---- +* Fine mesh geometry (default) +*---- + Geometry := GEO: :: CAR3D 3 1 1 + X- REFL X+ REFL Y- REFL Y+ SYME Z- REFL Z+ REFL + CELL FuelL LZC FuelR +*---- +* Fuel on the left +*---- + ::: FuelL := GEO: CARCELZ 8 3 4 4 + MESHX <<XLFmin>> <<XLFint>> <<XLFcent>> <<XLFmax>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rcool1>> <<rcool2>> + <<rcool3>> <<rcool4>> <<rpt>> + <<rgap>> <<rct>> <<rmod>> + OFFCENTER <<OffsetXL>> 0.0 0.0 + MIX <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBEZ 2 4 + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> SPLITR -3 1 + MIX <<mfL1>> <<ms>> + <<mfL1>> <<ms>> + <<mfL1>> <<ms>> + <<mfL1>> <<ms>> ; + ::: ROD2 := GEO: TUBEZ 2 4 + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> SPLITR -3 1 + MIX <<mfL2>> <<ms>> + <<mfL2>> <<ms>> + <<mfL2>> <<ms>> + <<mfL2>> <<ms>> ; + ::: ROD3 := GEO: TUBEZ 2 4 + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> SPLITR -3 1 + MIX <<mfL3>> <<ms>> + <<mfL3>> <<ms>> + <<mfL3>> <<ms>> + <<mfL3>> <<ms>> ; + ::: ROD4 := GEO: TUBEZ 2 4 + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> SPLITR -3 1 + MIX <<mfL4>> <<ms>> + <<mfL4>> <<ms>> + <<mfL4>> <<ms>> + <<mfL4>> <<ms>> ; + ; +*---- +* LZC +*---- + ::: LZC := GEO: CARCELY 2 2 4 4 + MESHX <<XLFmax>> 0.0 <<XRFmin>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<RinnGT>> <<RoutGT>> + MIX <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + CLUSTER S1 S2 S3 F1 F2 + ::: S1 := GEO: TUBEY 4 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZs1>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms11>> <<ms12>> <<ms13>> <<ms14>> + <<ms11>> <<ms12>> <<ms13>> <<ms14>> + <<ms11>> <<ms12>> <<ms13>> <<ms14>> + <<ms11>> <<ms12>> <<ms13>> <<ms14>> ; + ::: S2 := GEO: TUBEY 4 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZs2>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms21>> <<ms22>> <<ms23>> <<ms24>> + <<ms21>> <<ms22>> <<ms23>> <<ms24>> + <<ms21>> <<ms22>> <<ms23>> <<ms24>> + <<ms21>> <<ms22>> <<ms23>> <<ms24>> ; + ::: S3 := GEO: TUBEY 4 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZs3>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms31>> <<ms32>> <<ms33>> <<ms34>> + <<ms31>> <<ms32>> <<ms33>> <<ms34>> + <<ms31>> <<ms32>> <<ms33>> <<ms34>> + <<ms31>> <<ms32>> <<ms33>> <<ms34>> ; + ::: F1 := GEO: TUBEY 4 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZf1>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>> + MIX <<mf11>> <<mf12>> <<mf13>> <<mf14>> + <<mf11>> <<mf12>> <<mf13>> <<mf14>> + <<mf11>> <<mf12>> <<mf13>> <<mf14>> + <<mf11>> <<mf12>> <<mf13>> <<mf14>> ; + ::: F2 := GEO: TUBEY 4 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZf2>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>> + MIX <<mf21>> <<mf22>> <<mf23>> <<mf24>> + <<mf21>> <<mf22>> <<mf23>> <<mf24>> + <<mf21>> <<mf22>> <<mf23>> <<mf24>> + <<mf21>> <<mf22>> <<mf23>> <<mf24>> ; + ; +*---- +* Fuel on the right +*---- + ::: FuelR := GEO: CARCELZ 8 3 4 4 + MESHX <<XRFmin>> <<XRFcent>> <<XRFint>> <<XRFmax>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rcool1>> <<rcool2>> + <<rcool3>> <<rcool4>> <<rpt>> + <<rgap>> <<rct>> <<rmod>> + OFFCENTER <<OffsetXR>> 0.0 0.0 + MIX <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBEZ 2 4 + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> SPLITR -3 1 + MIX <<mfR1>> <<ms>> + <<mfR1>> <<ms>> + <<mfR1>> <<ms>> + <<mfR1>> <<ms>> ; + ::: ROD2 := GEO: TUBEZ 2 4 + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> SPLITR -3 1 + MIX <<mfR2>> <<ms>> + <<mfR2>> <<ms>> + <<mfR2>> <<ms>> + <<mfR2>> <<ms>> ; + ::: ROD3 := GEO: TUBEZ 2 4 + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> SPLITR -3 1 + MIX <<mfR3>> <<ms>> + <<mfR3>> <<ms>> + <<mfR3>> <<ms>> + <<mfR3>> <<ms>> ; + ::: ROD4 := GEO: TUBEZ 2 4 + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> SPLITR -3 1 + MIX <<mfR4>> <<ms>> + <<mfR4>> <<ms>> + <<mfR4>> <<ms>> + <<mfR4>> <<ms>> ; + ; + ; + EVALUATE NbAngles Density := 16 1000.0 ; + ENDIF ; +ENDIF ; +:: <<NbAngles>> <<Density>> ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/proc_LZC/GeoLZCM.c2m b/Dragon/data/proc_LZC/GeoLZCM.c2m new file mode 100755 index 0000000..740d351 --- /dev/null +++ b/Dragon/data/proc_LZC/GeoLZCM.c2m @@ -0,0 +1,1020 @@ +*DECK GeoLZCM.c2m +*---------- +* Name : GeoLZCM.c2m +* Mesh : DRAGON procedure +* Use : LZC geometry +* for G2 (cluster model) +* Author : G. Marleau +* Date : 2006/05/05 +* +* Procedure description: +* Geometry := GeoLZCM :: +* <<PrtLev>> <<Type>> <<Mesh>> +* <<NbAngles>> <<Density>> ; +* +* Procedure parameters: +* Geometry : XSM_FILE containing the LZC geometry +* for G2 +* PrtLev : Print Level +* 0 -> Minimum print +* 1 -> Standard print +* 10 -> Debug +* Type : Geometry type +* 'XY' -> 2-D projection of geometry in XY plane +* 'XZ' -> 2-D projection of geometry in XZ plane +* '3D' -> 3-D geometry (default value if Type +* is not known) +* Mesh : Geometry Mesh +* Coarse -> Corase mesh 3-D geometry +* Fine -> Fine mesh 3-D geometry (default) +* NbAngles : Suggested angular quadrature for tracking +* Density : Suggested line density for tracking +* +*---------- +* Define procedure parameters and options and read +* the procedure input data +*---- +PARAMETER Geometry :: EDIT 0 ::: XSM_FILE Geometry ; ; +STRING Type Mesh ; +INTEGER PrtLev NbAngles ; +REAL Density ; +:: >>PrtLev<< >>Type<< >>Mesh<< ; +*---- +* Local procedures and modules used in this procedure +*---- +PROCEDURE DimFuel DimCell DimLZC ; +MODULE GEO: END: ; +*---- +* Get overall cell dimensions using procedure DimCell +*---- +REAL rcool4 rpt rgap rct LPitch BLength ; +DimCell :: + >>rcool4<< >>rpt<< >>rgap<< >>rct<< >>LPitch<< >>BLength<< ; +IF PrtLev 9 > THEN + ECHO "Overall cell dimensions" ; + ECHO rcool4 rpt rgap rct LPitch BLength ; +ENDIF ; +*---- +* Get ring by ring description of fuel using procedure DimFuel +*---- +INTEGER Ring ; +INTEGER npinp1 npinp2 npinp3 npinp4 ; +REAL rradf1 rrads1 rpinp1 apinp1 ; +REAL rradf2 rrads2 rpinp2 apinp2 ; +REAL rradf3 rrads3 rpinp3 apinp3 ; +REAL rradf4 rrads4 rpinp4 apinp4 ; +EVALUATE Ring := 1 ; +DimFuel :: <<Ring>> + >>npinp1<< >>rpinp1<< >>apinp1<< >>rradf1<< >>rrads1<< ; +EVALUATE Ring := 2 ; +DimFuel :: <<Ring>> + >>npinp2<< >>rpinp2<< >>apinp2<< >>rradf2<< >>rrads2<< ; +EVALUATE Ring := 3 ; +DimFuel :: <<Ring>> + >>npinp3<< >>rpinp3<< >>apinp3<< >>rradf3<< >>rrads3<< ; +EVALUATE Ring := 4 ; +DimFuel :: <<Ring>> + >>npinp4<< >>rpinp4<< >>apinp4<< >>rradf4<< >>rrads4<< ; +IF PrtLev 9 > THEN + ECHO "Fuel rings 1 to 4" ; + ECHO npinp1 npinp2 npinp3 npinp4 ; + ECHO rradf1 rrads1 rpinp1 apinp1 ; + ECHO rradf2 rrads2 rpinp2 apinp2 ; + ECHO rradf3 rrads3 rpinp3 apinp3 ; + ECHO rradf4 rrads4 rpinp4 apinp4 ; +ENDIF ; +*---- +* Get LZC using procedure DimLZC +*---- +REAL RinnGT RoutGT + RinnBa RoutBa + RinnFe RoutFe + RinnBu RoutBu + RinnSc RoutSc + Rcent ; +DimLZC :: + >>RinnGT<< >>RoutGT<< + >>RinnBa<< >>RoutBa<< + >>RinnFe<< >>RoutFe<< + >>RinnBu<< >>RoutBu<< + >>RinnSc<< >>RoutSc<< + >>Rcent<< ; +IF PrtLev 9 > THEN + ECHO "LZC dimensions" ; + ECHO "Guide tube " RinnGT RoutGT ; + ECHO "Feeder tubes " RinnBa RoutBa RinnFe RoutFe ; + ECHO "Scavenger tubes " RinnBu RoutBu RinnSc RoutSc ; + ECHO "LZC location " Rcent ; +ENDIF ; +*---- +* Evaluate derived dimensions for supercell description +*---- +REAL HLPitch := LPitch 2.0 / ; +REAL HBLength := BLength 2.0 / ; +REAL LimCoCell := HLPitch rct - RoutGT + 2.0 / ; +REAL IntFCell := LPitch LimCoCell - ; +REAL IntYCell := HLPitch LimCoCell - ; +REAL OffsetXL OffsetXR := LimCoCell 2.0 / LimCoCell -2.0 / ; +* Left and right fuel cell X position +REAL XLFmin XLFint XLFcent XLFmax + XRFmin XRFint XRFcent XRFmax := + LPitch -1.0 * IntFCell -1.0 * HLPitch -1.0 * LimCoCell -1.0 * + LimCoCell IntFCell HLPitch LPitch ; +* Cell Y position +REAL YFmin YFintb YFcent YFintt YFmax := + XLFcent IntYCell -1.0 * 0.0 IntYCell XRFcent ; +* Z position +REAL ZFmin ZFintf ZFcent ZFintb ZFmax := + HBLength -1.0 * XLFmax 0.0 XRFmin HBLength ; +IF PrtLev 9 > THEN + ECHO "Derived dimensions for supercell" ; + ECHO HLPitch HBLength LimCoCell IntFCell OffsetXL OffsetXR ; + ECHO XLFmin XLFint XLFcent XLFmax ; + ECHO XRFmin XRFint XRFcent XRFmax ; + ECHO YFmin YFintb YFcent YFintt YFmax ; + ECHO ZFmin ZFintf ZFcent ZFintb ZFmax ; +ENDIF ; +*---- +* Derived dimensions for fuel +*---- +REAL eps1 := 0.1 ; +REAL rcool1 rcool2 rcool3 := + rrads1 eps1 + rpinp2 rrads2 - + 2.0 / + rpinp2 rrads2 + rpinp3 rrads3 - + 2.0 / + rpinp3 rrads3 + rpinp4 rrads4 - + 2.0 / ; +REAL rmod := HLPitch LimCoCell - eps1 - ; +IF PrtLev 9 > THEN + ECHO "Derived dimensions for fuel" ; + ECHO eps1 rcool1 rcool2 rcool3 rmod ; +ENDIF ; +*---- +* Define mixtures +*---- +INTEGER mcL1 mpt mgap mct mmod := 1 2 3 4 5 ; +INTEGER mfL1 mfL2 mfL3 mfL4 := 6 7 8 9 ; +INTEGER mfR1 mfR2 mfR3 mfR4 := 10 11 12 13 ; +INTEGER ms := 14 ; +INTEGER mfill mgt := 15 16 ; +INTEGER ms11 ms12 ms13 ms14 := + 17 18 19 20 ; +INTEGER ms21 ms22 ms23 ms24 := + ms11 4 + ms12 4 + ms13 4 + ms14 4 + ; +INTEGER ms31 ms32 ms33 ms34 := + ms11 8 + ms12 8 + ms13 8 + ms14 8 + ; +INTEGER mf11 mf12 mf13 mf14 := + ms11 12 + ms12 12 + ms13 12 + ms14 12 + ; +INTEGER mf21 mf22 mf23 mf24 := + ms11 16 + ms12 16 + ms13 16 + ms14 16 + ; +REAL PIO2 PIO6 := $Pi_R 2.0 / $Pi_R 6.0 / ; +REAL ApinZs1 ApinZs2 ApinZs3 + ApinZf1 ApinZf2 := + 0.0 PIO6 4.0 * PIO6 8.0 * + PIO6 2.0 * PIO6 6.0 * ; +REAL ApinXs1 ApinXs2 ApinXs3 + ApinXf1 ApinXf2 := + PIO2 ApinZs1 - PIO2 ApinZs2 - PIO2 ApinZs3 - + PIO2 ApinZf1 - PIO2 ApinZf2 - ; +ECHO ApinZs1 ApinZs2 ApinZs3 ApinZf1 ApinZf2 ; +IF PrtLev 9 > THEN + ECHO "Mixtures for coolant, PT, gap, CT and moderator" ; + ECHO mcL1 mpt mgap mct mmod ; + ECHO "Mixtures for left fuel in each ring and for sheathing" ; + ECHO mfL1 mfL2 mfL3 mfL4 ms ; + ECHO "Mixtures for right fuel in each ring and for sheathing" ; + ECHO mfR1 mfR2 mfR3 mfR4 ms ; + ECHO "Mixtures for scavenger tube 1" ; + ECHO ms11 ms12 ms13 ms14 ; + ECHO "Mixtures for scavenger tube 2" ; + ECHO ms21 ms22 ms23 ms24 ; + ECHO "Mixtures for scavenger tube 3" ; + ECHO ms31 ms32 ms33 ms34 ; + ECHO "Mixtures for feeder tube 1" ; + ECHO mf11 mf12 mf13 mf14 ; + ECHO "Mixtures for feeder tube 12" ; + ECHO mf11 mf22 mf23 mf24 ; +ENDIF ; +*---- +* Geometry +*---- +IF Type "XY" = THEN + IF Mesh "Coarse" = THEN +*---- +* Coarse mesh geometry +*---- + Geometry := GEO: :: CAR2D 3 1 + X- REFL X+ REFL Y- REFL Y+ SYME + CELL FuelL LZC FuelR +*---- +* Fuel on the left +*---- + ::: FuelL := GEO: CARCEL 7 1 1 + MESHX <<XLFmin>> <<XLFmax>> + MESHY <<YFmin>> <<YFmax>> + RADIUS 0.00000 <<rcool1>> <<rcool2>> + <<rcool3>> <<rcool4>> <<rpt>> + <<rgap>> <<rct>> + OFFCENTER <<OffsetXL>> 0.0 0.0 + MIX <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> + MIX <<mfL1>> <<ms>> ; + ::: ROD2 := GEO: TUBE 2 + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> + MIX <<mfL2>> <<ms>> ; + ::: ROD3 := GEO: TUBE 2 + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> + MIX <<mfL3>> <<ms>> ; + ::: ROD4 := GEO: TUBE 2 + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> + MIX <<mfL4>> <<ms>> ; + ; +*---- +* LZC +*---- + ::: LZC := GEO: CAR2D 1 1 + MESHX <<XLFmax>> <<XRFmin>> + MESHY <<YFmin>> <<YFmax>> + MIX <<mmod>> + ; +*---- +* Fuel on the right +*---- + ::: FuelR := GEO: CARCEL 7 1 1 + MESHX <<XRFmin>> <<XRFmax>> + MESHY <<YFmin>> <<YFmax>> + RADIUS 0.00000 <<rcool1>> <<rcool2>> + <<rcool3>> <<rcool4>> <<rpt>> + <<rgap>> <<rct>> + OFFCENTER <<OffsetXR>> 0.0 0.0 + MIX <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> + MIX <<mfR1>> <<ms>> ; + ::: ROD2 := GEO: TUBE 2 + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> + MIX <<mfR2>> <<ms>> ; + ::: ROD3 := GEO: TUBE 2 + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> + MIX <<mfR3>> <<ms>> ; + ::: ROD4 := GEO: TUBE 2 + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> + MIX <<mfR4>> <<ms>> ; + ; + ; + EVALUATE NbAngles Density := 15 30.0 ; + ELSE +*---- +* Fine mesh geometry (default) +*---- + Geometry := GEO: :: CAR2D 3 1 + X- REFL X+ REFL Y- REFL Y+ SYME + CELL FuelL LZC FuelR +*---- +* Fuel on the left +*---- + ::: FuelL := GEO: CARCEL 8 3 4 + MESHX <<XLFmin>> <<XLFint>> <<XLFcent>> <<XLFmax>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + RADIUS 0.00000 <<rcool1>> <<rcool2>> + <<rcool3>> <<rcool4>> <<rpt>> + <<rgap>> <<rct>> <<rmod>> + OFFCENTER <<OffsetXL>> 0.0 0.0 + MIX <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> + MIX <<mfL1>> <<ms>> ; + ::: ROD2 := GEO: TUBE 2 + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> + MIX <<mfL2>> <<ms>> ; + ::: ROD3 := GEO: TUBE 2 + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> + MIX <<mfL3>> <<ms>> ; + ::: ROD4 := GEO: TUBE 2 + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> + MIX <<mfL4>> <<ms>> ; + ; +*---- +* LZC +*---- + ::: LZC := GEO: CAR2D 2 4 + MESHX <<XLFmax>> 0.0 <<XRFmin>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + MIX <<mmod>> <<mmod>> + <<mmod>> <<mmod>> + <<mmod>> <<mmod>> + <<mmod>> <<mmod>> + ; +*---- +* Fuel on the right +*---- + ::: FuelR := GEO: CARCEL 8 3 4 + MESHX <<XRFmin>> <<XRFcent>> <<XRFint>> <<XRFmax>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + RADIUS 0.00000 <<rcool1>> <<rcool2>> + <<rcool3>> <<rcool4>> <<rpt>> + <<rgap>> <<rct>> <<rmod>> + OFFCENTER <<OffsetXR>> 0.0 0.0 + MIX <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> + MIX <<mfR1>> <<ms>> ; + ::: ROD2 := GEO: TUBE 2 + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> + MIX <<mfR2>> <<ms>> ; + ::: ROD3 := GEO: TUBE 2 + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> + MIX <<mfR3>> <<ms>> ; + ::: ROD4 := GEO: TUBE 2 + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> + MIX <<mfR4>> <<ms>> ; + ; + ; + EVALUATE NbAngles Density := 15 30.0 ; + ENDIF ; +ELSEIF Type "XZ" = THEN + IF Mesh "Coarse" = THEN +*---- +* Coarse mesh geometry +*---- + Geometry := GEO: :: CAR2D 3 1 + X- REFL X+ REFL Y- REFL Y+ REFL + CELL FuelL LZC FuelR +*---- +* Fuel on the left +*---- + ::: FuelL := GEO: CAR2D 1 1 + MESHX <<XLFmin>> <<XLFmax>> + MESHY <<ZFmin>> <<ZFmax>> + MIX <<mmod>> + ; +*---- +* LZC +*---- + ::: LZC := GEO: CARCEL 2 1 1 + MESHX <<XLFmax>> <<XRFmin>> + MESHY <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<RinnGT>> <<RoutGT>> + MIX <<mfill>> <<mgt>> <<mmod>> + CLUSTER S1 S2 S3 F1 F2 + ::: S1 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXs1>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms11>> <<ms12>> <<ms13>> <<ms14>> ; + ::: S2 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXs2>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms21>> <<ms22>> <<ms23>> <<ms24>> ; + ::: S3 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXs3>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms31>> <<ms32>> <<ms33>> <<ms34>> ; + ::: F1 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXf1>> + RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>> + MIX <<mf11>> <<mf12>> <<mf13>> <<mf14>> ; + ::: F2 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXf2>> + RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>> + MIX <<mf21>> <<mf22>> <<mf23>> <<mf24>> ; + ; +*---- +* Fuel on the right +*---- + ::: FuelR := GEO: CAR2D 1 1 + MESHX <<XRFmin>> <<XRFmax>> + MESHY <<ZFmin>> <<ZFmax>> + MIX <<mmod>> + ; + ; + EVALUATE NbAngles Density := 15 30.0 ; + ELSE +*---- +* Fine mesh geometry (default) +*---- + Geometry := GEO: :: CAR2D 3 1 + X- REFL X+ REFL Y- REFL Y+ REFL + CELL FuelL LZC FuelR +*---- +* Fuel on the left +*---- + ::: FuelL := GEO: CAR2D 3 4 + MESHX <<XLFmin>> <<XLFint>> <<XLFcent>> <<XLFmax>> + MESHY <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + MIX <<mmod>> <<mmod>> <<mmod>> + <<mmod>> <<mmod>> <<mmod>> + <<mmod>> <<mmod>> <<mmod>> + <<mmod>> <<mmod>> <<mmod>> + ; +*---- +* LZC +*---- + ::: LZC := GEO: CARCEL 2 2 4 + MESHX <<XLFmax>> 0.0 <<XRFmin>> + MESHY <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<RinnGT>> <<RoutGT>> + MIX <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + CLUSTER S1 S2 S3 F1 F2 + ::: S1 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXs1>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms11>> <<ms12>> <<ms13>> <<ms14>> ; + ::: S2 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXs2>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms21>> <<ms22>> <<ms23>> <<ms24>> ; + ::: S3 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXs3>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms31>> <<ms32>> <<ms33>> <<ms34>> ; + ::: F1 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXf1>> + RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>> + MIX <<mf11>> <<mf12>> <<mf13>> <<mf14>> ; + ::: F2 := GEO: TUBE 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinXf2>> + RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>> + MIX <<mf21>> <<mf22>> <<mf23>> <<mf24>> ; + ; +*---- +* Fuel on the right +*---- + ::: FuelR := GEO: CAR2D 3 4 + MESHX <<XRFmin>> <<XRFcent>> <<XRFint>> <<XRFmax>> + MESHY <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + MIX <<mmod>> <<mmod>> <<mmod>> + <<mmod>> <<mmod>> <<mmod>> + <<mmod>> <<mmod>> <<mmod>> + <<mmod>> <<mmod>> <<mmod>> + ; + ; + EVALUATE NbAngles Density := 15 30.0 ; + ENDIF ; +ELSE + IF Mesh "Coarse" = THEN +*---- +* Coarse mesh geometry +*---- + Geometry := GEO: :: CAR3D 3 1 1 + X- REFL X+ REFL Y- REFL Y+ SYME Z- REFL Z+ REFL + CELL FuelL LZC FuelR +*---- +* Fuel on the left +*---- + ::: FuelL := GEO: CARCELZ 7 1 1 1 + MESHX <<XLFmin>> <<XLFmax>> + MESHY <<YFmin>> <<YFmax>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rcool1>> <<rcool2>> + <<rcool3>> <<rcool4>> <<rpt>> + <<rgap>> <<rct>> + OFFCENTER <<OffsetXL>> 0.0 0.0 + MIX <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBEZ 2 1 + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> + MIX <<mfL1>> <<ms>> ; + ::: ROD2 := GEO: TUBEZ 2 1 + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> + MIX <<mfL2>> <<ms>> ; + ::: ROD3 := GEO: TUBEZ 2 1 + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> + MIX <<mfL3>> <<ms>> ; + ::: ROD4 := GEO: TUBEZ 2 1 + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> + MIX <<mfL4>> <<ms>> ; + ; +*---- +* LZC +*---- + ::: LZC := GEO: CARCELY 2 1 1 1 + MESHX <<XLFmax>> <<XRFmin>> + MESHY <<YFmin>> <<YFmax>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<RinnGT>> <<RoutGT>> + MIX <<mfill>> <<mgt>> <<mmod>> + CLUSTER S1 S2 S3 F1 F2 + ::: S1 := GEO: TUBEY 4 1 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZs1>> + MESHY <<YFmin>> <<YFmax>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms11>> <<ms12>> <<ms13>> <<ms14>> ; + ::: S2 := GEO: TUBEY 4 1 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZs2>> + MESHY <<YFmin>> <<YFmax>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms21>> <<ms22>> <<ms23>> <<ms24>> ; + ::: S3 := GEO: TUBEY 4 1 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZs3>> + MESHY <<YFmin>> <<YFmax>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms31>> <<ms32>> <<ms33>> <<ms34>> ; + ::: F1 := GEO: TUBEY 4 1 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZf1>> + MESHY <<YFmin>> <<YFmax>> + RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>> + MIX <<mf11>> <<mf12>> <<mf13>> <<mf14>> ; + ::: F2 := GEO: TUBEY 4 1 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZf2>> + MESHY <<YFmin>> <<YFmax>> + RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>> + MIX <<mf21>> <<mf22>> <<mf23>> <<mf24>> ; + ; +*---- +* Fuel on the right +*---- + ::: FuelR := GEO: CARCELZ 7 1 1 1 + MESHX <<XRFmin>> <<XRFmax>> + MESHY <<YFmin>> <<YFmax>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rcool1>> <<rcool2>> + <<rcool3>> <<rcool4>> <<rpt>> + <<rgap>> <<rct>> + OFFCENTER <<OffsetXR>> 0.0 0.0 + MIX <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBEZ 2 1 + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> + MIX <<mfR1>> <<ms>> ; + ::: ROD2 := GEO: TUBEZ 2 1 + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> + MIX <<mfR2>> <<ms>> ; + ::: ROD3 := GEO: TUBEZ 2 1 + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> + MIX <<mfR3>> <<ms>> ; + ::: ROD4 := GEO: TUBEZ 2 1 + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + MESHZ <<ZFmin>> <<ZFmax>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> + MIX <<mfR4>> <<ms>> ; + ; + ; + EVALUATE NbAngles Density := 8 50.0 ; + ELSE +*---- +* Fine mesh geometry (default) +*---- + Geometry := GEO: :: CAR3D 3 1 1 + X- REFL X+ REFL Y- REFL Y+ SYME Z- REFL Z+ REFL + CELL FuelL LZC FuelR +*---- +* Fuel on the left +*---- + ::: FuelL := GEO: CARCELZ 8 3 4 4 + MESHX <<XLFmin>> <<XLFint>> <<XLFcent>> <<XLFmax>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rcool1>> <<rcool2>> + <<rcool3>> <<rcool4>> <<rpt>> + <<rgap>> <<rct>> <<rmod>> + OFFCENTER <<OffsetXL>> 0.0 0.0 + MIX <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBEZ 2 4 + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> + MIX <<mfL1>> <<ms>> + <<mfL1>> <<ms>> + <<mfL1>> <<ms>> + <<mfL1>> <<ms>> ; + ::: ROD2 := GEO: TUBEZ 2 4 + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> + MIX <<mfL2>> <<ms>> + <<mfL2>> <<ms>> + <<mfL2>> <<ms>> + <<mfL2>> <<ms>> ; + ::: ROD3 := GEO: TUBEZ 2 4 + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> + MIX <<mfL3>> <<ms>> + <<mfL3>> <<ms>> + <<mfL3>> <<ms>> + <<mfL3>> <<ms>> ; + ::: ROD4 := GEO: TUBEZ 2 4 + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> + MIX <<mfL4>> <<ms>> + <<mfL4>> <<ms>> + <<mfL4>> <<ms>> + <<mfL4>> <<ms>> ; + ; +*---- +* LZC +*---- + ::: LZC := GEO: CARCELY 2 2 4 4 + MESHX <<XLFmax>> 0.0 <<XRFmin>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<RinnGT>> <<RoutGT>> + MIX <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + <<mfill>> <<mgt>> <<mmod>> + CLUSTER S1 S2 S3 F1 F2 + ::: S1 := GEO: TUBEY 4 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZs1>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms11>> <<ms12>> <<ms13>> <<ms14>> + <<ms11>> <<ms12>> <<ms13>> <<ms14>> + <<ms11>> <<ms12>> <<ms13>> <<ms14>> + <<ms11>> <<ms12>> <<ms13>> <<ms14>> ; + ::: S2 := GEO: TUBEY 4 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZs2>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms21>> <<ms22>> <<ms23>> <<ms24>> + <<ms21>> <<ms22>> <<ms23>> <<ms24>> + <<ms21>> <<ms22>> <<ms23>> <<ms24>> + <<ms21>> <<ms22>> <<ms23>> <<ms24>> ; + ::: S3 := GEO: TUBEY 4 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZs3>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + RADIUS 0.0 <<RinnBu>> <<RoutBu>> <<RinnSc>> <<RoutSc>> + MIX <<ms31>> <<ms32>> <<ms33>> <<ms34>> + <<ms31>> <<ms32>> <<ms33>> <<ms34>> + <<ms31>> <<ms32>> <<ms33>> <<ms34>> + <<ms31>> <<ms32>> <<ms33>> <<ms34>> ; + ::: F1 := GEO: TUBEY 4 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZf1>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>> + MIX <<mf11>> <<mf12>> <<mf13>> <<mf14>> + <<mf11>> <<mf12>> <<mf13>> <<mf14>> + <<mf11>> <<mf12>> <<mf13>> <<mf14>> + <<mf11>> <<mf12>> <<mf13>> <<mf14>> ; + ::: F2 := GEO: TUBEY 4 4 + NPIN 1 RPIN <<Rcent>> APIN <<ApinZf2>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + RADIUS 0.0 <<RinnBa>> <<RoutBa>> <<RinnFe>> <<RoutFe>> + MIX <<mf21>> <<mf22>> <<mf23>> <<mf24>> + <<mf21>> <<mf22>> <<mf23>> <<mf24>> + <<mf21>> <<mf22>> <<mf23>> <<mf24>> + <<mf21>> <<mf22>> <<mf23>> <<mf24>> ; + ; +*---- +* Fuel on the right +*---- + ::: FuelR := GEO: CARCELZ 8 3 4 4 + MESHX <<XRFmin>> <<XRFcent>> <<XRFint>> <<XRFmax>> + MESHY <<YFmin>> <<YFintb>> <<YFcent>> <<YFintt>> <<YFmax>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rcool1>> <<rcool2>> + <<rcool3>> <<rcool4>> <<rpt>> + <<rgap>> <<rct>> <<rmod>> + OFFCENTER <<OffsetXR>> 0.0 0.0 + MIX <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + <<mcL1>> <<mcL1>> <<mcL1>> <<mcL1>> + <<mpt>> <<mgap>> <<mct>> <<mmod>> <<mmod>> + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBEZ 2 4 + NPIN <<npinp1>> RPIN <<rpinp1>> APIN <<apinp1>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rradf1>> <<rrads1>> + MIX <<mfR1>> <<ms>> + <<mfR1>> <<ms>> + <<mfR1>> <<ms>> + <<mfR1>> <<ms>> ; + ::: ROD2 := GEO: TUBEZ 2 4 + NPIN <<npinp2>> RPIN <<rpinp2>> APIN <<apinp2>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rradf2>> <<rrads2>> + MIX <<mfR2>> <<ms>> + <<mfR2>> <<ms>> + <<mfR2>> <<ms>> + <<mfR2>> <<ms>> ; + ::: ROD3 := GEO: TUBEZ 2 4 + NPIN <<npinp3>> RPIN <<rpinp3>> APIN <<apinp3>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rradf3>> <<rrads3>> + MIX <<mfR3>> <<ms>> + <<mfR3>> <<ms>> + <<mfR3>> <<ms>> + <<mfR3>> <<ms>> ; + ::: ROD4 := GEO: TUBEZ 2 4 + NPIN <<npinp4>> RPIN <<rpinp4>> APIN <<apinp4>> + MESHZ <<ZFmin>> <<ZFintf>> <<ZFcent>> <<ZFintb>> <<ZFmax>> + RADIUS 0.00000 <<rradf4>> <<rrads4>> + MIX <<mfR4>> <<ms>> + <<mfR4>> <<ms>> + <<mfR4>> <<ms>> + <<mfR4>> <<ms>> ; + ; + ; + EVALUATE NbAngles Density := 16 100.0 ; + ENDIF ; +ENDIF ; +:: <<NbAngles>> <<Density>> ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/proc_LZC/MicLZC.c2m b/Dragon/data/proc_LZC/MicLZC.c2m new file mode 100644 index 0000000..81181b4 --- /dev/null +++ b/Dragon/data/proc_LZC/MicLZC.c2m @@ -0,0 +1,98 @@ +*DECK MicLZC.c2m +*---------- +* Name : MicLZC.c2m +* Type : DRAGON procedure +* Use : Create Gentilly-2 library for LZC analysis. +* (cluster model) +* Author : G. Marleau +* Date : 2006/03/01 +* +* Procedure description: +* MicroLib := MicLZC :: +* <<PrtLev>> <<LibSelect>> <<Position>> <<Type>> ; +* +* MicroLib : XSM_FILE containing the microlib where: +* Mix 1 = coolant at full density +* Mix 2 = pressure tube +* Mix 3 = helium for gap +* Mix 4 = Zr for calandrial tube and adjuster +* Mix 5 = moderator +* Mix 6 = left central fuel +* Mix 7 = left 6 pin fuel +* Mix 8 = left 12 pin fuel +* Mix 9 = left outer fuel +* Mix 10 = right central fuel +* Mix 11 = right 6 pin fuel +* Mix 12 = right 12 pin fuel +* Mix 13 = right outer fuel +* Mix 14 = fuel sheating +* Mix 15 = LZC zone to fill or empty +* Mix 16 = Guide tube +* Mix 17 = mixture in bubbler for scavenger 1 +* Mix 18 = mixture around bubbler for scavenger 1 +* Mix 19 = mixture in scavenger 1 +* Mix 20 = mixture around scavenger 1 +* Mix 21 = mixture in bubbler for scavenger 2 +* Mix 22 = mixture around bubbler for scavenger 2 +* Mix 23 = mixture in scavenger 2 +* Mix 24 = mixture around scavenger 2 +* Mix 25 = mixture in bubbler for scavenger 3 +* Mix 26 = mixture around bubbler for scavenger 3 +* Mix 27 = mixture in scavenger 3 +* Mix 28 = mixture around scavenger 3 +* Mix 29 = mixture in balance for feeder 1 +* Mix 30 = mixture around balance for feeder 1 +* Mix 31 = mixture in feeder 1 +* Mix 32 = mixture around feeder 1 +* Mix 33 = mixture in balance for feeder 2 +* Mix 34 = mixture around balance for feeder 2 +* Mix 35 = mixture in feeder 2 +* Mix 36 = mixture around feeder 2 +* PrtLev : Print Level +* 0 -> Minimum print +* 1 -> Standard print +* 10 -> Debug +* LibSelect : Selection of library with +* WU -> WLUP 69 groups WIMS-AECL library +* (LibName = IAEA69). +* WL -> WIMFRITH 69 groups WIMS-AECL library +* (LibName = WIMSLIB). +* E5 -> ENDF/B-V 89 groups WIMS-AECL library +* (LibName = E5WLIB). +* E6 -> ENDF/B-VI 89 groups WIMS-AECL library +* (LibName = E6MLIB). +* Position : Cobalt adjusters position +* NO -> Adjuster and guide tubes absent +* OUT -> Adjuster absent and guide tubes present +* IN -> Adjuster and guide tubes present +* (default) +* Type : +* +*---------- +* Define procedure parameters and options and read +* the procedure input data +*---- +PARAMETER MicroLib :: EDIT 0 ::: XSM_FILE MicroLib ; ; +INTEGER PrtLev ; +STRING Position LibSelect Type ; +:: >>PrtLev<< >>LibSelect<< >>Position<< >>Type<< ; +*---------- +* Define procedures +*---------- +PROCEDURE MicWUC MicWULZC ; +*PROCEDURE MicWU MicWLD ; +*PROCEDURE MicE5C MicE5D ; +*PROCEDURE MicE6C MicE6D ; +*---- +* Get fuel XS +*---- +IF LibSelect "WU" = THEN + MicroLib := MicWUC :: + <<PrtLev>> ; +*---- +* Get LZC cross sections +*---- + MicroLib := MicWULZC MicroLib :: + <<PrtLev>> <<Position>> <<Type>> ; +ENDIF ; +QUIT . diff --git a/Dragon/data/proc_LZC/MicWUC.c2m b/Dragon/data/proc_LZC/MicWUC.c2m new file mode 100644 index 0000000..d24bbec --- /dev/null +++ b/Dragon/data/proc_LZC/MicWUC.c2m @@ -0,0 +1,189 @@ +*DECK MicWUC.c2m +*---------- +* Name : MicWUC.c2m +* Type : DRAGON procedure +* Use : Create Gentilly-2 fuel microlib +* with cross sections from the WLUP based 69 +* groups (LibName = IAEA69). +* Author : G. Marleau +* Date : 2006/03/01 +* +* Procedure description: +* MicroLib := MicWUC :: +* <<PrtLev>> ; +* +* Procedure parameters: +* MicroLib : XSM_FILE containing the new microlib where: +* Mix 1 = coolant at full density +* Mix 2 = pressure tube +* Mix 3 = helium for gap and LZC +* Mix 4 = Zr for calandrial tube and adjuster +* Mix 5 = moderator +* Mix 6 = left central fuel +* Mix 7 = left 6 pin fuel +* Mix 8 = left 12 pin fuel +* Mix 9 = left outer fuel +* Mix 10 = right central fuel +* Mix 11 = right 6 pin fuel +* Mix 12 = right 12 pin fuel +* Mix 13 = right outer fuel +* Mix 14 = fuel sheating +* (see GeoLZC.c2m for mixture location in +* geometry) +* PrtLev : Print Level +* 0 -> Minimum print +* 1 -> Standard print +* 10 -> Debug +* +*---------- +* Define procedure parameters and options and read +* the procedure input data +*---- +PARAMETER MicroLib :: EDIT 0 ::: XSM_FILE MicroLib ; ; +INTEGER PrtLev ; +:: >>PrtLev<< ; +*---- +* Local modules used in this procedure +*---- +MODULE INFO: LIB: DELETE: ; +*---------- +* Define local fuel, coolant and moderator properties +* Fuel enrichment (%) = 0.71140; T = 941.26 K; RHO =10.4375 g/cc; +* Coolant D2 AT (%) = 99.222; T = 560.66 K; +* Moderator D2 AT (%) = 99.911; T = 345.66 K. +*---------- +REAL + Ftemp Fenri Fdens + Ctemp Cprty + Mtemp Mprty + Wprty := + 941.26 0.7114 10.4375 + 560.66 99.30 + 345.66 99.92 + 0.0 ; +REAL + FwgtU8 FwgtU5 FwgtO16 + CwgtH1 CwgtD2 CwgtO16 Cdens + MwgtH1 MwgtD2 MwgtO16 Mdens + WwgtH1 WwgtD2 WwgtO16 Wdens ; +REAL + BNZrNb BNZr BNZrII := + 0.00021 0.00031 0.00005962 ; +REAL + B10ZrNb B10Zr B10ZrII + B11ZrNb B11Zr B11ZrII := + 0.199 BNZrNb * 0.199 BNZr * 0.199 BNZrII * + 0.801 BNZrNb * 0.801 BNZr * 0.801 BNZrII * ; +*---------- +* Library and isotope names +*---------- +STRING LibName + IsoU5 IsoU8 IsoU6 IsoPu9 IsoO16 + IsoH1 IsoD2 IsoB10 IsoB11 IsoBn + IsoZr IsoFe IsoFeA IsoCr IsoNi + IsoNb IsoC IsoSi IsoMn IsoCu + IsoMo IsoAl IsoHe4 IsoN IsoTi + IsoCd IsoCo59 := + "IAEA69" + "2235" "8238" "236" "6239" "6016" + "3001" "3002" "10" "11" "1011" + "91" "2056" "2056" "52" "58" + "93" "2012" "29" "55" "3063" + "96" "27" "4" "14" "48" + "2113" "1059" ; +*---------- +* End region information +* endzrv = Volume of Zr in endregion +* endcoolv = Volume of coolant in endregion +* endgasv = Volume of gas (void) in endregion +* endtotv = Total volume of endregion +* endzrf = Volume fraction of Zr in endregion +* endcoolf = Volume fraction of coolant in endregion +* endgasf = Volume fraction of gas (void) in endregion +*---------- +REAL endzrv endcoolv endgasv := 42.40 77.43 8.65 ; +REAL endtotv := endzrv endcoolv endgasv + + ; +REAL endzrf endcoolf endgasf := + endzrv endtotv / endcoolv endtotv / endgasv endtotv / ; +*---------- +* Guide tube perforation +*---------- +REAL PerfM1 PerfM2 PerfM3 := 0.3566 0.02 0.005 ; +REAL PerfZ1 PerfZ2 PerfZ3 := + 1.0 PerfM1 - 1.0 PerfM2 - 1.0 PerfM3 - ; +*---------- +* COMPUTE WGT% FOR FUEL +* COMPUTE DENSITY AND WGT% FOR COOLANT AND MODERATOR +*---------- +INFO: :: + ENR: <<Fenri>> WGT% + LIB: WIMSD4 FIL: <<LibName>> + ISO: 3 <<IsoU5>> <<IsoU8>> <<IsoO16>> + CALC WGT% UO2 >>FwgtU5<< >>FwgtU8<< >>FwgtO16<< ; +INFO: :: + TMP: <<Ctemp>> K + PUR: <<Cprty>> WGT% + CALC DENS WATER >>Cdens<< + LIB: WIMSD4 FIL: <<LibName>> + ISO: 3 <<IsoH1>> <<IsoD2>> <<IsoO16>> + CALC WGT% D2O >>CwgtH1<< >>CwgtD2<< >>CwgtO16<< ; +INFO: :: + TMP: <<Mtemp>> K + PUR: <<Mprty>> WGT% + CALC DENS WATER >>Mdens<< + LIB: WIMSD4 FIL: <<LibName>> + ISO: 3 <<IsoH1>> <<IsoD2>> <<IsoO16>> + CALC WGT% D2O >>MwgtH1<< >>MwgtD2<< >>MwgtO16<< ; +INFO: :: + TMP: <<Mtemp>> K + PUR: <<Wprty>> WGT% + CALC DENS WATER >>Wdens<< + LIB: WIMSD4 FIL: <<LibName>> + ISO: 3 <<IsoH1>> <<IsoD2>> <<IsoO16>> + CALC WGT% D2O >>WwgtH1<< >>WwgtD2<< >>WwgtO16<< ; +*---------- +* Generate Microscopic cross section library +*---------- +MicroLib := LIB: :: + EDIT <<PrtLev>> + NMIX 14 CTRA WIMS + DEPL LIB: WIMSD4 FIL: <<LibName>> + MIXS LIB: WIMSD4 FIL: <<LibName>> + MIX 1 <<Ctemp>> <<Cdens>> + H1 = <<IsoH1>> <<CwgtH1>> D2 = <<IsoD2>> <<CwgtD2>> + O16 = <<IsoO16>> <<CwgtO16>> + MIX 2 <<Ctemp>> 6.57 + ZR = <<IsoZr>> 97.5 NB = <<IsoNb>> 2.5 +* BN = <<IsoBn>> <<BNZrNb>> + B10 = <<IsoB10>> <<B10ZrNb>> B11 = <<IsoB11>> <<B11ZrNb>> + MIX 3 <<Mtemp>> 0.0014 + HE4 = <<IsoHe4>> 100.0 + MIX 4 <<Mtemp>> 6.44 + ZR = <<IsoZr>> 99.71 FE = <<IsoFe>> 0.16 + CR = <<IsoCr>> 0.11 NI = <<IsoNi>> 0.06 +* BN = <<IsoBn>> <<BNZr>> + B10 = <<IsoB10>> <<B10Zr>> B11 = <<IsoB11>> <<B11Zr>> + MIX 5 <<Mtemp>> <<Mdens>> + H1 = <<IsoH1>> <<MwgtH1>> D2 = <<IsoD2>> <<MwgtD2>> + O16 = <<IsoO16>> <<MwgtO16>> + MIX 6 <<Ftemp>> <<Fdens>> + O16 = <<IsoO16>> <<FwgtO16>> U5 = <<IsoU5>> <<FwgtU5>> 1 + U8 = <<IsoU8>> <<FwgtU8>> 1 U6 = <<IsoU6>> 0.0 1 + PU9 = <<IsoPu9>> 0.0 1 + MIX 7 COMB 6 100.0 + MIX 8 COMB 6 100.0 + MIX 9 COMB 6 100.0 + MIX 10 <<Ftemp>> <<Fdens>> + O16 = <<IsoO16>> <<FwgtO16>> U5 = <<IsoU5>> <<FwgtU5>> 2 + U8 = <<IsoU8>> <<FwgtU8>> 2 U6 = <<IsoU6>> 0.0 2 + PU9 = <<IsoPu9>> 0.0 2 + MIX 11 COMB 10 100.0 + MIX 12 COMB 10 100.0 + MIX 13 COMB 10 100.0 + MIX 14 <<Ctemp>> 6.44 + ZR = <<IsoZr>> 99.71 FE = <<IsoFe>> 0.16 + CR = <<IsoCr>> 0.11 NI = <<IsoNi>> 0.06 +* BN = <<IsoBn>> <<BNZr>> + B10 = <<IsoB10>> <<B10Zr>> B11 = <<IsoB11>> <<B11Zr>> + ; +QUIT . diff --git a/Dragon/data/proc_LZC/MicWULZC.c2m b/Dragon/data/proc_LZC/MicWULZC.c2m new file mode 100644 index 0000000..0dcda38 --- /dev/null +++ b/Dragon/data/proc_LZC/MicWULZC.c2m @@ -0,0 +1,395 @@ +*DECK MicWULZC.c2m +*---------- +* Name : MicWULZC.c2m +* Type : DRAGON procedure +* Use : Add LZC mixtures to Gentilly-2 fuel microlib +* with cross sections from the WLUP based 69 +* groups (LibName = IAEA69). +* (cluster model) +* Author : G. Marleau +* Date : 2006/03/01 +* +* Procedure description: +* MicroLib := MicWULZC MicroLib :: +* <<PrtLev>> <<Position>> <<Type>> ; +* +* MicroLib : XSM_FILE containing the updated microlib where: +* Mix 15 = LZC zone to fill or empty +* Mix 16 = Guide tube +* Mix 17 = mixture in bubbler for scavenger 1 +* Mix 18 = mixture around bubbler for scavenger 1 +* Mix 19 = mixture in scavenger 1 +* Mix 20 = mixture around scavenger 1 +* Mix 21 = mixture in bubbler for scavenger 2 +* Mix 22 = mixture around bubbler for scavenger 2 +* Mix 23 = mixture in scavenger 2 +* Mix 24 = mixture around scavenger 2 +* Mix 25 = mixture in bubbler for scavenger 3 +* Mix 26 = mixture around bubbler for scavenger 3 +* Mix 27 = mixture in scavenger 3 +* Mix 28 = mixture around scavenger 3 +* Mix 29 = mixture in balance for feeder 1 +* Mix 30 = mixture around balance for feeder 1 +* Mix 31 = mixture in feeder 1 +* Mix 32 = mixture around feeder 1 +* Mix 33 = mixture in balance for feeder 2 +* Mix 34 = mixture around balance for feeder 2 +* Mix 35 = mixture in feeder 2 +* Mix 36 = mixture around feeder 2 +* PrtLev : Print Level +* 0 -> Minimum print +* 1 -> Standard print +* 10 -> Debug +* Position : LZC position +* NO -> LZC and guide tubes absent. +* OUT -> LZC empty. +* IN -> LZC filled. +* Type : LZC type +* 10 -> 1 scavenger tube 0 feeder tube (default) +* 21 -> 2 scavenger tubes 1 feeder tube +* 32 -> 3 scavenger tubes 2 feeder tubes +* +* +* +*---------- +* Define procedure parameters and options and read +* the procedure input data +*---- +PARAMETER MicroLib :: EDIT 0 ::: XSM_FILE MicroLib ; ; +INTEGER PrtLev ; +STRING Position Type ; +:: >>PrtLev<< >>Position<< >>Type<< ; +*---------- +* Define modules +*---------- +MODULE INFO: LIB: DELETE: ; +*---------- +* Library and isotope names +*---------- +REAL Mtemp Mprty := 345.66 99.92 ; +REAL MwgtH1 MwgtD2 MwgtO16 Mdens ; +REAL Wtemp Wprty := Mtemp 0.001 ; +REAL WwgtH1 WwgtD2 WwgtO16 Wdens ; +REAL + BNZrNb BNZr BNZrII := + 0.00021 0.00031 0.00005962 ; +REAL + B10ZrNb B10Zr B10ZrII + B11ZrNb B11Zr B11ZrII := + 0.199 BNZrNb * 0.199 BNZr * 0.199 BNZrII * + 0.801 BNZrNb * 0.801 BNZr * 0.801 BNZrII * ; +*---------- +* Library and isotope names +*---------- +STRING LibName + IsoU5 IsoU8 IsoU6 IsoPu9 IsoO16 + IsoH1 IsoD2 IsoB10 IsoB11 IsoBn + IsoZr IsoFe IsoFeA IsoCr IsoNi + IsoNb IsoC IsoSi IsoMn IsoCu + IsoMo IsoAl IsoHe4 IsoN IsoTi + IsoCd IsoCo59 := + "IAEA69" + "2235" "8238" "236" "6239" "6016" + "3001" "3002" "10" "11" "1011" + "91" "2056" "2056" "52" "58" + "93" "2012" "29" "55" "3063" + "96" "27" "4" "14" "48" + "2113" "1059" ; +INFO: :: + TMP: <<Mtemp>> K + PUR: <<Mprty>> WGT% + CALC DENS WATER >>Mdens<< + LIB: WIMSD4 FIL: <<LibName>> + ISO: 3 <<IsoH1>> <<IsoD2>> <<IsoO16>> + CALC WGT% D2O >>MwgtH1<< >>MwgtD2<< >>MwgtO16<< ; +INFO: :: + TMP: <<Wtemp>> K + PUR: <<Wprty>> WGT% + CALC DENS WATER >>Wdens<< + LIB: WIMSD4 FIL: <<LibName>> + ISO: 3 <<IsoH1>> <<IsoD2>> <<IsoO16>> + CALC WGT% D2O >>WwgtH1<< >>WwgtD2<< >>WwgtO16<< ; +*---------- +* Complete Microscopic cross section library +*---------- +IF Position "IN" = THEN + IF Type "32" = THEN + MicroLib := LIB: MicroLib :: + EDIT <<PrtLev>> + NMIX 36 CTRA WIMS + MIXS LIB: WIMSD4 FIL: <<LibName>> + MIX 15 <<Wtemp>> <<Wdens>> + H1 = <<IsoH1>> <<WwgtH1>> D2 = <<IsoD2>> <<WwgtD2>> + O16 = <<IsoO16>> <<WwgtO16>> + MIX 16 <<Mtemp>> 6.4003 + ZR = <<IsoZr>> 98.2083 FE = <<IsoFe>> 0.135 + CR = <<IsoCr>> 0.10 NI = <<IsoNi>> 0.055 +* BN = <<IsoBn>> <<BNZrII>> + B10 = <<IsoB10>> <<B10ZrII>> B11 = <<IsoB11>> <<B10ZrII>> +*Scavanger 1 + MIX 17 <<Mtemp>> 0.0014 + HE4 = <<IsoHe4>> 100.0 + MIX 18 COMB 16 100.0 + MIX 19 COMB 15 100.0 + MIX 20 COMB 16 100.0 +*Scavanger 2 + MIX 21 COMB 17 100.0 + MIX 22 COMB 18 100.0 + MIX 23 COMB 19 100.0 + MIX 24 COMB 20 100.0 +*Scavanger 3 + MIX 25 COMB 17 100.0 + MIX 26 COMB 18 100.0 + MIX 27 COMB 19 100.0 + MIX 28 COMB 20 100.0 +*Feeder 1 + MIX 29 COMB 17 100.0 + MIX 30 COMB 18 100.0 + MIX 31 COMB 19 100.0 + MIX 32 COMB 20 100.0 +*Feeder 2 + MIX 33 COMB 17 100.0 + MIX 34 COMB 18 100.0 + MIX 35 COMB 19 100.0 + MIX 36 COMB 20 100.0 + ; + ELSEIF Type "21" = THEN + MicroLib := LIB: MicroLib :: + EDIT <<PrtLev>> + NMIX 36 CTRA WIMS + MIXS LIB: WIMSD4 FIL: <<LibName>> + MIX 15 <<Wtemp>> <<Wdens>> + H1 = <<IsoH1>> <<WwgtH1>> D2 = <<IsoD2>> <<WwgtD2>> + O16 = <<IsoO16>> <<WwgtO16>> + MIX 16 <<Mtemp>> 6.4003 + ZR = <<IsoZr>> 98.2083 FE = <<IsoFe>> 0.135 + CR = <<IsoCr>> 0.10 NI = <<IsoNi>> 0.055 +* BN = <<IsoBn>> <<BNZrII>> + B10 = <<IsoB10>> <<B10ZrII>> B11 = <<IsoB11>> <<B10ZrII>> +*Scavanger 1 + MIX 17 <<Mtemp>> 0.0014 + HE4 = <<IsoHe4>> 100.0 + MIX 18 COMB 16 100.0 + MIX 19 COMB 15 100.0 + MIX 20 COMB 16 100.0 +*Scavanger 2 + MIX 21 COMB 17 100.0 + MIX 22 COMB 18 100.0 + MIX 23 COMB 19 100.0 + MIX 24 COMB 20 100.0 +*Scavanger 3 + MIX 25 COMB 15 100.0 + MIX 26 COMB 15 100.0 + MIX 27 COMB 15 100.0 + MIX 28 COMB 15 100.0 +*Feeder 1 + MIX 29 COMB 17 100.0 + MIX 30 COMB 18 100.0 + MIX 31 COMB 19 100.0 + MIX 32 COMB 20 100.0 +*Feeder 2 + MIX 33 COMB 15 100.0 + MIX 34 COMB 15 100.0 + MIX 35 COMB 15 100.0 + MIX 36 COMB 15 100.0 + ; + ELSE + MicroLib := LIB: MicroLib :: + EDIT <<PrtLev>> + NMIX 36 CTRA WIMS + MIXS LIB: WIMSD4 FIL: <<LibName>> + MIX 15 <<Wtemp>> <<Wdens>> + H1 = <<IsoH1>> <<WwgtH1>> D2 = <<IsoD2>> <<WwgtD2>> + O16 = <<IsoO16>> <<WwgtO16>> + MIX 16 <<Mtemp>> 6.4003 + ZR = <<IsoZr>> 98.2083 FE = <<IsoFe>> 0.135 + CR = <<IsoCr>> 0.10 NI = <<IsoNi>> 0.055 +* BN = <<IsoBn>> <<BNZrII>> + B10 = <<IsoB10>> <<B10ZrII>> B11 = <<IsoB11>> <<B10ZrII>> +*Scavanger 1 + MIX 17 <<Mtemp>> 0.0014 + HE4 = <<IsoHe4>> 100.0 + MIX 18 COMB 16 100.0 + MIX 19 COMB 15 100.0 + MIX 20 COMB 16 100.0 +*Scavanger 2 + MIX 21 COMB 15 100.0 + MIX 22 COMB 15 100.0 + MIX 23 COMB 15 100.0 + MIX 24 COMB 15 100.0 +*Scavanger 3 + MIX 25 COMB 15 100.0 + MIX 26 COMB 15 100.0 + MIX 27 COMB 15 100.0 + MIX 28 COMB 15 100.0 +*Feeder 1 + MIX 29 COMB 15 100.0 + MIX 30 COMB 15 100.0 + MIX 31 COMB 15 100.0 + MIX 32 COMB 15 100.0 +*Feeder 2 + MIX 33 COMB 15 100.0 + MIX 34 COMB 15 100.0 + MIX 35 COMB 15 100.0 + MIX 36 COMB 15 100.0 + ; + ENDIF ; +ELSEIF Position "OUT" = THEN + IF Type "32" = THEN + MicroLib := LIB: MicroLib :: + EDIT <<PrtLev>> + NMIX 36 CTRA WIMS + MIXS LIB: WIMSD4 FIL: <<LibName>> + MIX 15 <<Mtemp>> 0.0014 + HE4 = <<IsoHe4>> 100.0 + MIX 16 <<Mtemp>> 6.4003 + ZR = <<IsoZr>> 98.2083 FE = <<IsoFe>> 0.135 + CR = <<IsoCr>> 0.10 NI = <<IsoNi>> 0.055 +* BN = <<IsoBn>> <<BNZrII>> + B10 = <<IsoB10>> <<B10ZrII>> B11 = <<IsoB11>> <<B10ZrII>> +*Scavanger 1 + MIX 17 COMB 15 100.0 + MIX 18 COMB 16 100.0 + MIX 19 <<Wtemp>> <<Wdens>> + H1 = <<IsoH1>> <<WwgtH1>> D2 = <<IsoD2>> <<WwgtD2>> + O16 = <<IsoO16>> <<WwgtO16>> + MIX 20 COMB 16 100.0 +*Scavanger 2 + MIX 21 COMB 17 100.0 + MIX 22 COMB 18 100.0 + MIX 23 COMB 19 100.0 + MIX 24 COMB 20 100.0 +*Scavanger 3 + MIX 25 COMB 17 100.0 + MIX 26 COMB 18 100.0 + MIX 27 COMB 19 100.0 + MIX 28 COMB 20 100.0 +*Feeder 1 + MIX 29 COMB 17 100.0 + MIX 30 COMB 18 100.0 + MIX 31 COMB 19 100.0 + MIX 32 COMB 20 100.0 +*Feeder 2 + MIX 33 COMB 17 100.0 + MIX 34 COMB 18 100.0 + MIX 35 COMB 19 100.0 + MIX 36 COMB 20 100.0 + ; + ELSEIF Type "21" = THEN + MicroLib := LIB: MicroLib :: + EDIT <<PrtLev>> + NMIX 36 CTRA WIMS + MIXS LIB: WIMSD4 FIL: <<LibName>> + MIX 15 <<Mtemp>> 0.0014 + HE4 = <<IsoHe4>> 100.0 + MIX 16 <<Mtemp>> 6.4003 + ZR = <<IsoZr>> 98.2083 FE = <<IsoFe>> 0.135 + CR = <<IsoCr>> 0.10 NI = <<IsoNi>> 0.055 +* BN = <<IsoBn>> <<BNZrII>> + B10 = <<IsoB10>> <<B10ZrII>> B11 = <<IsoB11>> <<B10ZrII>> +*Scavanger 1 + MIX 17 COMB 15 100.0 + MIX 18 COMB 16 100.0 + MIX 19 <<Wtemp>> <<Wdens>> + H1 = <<IsoH1>> <<WwgtH1>> D2 = <<IsoD2>> <<WwgtD2>> + O16 = <<IsoO16>> <<WwgtO16>> + MIX 20 COMB 16 100.0 +*Scavanger 2 + MIX 21 COMB 17 100.0 + MIX 22 COMB 18 100.0 + MIX 23 COMB 19 100.0 + MIX 24 COMB 20 100.0 +*Scavanger 3 + MIX 25 COMB 15 100.0 + MIX 26 COMB 15 100.0 + MIX 27 COMB 15 100.0 + MIX 28 COMB 15 100.0 +*Feeder 1 + MIX 29 COMB 17 100.0 + MIX 30 COMB 18 100.0 + MIX 31 COMB 19 100.0 + MIX 32 COMB 20 100.0 +*Feeder 2 + MIX 33 COMB 15 100.0 + MIX 34 COMB 15 100.0 + MIX 35 COMB 15 100.0 + MIX 36 COMB 15 100.0 + ; + ELSE + MicroLib := LIB: MicroLib :: + EDIT <<PrtLev>> + NMIX 36 CTRA WIMS + MIXS LIB: WIMSD4 FIL: <<LibName>> + MIX 15 <<Mtemp>> 0.0014 + HE4 = <<IsoHe4>> 100.0 + MIX 16 <<Mtemp>> 6.4003 + ZR = <<IsoZr>> 98.2083 FE = <<IsoFe>> 0.135 + CR = <<IsoCr>> 0.10 NI = <<IsoNi>> 0.055 +* BN = <<IsoBn>> <<BNZrII>> + B10 = <<IsoB10>> <<B10ZrII>> B11 = <<IsoB11>> <<B10ZrII>> +*Scavanger 1 + MIX 17 COMB 15 100.0 + MIX 18 COMB 16 100.0 + MIX 19 <<Wtemp>> <<Wdens>> + H1 = <<IsoH1>> <<WwgtH1>> D2 = <<IsoD2>> <<WwgtD2>> + O16 = <<IsoO16>> <<WwgtO16>> + MIX 20 COMB 16 100.0 +*Scavanger 2 + MIX 21 COMB 15 100.0 + MIX 22 COMB 15 100.0 + MIX 23 COMB 15 100.0 + MIX 24 COMB 15 100.0 +*Scavanger 3 + MIX 25 COMB 15 100.0 + MIX 26 COMB 15 100.0 + MIX 27 COMB 15 100.0 + MIX 28 COMB 15 100.0 +*Feeder 1 + MIX 29 COMB 15 100.0 + MIX 30 COMB 15 100.0 + MIX 31 COMB 15 100.0 + MIX 32 COMB 15 100.0 +*Feeder 2 + MIX 33 COMB 15 100.0 + MIX 34 COMB 15 100.0 + MIX 35 COMB 15 100.0 + MIX 36 COMB 15 100.0 + ; + ENDIF ; +ELSE + MicroLib := LIB: MicroLib :: + EDIT <<PrtLev>> + NMIX 36 CTRA WIMS + MIXS LIB: WIMSD4 FIL: <<LibName>> + MIX 15 <<Mtemp>> <<Mdens>> + H1 = <<IsoH1>> <<MwgtH1>> D2 = <<IsoD2>> <<MwgtD2>> + O16 = <<IsoO16>> <<MwgtO16>> + MIX 16 COMB 15 100.0 +*Scavanger 1 + MIX 17 COMB 15 100.0 + MIX 18 COMB 15 100.0 + MIX 19 COMB 15 100.0 + MIX 20 COMB 15 100.0 +*Scavanger 2 + MIX 21 COMB 15 100.0 + MIX 22 COMB 15 100.0 + MIX 23 COMB 15 100.0 + MIX 24 COMB 15 100.0 +*Scavanger 3 + MIX 25 COMB 15 100.0 + MIX 26 COMB 15 100.0 + MIX 27 COMB 15 100.0 + MIX 28 COMB 15 100.0 +*Feeder 1 + MIX 29 COMB 15 100.0 + MIX 30 COMB 15 100.0 + MIX 31 COMB 15 100.0 + MIX 32 COMB 15 100.0 +*Feeder 2 + MIX 33 COMB 15 100.0 + MIX 34 COMB 15 100.0 + MIX 35 COMB 15 100.0 + MIX 36 COMB 15 100.0 + ; +ENDIF ; +QUIT . |