Subprogram const1 (kchi, pc1): pc1 array should be dimensioned to 9, and not 8 (well, the bug is now corrected...)
Subprogram diadic (kv, ki, pv, pw, pc1, pm): dimension of pc1 should be (9), not (8). Well, now the bug is corected.