k12 = k1local*k2local
I 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)
& - 4.0 _d 0*k1local*k2local*(1.0 _d 0-2.0 _d 0*gamm)
& (1.0 _d 0 + (k1local/hnew) + (k1local*k2local/(hnew*hnew)))