Initial implementation of spray capability for Manifold EOS

Source/Eos/Manifold.H

@@ -185,7 +185,7 @@ struct Manifold
   AMREX_FORCE_INLINE
   static void T2Cpi(const amrex::Real /*T*/, amrex::Real* /*Cpi[NUM_SPECIES]*/)
   {
-    // amrex::Abort("\nT2Cpi is not yet implemented for Maniold model");
+    amrex::Abort("\nT2Cpi is not yet implemented for Maniold model");
   }
 
   AMREX_GPU_HOST_DEVICE

Source/Spray/Drag.H

@@ -439,22 +439,15 @@ calculateSpraySource(
           amrex::Real X_norm = X_vapor[spf] / sumXVap;
           Ddiag[fspec] *= X_norm;
           rhoDtotal += Ddiag[fspec];
-          cp_fuel += Y_skin[fspec]; // * cp_n[fspec];
-
           sumYSkin += Y_skin[fspec];
         }
 
-#ifndef USE_MANIFOLD_EOS
         amrex::GpuArray<amrex::Real, NUM_SPECIES> Y_fuel_skin;
         Y_fuel_skin.fill(0.0);
         for (int spf = 0; spf < SPRAY_FUEL_NUM; ++spf) {
           Y_fuel_skin[fdat.indx[spf]] = Y_skin[fdat.indx[spf]] / sumYSkin;
         }
         eos.TY2Cp(T_skin, Y_fuel_skin.data(), cp_fuel);
-#else
-        // Manifold EOS: for now, just taking mixture Cp rather than Fuel Cp
-        eos.TY2Cp(T_skin, Y_skin.data(), cp_fuel);
-#endif
         cp_fuel *= SprayUnits::eng_conv;
 
         amrex::Real Sc_skin = mu_skin / rhoDtotal;