<html><body><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; color: #000000"><div><br></div><div>Hi Marco,<br></div><div><br data-mce-bogus="1"></div><div>I tried and it compiled . But before starting the time step it failed. <br data-mce-bogus="1"></div><div>Its a very simple case, the water flowing throgh a cube. So I have inlet, outlet, walls on upside and bottom and periodic conditions on side walls.<br data-mce-bogus="1"></div><div>I get this error:<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Matrix: 0 aa 6 6 6<br> 0 aa 8.33333E+01 -4.87842E+01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00<br> 0 aa -4.87842E+01 8.66667E+01 -4.87842E+01 7.11752E+00 -1.66667E+00 0.00000E+00<br> 0 aa 0.00000E+00 -4.87842E+01 8.33333E+01 -4.16667E+01 7.11752E+00 0.00000E+00<br> 0 aa 0.00000E+00 7.11752E+00 -4.16667E+01 8.33333E+01 -4.87842E+01 0.00000E+00<br> 0 aa 0.00000E+00 -1.66667E+00 7.11752E+00 -4.87842E+01 Infinity -Infinity<br> 0 aa 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 -Infinity Infinity<br><br>Matrix: 0 bb 6 6 6<br> 0 bb 4.16667E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00<br> 0 bb 0.00000E+00 1.66667E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00<br> 0 bb 0.00000E+00 0.00000E+00 4.16667E-02 0.00000E+00 0.00000E+00 0.00000E+00<br> 0 bb 0.00000E+00 0.00000E+00 0.00000E+00 4.16667E-02 0.00000E+00 0.00000E+00<br> 0 bb 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 8.33333E-03 0.00000E+00<br> 0 bb 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00<br><br>Matrix: 0 Aeig 6 6 6<br> 0 Aeig 8.33333E+01 -4.87842E+01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00<br> 0 Aeig -4.87842E+01 8.66667E+01 -4.87842E+01 7.11752E+00 -1.66667E+00 0.00000E+00<br> 0 Aeig 0.00000E+00 -4.87842E+01 8.33333E+01 -4.16667E+01 7.11752E+00 0.00000E+00<br> 0 Aeig 0.00000E+00 7.11752E+00 -4.16667E+01 8.33333E+01 -4.87842E+01 0.00000E+00<br> 0 Aeig 0.00000E+00 -1.66667E+00 7.11752E+00 -4.87842E+01 Infinity -Infinity<br> 0 Aeig 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 -Infinity Infinity<br><br>Matrix: 0 Deig 1 6 6<br> 0 Deig 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00<br> Error in generalev, info= 12 6 -6<br> <br>an error occured: dying ...<br><br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Could you please tell me why?<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Thank you<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Sijo GEORGE<br data-mce-bogus="1"></div><hr id="zwchr" data-marker="__DIVIDER__"><div data-marker="__HEADERS__"><b>De: </b>"nek5000-users" <nek5000-users@lists.mcs.anl.gov><br><b>À: </b>"nek5000-users" <nek5000-users@lists.mcs.anl.gov><br><b>Envoyé: </b>Jeudi 31 Mai 2018 14:52:40<br><b>Objet: </b>Re: [Nek5000-users] How to give NEK5000 boundary conditions to exodus mesh file<br></div><br><div data-marker="__QUOTED_TEXT__">I specify the boundary conditions in the usrdat subroutine using the following piece of code<div class=""><br class="">
</div>
<div class="">
<div class=""> do i=1,nelt</div>
<div class=""> do j=1,2*ndim</div>
<div class=""> if(bc(5,j,i,1).eq.1) then</div>
<div class=""> cbc(j,i,1)='mv '</div>
<div class=""> elseif(bc(5,j,i,1).eq.2) then</div>
<div class=""> cbc(j,i,1)='W '</div>
<div class=""> elseif(bc(5,j,i,1).eq.3) then</div>
<div class=""> cbc(j,i,1)='W '</div>
<div class=""> elseif(bc(5,j,i,1).eq.4) then</div>
<div class=""> cbc(j,i,1)='W '</div>
<div class=""> else</div>
<div class=""> cbc(j,i,1)='E '</div>
<div class=""> endif</div>
<div class=""> enddo</div>
<div class=""> enddo</div>
</div>
<div class=""><br class="">
</div>
<div class="">Marco</div>
<div><br class="">
<blockquote class="">
<div class="">On May 31, 2018, at 8:50 AM, <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="" target="_blank">
nek5000-users@lists.mcs.anl.gov</a> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div class="">
<div style="font-family: arial, helvetica, sans-serif; font-size: 12pt;" class="">
<div class="">Hi Neks,<br class="">
</div>
<div class=""><br class="">
</div>
<div class="">I did the conversion of a mesh from ICEM to exodus file then to re2 file. </div>
<div class=""><br class="">
</div>
<div class="">But I am bit confused on the boundary conditions. Because from ICEM before I export to exodus I give the B.Cs as side set numbers and after the conversion to re2 I wont be able to put the boundary conditions in NEK format.</div>
<div class=""><br class="">
</div>
<div class="">So could you please tell me how you put the BC in NEK5000 format to a mesh file with exodus format.<br class="">
</div>
<div class=""><br class="">
</div>
<div class="">Thank you<br class="">
</div>
</div>
</div>
_______________________________________________<br class="">
Nek5000-users mailing list<br class="">
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" class="" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br class="">
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br class="">
</div>
</blockquote>
</div>
<br class="">
<br>_______________________________________________<br>Nek5000-users mailing list<br>Nek5000-users@lists.mcs.anl.gov<br>https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br></div></div></body></html>