Thank you very much for your advice! I have changed the parameter LMAX2 to 12 (before was 5) in modules.F which leads to improvement in three orders: <n, K | n, K> = 0.999999992267 compared to the previous value of 0.999990071859 More hints are very welcome. Best regards, Martin