<div dir="ltr"><div><div><div><div>Thank you Klaus for getting back to me.<br><br></div>Looking  into CCSN/Grid_markRefineDerefine.F<wbr>90 (which is the closest to what I am looking for) I see that Grid_markRefineSpecialized is called right before Grid_markRefineDerefine returns. I made the call from the same location in my setup, and I still see no change. Following are the lines in my Grid_markRefineDerefin routine :<br><br>    <i>spec(1:3) = 0.0   ! setting origin coordinates, and of course spec was initially declared as a 4 dimensional array <br>    spec(4) = 8.0      ! limiting radius = 8 cm <br>    call Grid_markRefineSpecialized(<wbr>INRADIUS,4,spec,1)  ! All blocks below r = 8cm should be of refinement level = 1</i><br></div><br> </div>Unfortunately,  I get blocks below r = 8cm with refinement level > 1. Refer to attached pictures in which the AMR grid is displayed on pressure field.  <br>I am not sure whether the grid dimension and/or geometry makes any difference, but It seems not to work for me.<br><br></div><div><div>Best,<br><br></div><div>Ismael<br></div></div></div>