@@ -30,6 +30,7 @@ void wrapContinuousEnergyTable( python::module& module, python::module& ) {
   using ANDP = njoy::ACEtk::block::ANDP;
   using DLWP = njoy::ACEtk::block::DLWP;
   using YP = njoy::ACEtk::block::YP;
+  using UNR = njoy::ACEtk::block::UNR;
 
   // wrap views created by this table
 
@@ -52,14 +53,15 @@ void wrapContinuousEnergyTable( python::module& module, python::module& ) {
                   SIG, AND, DLW, std::optional< GPD >,
                   std::optional< MTRP >, std::optional< SIGP >,
                   std::optional< ANDP >, std::optional< DLWP >,
-                  std::optional< YP > >(),
+                  std::optional< YP >, std::optional< UNR > >(),
     python::arg( "z" ), python::arg( "a" ),
     python::arg( "header" ), python::arg( "esz" ),
     python::arg( "nu" ) = std::nullopt,
     python::arg( "mtr" ), python::arg( "lqr" ), python::arg( "tyr" ),
     python::arg( "sig" ), python::arg( "ang" ), python::arg( "dlw" ),
     python::arg( "gpd" ), python::arg( "mtrp" ), python::arg( "sigp" ),
     python::arg( "andp" ), python::arg( "dlwp" ), python::arg( "yp" ),
+    python::arg( "unr" ),
     "Initialise the table\n\n"
     "Arguments:\n"
     "    self      the table\n"
@@ -79,7 +81,8 @@ void wrapContinuousEnergyTable( python::module& module, python::module& ) {
     "    sigp      the photon production cross section data block\n"
     "    angp      the photon production angular distribution data block\n"
     "    dlwp      the photon production energy distribution data block\n"
-    "    yp        the photon multiplicity reaction number block"
+    "    yp        the photon multiplicity reaction number block\n"
+    "    unr       the unresolved resonance probability table block"
   )
   .def_property_readonly(
 
 
@@ -52,6 +52,7 @@ ContinuousEnergyTable& operator=( ContinuousEnergyTable&& base ) {
  *  @param[in] angp      the photon production angular distribution data block
  *  @param[in] dlwp      the photon production energy distribution data block
  *  @param[in] yp        the photon multiplicity reaction number block
+ *  @param[in] unr       the unresolved resonance probability table block
  */
 ContinuousEnergyTable( unsigned int z, unsigned int a, Header header,
                        block::ESZ esz, std::optional< block::NU > nu,
@@ -62,7 +63,8 @@ ContinuousEnergyTable( unsigned int z, unsigned int a, Header header,
                        std::optional< block::SIGP > sigp,
                        std::optional< block::ANDP > andp,
                        std::optional< block::DLWP > dlwp,
-                       std::optional< block::YP > yp ) :
+                       std::optional< block::YP > yp,
+                       std::optional< block::UNR > unr ) :
   ContinuousEnergyTable(
       Table( std::move( header ),
              generateData( z, a, std::move( esz ), std::move( nu ),
@@ -71,4 +73,5 @@ ContinuousEnergyTable( unsigned int z, unsigned int a, Header header,
                            std::move( ang ), std::move( dlw ),
                            std::move( gpd ), std::move( mtrp ),
                            std::move( sigp ), std::move( andp ),
-                           std::move( dlwp ), std::move( yp ) ) ) ) {}
+                           std::move( dlwp ), std::move( yp ),
+                           std::move( unr ) ) ) ) {}
@@ -7,7 +7,8 @@ Data generateData( unsigned int z, unsigned int a,
                    std::optional< block::SIGP >&& sigp,
                    std::optional< block::ANDP >&& andp,
                    std::optional< block::DLWP >&& dlwp,
-                   std::optional< block::YP >&& yp ) {
+                   std::optional< block::YP >&& yp,
+                   std::optional< block::UNR >&& unr ) {
 
   std::array< int32_t, 16 > iz;
   std::array< double, 16 > aw;
@@ -69,7 +70,7 @@ Data generateData( unsigned int z, unsigned int a,
   if ( nu ) {
 
     jxs[1] = xss.size() + 1;
-    xss.insert( xss.end(), nu.value().begin(), nu.value().end() );
+    xss.insert( xss.end(), nu->begin(), nu->end() );
   }
   jxs[2] = xss.size() + 1;
   xss.insert( xss.end(), mtr.begin(), mtr.end() );
@@ -101,12 +102,19 @@ Data generateData( unsigned int z, unsigned int a,
     jxs[17] = xss.size() + 1;
     jxs[18] = xss.size() + ntrp + 1;
     xss.insert( xss.end(), dlwp->begin(), dlwp->end() );
+    jxs[19] = xss.size() + 1;
+    xss.insert( xss.end(), yp->begin(), yp->end() );
   }
   if ( mtr.hasReactionNumber( 18 ) ) {
 
     jxs[20] = jxs[8] + sig.crossSectionLocator( mtr.index( 18 ) ) - 1;
   }
   jxs[21] = xss.size() + 1;
+  if ( unr ) {
+
+    jxs[22] = xss.size() + 1;
+    xss.insert( xss.end(), unr->begin(), unr->end() );
+  }
 
   // set the nxs values for the continuous energy table
   // NXS(1) = length