def METHOD_NAME(inps):
    print('\n-------------------READ INPUTS -------------------')
    print(f'Read metadata from file: {inps.file}')
    atr = readfile.read_attribute(inps.file)
    inps.pix_box, inps.geo_box = subset.subset_input_dict2box(vars(inps), atr)
    inps.outfile = inps.outfile if inps.outfile else atr['PROJECT_NAME']
    # date1/2 and dset
    ftype = atr['FILE_TYPE']
    if ftype in ['timeseries', 'HDFEOS']:
        date1, date2 = inps.dset.split('_')
        inps.dset = date2
    elif ftype == 'ifgramStack':
        date1, date2 = inps.dset.split('-')[1].split('_')
    else:
        # velocity, unw
        date1, date2 = ptime.yyyymmdd(atr['DATE12'].replace('_','-').split('-'))
        if inps.dset.startswith('step'):
            date1 = date2 = inps.dset.split('step')[-1]
    print(f'InSAR start / end date: {date1} / {date2}')
    ## read data
    print(f'Read {inps.dset} from file: {inps.file}')
    dis, atr = readfile.read(inps.file, datasetName=inps.dset, box=inps.pix_box)
    if ftype == 'timeseries':
        print(f'Read {date1} from file: {inps.file}')
        dis -= readfile.read(inps.file, datasetName=date1, box=inps.pix_box)[0]
    # convert radians to meters
    if atr['UNIT'] == 'radian':
        phase2range = float(atr['WAVELENGTH']) / (-4 * np.pi)
        dis *= phase2range
    # mask
    if inps.mask_file is not None:
        mask = readfile.read(inps.mask_file, box=inps.pix_box)[0]
        print(f'Set data to NaN for pixels with zero value in file: {inps.mask_file}')
        dis[mask==0] = np.nan
    # read geometry incidence / azimuth angle
    print(f'\nread incidence / azimuth angle from file: {inps.geom_file}')
    inc_angle = readfile.read(inps.geom_file, datasetName='incidenceAngle', box=inps.pix_box)[0]
    az_angle  = readfile.read(inps.geom_file, datasetName='azimuthAngle',   box=inps.pix_box)[0]
    print(f'Mean LOS incidence angle: {np.nanmean(inc_angle):.2f}°')
    print(f'Mean LOS azimuth   angle: {np.nanmean(az_angle):.2f}°')
    # update attributes
    if inps.subset_lat is not None or inps.subset_x is not None:
        atr = attr.update_attribute4subset(atr, inps.pix_box)
    ## create kite container
    create_kite_container(
        dis,
        atr,
        date1,
        date2,
        inc_angle,
        az_angle,
        out_file=inps.outfile,
    )
    return