Why I cannot load UStoredepth corretly

[pipipercy]

Hi, i have confronted an issue when run the script (https://github.com/openstreams/wflow/tree/master/notebooks/wflow-test-sbm.ipynb). in [7], i have gotten the map dislike the map in the demo, as follows,

which is different from the figure as in the notebook like

and i have tried to directly import this data using

import pcraster as pcr

demfilename = "instate/UStoreDepth.map"
dem = pcr.readmap(demfilename)

and got the map as follows,

i am confused about the difference and have not find the dir in the original code source.

Could you please help me figure it out? Thank you very much!

[verseve]

Hi, I made a couple of changes to the notebook so it is line with the current code (these haven't been update for quite a while), and should work now.

UstoreDepth is not a state anymore since we started to support multiple soil layers, this is now UStoreLayerDepth (a list). UstoreDepth is still used in the code, and is the sum of UStoreLayerDepth (amount of available water in the unsaturated zone). Also, FirstZoneCapacity has been replaced by SatWaterDepth.