...

Examples/RMI/RMI_hybrid.py

@@ -123,10 +123,10 @@ def calc_maxvelo(simu, v):
 
 ctime = MPI.Wtime()
 # Domain
-box = hysop.Box(length=[2., 2., 4.], origin=[-1., -1., -2.],
+box = hysop.Box(length=[2., 2., 2.], origin=[-1., -1., -1.],
                 proc_tasks=PROC_TASKS)
-bc_b = SubBox(length=[2., 2., 0.1], origin=[-1., -1., -2.], parent=box)
-bc_t = SubBox(length=[2., 2., 0.1], origin=[-1., -1., 1.9], parent=box)
+bc_b = SubBox(length=[2., 2., 0.1], origin=[-1., -1., -1.], parent=box)
+bc_t = SubBox(length=[2., 2., 0.1], origin=[-1., -1., 0.9], parent=box)
 mpi_params = MPIParams(comm=box.comm_task, task_id=PROC_TASKS[main_rank])
 mpi_params_S = MPIParams(comm=box.comm_task, task_id=TASK_SCALAR)
 mpi_params_UW = MPIParams(comm=box.comm_task, task_id=TASK_UW)

HySoP/hysop/operator/discrete/baroclinic.py

@@ -136,7 +136,7 @@ class Baroclinic(DiscreteOperator):
             self._tempGrad[d][iCompute] *= self.viscosity
             self._result[d][iCompute] -= self._tempGrad[d][iCompute]
         # gravity term : result = result - g
-        self._result[2][iCompute] -= 9.81
+        self._result[2][iCompute] -= 1.  # 9.81
 
         # baroclinicTerm = -(gradRho/rho) x (gradP/rho)
         self._tempGrad = self._gradOp(self.density[0:1], self._tempGrad)