#include "solution.hpp"

Solution::Solution(const Geometry& geometry){
  size_t nX=geometry.nX;
  P=new double*[nX];
  for(size_t ix=0;ix<nX;++ix){
    P[ix]=new double[geometry.nZ[ix]];
  }
  hydr=new double[nX];
  hov=new double[nX];
  hsat=new double[nX];
  l=new double[nX];
  Pl=new double[nX];
}