Charbonnier diffusivity function added