bug fixed in mk_crystal_unit.F.