<div><span style="color:rgb(0,0,0);font-family:Helvetica;font-size:medium">The fact is I do use a uniform cartesian grid, but the data structure in the code is in "AIJ" instead of the "IJK" format.</span><br style="color:rgb(0,0,0);font-family:Helvetica;font-size:medium">
<span style="color:rgb(0,0,0);font-family:Helvetica;font-size:medium">This was done because the ghost cells (from the immersed boundary method) were updated from multiple directions and this cannot be handled using the "ijk" format.</span><br style="color:rgb(0,0,0);font-family:Helvetica;font-size:medium">
<span style="color:rgb(0,0,0);font-family:Helvetica;font-size:medium">The computational grid is a staggered grid with velocities at the cell edges and pressure at cell centres.</span><br style="color:rgb(0,0,0);font-family:Helvetica;font-size:medium">
<span style="color:rgb(0,0,0);font-family:Helvetica;font-size:medium">I guess , that puts the DMDA out of the solution then.</span><br style="color:rgb(0,0,0);font-family:Helvetica;font-size:medium"><br style="color:rgb(0,0,0);font-family:Helvetica;font-size:medium">
<span style="color:rgb(0,0,0);font-family:Helvetica;font-size:medium">So, though I am on a structured uniform cartesian grid (planning on local mesh refinement in the future), the boundary treatment doesn't allow me to use the procedure that would be normally followed.</span><br style="color:rgb(0,0,0);font-family:Helvetica;font-size:medium">
<span style="color:rgb(0,0,0);font-family:Helvetica;font-size:medium">Am I right in this assessment? So then I wait for the manual on DMComplex and get pets-dev.</span><br style="color:rgb(0,0,0);font-family:Helvetica;font-size:medium">
<span style="color:rgb(0,0,0);font-family:Helvetica;font-size:medium">What do you suggest?</span><br style="color:rgb(0,0,0);font-family:Helvetica;font-size:medium"><span style="color:rgb(0,0,0);font-family:Helvetica;font-size:medium">Thanks a lot!</span><br>
</div>-- <br><div dir="ltr">Arun Kamath <div>PhD Stipendiat</div><div>NTNU, Trondheim</div></div><br>