ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c---> ior=order of SL remapping accuracy/2;
c---> only even order schemes are considered
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

      parameter (ior=1)

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c---> ihlag is a SL halo size; ihlag .ge. max(C, ior + 1, 3)
c---> if ior=1,2 => ihlag=3
c---> if ior=3   => ihlag=4
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

      parameter (ihlag=3)