parse integrated boundary flux instead of altitude

tduretz · Feb 20, 2024 · 728da6a · 728da6a
1 parent d899378
commit 728da6a
Show file tree

Hide file tree

Showing 2 changed files with 11 additions and 3 deletions.
diff --git a/MDLIB/Setup.c b/MDLIB/Setup.c
@@ -391,6 +391,10 @@ void SetBCs(SetBCs_ff setBCs, MdoodzInput *instance, grid *mesh, surface *topo)
     free(newBoundary);
   }
 
+  // override by integrated velocity values
+  instance->topo_height->west = VxWestSum*mesh->dz;
+  instance->topo_height->east = VxEastSum*mesh->dz;
+
   /* --------------------------------------------------------------------------------------------------------*/
   /* Set the BCs for Vz on all grid levels */
   /* Type  0: Dirichlet point that matches the physical boundary (Vx:

diff --git a/SETS/RiverTomCrouteinf.c b/SETS/RiverTomCrouteinf.c
@@ -129,9 +129,13 @@ SetBC SetBCVz(MdoodzInput *instance, POSITION position, Coordinates coordinates)
   const double surf_Winc = instance->model.surf_Winc;
   const double surf_Vinc = instance->model.surf_Vinc;
   const double Vz_corr   = surf_Winc*surf_Vinc / Lx;
-  const double hW        = instance->topo_height->west - instance->model.zmin;
-  const double hE        = instance->topo_height->east - instance->model.zmin;
-  const double VzS       = -0.5 * V_tot * (hW + hE) / Lx + Vz_corr;
+  // const double hW        = instance->topo_height->west - instance->model.zmin;
+  // const double hE        = instance->topo_height->east - instance->model.zmin;
+  // const double VzS       = -0.5 * V_tot * (hW + hE) / Lx + Vz_corr;
+  const double VxW_H     = instance->topo_height->west;
+  const double VxE_H     = instance->topo_height->east;
+  const double VzS       = -(fabs(VxW_H) + fabs(VxE_H))/Lx + Vz_corr;
+
 
   // Here we need the total inflow flux: ncell_W * dz * VxW + ncell_E * dz * VxE
   // ncell_W is the number of Vx points for which we apply a BC value on the west