void Parameter<EvalT, TRAITS>::
evaluateFields(typename TRAITS::EvalData workset)
{ 
  //std::cout << "ROGER ParamValue = " << param->getValue() << std::endl;

  for (index_t cell = 0; cell < workset.num_cells; ++cell) {
    for (typename PHX::MDField<ScalarT, Cell, Point>::size_type pt = 0;
	 pt < target_field.dimension(1); ++pt) {
      target_field(cell,pt) = param->getValue();
    }
  }

}
void MultiVariateParameter<EvalT, TRAITS>::
evaluateFields(typename TRAITS::EvalData workset)
{
  ScalarT sum = 0;
  const int num_param = param.size();
  for (int i=0; i<num_param; ++i)
    sum += param[i]->getValue();

  for (std::size_t cell = 0; cell < workset.num_cells; ++cell) {
    for (typename PHX::MDField<ScalarT, Cell, Point>::size_type pt = 0;
         pt < target_field.dimension(1); ++pt) {
      target_field(cell,pt) = sum;
    }
  }

}