def METHOD_NAME(options):
    """
    Create a `Flow` for phonon calculations including
    eps_inf, Born effective charges and  dynamical quadrupoles
    """
    # Working directory (default is the name of the script with '.py' removed and "run_" replaced by "flow_")
    if not options.workdir:
        options.workdir = os.path.basename(sys.argv[0]).replace(".py", "").replace("run_", "flow_")
    pseudos = abidata.pseudos("31-Ga.LDA.fhi", "15-P.LDA.fhi")
    # Initialize the structure from an abinit string
    # Other approaches (e.g. from cif files) are available as well.
    structure = abilab.Structure.from_abistring("""
 acell 1.0 1.0 1.0
 rprim 0.0 5.023016083002 5.023016083002
       5.023016083002 0.0 5.023016083002
       5.023016083002 5.023016083002 0.0
 # Definition of the atom types and positions
 # ******************************************
 ntypat 2
 znucl  31 15
 natom 2
 typat 1 2
 xred
         0.0000000      0.0000000      0.0000000
         0.2500000      0.2500000      0.2500000