FATES not recognizing changes to zbedrock

[KatieMurenbeeld]

Hi FATES folks. I’m helping my lab group to set up some experiments using FATES where we want to change the depth to bedrock. After playing around, I’ve noticed that when I change zbedrock in the surface data file (generated using the subset_data script from ctsm/tools/site_and_regional) CLM without FATES will recognize the modified depth to bedrock, but FATES does not. FATES doesn’t recognize zbedrock at all. Is this correct? Is there a way to modify depth to bedrock for FATES?

For the image below the FATES and CLM cases were run the same expect for use_fates = true or false.

Let me know if you need more information.

[ekluzek]

@KatieMurenbeeld there's a variable called use_bedrock which by default is off when FATES is on. I'm not sure this is required. But, likely if you explicitly add

use_bedrock=TRUE

in your user_nl_clm you can at least see if it will run. Since, bedrock is a below ground feature and FATES is above ground, I would think it has a chance at working. But, there must be some reason we turned bedrock off for FATES.

[KatieMurenbeeld]

@ekluzek I reran the simulations with use_bedrock = .true. in the user_nl_clm file, and FATES did recognize the depth to bedrock changes.

[wwieder]

@glemieux is this draft PR still something that can be brought in?

[wwieder]

@adrifoster should we mention this in the technote related to FATES?

[glemieux]

@wwieder I think it needs more work. It looked like an easy fix, but we ran into unexpected issues when we tested it. That said, I can bring it up to date and characterize the scope of work at least.