<div dir="ltr">Hi Emmanuel,<br><br>1. The Circuits model changed a decent bit from FLASH4.7.1 to FLASH4.8<div>2. I don't know that 2D cartesian is supported for CESZAR currently. Only 1D unit tests for CESZAR exist that I've seen.</div><div><br></div><div>General-ish debugging help:</div><div> - navigate to your object directory</div><div> - grep (-i to be case insensitive because Fortran is case insensitive) the abort message (or a suitable subset, in this case "getCellFaceArea")</div><div><br></div><div>In your case, you will see that gr_getCellFaceArea.F90 expects to only compute face areas for the K-face if the problem is 3D (whereas you're doing 2D so it aborts).</div><div><br></div><div>Next general step: find what is calling the offending subroutine<br> - grep -in "call gr_getCellFaceArea" *.F90</div><div> - grep -in "call Grid_getBlkData" *.F90</div><div><br></div><div>-> Circuit_fluxCalc.F90<br></div><div><br></div><div>copy to your problem directory (probably CurrentSheetCESZAR) and modify based on your desired problem (looks like you'd want to compute fluxes in y not z).</div><div><br><div><div><div>Best wishes,<br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr">--------<div>Ryan</div></div></div></div></div></div></div><br></div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Jul 10, 2024 at 11:43 AM Emanuel Escobar <<a href="mailto:emanuel.escobar.lavin@gmail.com">emanuel.escobar.lavin@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Dear Flash community, I hope this message finds everyone well. I am writing about a problem I can't fix now. I am trying to simulate a modification of the CurrentSheet problem. In my case, I have added boundary conditions as shown in the file ProblemScheme.jpg and changed the variable condition of the Magnetic field by a constant value. At this point, the changes are all working well.</div><div>However, when I try to add source terms associated with a current model (CESZAR in my case) the simulation is aborted (abort_message: getCellFaceArea : this face doesn't exit). The error messages are shown in the files Screenshot_from_2024-07-10 14-33-30.png and currentSheet_mhd_2d.log.</div><div>Please, if anyone can give me some idea of how I can fix and understand this problem I would really appreciated.</div><div>Best regards<br></div><div><div><div><div dir="ltr" class="gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Emanuel Escobar Lavín<br>B.Eng., Aerospace Engineer<br>
<div>M.Sc.Eng., major in Mechanics<br></div><div>Phone: +56 (2) 29761906 - +56 (9) 72133822<br></div> </div></div></div></div></div></div></div></div></div></div>
_______________________________________________<br>
flash-users mailing list<br>
<a href="mailto:flash-users@flash.rochester.edu" target="_blank">flash-users@flash.rochester.edu</a><br>
<br>
For list info, including unsubscribe:<br>
<a href="https://flash.rochester.edu/mailman/listinfo/flash-users" rel="noreferrer" target="_blank">https://flash.rochester.edu/mailman/listinfo/flash-users</a><br>
</blockquote></div>