4.4. VERA_3A benchmark

!===============================================================================
!=== material
!===============================================================================
material:

  library: micro
    mg_lib  = default
    ce_lib  = default
  /library

  mat: fuel_3.1
    id = 922340; density = 6.11874e-6
    id = 922350; density = 7.18132e-4
    id = 922360; density = 3.29861E-6
    id = 922380; density = 2.21546e-2
    id = 80160 ; density = 4.57642e-2
  /mat

  mat: moderator_3a
    id = 80160;  density = 2.48112E-02
    id = 10010;  density = 4.96224E-02
    id = 50100;  density = 1.07070E-05
    id = 50110;  density = 4.30971E-05
  /mat

  mat: helium
    id = 20040; density = 2.68714e-05
  /mat

  mat: zircaloy-4
    id = 240500; density = 3.30121E-06
    id = 240520; density = 6.36606E-05
    id = 240530; density = 7.21860E-06
    id = 240540; density = 1.79686E-06
    id = 260540; density = 8.68307E-06
    id = 260560; density = 1.36306E-04
    id = 260570; density = 3.14789E-06
    id = 260580; density = 4.18926E-07
    id = 400900; density = 2.18865E-02
    id = 400910; density = 4.77292E-03
    id = 400920; density = 7.29551E-03
    id = 400940; density = 7.39335E-03
    id = 400960; density = 1.19110E-03
    id = 501120; density = 4.68066E-06
    id = 501140; density = 3.18478E-06
    id = 501150; density = 1.64064E-06
    id = 501160; density = 7.01616E-05
    id = 501170; density = 3.70592E-05
    id = 501180; density = 1.16872E-04
    id = 501190; density = 4.14504E-05
    id = 501200; density = 1.57212E-04
    id = 501220; density = 2.23417E-05
    id = 501240; density = 2.79392E-05
    id = 721740; density = 3.54138E-09
    id = 721760; density = 1.16423E-07
    id = 721770; density = 4.11686E-07
    id = 721780; density = 6.03806E-07
    id = 721790; density = 3.01460E-07
    id = 721800; density = 7.76449E-07
  /mat

  mat: inconel
    id = 140280; density = 4.04885E-03
    id = 140290; density = 2.05685E-04
    id = 140300; density = 1.35748E-04
    id = 220460; density = 2.12518E-04
    id = 220470; density = 1.91652E-04
    id = 220480; density = 1.89901E-03
    id = 220490; density = 1.39360E-04
    id = 220500; density = 1.33435E-04
    id = 240500; density = 6.18222E-04
    id = 240520; density = 1.19218E-02
    id = 240530; density = 1.35184E-03
    id = 240540; density = 3.36501E-04
    id = 260540; density = 3.61353E-04
    id = 260560; density = 5.67247E-03
    id = 260570; density = 1.31002E-04
    id = 260580; density = 1.74340E-05
    id = 280580; density = 4.17608E-02
    id = 280600; density = 1.60862E-02
    id = 280610; density = 6.99255E-04
    id = 280620; density = 2.22953E-03
    id = 280640; density = 5.67796E-04
  /mat

  mat: top_nozzle
    id =  10010; density = 4.01211E-02
    id =  50100; density = 8.65222E-06
    id =  50110; density = 3.48263E-05
    id =  60000; density = 6.14459E-05
    id =  80160; density = 2.00606E-02
    id = 140280; density = 3.02920E-04
    id = 140290; density = 1.53886E-05
    id = 140300; density = 1.01561E-05
    id = 150310; density = 1.34026E-05
    id = 240500; density = 1.46468E-04
    id = 240520; density = 2.82449E-03
    id = 240530; density = 3.20275E-04
    id = 240540; density = 7.97232E-05
    id = 250550; density = 3.35836E-04
    id = 260540; density = 6.60188E-04
    id = 260560; density = 1.03635E-02
    id = 260570; density = 2.39339E-04
    id = 260580; density = 3.18517E-05
    id = 280580; density = 1.01650E-03
    id = 280600; density = 3.91552E-04
    id = 280610; density = 1.70205E-05
    id = 280620; density = 5.42688E-05
    id = 280640; density = 1.38207E-05
  /mat

  mat: bottom_nozzle
    id =  10010; density = 3.57661E-02
    id =  50100; density = 7.70514E-06
    id =  50110; density = 3.10142E-05
    id =  60000; density = 8.96008E-05
    id =  80160; density = 1.78830E-02
    id = 140280; density = 4.41720E-04
    id = 140290; density = 2.24397E-05
    id = 140300; density = 1.48097E-05
    id = 150310; density = 1.95438E-05
    id = 240500; density = 2.13581E-04
    id = 240520; density = 4.11869E-03
    id = 240530; density = 4.67027E-04
    id = 240540; density = 1.16253E-04
    id = 250550; density = 4.89719E-04
    id = 260540; density = 9.62690E-04
    id = 260560; density = 1.51122E-02
    id = 260570; density = 3.49006E-04
    id = 260580; density = 4.64463E-05
    id = 280580; density = 1.48226E-03
    id = 280600; density = 5.70964E-04
    id = 280610; density = 2.48194E-05
    id = 280620; density = 7.91351E-05
    id = 280640; density = 2.01534E-05
  /mat

  mat: plates
    id =  10010; density = 2.48112E-02
    id =  50100; density = 5.33040E-06
    id =  50110; density = 2.14555E-05
    id =  60000; density = 1.60447E-04
    id =  80160; density = 1.24056E-02
    id = 140280; density = 7.90985E-04
    id = 140290; density = 4.01826E-05
    id = 140300; density = 2.65197E-05
    id = 150310; density = 3.49969E-05
    id = 240500; density = 3.82458E-04
    id = 240520; density = 7.37532E-03
    id = 240530; density = 8.36302E-04
    id = 240540; density = 2.08173E-04
    id = 250550; density = 8.76936E-04
    id = 260540; density = 1.72388E-03
    id = 260560; density = 2.70613E-02
    id = 260570; density = 6.24963E-04
    id = 260580; density = 8.31710E-05
    id = 280580; density = 2.65427E-03
    id = 280600; density = 1.02242E-03
    id = 280610; density = 4.44439E-05
    id = 280620; density = 1.41707E-04
    id = 280640; density = 3.60885E-05
  /mat
