Extracting depth to water table

[MABOSSANI AKARA]

Hello All,

I am using Parflow to simulate temporal changes in water table.

I am using an orthogonal grid with solids to define the active model domain.

I was able to run the model, but I am having issues extracting the depth to water table. The saturation file shows that the position of the water table is correct. 

But when I try extracting the depth to water using  "hydro.calculate_water_table_depth(pressure, saturation, dz)", the output does not look right.

I have attached the plotting ('plot_wt_og.py.txt') scripts and will appreciate any suggestions.