<div>Thank you Jed,<br></div><div><br></div><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'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>
<div><br></div><div>Cheers,</div><div>Ilyas</div><br><div class="gmail_quote">2011/4/28 Jed Brown <span dir="ltr"><<a href="mailto:jed@59a2.org">jed@59a2.org</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div class="im"><div class="gmail_quote">On Thu, Apr 28, 2011 at 15:41, ilyas ilyas <span dir="ltr"><<a href="mailto:ilyascfd@gmail.com" target="_blank">ilyascfd@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Would you please explain more the first option?</blockquote></div><br></div><div>You create a system with larger block size. If you have two node-centered values and one cell-centered, you would use a block size of 3. Each center would be associated with the node to its lower left, for example. There will be a fringe of "cell-centers" that extend out of domain on the right and top, you use x = 0 for these equations and the Jacobian for those equations will just have a 1 on the diagonal.</div>
</blockquote></div><br>