40   40   40         : nrg, ntg, npg ; Grid numbers in (r,theta,phi)
   20                   : nrgin         ; Grid numbers in  rgin < r < rgmid
   10                   : nlg           ; Number of multipoles l
 0.000e+00              : rgin          ; Computational domain in r.
 2.000e+00              : rgmid
 1.000e+01              : rgout  


  read(2,'(4i5)') nrg, ntg, npg
  read(2,'(3i5)') nrgin
  read(2,'(3i5)') nlg
  read(2,'(1e10.3)') rgin
  read(2,'(1e10.3)') rgmid
  read(2,'(1e10.3)') rgout