exf_BulkRhn = (exf_one - stab)*cstanton_1 + stab*cstanton_2
exf_BulkRhn = 0.0
_RL function exf_BulkRhn(
_RL exf_BulkRhn
tmpbulk= exf_BulkRhn(stable)