If nrootlayer is greater than nlayer, most of the wat*.dat files are populated with zeros. No warning is given that nrootlayer is greater than nlayer, nor is there a mention in the documentation that these values should make sense together.
1.) Edit the documention to make this explicitly clear
2.) Include a routine that compares nrootlayer to nlayer. If it is larger, set nrootlayer equal to nlayer
pseudocode:
if nrootlayer>nlayer: nrootlater=nlayer
If nrootlayer is greater than nlayer, most of the wat*.dat files are populated with zeros. No warning is given that nrootlayer is greater than nlayer, nor is there a mention in the documentation that these values should make sense together.
1.) Edit the documention to make this explicitly clear
2.) Include a routine that compares nrootlayer to nlayer. If it is larger, set nrootlayer equal to nlayer
pseudocode:
if nrootlayer>nlayer: nrootlater=nlayer