Create calculate_kernel.f90