What is going on: It looks like the layer homogeneity is causing a

NoneType Blow up after 8000 timesteps about lgar-py HOT 3 CLOSED

taddyb commented on August 17, 2024

NoneType Blow up after 8000 timesteps

from lgar-py.

Comments (3)

taddyb commented on August 17, 2024

The code works for the whole month December for 2001 if run by itself. It looks to be a stability issue

from lgar-py.

taddyb commented on August 17, 2024

Here is the self.top_layer.print() that showed me there was an issue:

self.print()
[2023-07-17 14:30:01,932][models.physics.layers.Layer][INFO] - Printing the top_layer for rank: 4
[2023-07-17 14:30:01,933][models.physics.layers.Layer][INFO] - [     Depth       Theta   Layer_num   dzdt     k_cm_hr    psi   ]
[2023-07-17 14:30:01,933][models.physics.layers.WettingFront][INFO] - [10.0000, 0.1561597759, 0, 0.000000, 0.000000, 2011.4011]
[2023-07-17 14:30:01,934][models.physics.layers.WettingFront][INFO] - [25.0000, 0.1561597759, 1, 0.000000, 0.000000, 2011.4011]
[2023-07-17 14:30:01,934][models.physics.layers.WettingFront][INFO] - [109.2912, 0.1561597759, 2, 0.000000, 0.000000, 2011.4011]

If you look at the theta values, they are all the same. This means that inserting water will affect the whole column, which isn't what happens.

A proposed solution is to create a NN to estimate Van Genuchten params from soil %. This can give heterogeneity for the whole column!

from lgar-py.

taddyb commented on August 17, 2024

Closing this for now. Not sure how to address in the scope of this project

from lgar-py.

Recommend Projects