IF (halo_dat(ihalo)(1:5) .EQ. 'Hn ') THEN
ELSEIF (halo_dat(ihalo)(1:5) .EQ. 'Hy ') THEN
& halo_dat(CW_MAX_LOC)
data halo_dat /
& '_', halo_dat(ihalo)(1:2), '__'