Read Wbar from table (#619)

* 1. Added a new index for reading molecular mass from table. 2. Updated Y2Wbar function to interpolate Wbar value. To be used with PR#9 of cmlm

* Removed fixing Wbar to -1

Author: SreejithNREL

Source/Eos/EosParams.H

@@ -43,6 +43,7 @@ struct EosParm<Manifold>
   density_lookup_type dens_lookup;
   bool compute_temperature{false};
   int idx_density{0};
+  int idx_Wbar{0};
   int idx_T{0};
   int idx_Wdot[MANIFOLD_DIM];
   int is_variance_of[MANIFOLD_DIM];
@@ -146,6 +147,11 @@ struct InitParm<eos::EosParm<eos::Manifold>>
 
     // Get important indices
     parm_in->m_h_parm.idx_T = get_var_index("T", h_manf_data_in);
+#ifdef PELE_USE_SPRAY
+    parm_in->m_h_parm.idx_Wbar = get_var_index("WBAR", h_manf_data_in);
+#else
+    parm_in->m_h_parm.idx_Wbar = get_var_index("WBAR", h_manf_data_in, false);
+#endif
 
     // For manifold table parameter source terms, assume if index not found,
     // source term is 0 For neural net, require a definition to be supplied for

Source/Eos/Manifold.H

@@ -78,9 +78,13 @@ struct Manifold
 
   AMREX_GPU_HOST_DEVICE
   AMREX_FORCE_INLINE
-  static void Y2WBAR(const amrex::Real* /*Y*/, amrex::Real& /*WBAR*/)
+  void Y2WBAR(const amrex::Real* Y, amrex::Real& WBAR)
   {
-    amrex::Error("Y2WBAR is not yet implemented for Manifold EOS");
+    if (eosparm->idx_Wbar >= 0) {
+      manfunc.get_func()->get_value(eosparm->idx_Wbar, Y, WBAR);
+    } else {
+      amrex::Error("Wbar not in table");
+    }
   }
 
   AMREX_GPU_HOST_DEVICE