import gspy_interface_robust

print('Groundwater response progression:')
for i in range(5):
    result = gspy_interface_robust.run_modflow_step(0.075)
    print(f'Call {i+3}: {result[0]:.4f} m³/s')

print('\nTesting response to different recharge rates:')
for rate in [0.0, 0.05, 0.1, 0.15]:
    result = gspy_interface_robust.run_modflow_step(rate)
    print(f'Recharge {rate:.2f} m/day → {result[0]:.4f} m³/s')