def METHOD_NAME(ngkpt):
    """
    This function constructs the input file for the GS calculation:
    """
    structure = dict(
     angdeg= 3*[60.0],
     acell=3* [ 7.1992351952],
     natom=2,
     ntypat=2,
     typat=[1, 2],
     znucl=[31, 15],
     xred=[
        0.0000000000, 0.0000000000, 0.0000000000,
        0.2500000000, 0.2500000000, 0.2500000000,
    ])
    #pseudos = ["Ga.psp8", "P.psp8"]
    pseudos = abidata.pseudos("Ga.oncvpsp", "P.psp8")
    gs_inp = abilab.AbinitInput(structure, pseudos=pseudos)
    gs_inp.set_vars(
        nband=8,
        #nband=13,
        ecut=30.0,
        #ecut=24.0,
        ngkpt=ngkpt,
        nshiftk=1,
        shiftk=[0, 0, 0],
        tolvrs=1.0e-10,
        nstep=150,
    )
    return gs_inp