<div dir="ltr"><div>Dear Friends,</div><div>I am experimenting with SodStep in cylindrical geometry. I run the following setup</div><div>./setup SodStep -auto -2d -geometry=cylindrical <br></div><div>When I change initial conditions and domain to ones shown in without.pdf file <br></div><div>(sim_rhoLeft     = 0.125<br><div>sim_rhoRight        = 1.<br>sim_pLeft         = 0.1<br>sim_pRight               = 1.<br>sim_uLeft         = 0.<br>sim_uRight                = 0.) </div>everything works just fine. However when I further change initial conditions adding following lines</div><div>          if ((yy < 0.188).or.(yy > 0.812)) then            <br>              rhoZone = sim_rhoLeft<br>              presZone = sim_pLeft<br>            endif   <br>to  Simulation_initBlock.F90 (shown in with.pdf) simulation starts to work fine but after first mesh refinement I receive "inconsistent grid state" and "missing neighbor" errors after which code aborts (please see sod.log). Any ideas on what causes this problem and what can be done would be greatly appreciated.</div><div>Many thanks!</div><div>MAksm Kozlov<br></div></div>