I'm trying to reproduce Medlyn gs values calculated by CLM by using the constituent variables applied to equation 9.1 of the CLM5 technote, but my values are coming out different to the gs values written by CLM (by low ten's to several thousand). Should I expect eq 9.1 to give me the same value as CLM?
I note in the code there is an equivalent check for Ball-Berry gs values, from PhotosynthesisMod:
I note in the code there is an equivalent check for Ball-Berry gs values, from PhotosynthesisMod:
These checks pass without a problem but there's not an equivalent check like this for Medlyn gs (which is what I'm trying to do myself). I'm not extremely familiar with these models and equations so hopefully not missing something obvious, any insight would be appreciated.gs_mol_err = mbb(p)*max(an(p,iv), 0._r8)*hs/cs*forc_pbot(c) + bbb(p)
if (abs(gs_mol(p,iv)-gs_mol_err) > 1.e-01_r8) then
write (iulog,*) 'Ball-Berry error check - stomatal conductance error:'
write (iulog,*) gs_mol(p,iv), gs_mol_err
end if