<div class="gmail_quote">On Sat, Apr 30, 2011 at 13:17, ilyas ilyas <span dir="ltr">&lt;<a href="mailto:ilyascfd@gmail.com">ilyascfd@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div>I guess the thing that makes it complicated application of boundary conditions.</div><div>Since XYZGHOSTED is not available in fortran,</div><div>I am using XYZPERIODIC in order to implement bc&#39;s, as it is done in ex31.c in SNES</div>

<div>Combining  XYZPERIODIC with larger block size is still not much clear for me.</div><div><br></div><div>By the way, what is the current status of DMDACreate3D </div><div>and other DMDA routines and their fortran support ?</div>

<div>According to ex11f90.F and ex22f.F in the manual page of DMDACreate3d ,</div><div>these routine(s) provide ghost cell support for fortran in 3D ?</div><div>If ghost cell support is available, implementation would be relatively easy</div>
</blockquote></div><br><div>There is Fortran support for periodic boundaries. Also, petsc-dev has support for ghost cells even when they do not imply any wrapping, see DMDA_BOUNDARY_GHOST.</div><div><br></div><div><a href="http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-dev/docs/manualpages/DM/DMDACreate3d.html">http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-dev/docs/manualpages/DM/DMDACreate3d.html</a></div>