def METHOD_NAME(nframes, atom_types):
    natoms = len(atom_types)
    data = {}
    data["coords"] = np.random.random([nframes, natoms, 3])
    data["forces"] = np.random.random([nframes, natoms, 3])
    data["cells"] = np.random.random([nframes, 9])
    data["energies"] = np.random.random([nframes, 1])
    types = list(set(atom_types))
    types.sort()
    data["atom_names"] = []
    data["atom_numbs"] = []
    for ii in range(len(types)):
        data["atom_names"].append("TYPE_%d" % ii)
        data["atom_numbs"].append(np.sum(atom_types == ii))
    data["atom_types"] = np.array(atom_types, dtype=int)
    return data