def dqdsst_adjustment(file: str, group_files: str, cdo, options, verbose,
    **kwargs):
    """
        Compute the kinematic surface net heat flux sensitivity to the sea surface temperature: dQdSST
        Q_model ~ Q + dQdSST * (T_model - SST)
        dQdSST = - 4 * eps * stef * T^3 - rho_atm * Cp * CH * U - rho_atm * CE * L * U * 2353 * ln(10 * q_s / T^2)
        SST: sea surface temperature
        t_air: sea surface atmospheric temperature
        rho_air: atmospheric density
        u_air: wind speed
        humidity: sea level specific humidity
        dQdSST: kinematic surface net heat flux sensitivity to the sea surface temperature (Watts meter-2 Celsius-1)
        :param preprocessor: PreProcessor object, contains a lot of necessary adjustments
        :param file: file which contains all necessary parameters to calculate dQdSST
        :param group_files: list of all files which contain necessary parameters to calculate dQdSST
        :param flags: any extra info necessary
    """
    if verbose:
        print('Calculating dQdSST')
    Cp = str(1004.8)
    Ch = str(0.00066)
    Ce = str(0.00115)
    stef = str(5.6697e-08)
    sstk_formula = '_sst = SST + 273.15;'
    l = '(2.008E6 - 2.3E3 * t_air)'
    q1 = '(-4.0 * ' + stef + '* _sst * _sst * _sst)'
    q2 = '(-rho_air * ' + Cp + ' * ' + Ch + '* u_air)'
    dqsdt = (
        '(5417.9827238149900949792936444282531738281250 * humidity / (_sst * _sst))'
        )
    q3 = '(-rho_air * ' + Ce + ' * ' + l + ' * u_air * ' + dqsdt + ' )'
    dqdsst_formula = ("'" + sstk_formula + 'dQdSST=' + q1 + ' + ' + q2 +
        ' + ' + q3 + "'")
    temp = cdo.setattribute(
        "dQdSST@long_name='surface net heat flux sensitivity to SST',dQdSST@unit='Watts meter-2 Celsius-1',dQdSST@files='"
         + group_files + "'", input='-delname,sat,rho_atm,U,qsea -aexpr,' +
        dqdsst_formula + ' ' + file, options=options)
    cdo.copy(input=temp, output=file, options=options)
    if verbose:
        print('dQdSST has been added.')