/material

Geometry:
!===============================================================================
!=== Pin_geometry
!===============================================================================
Pin_geometry:
   pin_geom_type: pg_fuel
     size_x           =  1.26
     size_y           =  1.26
     radii            =  0.4096 0.418 0.475
     r_mesh           =  3 1 1 1
     azi_mesh         =  8 8 8 8
   /pin_geom_type

   pin_geom_type: pg_gt
     size_x           =  1.26
     size_y           =  1.26
     radii            =  0.561 0.602
     r_mesh           =  3 1 1
     azi_mesh         =  8 8 8
   /pin_geom_type

   pin_geom_type: rectangle
     size_x           =  1.26
     size_y           =  1.26
     x_mesh           =  8
     y_mesh           =  8
   /pin_geom_type

   pin_geom_type: pg_plate_nozzle
     size_x           =  1.26
     size_y           =  1.26
     radii            =  0.602
     r_mesh           =  3 1
     azi_mesh         =  8 8
   /pin_geom_type

   pin_geom_type:  pg_plug
     size_x           =  1.26
     size_y           =  1.26
     radii            =  0.475
     r_mesh           =  3 1
     azi_mesh         =  8 8
   /pin_geom_type

   pin_geom_type: plenum
     size_x           =  1.26
     size_y           =  1.26
     radii            =  0.418 0.475
     r_mesh           =  3 1 1
     azi_mesh         =  8 8 8
   /pin_geom_type
/ Pin_geometry

