add
    call fftw_free(c%C_bar_DD_k_ptr)
    call fftw_free(c%C_tilde_short_k_ptr)
    deallocate(c%C_tilde_k)
matrices are put in coordinates_module.F