[WIP] Add NVP (Non-Vascular Plant) Layer in CTSM

[huitang-earth]

Description:

This PR introduces a physically explicit representation of non-vascular plants (NVP; e.g., moss and lichen) in CTSM–FATES and its biophysical effect by assigning them to a dedicated vertical layer in CTSM, and a NVP PFT in FATES. The core architectural decision is to use layer 0 of the soil-snow layer structure to accommodate NVP and its growth (i.e., changes in height) when it is present. This layer integrates consistently into CLM’s radiation (under snow), thermal, and hydrology schemes, while coupling to FATES (NGEET/fates#1556) for canopy radiation (without snow) and photosynthesis (see the following call graph).
Note: The implementation is not fully ready yet and is still open for discussion.

Key implementation details in CLM:

• Adds NVP as an explicit surface layer at index 0, controlled by col%nvp_layer_active and col%jbot_sno, with geometry (dz, z, zi) updated each daily timestep via UpdateNVPLayer.
• Extends snow radiation to include NVP absorption at layer 0 with NVP specific optical parameters (when under snow).
• Integrates NVP into the thermal solver with its own thermal properties.
• Adds NVP specific hydraulic properties, water retention curve, and water balance subroutine NVPWaterBalance_Column.
• Adds interface with FATES for exchanging NVP structure properties, radiation flux, temperature and water status.

Collaborators:

@rosiealice

Expectation of Answer Changes:

This PR is expected to produce answer changes when NVP is active and NVP layer is present:

• Surface energy balance modified due to new top boundary condition at layer 0 (CLM, FATES)
• Soil temperature profiles altered due to NVP thermal properties (CLM)
• Hydrology affected by additional water storage and evaporation pathway (CLM)
• Carbon fluxes (GPP/NPP) changed due to NVP photosynthesis and modified PAR environment (FATES)

No answer changes are expected when NVP is disabled.

Description of generative AI usage

Claude code has been used for porting changes into the new CTSM version and designing the overall structure. The call graph is also generated by Claude code.

Any User Interface Changes (namelist or namelist defaults changes): Yes.

Does this create a need to change or add documentation? Yes

[samsrabin]

I will change this to be frozen only if temperature is ≤ freezing. That would be consistent with WaterStateType%InitCold().

Some text in the subroutine description suggests this is handled elsewhere, but it's cryptic, and I don't understand it. It would be easier to understand if this subroutine chose liquid or frozen itself. (Make sure to update the subroutine description!)

[samsrabin]

average equally over all patches in the site

Why? I would have thought these should be weighted by the area of the column on each patch.

[samsrabin]

Shouldn't the case of small nvp_frac_sum be handled?

[samsrabin]

It looks like you're never passing the nvp-specific stomatal_model to FATES:

if (trim(fates_stomatal_model) == 'ballberry1987') then
   pass_stomatal_model = 1
else if (trim(fates_stomatal_model) == 'medlyn2011') then
   pass_stomatal_model = 2
end if
call set_fates_ctrlparms('stomatal_model',ival=pass_stomatal_model)

Is that intentional? There are, after all, parts of the FATES PR (e.g.) that look for the nvp_model value (3).