def METHOD_NAME(parameter=None, level=None, **kwargs):
    """
    The latest composited and interpolated GFS data from S3
    Arguments
    ---------
    parameter : str
        parameter, e.g. 'SOIM'.
    level : str
        source depth, e.g. "0-10 m DPTH"
    Returns
    -------
    node : GFS
        GFS node with the latest forecast data available for the given parameter and level.
    """
    s3 = s3fs.S3FileSystem(anon=True)
    # get latest date
    prefix = "%s/%s/%s/" % (BUCKET, parameter, level)
    dates = [path.replace(prefix, "") for path in s3.ls(prefix)]
    if not dates:
        raise RuntimeError("No data found at '%s'" % prefix)
    date = max(dates)
    # get latest hour
    prefix = "%s/%s/%s/%s/" % (BUCKET, parameter, level, date)
    hours = [path.replace(prefix, "") for path in s3.ls(prefix)]
    if not hours:
        raise RuntimeError("No data found at '%s'" % prefix)
    hour = max(hours)
    # node
    return GFS(parameter=parameter, level=level, date=date, hour=hour, **kwargs)