HorizontalProblem hydr computation

src/qt/output/output_view.cpp

@@ -157,14 +157,18 @@ QtOutputView::setColor(size_t ix,size_t iz){
 
 void
 QtOutputView::drawVector(size_t ix,size_t iz,double u,double v){
+  double h=height();
+  double w=width();
+
   double x1=ix*dX;
   double y1=solver->get_geometry().Z[ix][iz];
   double x2=x1+u;
   double y2=y1+v;
-  double x3=x2+(-0.866*u+0.5*v)*scale_arrow_head*lX;
-  double y3=y2+(-0.5*u-0.866*v)*scale_arrow_head*factor;
-  double x4=x2+(-0.866*u-0.5*v)*scale_arrow_head*lX;
-  double y4=y2+(0.5*u-0.866*v)*scale_arrow_head*factor;
+  double x3=x2+(-0.866*u/lX*w+0.5*v/factor*h)*scale_arrow_head*lX/w;//*lX;
+  double y3=y2+(-0.5*u/lX*w-0.866*v/factor*h)*scale_arrow_head*factor/h;
+
+  double x4=x2+(-0.866*u/lX*w-0.5*v/factor*h)*scale_arrow_head*lX/w;//*lX;
+  double y4=y2+(0.5*u/lX*w-0.866*v/factor*h)*scale_arrow_head*factor/h;
 
   glColor3f(1,1,0);
   glBegin(GL_LINE_STRIP);

src/qt/output/output_view.hpp

@@ -14,7 +14,7 @@ private:
   static constexpr size_t const arrow_step_x=3;
   static constexpr size_t const arrow_step_z=7;
   static constexpr double const scale_arrow=12000;
-  static constexpr double const scale_arrow_head=0.1;
+  static constexpr double const scale_arrow_head=0.3;
 
   size_t nX;
   double dX,lX;

todo.txt

@@ -2,3 +2,7 @@
 HorizontalProblem::run
 Overland::run
 indice_x_Richards
+
+
+[Limitation]
+Ici k0 est constant -> impact sur horizontal problem