91   integer ipp2, idij, idlj, idot, ipph;  in passfb()  local
204     idij = 0;  in passfb()
206       idij += 2;  in passfb()
208         idij += 2;  in passfb()
210 …       c1_ref(i - 1, k, j) = wa[idij - 1] * ch_ref(i - 1, k, j) - fsign*wa[idij] * ch_ref(i, k, j);  in passfb()
211           c1_ref(i, k, j) = wa[idij - 1] * ch_ref(i, k, j) + fsign*wa[idij] * ch_ref(i - 1, k, j);  in passfb()
221       idij = idj;  in passfb()
223         idij += 2;  in passfb()
224         c1_ref(i - 1, k, j) = wa[idij - 1] * ch_ref(i - 1, k, j) - fsign*wa[idij] * ch_ref(i, k, j);  in passfb()
225         c1_ref(i, k, j) = wa[idij - 1] * ch_ref(i, k, j) + fsign*wa[idij] * ch_ref(i - 1, k, j);  in passfb()
1044   integer idp2, ipp2, idij, ipph;  in radbg()  local
1204       idij = is;  in radbg()
1206         idij += 2;  in radbg()
1208           c1_ref(i - 1, k, j) = wa[idij - 1] * ch_ref(i - 1, k, j)   in radbg()
1209             - wa[idij] * ch_ref(i, k, j);  in radbg()
1210           c1_ref(i, k, j) = wa[idij - 1] * ch_ref(i, k, j) + wa[idij] * ch_ref(i - 1, k, j);  in radbg()
1219         idij = is;  in radbg()
1221           idij += 2;  in radbg()
1222           c1_ref(i - 1, k, j) = wa[idij - 1] * ch_ref(i - 1, k, j)   in radbg()
1223             - wa[idij] * ch_ref(i, k, j);  in radbg()
1224           c1_ref(i, k, j) = wa[idij - 1] * ch_ref(i, k, j) + wa[idij] * ch_ref(i - 1, k, j);  in radbg()
1560   integer idp2, ipp2, idij, ipph;  in radfg()  local
1607         idij = is;  in radfg()
1609           idij += 2;  in radfg()
1611             ch_ref(i - 1, k, j) = wa[idij - 1] * c1_ref(i - 1, k, j)  in radfg()
1612               + wa[idij] * c1_ref(i, k, j);  in radfg()
1613             ch_ref(i, k, j) = wa[idij - 1] * c1_ref(i, k, j) - wa[  in radfg()
1614                                                                   idij] * c1_ref(i - 1, k, j);  in radfg()
1623           idij = is;  in radfg()
1625             idij += 2;  in radfg()
1626             ch_ref(i - 1, k, j) = wa[idij - 1] * c1_ref(i - 1, k, j)   in radfg()
1627               + wa[idij] * c1_ref(i, k, j);  in radfg()
1628             ch_ref(i, k, j) = wa[idij - 1] * c1_ref(i, k, j) - wa[  in radfg()
1629                                                                   idij] * c1_ref(i - 1, k, j);  in radfg()