How to distribute the CN column state vector into vertical state vector?

[jinmuluo]

Hi CLM Community,

CLM uses multiple variables to represent one species's different state in model; Like NH4, it has "nh4_col(c)" and "nh4_vr_col(c, d)". I know the column state vector is integrated from the vertical state vector like that: nh4_col(c) = sum(nh4_vr_col(c, :), dim=2); My question is that does CLM provides any way to distribute the nh4_col(c) into nh4_vr_col(c, d)?

Best,

Jinmu

 

[slevis]

Dear @jinmuluo,

First the model calculates nh4_vr_col(c, d) in every soil layer, d, and then it integrates this data to the whole column (as you said). So there is no need to use the column data to recalculate the distributed soil-layer data. I hope this explanation helps.

Sam Levis

 

[jinmuluo]

Dear @jinmuluo,

First the model calculates nh4_vr_col(c, d) in every soil layer, d, and then it integrates this data to the whole column (as you said). So there is no need to use the column data to recalculate the distributed soil-layer data. I hope this explanation helps.

Sam Levis

Hi Sam,

Thanks so much for your help! I do understand there is no need to redistribute the column data to soil layer data. But one part of my code needs to calculate something in whole column then subtract the column value, in this time I need to distribute this subtracted value into per layer. I used the "ndep_profile" to add this subtracted value into per layer like the addition of deposition flux, but some layers show negative value then model force me out. So, I want find a way to redistribute the new column value (positive) into per layer not add the subtracted value which turn some layer into negative, do you have any other suggestion? Thanks so much!

Best,

Jinmu

 

[slevis]

I'm sorry, I have no suggestions.

 

[jinmuluo]

I'm sorry, I have no suggestions.

Thank you anyway!