!===============================================================================
!=== Pin
!===============================================================================
Pin:
   pin_type: pin_3a
     name             =  pin
     int_ID           =  1
     geom             =  pg_fuel
     mat              =  fuel_3.1  helium  zircaloy-4  moderator_3a
     temperature      =  600   600  600   600
   /pin_type

   pin_type: pin_gt
     name             =  pin_gt
     int_ID           =  2
     geom             =  pg_gt
     mat              =  moderator_3a  zircaloy-4  moderator_3a
     temperature      =  600   600  600
   /pin_type

   pin_type: rectangle
     name             =  rectangle
     int_ID           =  13
     geom             =  rectangle
     mat              =  moderator_3a
     temperature      =  600
   /pin_type

   pin_type: pin_refl
     name             =  pin_refl
     int_ID           =  3
     geom             =  pg_fuel
     mat              =  moderator_3a  moderator_3a  moderator_3a  moderator_3a
     temperature      =  600   600  600   600
   /pin_type

   pin_type: pin_plate_cylinder
     name             =  pin_plate_cylinder
     int_ID           =  4
     geom             =  pg_plate_nozzle
     mat              =  plates plates
     temperature      =  600 600
   /pin_type

   pin_type: pin_plate_rectangle
     name             =  pin_plate_rectangle
     int_ID           =  14
     geom             =  rectangle
     mat              =  plates
     temperature      =  600
   /pin_type

   pin_type: pin_plate2
     name             =  pin_plate2
     int_ID           =  5
     geom             =  rectangle
     mat              =  moderator_3a
     temperature      =  600
   /pin_type

   pin_type: bottom_nozzle_rect
     name             =  bottom_nozzle_rect
     int_ID           =  6
     geom             =  rectangle
     mat              =  bottom_nozzle
     temperature      =  600
   /pin_type

   pin_type: bottom_nozzle_cylinder
     name             =  bottom_nozzle_cylinder
     int_ID           =  16
     geom             =  pg_plate_nozzle
     mat              =  bottom_nozzle bottom_nozzle
     temperature      =  600 600
   /pin_type

   pin_type: pin_nozzle2
     name             =  pin_nozzle2
     int_ID           =  7
     geom             =  rectangle
     mat              =  moderator_3a
     temperature      =  600
   /pin_type

   pin_type: pin_plug
     name             =  pin_plug
     int_ID           =  8
     geom             =  pg_plug
     mat              =  zircaloy-4 moderator_3a
     temperature      =  600 600
   /pin_type

   pin_type: pin_plenum
     name             =  pin_plenum
     int_ID           =  9
     geom             =  plenum
     mat              =  helium  zircaloy-4  moderator_3a
     temperature      =  600   600  600
   /pin_type

    pin_type: top_nozzle
     name             =  top_nozzle
     int_ID           =  10
     geom             =  pg_plate_nozzle
     mat              =  top_nozzle top_nozzle
     temperature      =  600 600
   /pin_type

/Pin

