Copilot comments

Author: dkachuma

inputFiles/singlePhaseWell/singlePhaseWell.ats

@@ -60,7 +60,7 @@ decks = [
         restartcheck_params=RestartcheckParameters(**restartcheck_params)),
     TestDeck(
         name="thermal_compressible_single_phase_wells_1d",
-        description='Thermal well test (1D, thermal, 1 wells)',
+        description='Thermal well test (1D, thermal, 1 well)',
         partitions=((1, 1, 1),),
         restart_step=10,
         check_step=10,

src/coreComponents/constitutive/thermalConductivity/SinglePhaseThermalConductivity.cpp

@@ -138,9 +138,13 @@ void SinglePhaseThermalConductivity::initializeState( localIndex const size ) co
   {
     for( localIndex q = 0; q < numQuad; ++q )
     {
-      referenceThermalConductivity[ei][q][0] = defaultThermalConductivityComponents[0];
-      referenceThermalConductivity[ei][q][1] = defaultThermalConductivityComponents[1];
-      referenceThermalConductivity[ei][q][2] = defaultThermalConductivityComponents[2];
+      for( integer dim=0; dim < 3; ++dim )
+      {
+        if( referenceThermalConductivity[ei][q][dim] < 0 )
+        {
+          referenceThermalConductivity[ei][q][dim] = defaultThermalConductivityComponents[dim];
+        }
+      }
     }
   } );
 }

src/coreComponents/constitutive/thermalConductivity/SinglePhaseThermalConductivityBase.hpp

@@ -55,7 +55,9 @@ class SinglePhaseThermalConductivityBaseUpdate
 
   /**
    * @brief Constructor for the class performing the thermal conductivity updates
+   * @param referenceThermalConductivity the array of cell-wise reference conductivities in the subregion
    * @param effectiveConductivity the array of cell-wise effective conductivities in the subregion
+   * @param dEffectiveConductivity_dT the array of cell-wise derivatives of the effective conductivities w.r.t. temperature in the subregion
    */
   SinglePhaseThermalConductivityBaseUpdate( arrayView3d< real64 > const & referenceThermalConductivity,
                                             arrayView3d< real64 > const & effectiveConductivity,

src/coreComponents/physicsSolvers/surfaceGeneration/EmbeddedSurfaceGenerator.cpp

@@ -288,7 +288,7 @@ real64 EmbeddedSurfaceGenerator::solverStep( real64 const & GEOS_UNUSED_PARAM( t
     string const thermalCondModelName = getConstitutiveName< SinglePhaseThermalConductivityBase >( fractureSubRegion );
     if( !thermalCondModelName.empty() )
     {
-      // if a thermal conductivity model exists we need to set the intial value to something meaningful
+      // if a thermal conductivity model exists we need to set the initial value to something meaningful
       SinglePhaseThermalConductivityBase & thermalCondModel = getConstitutiveModel< SinglePhaseThermalConductivityBase >( fractureSubRegion, thermalCondModelName );
       thermalCondModel.initializeState();
     }