<div dir="ltr">Hi all! <div>Continuing with Greg's question, I have the following doubt. </div><div>Implementing boundary conditions for the magnetic field using the unsplit staggered mesh, </div><div>in the boundary region: </div><div>Is essential to define the field either just in the cell center or just in the cell faces, or in both places is necesary? </div><div><br></div><div>Thank you!</div><div><br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr">El vie., 24 ago. 2018 a las 11:59, Ernesto Zurbriggen (<<a href="mailto:ezurbriggen@unc.edu.ar">ezurbriggen@unc.edu.ar</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Greg! <div>I'm not in the details, but remember that in cylindrical geometry (IAXIS,JAXIS,KAXIS)=(r,z,theta). </div><div><br></div><div>>>I am using the Unsplit Staggered Mesh MHD solver, and have included a modified Grid_bcApplyToRegion.F90 that alters the standard "reflect" condition on the HIGH face of IAXIS to: <br></div><div>So, if you would like to modify the BC on the HIGH face of IAXIS, </div><div><br></div><div>>> if(ivar == MAGZ_VAR) regionData(k-i,1:je,1:ke,ivar)= 1.0 <br></div><div>maybe you should put ivar == MAGX_VAR instead. </div><div><br></div><div>Best!</div><div><br></div><div>Ernesto </div><br class="m_-2940089929125774333gmail-Apple-interchange-newline"></div><br><div class="gmail_quote"><div dir="ltr">El vie., 24 ago. 2018 a las 5:48, Boyle, Gregory (<<a href="mailto:gregory.boyle@desy.de" target="_blank">gregory.boyle@desy.de</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear FLASH community,<br>
<br>
I am attempting to modify the the 2D NohCylindrical test problem (32.2.4) to include a permanent B_{\theta} condition at rightmost r boundary. This is a first step towards including the effect of a longitudinal current.<br>
<br>
I am using the Unsplit Staggered Mesh MHD solver, and have included a modified Grid_bcApplyToRegion.F90 that alters the standard "reflect" condition on the HIGH face of IAXIS to:<br>
<br>
k = 2*guard+1<br>
if(isFace)k=k+1<br>
do i = 1,guard<br>
   regionData(k-i,1:je,1:ke,ivar)= regionData(i,1:je,1:ke,ivar)*sign<br>
   if(ivar == MAGZ_VAR) regionData(k-i,1:je,1:ke,ivar)= 1.0<br>
end do<br>
<br>
However, when I look at the output it seems the condition of setting MAGZ_VAR to 1.0 is being ignored.<br>
<br>
Any suggestions would be greatly appreciated!<br>
<br>
Cheers,<br>
Greg.<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="m_-2940089929125774333gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font face="monospace, monospace" style="background-color:rgb(234,209,220)"><i><font color="#000000">Dr. Ernesto Zurbriggen</font></i></font><div><font face="monospace, monospace"><i style="background-color:rgb(234,209,220)"><br></i></font><div><span style="background-color:rgb(234,209,220)"><font face="monospace, monospace"><i><font color="#444444">Consejo Nacional de Investigaciones Científicas y Técnicas </font></i></font><i style="font-size:12.8px;font-family:monospace,monospace"><font color="#274e13">(CONICET)</font></i><i style="font-family:monospace,monospace;font-size:12.8px"><font color="#444444">, Instituto de Astronomía Teórica y Experimental </font><font color="#274e13">(IATE)</font><font color="#444444">.</font></i></span></div><div><i style="font-size:12.8px;color:rgb(68,68,68);font-family:monospace,monospace;background-color:rgb(234,209,220)"><br></i></div><div><span style="background-color:rgb(234,209,220)"><i style="font-size:12.8px;color:rgb(68,68,68);font-family:monospace,monospace">Universidad Nacional de Córdoba </i><i style="font-size:12.8px;font-family:monospace,monospace"><font color="#274e13">(UNC)</font></i><i style="font-size:12.8px;color:rgb(68,68,68);font-family:monospace,monospace">, </i><font face="monospace, monospace" style="font-size:12.8px"><i><font color="#444444">Observatorio Astronómico de Córdoba </font></i></font><i style="font-size:12.8px;font-family:monospace,monospace"><font color="#274e13">(OAC)</font></i><i style="font-size:12.8px;color:rgb(68,68,68);font-family:monospace,monospace">.</i></span></div></div><div><i style="color:rgb(68,68,68);font-family:monospace,monospace;font-size:12.8px;background-color:rgb(234,209,220)"><br></i></div><div><i style="color:rgb(68,68,68);font-family:monospace,monospace;font-size:12.8px;background-color:rgb(234,209,220)">Teléfono: +54 0351 4331064-5, interno 222.</i></div><div><span style="font-size:12.8px"><i><span style="background-color:rgb(234,209,220)"><font color="#0000ff"><u>Laprida 854, Bº Observatorio (CP 5000), Córdoba, Argentina</u></font><font color="#444444">.</font></span><br></i></span><font color="#444444"><br><br></font></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font face="monospace, monospace" style="background-color:rgb(234,209,220)"><i><font color="#000000">Dr. Ernesto Zurbriggen</font></i></font><div><font face="monospace, monospace"><i style="background-color:rgb(234,209,220)"><br></i></font><div><span style="background-color:rgb(234,209,220)"><font face="monospace, monospace"><i><font color="#444444">Consejo Nacional de Investigaciones Científicas y Técnicas </font></i></font><i style="font-size:12.8px;font-family:monospace,monospace"><font color="#274e13">(CONICET)</font></i><i style="font-family:monospace,monospace;font-size:12.8px"><font color="#444444">, Instituto de Astronomía Teórica y Experimental </font><font color="#274e13">(IATE)</font><font color="#444444">.</font></i></span></div><div><i style="font-size:12.8px;color:rgb(68,68,68);font-family:monospace,monospace;background-color:rgb(234,209,220)"><br></i></div><div><span style="background-color:rgb(234,209,220)"><i style="font-size:12.8px;color:rgb(68,68,68);font-family:monospace,monospace">Universidad Nacional de Córdoba </i><i style="font-size:12.8px;font-family:monospace,monospace"><font color="#274e13">(UNC)</font></i><i style="font-size:12.8px;color:rgb(68,68,68);font-family:monospace,monospace">, </i><font face="monospace, monospace" style="font-size:12.8px"><i><font color="#444444">Observatorio Astronómico de Córdoba </font></i></font><i style="font-size:12.8px;font-family:monospace,monospace"><font color="#274e13">(OAC)</font></i><i style="font-size:12.8px;color:rgb(68,68,68);font-family:monospace,monospace">.</i></span></div></div><div><i style="color:rgb(68,68,68);font-family:monospace,monospace;font-size:12.8px;background-color:rgb(234,209,220)"><br></i></div><div><i style="color:rgb(68,68,68);font-family:monospace,monospace;font-size:12.8px;background-color:rgb(234,209,220)">Teléfono: +54 0351 4331064-5, interno 222.</i></div><div><span style="font-size:12.8px"><i><span style="background-color:rgb(234,209,220)"><font color="#0000ff"><u>Laprida 854, Bº Observatorio (CP 5000), Córdoba, Argentina</u></font><font color="#444444">.</font></span><br></i></span><font color="#444444"><br><br></font></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>