!===============================================================================
!=== Lattice_geometry
!===============================================================================
Lattice_geometry
   lattice_type: reflectors
    pinmap =
        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
        3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
    /lattice_type

   lattice_type: plates
    pinmap =
        4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4
        4 14 14 14 14 14 14 14 14 14 14 14 14 14 14 14  4
        4 14 14 14 14 14 14 14 14 14 14 14 14 14 14 14  4
        4 14 14 14 14 14 14 14 14 14 14 14 14 14 14 14  4
        4 14 14 14 14 14 14 14 14 14 14 14 14 14 14 14  4
        4 14 14 14 14 14 14 14 14 14 14 14 14 14 14 14  4
        4 14 14 14 14 14 14 14 14 14 14 14 14 14 14 14  4
        4 14 14 14 14 14 14 14 14 14 14 14 14 14 14 14  4
        4 14 14 14 14 14 14 14 14 14 14 14 14 14 14 14  4
        4 14 14 14 14 14 14 14 14 14 14 14 14 14 14 14  4
        4 14 14 14 14 14 14 14 14 14 14 14 14 14 14 14  4
        4 14 14 14 14 14 14 14 14 14 14 14 14 14 14 14  4
        4 14 14 14 14 14 14 14 14 14 14 14 14 14 14 14  4
        4 14 14 14 14 14 14 14 14 14 14 14 14 14 14 14  4
        4 14 14 14 14 14 14 14 14 14 14 14 14 14 14 14  4
        4 14 14 14 14 14 14 14 14 14 14 14 14 14 14 14  4
        4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4
    /lattice_type

   lattice_type: bottom_nozzle
    pinmap =
        16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16
        16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16
        16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16
        16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16
        16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16
        16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16
        16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16
        16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16
        16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16
        16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16
        16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16
        16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16
        16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16
        16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16
        16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16
        16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16
        16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16
    /lattice_type

   lattice_type: plugs
    pinmap =
        8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8
        8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8
        8 8 8 8 8 2 8 8 2 8 8 2 8 8 8 8 8
        8 8 8 2 8 8 8 8 8 8 8 8 8 2 8 8 8
        8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8
        8 8 2 8 8 2 8 8 2 8 8 2 8 8 2 8 8
        8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8
        8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8
        8 8 2 8 8 2 8 8 2 8 8 2 8 8 2 8 8
        8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8
        8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8
        8 8 2 8 8 2 8 8 2 8 8 2 8 8 2 8 8
        8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8
        8 8 8 2 8 8 8 8 8 8 8 8 8 2 8 8 8
        8 8 8 8 8 2 8 8 2 8 8 2 8 8 8 8 8
        8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8
        8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8
    /lattice_type

   lattice_type: fuel
    pinmap =
        1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
        1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
        1 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 1
        1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
        1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
        1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
        1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
        1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
        1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
        1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
        1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
        1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
        1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
        1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
        1 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 1
        1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
        1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
    /lattice_type

   lattice_type: plenum
    pinmap =
        9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9
        9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9
        9 9 9 9 9 2 9 9 2 9 9 2 9 9 9 9 9
        9 9 9 2 9 9 9 9 9 9 9 9 9 2 9 9 9
        9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9
        9 9 2 9 9 2 9 9 2 9 9 2 9 9 2 9 9
        9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9
        9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9
        9 9 2 9 9 2 9 9 2 9 9 2 9 9 2 9 9
        9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9
        9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9
        9 9 2 9 9 2 9 9 2 9 9 2 9 9 2 9 9
        9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9
        9 9 9 2 9 9 9 9 9 9 9 9 9 2 9 9 9
        9 9 9 9 9 2 9 9 2 9 9 2 9 9 9 9 9
        9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9
        9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9
    /lattice_type

   lattice_type: top_nozzle
    pinmap =
        10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
        10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
        10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
        10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
        10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
        10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
        10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
        10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
        10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
        10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
        10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
        10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
        10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
        10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
        10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
        10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
        10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
    /lattice_type

   lattice_type: gap
    pinmap =
        3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3
        3 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13  3
        3 13 13 13 13 2  13 13 2  13 13 2  13 13 13 13  3
        3 13 13 2  13 13 13 13 13 13 13 13 13 2  13 13  3
        3 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13  3
        3 13 2  13 13 2  13 13 2  13 13 2  13 13 2  13  3
        3 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13  3
        3 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13  3
        3 13 2  13 13 2  13 13 2  13 13 2  13 13 2  13  3
        3 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13  3
        3 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13  3
        3 13 2  13 13 2  13 13 2  13 13 2  13 13 2  13  3
        3 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13  3
        3 13 13 2  13 13 13 13 13 13 13 13 13 2  13 13  3
        3 13 13 13 13 2  13 13 2  13 13 2  13 13 13 13  3
        3 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13  3
        3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3
    /lattice_type

/Lattice_geometry

