RichardsFastSlow/src/kernel/overland.cpp

#include "overland.hpp"
#include "debug.hpp"

namespace Kernel{

  void Overland::run(){
    if(Debug::level>0) cout<<"   [Overland::run] start"<<endl;
    for(size_t ix=0;ix<geometry->nX;++ix){
      hov[ix]=in_hov[ix];
      Psoil[ix]=previous_P[ix][geometry->nZ[ix]-1];
    }
    error=0;
    if(Debug::level>0) cout<<"   [Overland::run] end"<<endl;
  }
}