if(reader->getNumberOfQuantities()!=para->getParD(level)->precursorBC.numberOfQuantities)throwstd::exception("Number of quantities in reader and number of quantities needed for precursor don't match!");
if(reader->getNumberOfQuantities()!=para->getParD(level)->precursorBC.numberOfQuantities)throwstd::runtime_error("Number of quantities in reader and number of quantities needed for precursor don't match!");