Dear all, It looks like a bug in espresso-3.1/PH/clinear.f90 and espresso-3.1.1/PH/clinear.f90: from line 14: nkr = nk2*nk2*nk3 I think it should be: nkr = nk1*nk2*nk3 Is it right? Best Regards