How to find out grid number of failure cells?

[Qing-sudo]

Hi all,

I'm running a mf6 model using flopy. The model has reported 'convergence failure' error (see the error message below). I wonder how can I know at which specific cell(s) the model is not converged?

I tried some functions like 'get_failure_cells()' that might be available in mf05 or other versions, but there seems no direct function in FloPy (in terms of mf6 version) to achieve this purpose.

The error message:

Writing simulation list file: mfsim.lst
Using Simulation name file: mfsim.nam

Solving:  Stress period:     1    Time step:     1

 Run end date and time (yyyy/mm/dd hh:mm:ss): 2023/06/01 20:15:50
 Elapsed run time: 22.891 Seconds


ERROR REPORT:

 1.  Simulation convergence failure occurred 1 time(s).
Premature termination of simulation.

The code I tried to debug my model:

model = flopy.mf6.MFSimulation.load(sim_ws='modelpath')
dis = model.get_model('modelname').get_package('dis')
rows, cols = dis.get_failure_cells()

Would be grateful to have some help here.