!===============================================================================
!=== Assembly geometry
!===============================================================================

Assembly_geometry
   assembly_type: A1
      grid = yes
      gap = 0.08
      gap_mat = moderator_3a
      axial:
        height = 10.000; lattice = reflectors
        height = 10.000; lattice = reflectors
        height = 10.000; lattice = reflectors
        height = 5.0000; lattice = plates
        height = 6.0530; lattice = bottom_nozzle
        height = 4.2280; lattice = gap
        height = 1.6700; lattice = plugs
        height = 12.0771; lattice = fuel
        height = 8.2112; lattice = fuel
        height = 8.2111; lattice = fuel
        height = 8.2112; lattice = fuel
        height = 8.2111; lattice = fuel
        height = 8.2112; lattice = fuel
        height = 8.2111; lattice = fuel
        height = 11.875; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 11.875; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 11.875; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 11.875; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 11.875; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 8.0650; lattice = fuel
        height = 11.7312; lattice = fuel
        height = 7.9212; lattice = fuel
        height = 7.9212; lattice = fuel
        height = 7.9212; lattice = fuel
        height = 7.9212; lattice = fuel
        height = 8.5560; lattice = plenum
        height = 7.4440; lattice = plenum
        height = 1.6700; lattice = plugs
        height = 2.1200; lattice = gap
        height = 8.8270; lattice = top_nozzle
        height = 7.6000; lattice = plates
        height = 10.000; lattice = reflectors
        height = 10.000; lattice = reflectors
        height = 10.000; lattice = reflectors
      /axial
   /assembly_type
/Assembly_geometry

!===============================================================================
!=== Grid geometry
!===============================================================================
Grid_geometry
    grid_type: inter
      height = 3.81
      type = nmv
      inner_strap:
        mat = zircaloy-4; thickness = 0.04902  !mass = 875g
      /inner_strap
      outer_strap:
        mat = moderator_3a;
        thickness = 0.02  !mass = 0
      /outer_strap
    /grid_type

    grid_type: end
      height = 3.866
      type = mv
      inner_strap:
        mat = inconel; thickness = 0.0449  !mass = 1017g
      /inner_strap
      outer_strap:
        mat = moderator_3a;
        thickness = 0.02  !mass = 0
      /outer_strap
    /grid_type

    grid_map:
      grid = end;   position = 48.884   ! refer to lower bottom reflector
      grid = inter; position = 110.2
      grid = inter; position = 162.4
      grid = inter; position = 214.6
      grid = inter; position = 266.8
      grid = inter; position = 319.0
      grid = inter; position = 371.2
      grid = end;   position = 423.2
    /grid_map

    configure = explicit
/Grid_geometry

Core_geometry
    symmetric = quarter
    nx_full = 1
    ny_full = 1
    bc_cond:
       south = 1.0, east = 1.0, north = 1.0
       west = 1.0, top = 0.0, bottom = 0.0
    /bc_cond
    layout = 1
    asy_map = A1
/Core_geometry

/Geometry

!===============================================================================
!=== options
!===============================================================================
options:
  solver
    axial_solver = sn
    radial_solver = moc
    ! 66 axial layers
    parallel = 1 12 1 1 1 !energy axial radial angle thread
    scatt_meth = p1
  /solver

  cmfd
    is_cmfd = yes
    relaxation = 1.0
    solver = sweep
  /cmfd

  ray
    spacing = 0.03
    quad_type = chebyshev_yamamoto
    azimuthal = 8
    polar = 3
    modular_type = quarter_assembly
  /ray

  iteration
    n_inner = 1
    n_outer = 200
    k_criteria = 1.0e-6
    fis_criteria = 1.0e-5
  /iteration

  resonance
    method = global_local
    interfer_method = pmrn
    with_energy_sph = yes
    n_fsr_fuel = 5
    n_fsr_mod = 1
    ssxs_store_scheme = shielded_xs
  /resonance

/options