Extracting heads from Modflow binary files using flopy

[zargol7419]

Hi everybody
I wanna extract heads derived from running an existing modflow model by using flopy package. I have executed different codes but I haven't received my expected answer. It seems that expected heads are not calculated before or I have executed a wrong code. I can't understand what is the problem.
I want to save all calculated heads in an array and then use them for another calculation. Could you please provide me with a solution for this problem.
Here is my executed codes:
First one:

Second one:

Thanks in advance for your answers.

[jdhughes-usgs]

Your first example has extracted the heads for all times. The data has dimensions [output_times, nlay, nrow, ncol], where output_times is the number of times you have output binary head data (specified in the oc file). It looks like you are using MODFLOW-2000, so the -999 values are either the HNOFLO or HDRY value you have specified in your input files. You can retrieve the head for individual times using .get_data().

[zargol7419]

Thank you so much. I understood what should I do. I really appreciate your cooperation.