RockCompressibility::RockCompressibility(const EclipseGridParser& deck) : pref_(0.0), rock_comp_(0.0) { if (deck.hasField("ROCKTAB")) { const table_t& rt = deck.getROCKTAB().rocktab_; int n = rt[0][0].size(); p_.resize(n); poromult_.resize(n); for (int i = 0; i < n; ++i) { p_[i] = rt[0][0][i]; poromult_[i] = rt[0][1][i]; } } else if (deck.hasField("ROCK")) { const ROCK& r = deck.getROCK(); pref_ = r.rock_compressibilities_[0][0]; rock_comp_ = r.rock_compressibilities_[0][1]; } else { std::cout << "**** warning: no rock compressibility data found in deck (ROCK or ROCKTAB)." << std::endl; } }