<div dir="ltr"><div dir="ltr"><img class="gmail-ajT" src="https://ssl.gstatic.com/ui/v1/icons/mail/images/cleardot.gif">Hi Aodhan! <div><br></div><div>Defining the initial magnetic field (bx,by,bz) in the way below for a 2D-cartesian mesh at least worked for me. </div><div><br></div><div>I hope that it can help you!</div><div><br></div><div><br></div><div><font face="monospace">! Define face quantities<br>! Cell face-centred variables for the StaggeredMesh scheme<br></font></div><font face="monospace"> do k=blkGC(LOW,KAXIS),blkGC(HIGH,KAXIS)<br> do j=blkGC(LOW,JAXIS),blkGC(HIGH,JAXIS)<br> do i=blkGC(LOW,IAXIS),blkGC(HIGH,IAXIS)<br> . </font><div><font face="monospace"> . </font></div><div><font face="monospace"> . <br> Sfacex(MAG_FACE_VAR,i,j,k)=Bx<br> Sfacey(MAG_FACE_VAR,i,j,k)=By<br><br> !*** Campo magnético Bz<br> ! Just possible cell-centred<br> Sfacez(MAG_FACE_VAR,i,j,k)=Bz<br><br> enddo<br> enddo<br> enddo<br></font><div><font face="monospace"> </font></div><div><span style="font-family:monospace">! Define centred quantities</span></div><div><font face="monospace"> do k = blkGC(LOW,KAXIS),blkGC(HIGH,KAXIS)<br> do j = blkGC(LOW,JAXIS),blkGC(HIGH,JAXIS)<br> do i = blkGC(LOW,IAXIS),blkGC(HIGH,IAXIS)<br> . </font></div><div><font face="monospace"> . </font></div><div><font face="monospace"> . <br> S(MAGX_VAR,i,j,k)=0.5*(Sfacex(MAG_FACE_VAR,i,j,k)+Sfacex(MAG_FACE_VAR,i+1, j , k ))<br> S(MAGY_VAR,i,j,k)=0.5*(Sfacey(MAG_FACE_VAR,i,j,k)+Sfacey(MAG_FACE_VAR, i ,j+1, k ))<br> S(MAGZ_VAR,i,j,k)= Sfacez(MAG_FACE_VAR,i,j,k) <br> . </font></div><div><font face="monospace"> . <br> enddo<br> enddo<br> enddo</font><br></div><div><br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El vie, 22 jul 2022 a las 8:07, Aodhan McIlvenny (<<a href="mailto:amcilvenny01@qub.ac.uk">amcilvenny01@qub.ac.uk</a>>) escribió:<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 lang="EN-GB" style="overflow-wrap: break-word;">
<div class="gmail-m_-3459943047341250994WordSection1">
<p class="MsoNormal"><span style="color:black">Hello all,<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:black"> <u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:black">I’ve set up an MHD simulation with the LaserSlab set up with a magnetic field by combining a few examples. The code works as expected when defining a magnetic field in the plane of the
2D simulation (cartesian or cylindrical) but computes a 0 or negative dt when defining it in the orthogonal direction (z) which is not resolved. Am I missing something obvious to run the simulation with a uniform magnetic field transverse to the simulation
plane? <u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:black"> <u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:black">Thanks in advance<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:black"> <u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:black">Best wishes,<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:black">Aodhan<u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
</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><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr" style="color:inherit;font-family:inherit;font-style:inherit;font-variant:inherit;font-weight:inherit;margin:0px;padding:0px;border:0px;font-stretch:inherit;line-height:inherit;vertical-align:baseline"><font style="background-color:rgb(234,209,220)" face="monospace, monospace"><font color="#000000"><span style="font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;margin:0px;padding:0px;border:0px;vertical-align:baseline;color:inherit"><div style="font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;margin:0px;padding:0px;border:0px;vertical-align:baseline;color:inherit"><i style="font-family:"Courier New",monospace;background-color:rgb(255,255,255);color:rgb(68,68,68);line-height:normal">Dr. Ernesto Zurbriggen, <br></i></div><div style="font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;margin:0px;padding:0px;border:0px;vertical-align:baseline;color:inherit"><i style="font-family:"Courier New",monospace;background-color:rgb(255,255,255);color:rgb(68,68,68);line-height:normal">Investigador Asistente (CONICET)</i><i style="color:rgb(68,68,68)"><br></i></div></span></font></font><font style="background-color:rgb(234,209,220)" face="monospace, monospace"><font color="#000000"><span style="font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;margin:0px;padding:0px;border:0px;vertical-align:baseline;color:inherit"><div style="font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;margin:0px;padding:0px;border:0px;vertical-align:baseline;color:inherit"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;vertical-align:baseline;color:inherit;background-color:rgb(255,255,255)"><i style="color:rgb(68,68,68)"></i></span></div></span></font></font></div><div dir="ltr" style="color:inherit;font-family:inherit;font-style:inherit;font-variant:inherit;font-weight:inherit;margin:0px;padding:0px;border:0px;font-stretch:inherit;line-height:inherit;vertical-align:baseline"><div style="margin:0px;padding:0px;border:0px;font-stretch:inherit;line-height:inherit;font-family:Calibri,Arial,Helvetica,sans-serif;vertical-align:baseline;color:rgb(0,0,0)"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:normal;font-family:"Courier New",monospace;vertical-align:baseline;color:inherit"><i style="color:rgb(68,68,68);line-height:normal"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;vertical-align:baseline;color:inherit"><i style="font-family:monospace,monospace;background-color:rgb(234,209,220)"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;vertical-align:baseline;color:inherit"><div style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;vertical-align:baseline;color:inherit"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:normal;font-family:"Courier New",monospace;vertical-align:baseline;color:inherit;background-color:rgb(255,255,255)">Correo:</span><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:normal;font-family:inherit;vertical-align:baseline;color:inherit;background-color:rgb(255,255,255)"><i style="line-height:normal"> <a href="mailto:ernesto.zurbriggen@craam.mackenzie.br" target="_blank">ezurbriggen@</a><a href="http://unc.edu.ar" target="_blank">unc.edu.ar</a></i></span><i><br></i></div><div style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;vertical-align:baseline;color:inherit"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:normal;font-family:inherit;vertical-align:baseline;color:inherit;background-color:rgb(255,255,255)"><i style="line-height:normal">Teléfono<font color="#888888"> </font></i></span><i><font color="#888888"><span style="color:rgb(102,102,102)"><span style="font-size:12.8px">IP: (+54) 351-535 3776</span></span></font> (int. 75601)<br></i></div></span></i></span><br></i></span></div><div style="margin:0px;padding:0px;border:0px;font-stretch:inherit;line-height:inherit;font-family:Calibri,Arial,Helvetica,sans-serif;vertical-align:baseline;color:rgb(0,0,0)"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:normal;font-family:"Courier New",monospace;vertical-align:baseline;color:inherit"><i style="color:rgb(68,68,68);line-height:normal"><i style="font-weight:inherit;line-height:normal">Instituto de Astronomía Teórica y Experimental (IATE-CONICET), </i></i></span></div><div style="margin:0px;padding:0px;border:0px;font-stretch:inherit;line-height:inherit;font-family:Calibri,Arial,Helvetica,sans-serif;vertical-align:baseline;color:rgb(0,0,0)"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:normal;font-family:"Courier New",monospace;vertical-align:baseline;color:inherit"><i style="color:rgb(68,68,68);line-height:normal"><i style="font-weight:inherit;line-height:normal">calle Francisco Laprida 922</i></i></span><i style="font-family:"Courier New",monospace;font-weight:inherit;color:rgb(68,68,68);line-height:normal">, CP 5000, Córdoba, Argentina. </i></div></div></div><div dir="ltr"><div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>