I k1local,k2local,
k12 = k1local*k2local
& + k1local*k2local)
_RL k1local, k2local
& - 4.0*k1local*k2local*(1.0-2.0*gamm)
& (1.0 + (k1local/hnew) + (k1local*k2local/(hnew*hnew)))
co3local = k1local*k2local*diclocal /
& (hnew*hnew + k1local*hnew + k1local*k2local)