foxesscloud 2.5.8__py3-none-any.whl → 2.6.0__py3-none-any.whl

foxesscloud/foxesscloud.py

@@ -1,7 +1,7 @@
 ##################################################################################################
 """
 Module:   Fox ESS Cloud
-Updated:  01 October 2024
+Updated:  05 October 2024
 By:       Tony Matthews
 """
 ##################################################################################################
@@ -10,7 +10,7 @@ By:       Tony Matthews
 # ALL RIGHTS ARE RESERVED © Tony Matthews 2023
 ##################################################################################################
 
-version = "1.7.0"
+version = "1.7.2"
 print(f"FoxESS-Cloud version {version}")
 
 debug_setting = 1
@@ -576,10 +576,28 @@ def get_firmware():
 
 battery = None
 battery_settings = None
-residual_handling = 1 # set to 2 if Residual returns current capacity
 
-def get_battery(info=0):
-    global token, device_id, battery, debug_setting, messages
+# 1 = Residual Energy, 2 = Residual Capacity
+residual_handling = 1
+
+# charge rates based on residual_handling
+battery_params = {
+#    cell temp    -5  0   5  10  15  20  25  30  35  40  45  50 55
+#    bms temp      5 10  15  20  25  30  35  40  45  50  55  60 65 
+    1: {'table': [ 0, 2, 10, 15, 25, 50, 50, 50, 50, 50, 30, 20, 0],
+        'step': 5,
+        'offset': 5,
+        'charge_loss': 0.975,
+        'discharge_loss': 0.975},
+    2: {'table': [ 0, 2, 10, 10, 15, 15, 25, 50, 50, 50, 30, 20, 0],
+        'step': 5,
+        'offset': 5,
+        'charge_loss': 1.040,
+        'discharge_loss': 0.975},
+}
+
+def get_battery(info=1):
+    global token, device_id, battery, debug_setting, messages, residual_handling, battery_params
     if get_device() is None:
         return None
     output(f"getting battery", 2)
@@ -593,14 +611,11 @@ def get_battery(info=0):
         errno = response.json().get('errno')
         output(f"** get_battery(), no result data, {errno_message(errno)}")
         return None
+    saved_info = battery['info'] if battery is not None and battery.get('info') is not None else None
     battery = result
-    if battery.get('residual') is not None:
-        battery['residual'] /=1000
-    if residual_handling == 2:
-        capacity = battery.get('residual')
-        soc = battery.get('soc')
-        battery['residual'] = capacity * soc / 100 if capacity is not None and soc is not None else capacity
-    if info == 1:
+    if saved_info is not None:
+        battery['info'] = saved_info
+    elif info == 1:
         response = signed_get(path="/generic/v0/device/battery/list", params=params)
         if response.status_code != 200:
             output(f"** get_battery().info got response code {response.status_code}: {response.reason}")
@@ -611,6 +626,26 @@ def get_battery(info=0):
                 output(f"** get_battery().info, no result data, {errno_message(errno)}")
             else:
                 battery['info'] = result['batteries']
+                if battery['info'][0]['masterVersion'] >= '1.014':
+                    residual_handling = 2
+    if battery.get('residual') is not None:
+        battery['residual'] /= 1000
+    if residual_handling == 2:
+        capacity = battery.get('residual')
+        soc = battery.get('soc')
+        residual = capacity * soc / 100 if capacity is not None and soc is not None else capacity
+    else:
+        residual = battery.get('residual')
+        soc = battery.get('soc')
+        capacity = residual / soc * 100 if residual is not None and soc is not None and soc > 0 else None 
+    battery['capacity'] = round(capacity, 3)
+    battery['residual'] = round(residual, 3)
+    battery['charge_rate'] = 50
+    params = battery_params[residual_handling]
+    battery['charge_loss'] = params['charge_loss']
+    battery['discharge_loss'] = params['discharge_loss']
+    if battery.get('temperature') is not None:
+        battery['charge_rate'] = params['table'][int((battery['temperature'] - params['offset']) / params['step'])]
     return battery
 
 ##################################################################################################
@@ -2543,7 +2578,7 @@ def forecast_value_timed(forecast, today, tomorrow, base_hour, run_time, time_of
     profile = []
     h = base_hour - time_offset
     while h < 0:
-        profile.append(0.0)
+        profile.append(None)
         h += 1 / steps_per_hour
     while h < 48:
         day = today if h < 24 else tomorrow
@@ -2553,10 +2588,10 @@ def forecast_value_timed(forecast, today, tomorrow, base_hour, run_time, time_of
             value = forecast.daily[day]['hourly'].get(int(h % 24))
         else:
             value = forecast.daily[day]['pt30'].get(hours_time(int(h * 2) / 2))
-        profile.append(c_float(value))
+        profile.append(value)
         h += 1 / steps_per_hour
     while len(profile) < run_time:
-        profile.append(0.0)
+        profile.append(None)
     return profile[:run_time]
 
 # build the timed work mode profile from the tariff strategy:
@@ -2592,23 +2627,23 @@ def strategy_timed(timed_mode, base_hour, run_time, min_soc=10, max_soc=100, cur
     return work_mode_timed
 
 # build the timed battery residual from the charge / discharge, work mode and min_soc
+# all power values are as measured at the inverter battery connection
 def battery_timed(work_mode_timed, kwh_current, capacity, time_to_next, kwh_min=None, reserve_drain=None):
     global charge_config, steps_per_hour
-    bat_timed = []
     allowed_drain = charge_config['allowed_drain'] if charge_config.get('allowed_drain') is not None else 4
     bms_loss = (charge_config['bms_power'] / 1000 if charge_config.get('bms_power') is not None else 0.05)
     charge_loss = charge_config['charge_loss']
+    discharge_loss = charge_config['discharge_loss']
     charge_limit = charge_config['charge_limit']
     float_charge = charge_config['float_charge']
     for i in range(0, len(work_mode_timed)):
-        bat_timed.append(kwh_current)
         w = work_mode_timed[i]
         w['kwh'] = kwh_current
         max_now = w['max_soc'] * capacity / 100
         if kwh_current < max_now and w['charge'] > 0.0:
             kwh_current += min([w['charge'], charge_limit - w['pv']]) * charge_loss / steps_per_hour
             kwh_current = max_now if kwh_current > max_now else kwh_current
-        kwh_current += (w['pv'] - w['discharge']) / charge_loss / steps_per_hour
+        kwh_current += (w['pv'] * charge_loss - w['discharge'] / discharge_loss) / steps_per_hour
         if kwh_current > capacity:
             # battery is full
             kwh_current = capacity
@@ -2629,16 +2664,16 @@ def battery_timed(work_mode_timed, kwh_current, capacity, time_to_next, kwh_min=
             reserve_drain = reserve_now
         if kwh_min is not None and kwh_current < kwh_min and i >= time_to_next:       # track minimum without charge
             kwh_min = kwh_current
-    return (bat_timed, kwh_min)
+    return ([work_mode_timed[i]['kwh'] for i in range(0, len(work_mode_timed))], kwh_min)
 
+# use work_mode_timed to generate time periods for the inverter schedule
 def charge_periods(work_mode_timed, base_hour, min_soc, capacity):
     global steps_per_hour
-    output(f"\nConfiguring schedule:",1)
     strategy = []
     start = base_hour
-    periods = []
+    times = []
     for t in range(0, min([24 * steps_per_hour, len(work_mode_timed)])):
-        period = periods[0] if len(periods) > 0 else work_mode_timed[0]
+        period = times[0] if len(times) > 0 else work_mode_timed[0]
         next_period = work_mode_timed[t]
         h = base_hour + t / steps_per_hour
         if h == 24 or period['mode'] != next_period['mode'] or period['hold'] != next_period['hold']:
@@ -2650,15 +2685,19 @@ def charge_periods(work_mode_timed, base_hour, min_soc, capacity):
                 s['max_soc'] = period.get('max_soc')
             elif period['mode'] == 'SelfUse' and period['hold'] == 1:
                 s['min_soc'] = min([int(period['kwh'] / capacity * 100 + 0.5), 100])
-                for p in periods:
+                s['end'] = (start + 1 / steps_per_hour) % 24
+                for p in times:
                     p['min_soc'] = s['min_soc']
             if s['mode'] != 'SelfUse' or s['min_soc'] != min_soc:
                 strategy.append(s)
             start = h
-            periods = []
-        periods.append(work_mode_timed[t])
-    if len(strategy) > 0 and strategy[-1]['min_soc'] != min_soc:
+            times = []
+        times.append(work_mode_timed[t])
+    if len(strategy) == 0:
+        return []
+    if strategy[-1]['min_soc'] != min_soc:
         strategy.append({'start': start %24, 'end': (start + 1 / steps_per_hour) % 24, 'mode': 'SelfUse', 'min_soc': min_soc})
+    output(f"\nConfiguring schedule:",1)
     periods = []
     for s in strategy:
         periods.append(set_period(segment = s, quiet=0))
@@ -2682,9 +2721,9 @@ charge_config = {
     'charge_current': None,           # max battery charge current setting in A
     'discharge_current': None,        # max battery discharge current setting in A
     'export_limit': None,             # maximum export power in kW
-    'discharge_loss': 0.97,           # loss converting battery discharge power to grid power
-    'pv_loss': 0.95,                  # loss converting PV power to battery charge power
-    'grid_loss': 0.975,               # loss converting grid power to battery charge power
+    'dc_ac_loss': 0.970,              # loss converting battery DC power to AC grid power
+    'pv_loss': 0.950,                 # loss converting PV power to DC battery charge power
+    'ac_dc_loss': 0.960,              # loss converting AC grid power to DC battery charge power
     'inverter_power': 101,            # Inverter power consumption in W
     'bms_power': 50,                  # BMS power consumption in W
     'force_charge_power': 5.00,       # charge power in kW when using force charge
@@ -2705,11 +2744,8 @@ charge_config = {
     'special_contingency': 33,        # contingency for special days when consumption might be higher
     'special_days': ['12-25', '12-26', '01-01'],
     'full_charge': None,              # day of month (1-28) to do full charge, or 'daily' or 'Mon', 'Tue' etc
-    'derate_temp': 22,                # battery temperature where cold derating starts to be applied
-    'derate_step': 5,                 # scale for derating factors in C
-    'derating': [24, 15, 10, 2],      # max charge current e.g. 5C step = 22C, 17C, 12C, 7C
     'data_wrap': 6,                   # data items to show per line
-    'target_soc': None,               # set the target SoC for charging
+    'target_soc': None,               # the target SoC for charging (over-rides calculated value)
     'shading': {                      # effect of shading on Solcast / forecast.solar
         'solcast': {'adjust': 0.95, 'am_delay': 1.0, 'am_loss': 0.2, 'pm_delay': 1.0, 'pm_loss': 0.2},
         'solar':   {'adjust': 1.20, 'am_delay': 1.0, 'am_loss': 0.2, 'pm_delay': 1.0, 'pm_loss': 0.2}
@@ -2732,7 +2768,7 @@ charge_needed_app_key = "awcr5gro2v13oher3v1qu6hwnovp28"
 def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=None, show_plot=None, run_after=None, reload=2,
         forecast_times=None, force_charge=0, test_time=None, test_soc=None, test_charge=None, **settings):
     global device, seasonality, solcast_api_key, debug_setting, tariff, solar_arrays, legend_location, time_shift, charge_needed_app_key
-    global timed_strategy, steps_per_hour, base_time, storage
+    global timed_strategy, steps_per_hour, base_time, storage, battery
     print(f"\n---------------- charge_needed ----------------")
     # validate parameters
     args = locals()
@@ -2832,8 +2868,22 @@ def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=N
         output(f"base_hour = {base_hour}, hour_adjustment = {hour_adjustment}, change_hour = {change_hour}, time_change = {time_change}")
         output(f"time_to_start = {time_to_start}, run_time = {run_time}, charge_today = {charge_today}")
         output(f"time_to_next = {time_to_next}, full_charge = {full_charge}")
+    if test_soc is not None:
+        current_soc = test_soc
+        capacity = 14.54
+        residual = test_soc * capacity / 100
+        bat_volt = 317.4
+        bat_power = 0.0
+        temperature = 30
+        bms_charge_current = 25
+        charge_loss = 1.040
+        discharge_loss = 0.974
+        bat_current = 0.0
+        device_power = 6.0
+        device_current = 35
+        model = 'H1-6.0-E'
+    else:
     # get device and battery info from inverter
-    if test_soc is None:
         get_battery()
         if battery is None or battery['status'] != 1:
             output(f"\nBattery status is not available")
@@ -2844,27 +2894,16 @@ def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=N
         bat_current = battery['current']
         temperature = battery['temperature']
         residual = battery['residual']
-        if charge_config.get('capacity') is not None:
-            capacity = charge_config['capacity']
-        elif residual is not None and residual > 0.2 and current_soc is not None and current_soc > 1:
-            capacity = residual * 100 / current_soc
-        else:
+        capacity = charge_config['capacity'] if charge_config.get('capacity') is not None else battery.get('capacity')
+        if capacity is None:
             output(f"Battery capacity could not be estimated. Please add the parameter 'capacity=xx' in kWh")
             return None
+        bms_charge_current = battery.get('charge_rate')
+        charge_loss = battery['charge_loss'] if battery.get('charge_loss') is not None else 1.0
+        discharge_loss = battery['discharge_loss'] if battery.get('discharge_loss') is not None else 1.0
         device_power = device.get('power')
         device_current = device.get('max_charge_current')
         model = device.get('deviceType')
-    else:
-        current_soc = test_soc
-        capacity = 14.54
-        residual = test_soc * capacity / 100
-        bat_volt = 315.4
-        bat_power = 0.0
-        temperature = 30
-        bat_current = 0.0
-        device_power = 6.0
-        device_current = 25
-        model = 'H1-6.0-E'
     min_soc = charge_config['min_soc'] if charge_config['min_soc'] is not None else 10
     max_soc = charge_config['max_soc'] if charge_config['max_soc'] is not None else 100
     volt_curve = charge_config['volt_curve']
@@ -2882,62 +2921,52 @@ def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=N
     output(f"  Min SoC:     {min_soc}% ({reserve:.2f}kWh)")
     output(f"  Current SoC: {current_soc}%")
     output(f"  Max SoC:     {max_soc}% ({capacity * max_soc / 100:.2f}kWh)")
+    output(f"  Charge Rate: {bms_charge_current:.1f}A")
     output(f"  Temperature: {temperature:.1f}°C")
     output(f"  Resistance:  {bat_resistance:.2f} ohms")
     output(f"  Nominal OCV: {bat_ocv:.1f}V at {nominal_soc}% SoC")
-    # charge times are derated based on temperature
+    # charge current may be derated based on temperature
     charge_current = device_current if charge_config['charge_current'] is None else charge_config['charge_current']
-    derate_temp = charge_config['derate_temp']
-    if temperature > 36:
-        output(f"\nHigh battery temperature may affect the charge rate")
-    elif round(temperature, 0) <= derate_temp:
-        output(f"\nLow battery temperature may affect the charge rate")
-        derating = charge_config['derating']
-        derate_step = charge_config['derate_step']
-        i = int((derate_temp - temperature) / (derate_step if derate_step is not None and derate_step > 0 else 1))
-        if derating is not None and type(derating) is list and i < len(derating):
-            derated_current = derating[i]
-            if derated_current < charge_current:
-                output(f"  Charge current reduced from {charge_current:.0f}A to {derated_current:.0f}A" )
-                charge_current = derated_current
-        else:
-            bat_hold = 2
-            output(f"  Full charge set")
+    if charge_current > bms_charge_current:
+        charge_current = bms_charge_current
     # inverter losses
     inverter_power = charge_config['inverter_power'] if charge_config['inverter_power'] is not None else round(device_power, 0) * 25
     operating_loss = inverter_power / 1000
     bms_power = charge_config['bms_power']
     bms_loss = bms_power / 1000
     # work out charge limit, power and losses. Max power going to the battery after ac conversion losses
+    ac_dc_loss = charge_config['ac_dc_loss']
     charge_limit = min([charge_current * (bat_ocv + charge_current * bat_resistance) / 1000, max([6, device_power])])
     if charge_limit < 0.1:
         output(f"** charge_current is too low ({charge_current:.1f}A)")
-    charge_loss = 1.0 - charge_limit * 1000 * bat_resistance / bat_ocv ** 2
     force_charge_power = charge_config['force_charge_power'] if timed_mode > 1 and charge_config.get('force_charge_power') is not None else 100
-    grid_loss = charge_config['grid_loss']
-    charge_power = min([(device_power - operating_loss) * grid_loss, force_charge_power * grid_loss, charge_limit])
+    charge_power = min([(device_power - operating_loss) * ac_dc_loss, force_charge_power * ac_dc_loss, charge_limit])
     float_charge = (charge_config['float_current'] if charge_config.get('float_current') is not None else 4) * bat_ocv / 1000
-    charge_config['charge_loss'] = charge_loss
-    charge_config['charge_limit'] = charge_limit
-    charge_config['charge_power'] = charge_power
-    charge_config['float_charge'] = float_charge
+    pv_loss = charge_config['pv_loss']
     # work out discharge limit = max power coming from the battery before ac conversion losses
-    discharge_loss = charge_config['discharge_loss']
-    discharge_limit = device_power / discharge_loss
+    dc_ac_loss = charge_config['dc_ac_loss']
+    discharge_limit = device_power / dc_ac_loss
     discharge_current = device_current if charge_config['discharge_current'] is None else charge_config['discharge_current']
     discharge_power = discharge_current * bat_ocv / 1000
     discharge_limit = discharge_power if discharge_power < discharge_limit else discharge_limit
     # charging happens if generation exceeds export limit in feedin work mode
     export_power = device_power if charge_config['export_limit'] is None else charge_config['export_limit']
-    export_limit = export_power / discharge_loss
+    export_limit = export_power / dc_ac_loss
     current_mode = get_work_mode()
+    # set parameters for battery_timed()
+    charge_config['charge_limit'] = charge_limit
+    charge_config['charge_power'] = charge_power
+    charge_config['float_charge'] = float_charge
+    charge_config['charge_loss'] = charge_loss
+    charge_config['discharge_loss'] = discharge_loss
+    # display what we have
     output(f"\ncharge_config = {json.dumps(charge_config, indent=2)}", 3)
     output(f"\nDevice Info:")
     output(f"  Model:     {model}")
     output(f"  Rating:    {device_power:.2f}kW")
     output(f"  Export:    {export_power:.2f}kW")
-    output(f"  Charge:    {charge_current:.1f}A, {charge_power:.2f}kW, {charge_loss * 100:.1f}% efficient")
-    output(f"  Discharge: {discharge_current:.1f}A, {discharge_limit:.2f}kW, {discharge_loss * 100:.1f}% efficient")
+    output(f"  Charge:    {charge_current:.1f}A, {charge_power:.2f}kW, {ac_dc_loss * 100:.1f}% efficient")
+    output(f"  Discharge: {discharge_current:.1f}A, {discharge_limit:.2f}kW, {dc_ac_loss * 100:.1f}% efficient")
     output(f"  Inverter:  {inverter_power:.0f}W power consumption")
     output(f"  BMS:       {bms_power:.0f}W power consumption")
     if current_mode is not None:
@@ -3040,8 +3069,8 @@ def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=N
             output(f"\nSettings will not be updated when forecast is not available")
             update_settings = 0
     # produce time lines for charge, discharge and work mode
-    charge_timed = [min([charge_limit, x * charge_config['pv_loss']]) for x in generation_timed]
-    discharge_timed = [min([discharge_limit, x / discharge_loss]) + bms_loss for x in consumption_timed]
+    charge_timed = [min([charge_limit, c_float(x) * pv_loss]) for x in generation_timed]
+    discharge_timed = [min([discharge_limit, c_float(x) / dc_ac_loss]) + bms_loss for x in consumption_timed]
     work_mode_timed = strategy_timed(timed_mode, base_hour, run_time, min_soc=min_soc, max_soc=max_soc, current_mode=current_mode)
     for i in range(0, len(work_mode_timed)):
         # get work mode
@@ -3052,9 +3081,9 @@ def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=N
             discharge_timed[i] = discharge_timed[i] * (1.0 - duration)
             work_mode_timed[i]['charge'] = charge_power * duration
         elif timed_mode > 0 and work_mode == 'ForceDischarge':
-            fdpwr = work_mode_timed[i]['fdpwr'] / discharge_loss / 1000
-            fdpwr = min([discharge_limit, export_limit + discharge_timed[i] * charge_loss, fdpwr])
-            discharge_timed[i] = fdpwr * duration / charge_loss + discharge_timed[i] * (1.0 - duration) - charge_timed[i] * duration
+            fdpwr = work_mode_timed[i]['fdpwr'] / dc_ac_loss / 1000
+            fdpwr = min([discharge_limit, export_limit + discharge_timed[i], fdpwr])
+            discharge_timed[i] = fdpwr * duration + discharge_timed[i] * (1.0 - duration) - charge_timed[i] * duration
         elif bat_hold > 0 and i >= int(time_to_start) and i < int(time_to_end):
             discharge_timed[i] = bms_loss
             if timed_mode > 1:
@@ -3099,13 +3128,13 @@ def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=N
         end_timed = time_to_end
         end_soc = int(end_residual / capacity * 100 + 0.5)
     else:
-        if test_charge is None:
-            output(f"\nCharge needed: {kwh_needed:.2f}kWh")
-            charge_message = "with charge added"
-        output(f"  SoC now:     {current_soc:.0f}% at {hours_time(hour_now)} on {today}")
         # work out time to add kwh_needed to battery
         charge_rate = charge_power * charge_loss
         hours = kwh_needed / charge_rate
+        if test_charge is None:
+            output(f"\nCharge needed: {kwh_needed:.2f}kWh ({hours_time(hours)})")
+            charge_message = "with charge added"
+        output(f"  SoC now:     {current_soc:.0f}% at {hours_time(hour_now)} on {today}")
         # check if charge time exceeded or charge needed exceeds capacity
         hours_to_full = (capacity - start_residual) / charge_rate
         if hours > charge_time:
@@ -3146,7 +3175,7 @@ def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=N
             if i >= start_timed and i < end_timed:
                 work_mode_timed[i]['mode'] = 'ForceCharge'
                 work_mode_timed[i]['charge'] = charge_power * t
-                work_mode_timed[i]['max_soc'] = target_soc if target_soc is not None else 100
+                work_mode_timed[i]['max_soc'] = target_soc if target_soc is not None else max_soc
                 work_mode_timed[i]['discharge'] *= (1-t)
     # rebuild the battery residual with the charge added and min_soc
     (bat_timed, x) = battery_timed(work_mode_timed, kwh_current, capacity, time_to_next)
@@ -3218,7 +3247,6 @@ def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=N
         data['capacity'] = capacity
         data['config'] = charge_config
         data['time'] = time_line
-        data['bat'] = bat_timed
         data['work_mode'] = work_mode_timed
         data['generation'] = generation_timed
         data['consumption'] = consumption_timed
@@ -3233,10 +3261,14 @@ def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=N
 # CHARGE_COMPARE - load saved data and compare with actual
 ##################################################################################################
 
-def charge_compare(save=None, v=None, show_data=1, show_plot=3):
+def charge_compare(save=None, v=None, show_data=1, show_plot=3, d=None):
     global charge_config, storage
+    now = datetime.now() if d is None else datetime.strptime(d, '%Y-%m-%d %H:%M')
+    yesterday = datetime.strftime(datetime.date(now - timedelta(days=1)), '%Y-%m-%d')
     if save is None and charge_config.get('save') is not None:
-        save = charge_config.get('save').replace('###', base_time.replace(' ', 'T'))
+        save = charge_config.get('save').replace('###', yesterday)
+        if not os.path.exists(storage + save):
+            save = None
     if save is None:
         print(f"** charge_compare(): please provide a saved file to load")
         return
@@ -3253,10 +3285,10 @@ def charge_compare(save=None, v=None, show_data=1, show_plot=3):
     steps_per_hour = data.get('steps')
     capacity = data.get('capacity')
     time_line = data.get('time')
-    bat_timed = data.get('bat')
     generation_timed = data.get('generation')
     consumption_timed = data.get('consumption')
     work_mode_timed = data.get('work_mode')
+    bat_timed = data['bat'] if data.get('bat') is not None else [work_mode_timed[t]['kwh'] for t in range(0, len(work_mode_timed))]
     run_time = len(time_line)
     base_hour = int(time_hours(base_time[11:16]))
     start_day = base_time[:10]
@@ -3298,7 +3330,7 @@ def charge_compare(save=None, v=None, show_data=1, show_plot=3):
         s = f"\nBattery Energy kWh:" if show_data == 2 else f"\nBattery SoC:"
         h = base_hour
         t = 0
-        while t < len(time_line) and bat_timed[t] is not None:
+        while t < len(time_line) and bat_timed[t] is not None and plots['SoC'][t] is not None:
             col = h % data_wrap
             s += f"\n  {hours_time(time_line[t])}" if t == 0 or col == 0 else ""
             s += f" {plots['SoC'][t]:5.2f}" if show_data == 2 else f"  {plots['SoC'][t] / capacity * 100:3.0f}%"
@@ -3380,14 +3412,14 @@ battery_info_app_key = "aug938dqt5cbqhvq69ixc4v39q6wtw"
 
 # show information about the current state of the batteries
 def battery_info(log=0, plot=1, count=None, info=1):
-    global debug_setting, battery_info_app_key
+    global debug_setting, battery_info_app_key, residual_handling
     output_spool(battery_info_app_key)
     bat = get_battery(info=info)
     if bat is None:
         output_close()
         return None
     nbat = None
-    if bat.get('info') is not None:
+    if info == 1 and bat.get('info') is not None:
         for b in bat['info']:
             output(f"\nSN {b['masterSN']}, {b['masterBatType']}, Version {b['masterVersion']} (BMS)")
             nbat = 0
@@ -3400,7 +3432,7 @@ def battery_info(log=0, plot=1, count=None, info=1):
     bat_current = bat['current']
     bat_power = bat['power']
     bms_temperature = bat['temperature']
-    capacity = residual / current_soc * 100
+    capacity = bat['capacity']
     cell_volts = get_cell_volts()
     if cell_volts is None:
         output_close()
@@ -3455,6 +3487,7 @@ def battery_info(log=0, plot=1, count=None, info=1):
     output(f"Cell Volts:          {avg(cell_volts):.3f}V average, {max(cell_volts):.3f}V maximum, {min(cell_volts):.3f}V minimum")
     output(f"Cell Imbalance:      {imbalance(cell_volts):.2f}%:")
     output(f"BMS Temperature:     {bms_temperature:.1f}°C")
+    output(f"BMS Charge Rate:     {bat.get('charge_rate'):.1f}A (estimated)")
     output(f"Battery Temperature: {avg(cell_temps):.1f}°C average, {max(cell_temps):.1f}°C maximum, {min(cell_temps):.1f}°C minimum")
     output(f"\nInfo by battery:")
     for i in range(0, nbat):

foxesscloud/openapi.py

@@ -1,7 +1,7 @@
 ##################################################################################################
 """
 Module:   Fox ESS Cloud using Open API
-Updated:  01 October 2024
+Updated:  05 October 2024
 By:       Tony Matthews
 """
 ##################################################################################################
@@ -10,7 +10,7 @@ By:       Tony Matthews
 # ALL RIGHTS ARE RESERVED © Tony Matthews 2024
 ##################################################################################################
 
-version = "2.5.8"
+version = "2.6.0"
 print(f"FoxESS-Cloud Open API version {version}")
 
 debug_setting = 1
@@ -541,10 +541,28 @@ battery = None
 battery_settings = None
 battery_vars = ['SoC', 'invBatVolt', 'invBatCurrent', 'invBatPower', 'batTemperature', 'ResidualEnergy' ]
 battery_data = ['soc', 'volt', 'current', 'power', 'temperature', 'residual']
-residual_handling = 1 # set to 2 if Residual returns current capacity
+
+# 1 = returns Residual Energy. 2 = resturns Residual Capacity
+residual_handling = 1
+
+# charge rates based on residual_handling
+battery_params = {
+#    cell temp    -5  0   5  10  15  20  25  30  35  40  45  50 55
+#    bms temp      5 10  15  20  25  30  35  40  45  50  55  60 65 
+    1: {'table': [ 0, 2, 10, 15, 25, 50, 50, 50, 50, 50, 30, 20, 0],
+        'step': 5,
+        'offset': 5,
+        'charge_loss': 0.975,
+        'discharge_loss': 0.975},
+    2: {'table': [ 0, 2, 10, 10, 15, 15, 25, 50, 50, 50, 30, 20, 0],
+        'step': 5,
+        'offset': 5,
+        'charge_loss': 1.040,
+        'discharge_loss': 0.975},
+}
 
 def get_battery(v = None, info=0):
-    global device_sn, battery, debug_setting, residual_handling
+    global device_sn, battery, debug_setting, residual_handling, battery_params
     if get_device() is None:
         return None
     output(f"getting battery", 2)
@@ -558,10 +576,20 @@ def get_battery(v = None, info=0):
     if residual_handling == 2:
         capacity = battery.get('residual')
         soc = battery.get('soc')
-        battery['residual'] = capacity * soc / 100 if capacity is not None and soc is not None else capacity
-    if info == 1:
-        output(f"** get_battery(): info is not available via Open API")
+        residual = capacity * soc / 100 if capacity is not None and soc is not None else capacity
+    else:
+        residual = battery.get('residual')
+        soc = battery.get('soc')
+        capacity = residual / soc * 100 if residual is not None and soc is not None and soc > 0 else None
+    battery['capacity'] = round(capacity, 3)
+    battery['residual'] = round(residual, 3)
     battery['status'] = 1
+    battery['charge_rate'] = 50
+    params = battery_params[residual_handling]
+    battery['charge_loss'] = params['charge_loss']
+    battery['discharge_loss'] = params['discharge_loss']
+    if battery.get('temperature') is not None:
+        battery['charge_rate'] = params['table'][int((battery['temperature'] - params['offset']) / params['step'])]
     return battery
 
 ##################################################################################################
@@ -2406,7 +2434,7 @@ def forecast_value_timed(forecast, today, tomorrow, base_hour, run_time, time_of
     profile = []
     h = base_hour - time_offset
     while h < 0:
-        profile.append(0.0)
+        profile.append(None)
         h += 1 / steps_per_hour
     while h < 48:
         day = today if h < 24 else tomorrow
@@ -2416,10 +2444,10 @@ def forecast_value_timed(forecast, today, tomorrow, base_hour, run_time, time_of
             value = forecast.daily[day]['hourly'].get(int(h % 24))
         else:
             value = forecast.daily[day]['pt30'].get(hours_time(int(h * 2) / 2))
-        profile.append(c_float(value))
+        profile.append(value)
         h += 1 / steps_per_hour
     while len(profile) < run_time:
-        profile.append(0.0)
+        profile.append(None)
     return profile[:run_time]
 
 # build the timed work mode profile from the tariff strategy:
@@ -2455,23 +2483,23 @@ def strategy_timed(timed_mode, base_hour, run_time, min_soc=10, max_soc=100, cur
     return work_mode_timed
 
 # build the timed battery residual from the charge / discharge, work mode and min_soc
+# note: all power values are as measured at the inverter battery connection
 def battery_timed(work_mode_timed, kwh_current, capacity, time_to_next, kwh_min=None, reserve_drain=None):
     global charge_config, steps_per_hour
-    bat_timed = []
     allowed_drain = charge_config['allowed_drain'] if charge_config.get('allowed_drain') is not None else 4
     bms_loss = (charge_config['bms_power'] / 1000 if charge_config.get('bms_power') is not None else 0.05)
     charge_loss = charge_config['charge_loss']
+    discharge_loss = charge_config['discharge_loss']
     charge_limit = charge_config['charge_limit']
     float_charge = charge_config['float_charge']
     for i in range(0, len(work_mode_timed)):
-        bat_timed.append(kwh_current)
         w = work_mode_timed[i]
         w['kwh'] = kwh_current
         max_now = w['max_soc'] * capacity / 100
         if kwh_current < max_now and w['charge'] > 0.0:
             kwh_current += min([w['charge'], charge_limit - w['pv']]) * charge_loss / steps_per_hour
             kwh_current = max_now if kwh_current > max_now else kwh_current
-        kwh_current += (w['pv'] - w['discharge']) / charge_loss / steps_per_hour
+        kwh_current += (w['pv'] * charge_loss - w['discharge'] / discharge_loss) / steps_per_hour
         if kwh_current > capacity:
             # battery is full
             kwh_current = capacity
@@ -2492,16 +2520,16 @@ def battery_timed(work_mode_timed, kwh_current, capacity, time_to_next, kwh_min=
             reserve_drain = reserve_now
         if kwh_min is not None and kwh_current < kwh_min and i >= time_to_next:       # track minimum without charge
             kwh_min = kwh_current
-    return (bat_timed, kwh_min)
+    return ([work_mode_timed[i]['kwh'] for i in range(0, len(work_mode_timed))], kwh_min)
 
+# use work_mode_timed to generate time periods for the inverter schedule
 def charge_periods(work_mode_timed, base_hour, min_soc, capacity):
     global steps_per_hour
-    output(f"\nConfiguring schedule:",1)
     strategy = []
     start = base_hour
-    periods = []
+    times = []
     for t in range(0, min([24 * steps_per_hour, len(work_mode_timed)])):
-        period = periods[0] if len(periods) > 0 else work_mode_timed[0]
+        period = times[0] if len(times) > 0 else work_mode_timed[0]
         next_period = work_mode_timed[t]
         h = base_hour + t / steps_per_hour
         if h == 24 or period['mode'] != next_period['mode'] or period['hold'] != next_period['hold']:
@@ -2513,15 +2541,19 @@ def charge_periods(work_mode_timed, base_hour, min_soc, capacity):
                 s['max_soc'] = period.get('max_soc')
             elif period['mode'] == 'SelfUse' and period['hold'] == 1:
                 s['min_soc'] = min([int(period['kwh'] / capacity * 100 + 0.5), 100])
-                for p in periods:
+                s['end'] = (start + 1 / steps_per_hour) % 24
+                for p in times:
                     p['min_soc'] = s['min_soc']
             if s['mode'] != 'SelfUse' or s['min_soc'] != min_soc:
                 strategy.append(s)
             start = h
-            periods = []
-        periods.append(work_mode_timed[t])
-    if len(strategy) > 0 and strategy[-1]['min_soc'] != min_soc:
+            times = []
+        times.append(work_mode_timed[t])
+    if len(strategy) == 0:
+        return []
+    if strategy[-1]['min_soc'] != min_soc:
         strategy.append({'start': start %24, 'end': (start + 1 / steps_per_hour) % 24, 'mode': 'SelfUse', 'min_soc': min_soc})
+    output(f"\nConfiguring schedule:",1)
     periods = []
     for s in strategy:
         periods.append(set_period(segment = s, quiet=0))
@@ -2545,9 +2577,9 @@ charge_config = {
     'charge_current': None,           # max battery charge current setting in A
     'discharge_current': None,        # max battery discharge current setting in A
     'export_limit': None,             # maximum export power in kW
-    'discharge_loss': 0.97,           # loss converting battery discharge power to grid power
-    'pv_loss': 0.95,                  # loss converting PV power to battery charge power
-    'grid_loss': 0.975,               # loss converting grid power to battery charge power
+    'dc_ac_loss': 0.97,               # loss converting battery DC power to AC grid power
+    'pv_loss': 0.95,                  # loss converting PV power to DC battery charge power
+    'ac_dc_loss': 0.962,              # loss converting AC grid power to DC battery charge power
     'inverter_power': 101,            # Inverter power consumption in W
     'bms_power': 50,                  # BMS power consumption in W
     'force_charge_power': 5.00,       # charge power in kW when using force charge
@@ -2568,11 +2600,8 @@ charge_config = {
     'special_contingency': 33,        # contingency for special days when consumption might be higher
     'special_days': ['12-25', '12-26', '01-01'],
     'full_charge': None,              # day of month (1-28) to do full charge, or 'daily' or 'Mon', 'Tue' etc
-    'derate_temp': 22,                # battery temperature where cold derating starts to be applied
-    'derate_step': 5,                 # scale for derating factors in C
-    'derating': [24, 15, 10, 2],      # max charge current e.g. 5C step = 22C, 17C, 12C, 7C
     'data_wrap': 6,                   # data items to show per line
-    'target_soc': None,               # set the target SoC for charging
+    'target_soc': None,               # the target SoC for charging (over-rides calculated value)
     'shading': {                      # effect of shading on Solcast / forecast.solar
         'solcast': {'adjust': 0.95, 'am_delay': 1.0, 'am_loss': 0.2, 'pm_delay': 1.0, 'pm_loss': 0.2},
         'solar':   {'adjust': 1.20, 'am_delay': 1.0, 'am_loss': 0.2, 'pm_delay': 1.0, 'pm_loss': 0.2}
@@ -2595,7 +2624,7 @@ charge_needed_app_key = "awcr5gro2v13oher3v1qu6hwnovp28"
 def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=None, show_plot=None, run_after=None, reload=2,
         forecast_times=None, force_charge=0, test_time=None, test_soc=None, test_charge=None, **settings):
     global device, seasonality, solcast_api_key, debug_setting, tariff, solar_arrays, legend_location, time_shift, charge_needed_app_key
-    global timed_strategy, steps_per_hour, base_time, storage
+    global timed_strategy, steps_per_hour, base_time, storage, battery
     print(f"\n---------------- charge_needed ----------------")
     # validate parameters
     args = locals()
@@ -2695,8 +2724,22 @@ def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=N
         output(f"base_hour = {base_hour}, hour_adjustment = {hour_adjustment}, change_hour = {change_hour}, time_change = {time_change}")
         output(f"time_to_start = {time_to_start}, run_time = {run_time}, charge_today = {charge_today}")
         output(f"time_to_next = {time_to_next}, full_charge = {full_charge}")
+    if test_soc is not None:
+        current_soc = test_soc
+        capacity = 14.54
+        residual = test_soc * capacity / 100
+        bat_volt = 317.4
+        bat_power = 0.0
+        temperature = 30
+        bms_charge_current = 25
+        charge_loss = 1.040
+        discharge_loss = 0.974
+        bat_current = 0.0
+        device_power = 6.0
+        device_current = 35
+        model = 'H1-6.0-E'
+    else:
     # get device and battery info from inverter
-    if test_soc is None:
         get_battery()
         if battery is None or battery['status'] != 1:
             output(f"\nBattery status is not available")
@@ -2707,27 +2750,16 @@ def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=N
         bat_current = battery['current']
         temperature = battery['temperature']
         residual = battery['residual']
-        if charge_config.get('capacity') is not None:
-            capacity = charge_config['capacity']
-        elif residual is not None and residual > 0.2 and current_soc is not None and current_soc > 1:
-            capacity = residual * 100 / current_soc
-        else:
+        capacity = charge_config['capacity'] if charge_config.get('capacity') is not None else battery.get('capacity')
+        if capacity is None:
             output(f"Battery capacity could not be estimated. Please add the parameter 'capacity=xx' in kWh")
             return None
+        bms_charge_current = battery.get('charge_rate')
+        charge_loss = battery['charge_loss'] if battery.get('charge_loss') is not None else 1.0
+        discharge_loss = battery['discharge_loss'] if battery.get('discharge_loss') is not None else 1.0
         device_power = device.get('power')
         device_current = device.get('max_charge_current')
         model = device.get('deviceType')
-    else:
-        current_soc = test_soc
-        capacity = 14.54
-        residual = test_soc * capacity / 100
-        bat_volt = 315.4
-        bat_power = 0.0
-        temperature = 30
-        bat_current = 0.0
-        device_power = 6.0
-        device_current = 25
-        model = 'H1-6.0-E'
     min_soc = charge_config['min_soc'] if charge_config['min_soc'] is not None else 10
     max_soc = charge_config['max_soc'] if charge_config['max_soc'] is not None else 100
     volt_curve = charge_config['volt_curve']
@@ -2746,61 +2778,51 @@ def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=N
     output(f"  Current SoC: {current_soc}%")
     output(f"  Max SoC:     {max_soc}% ({capacity * max_soc / 100:.2f}kWh)")
     output(f"  Temperature: {temperature:.1f}°C")
+    output(f"  Charge Rate: {bms_charge_current:.1f}A")
     output(f"  Resistance:  {bat_resistance:.2f} ohms")
     output(f"  Nominal OCV: {bat_ocv:.1f}V at {nominal_soc}% SoC")
-    # charge times are derated based on temperature
+    # charge current may be derated based on temperature
     charge_current = device_current if charge_config['charge_current'] is None else charge_config['charge_current']
-    derate_temp = charge_config['derate_temp']
-    if temperature > 36:
-        output(f"\nHigh battery temperature may affect the charge rate")
-    elif round(temperature, 0) <= derate_temp:
-        output(f"\nLow battery temperature may affect the charge rate")
-        derating = charge_config['derating']
-        derate_step = charge_config['derate_step']
-        i = int((derate_temp - temperature) / (derate_step if derate_step is not None and derate_step > 0 else 1))
-        if derating is not None and type(derating) is list and i < len(derating):
-            derated_current = derating[i]
-            if derated_current < charge_current:
-                output(f"  Charge current reduced from {charge_current:.0f}A to {derated_current:.0f}A" )
-                charge_current = derated_current
-        else:
-            bat_hold = 2
-            output(f"  Full charge set")
+    if charge_current > bms_charge_current:
+        charge_current = bms_charge_current
     # inverter losses
     inverter_power = charge_config['inverter_power'] if charge_config['inverter_power'] is not None else round(device_power, 0) * 25
     operating_loss = inverter_power / 1000
     bms_power = charge_config['bms_power']
     bms_loss = bms_power / 1000
     # work out charge limit, power and losses. Max power going to the battery after ac conversion losses
+    ac_dc_loss = charge_config['ac_dc_loss']
     charge_limit = min([charge_current * (bat_ocv + charge_current * bat_resistance) / 1000, max([6, device_power])])
     if charge_limit < 0.1:
         output(f"** charge_current is too low ({charge_current:.1f}A)")
-    charge_loss = 1.0 - charge_limit * 1000 * bat_resistance / bat_ocv ** 2
     force_charge_power = charge_config['force_charge_power'] if timed_mode > 1 and charge_config.get('force_charge_power') is not None else 100
-    grid_loss = charge_config['grid_loss']
-    charge_power = min([(device_power - operating_loss) * grid_loss, force_charge_power * grid_loss, charge_limit])
+    charge_power = min([(device_power - operating_loss) * ac_dc_loss, force_charge_power * ac_dc_loss, charge_limit])
     float_charge = (charge_config['float_current'] if charge_config.get('float_current') is not None else 4) * bat_ocv / 1000
-    charge_config['charge_loss'] = charge_loss
-    charge_config['charge_limit'] = charge_limit
-    charge_config['charge_power'] = charge_power
-    charge_config['float_charge'] = float_charge
+    pv_loss = charge_config['pv_loss']
     # work out discharge limit = max power coming from the battery before ac conversion losses
-    discharge_loss = charge_config['discharge_loss']
-    discharge_limit = device_power / discharge_loss
+    dc_ac_loss = charge_config['dc_ac_loss']
+    discharge_limit = device_power / dc_ac_loss
     discharge_current = device_current if charge_config['discharge_current'] is None else charge_config['discharge_current']
     discharge_power = discharge_current * bat_ocv / 1000
     discharge_limit = discharge_power if discharge_power < discharge_limit else discharge_limit
     # charging happens if generation exceeds export limit in feedin work mode
     export_power = device_power if charge_config['export_limit'] is None else charge_config['export_limit']
-    export_limit = export_power / discharge_loss
+    export_limit = export_power / dc_ac_loss
     current_mode = get_work_mode()
+    # set parameters for battery_timed()
+    charge_config['charge_limit'] = charge_limit
+    charge_config['charge_power'] = charge_power
+    charge_config['float_charge'] = float_charge
+    charge_config['charge_loss'] = charge_loss
+    charge_config['discharge_loss'] = discharge_loss
+    # display what we have
     output(f"\ncharge_config = {json.dumps(charge_config, indent=2)}", 3)
     output(f"\nDevice Info:")
     output(f"  Model:     {model}")
     output(f"  Rating:    {device_power:.2f}kW")
     output(f"  Export:    {export_power:.2f}kW")
-    output(f"  Charge:    {charge_current:.1f}A, {charge_power:.2f}kW, {charge_loss * 100:.1f}% efficient")
-    output(f"  Discharge: {discharge_current:.1f}A, {discharge_limit:.2f}kW, {discharge_loss * 100:.1f}% efficient")
+    output(f"  Charge:    {charge_current:.1f}A, {charge_power:.2f}kW, {ac_dc_loss * 100:.1f}% efficient")
+    output(f"  Discharge: {discharge_current:.1f}A, {discharge_limit:.2f}kW, {dc_ac_loss * 100:.1f}% efficient")
     output(f"  Inverter:  {inverter_power:.0f}W power consumption")
     output(f"  BMS:       {bms_power:.0f}W power consumption")
     if current_mode is not None:
@@ -2902,8 +2924,8 @@ def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=N
             output(f"\nSettings will not be updated when forecast is not available")
             update_settings = 0
     # produce time lines for charge, discharge and work mode
-    charge_timed = [min([charge_limit, x * charge_config['pv_loss']]) for x in generation_timed]
-    discharge_timed = [min([discharge_limit, x / discharge_loss]) + bms_loss for x in consumption_timed]
+    charge_timed = [min([charge_limit, c_float(x) * pv_loss]) for x in generation_timed]
+    discharge_timed = [min([discharge_limit, c_float(x) / dc_ac_loss]) + bms_loss for x in consumption_timed]
     work_mode_timed = strategy_timed(timed_mode, base_hour, run_time, min_soc=min_soc, max_soc=max_soc, current_mode=current_mode)
     for i in range(0, len(work_mode_timed)):
         # get work mode
@@ -2914,9 +2936,9 @@ def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=N
             discharge_timed[i] = discharge_timed[i] * (1.0 - duration)
             work_mode_timed[i]['charge'] = charge_power * duration
         elif timed_mode > 0 and work_mode == 'ForceDischarge':
-            fdpwr = work_mode_timed[i]['fdpwr'] / discharge_loss / 1000
-            fdpwr = min([discharge_limit, export_limit + discharge_timed[i] * charge_loss, fdpwr])
-            discharge_timed[i] = fdpwr * duration / charge_loss + discharge_timed[i] * (1.0 - duration) - charge_timed[i] * duration
+            fdpwr = work_mode_timed[i]['fdpwr'] / dc_ac_loss / 1000
+            fdpwr = min([discharge_limit, export_limit + discharge_timed[i], fdpwr])
+            discharge_timed[i] = fdpwr * duration + discharge_timed[i] * (1.0 - duration) - charge_timed[i] * duration
         elif bat_hold > 0 and i >= int(time_to_start) and i < int(time_to_end):
             discharge_timed[i] = bms_loss
             if timed_mode > 1:
@@ -2961,13 +2983,13 @@ def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=N
         end_timed = time_to_end
         end_soc = int(end_residual / capacity * 100 + 0.5)
     else:
-        if test_charge is None:
-            output(f"\nCharge needed: {kwh_needed:.2f}kWh:")
-            charge_message = "with charge added"
-        output(f"  SoC now:     {current_soc:.0f}% at {hours_time(hour_now)} on {today}")
         # work out time to add kwh_needed to battery
         charge_rate = charge_power * charge_loss
         hours = kwh_needed / charge_rate
+        if test_charge is None:
+            output(f"\nCharge needed: {kwh_needed:.2f}kWh ({hours_time(hours)})")
+            charge_message = "with charge added"
+        output(f"  SoC now:     {current_soc:.0f}% at {hours_time(hour_now)} on {today}")
         # check if charge time exceeded or charge needed exceeds capacity
         hours_to_full = (capacity - start_residual) / charge_rate
         if hours > charge_time:
@@ -3008,7 +3030,7 @@ def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=N
             if i >= start_timed and i < end_timed:
                 work_mode_timed[i]['mode'] = 'ForceCharge'
                 work_mode_timed[i]['charge'] = charge_power * t
-                work_mode_timed[i]['max_soc'] = target_soc if target_soc is not None else 100
+                work_mode_timed[i]['max_soc'] = target_soc if target_soc is not None else max_soc
                 work_mode_timed[i]['discharge'] *= (1-t)
     # rebuild the battery residual with any charge added and min_soc
     (bat_timed, x) = battery_timed(work_mode_timed, kwh_current, capacity, time_to_next)
@@ -3080,7 +3102,6 @@ def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=N
         data['capacity'] = capacity
         data['config'] = charge_config
         data['time'] = time_line
-        data['bat'] = bat_timed
         data['work_mode'] = work_mode_timed
         data['generation'] = generation_timed
         data['consumption'] = consumption_timed
@@ -3096,8 +3117,12 @@ def charge_needed(forecast=None, update_settings=0, timed_mode=None, show_data=N
 
 def charge_compare(save=None, v=None, show_data=1, show_plot=3):
     global charge_config, storage
+    now = datetime.now() if d is None else datetime.strptime(d, '%Y-%m-%d %H:%M')
+    yesterday = datetime.strftime(datetime.date(now - timedelta(days=1)), '%Y-%m-%d')
     if save is None and charge_config.get('save') is not None:
-        save = charge_config.get('save').replace('###', base_time.replace(' ', 'T'))
+        save = charge_config.get('save').replace('###', yesterday)
+        if not os.path.exists(storage + save):
+            save = None
     if save is None:
         print(f"** charge_compare(): please provide a saved file to load")
         return
@@ -3114,10 +3139,10 @@ def charge_compare(save=None, v=None, show_data=1, show_plot=3):
     steps_per_hour = data.get('steps')
     capacity = data.get('capacity')
     time_line = data.get('time')
-    bat_timed = data.get('bat')
     generation_timed = data.get('generation')
     consumption_timed = data.get('consumption')
     work_mode_timed = data.get('work_mode')
+    bat_timed = data['bat'] if data.get('bat') is not None else [work_mode_timed[t]['kwh'] for t in range(0, len(work_mode_timed))]
     run_time = len(time_line)
     base_hour = int(time_hours(base_time[11:16]))
     start_day = base_time[:10]
@@ -3159,7 +3184,7 @@ def charge_compare(save=None, v=None, show_data=1, show_plot=3):
         s = f"\nBattery Energy kWh:" if show_data == 2 else f"\nBattery SoC:"
         h = base_hour
         t = 0
-        while t < len(time_line) and bat_timed[t] is not None:
+        while t < len(time_line) and bat_timed[t] is not None and plots['SoC'][t] is not None:
             col = h % data_wrap
             s += f"\n  {hours_time(time_line[t])}" if t == 0 or col == 0 else ""
             s += f" {plots['SoC'][t]:5.2f}" if show_data == 2 else f"  {plots['SoC'][t] / capacity * 100:3.0f}%"
@@ -3240,14 +3265,14 @@ battery_info_app_key = "aug938dqt5cbqhvq69ixc4v39q6wtw"
 
 # show information about the current state of the batteries
 def battery_info(log=0, plot=1, count=None, info=1):
-    global debug_setting, battery_info_app_key
+    global debug_setting, battery_info_app_key, residual_handling
     output_spool(battery_info_app_key)
     bat = get_battery(info=info)
     if bat is None:
         output_close()
         return None
     nbat = None
-    if bat.get('info') is not None:
+    if info == 1 and bat.get('info') is not None:
         for b in bat['info']:
             output(f"\nSN {b['masterSN']}, {b['masterBatType']}, Version {b['masterVersion']} (BMS)")
             nbat = 0
@@ -3260,7 +3285,7 @@ def battery_info(log=0, plot=1, count=None, info=1):
     bat_current = bat['current']
     bat_power = bat['power']
     bms_temperature = bat['temperature']
-    capacity = residual / current_soc * 100
+    capacity = bat['capacity']
     cell_volts = get_cell_volts()
     if cell_volts is None:
         output_close()
@@ -3315,6 +3340,7 @@ def battery_info(log=0, plot=1, count=None, info=1):
     output(f"Cell Volts:          {avg(cell_volts):.3f}V average, {max(cell_volts):.3f}V maximum, {min(cell_volts):.3f}V minimum")
     output(f"Cell Imbalance:      {imbalance(cell_volts):.2f}%:")
     output(f"BMS Temperature:     {bms_temperature:.1f}°C")
+    output(f"BMS Charge Rate:     {bat.get('charge_rate'):.1f}A (estimated)")
     output(f"Battery Temperature: {avg(cell_temps):.1f}°C average, {max(cell_temps):.1f}°C maximum, {min(cell_temps):.1f}°C minimum")
     output(f"\nInfo by battery:")
     for i in range(0, nbat):

{foxesscloud-2.5.8.dist-info → foxesscloud-2.6.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: foxesscloud
-Version: 2.5.8
+Version: 2.6.0
 Summary: library for accessing Fox ESS cloud data using Open API
 Author-email: Tony Matthews <tony@quasair.co.uk>
 Project-URL: Homepage, https://github.com/TonyM1958/FoxESS-Cloud
@@ -116,7 +116,12 @@ Each of these calls will return a dictionary or list containing the relevant inf
 
 get_generation() will return the latest generation information for the device. The results are also stored in f.device as 'generationToday', 'generationMonth' and 'generationTotal'.
 
-get_battery() returns the current battery status, including 'soc', 'volt', 'current', 'power', 'temperature' and 'residual'. The result also updates f.battery.
+get_battery() returns the current battery status, including 'soc', 'volt', 'current', 'power', 'temperature' and 'residual'. The result also updates f.battery. Additional battery attributes include:
++ 'info': a list of BMS and battery serial numbers and firmware versions
++ 'capacity': the estimated battery capacity, derrived from 'residual' and 'soc'
++ 'charge_rate': the estimated BMS charge rate available, based on the current 'temperature' of the BMS
++ 'charge_loss': the ratio of the kWh added to the battery for each kWh applied during charging
++ 'discharge_loss': the ratio of the kWh available for each kWh removed from the battery during during discharging
 
 get_settings() will return the battery settings and is equivalent to get_charge() and get_min(). The results are stored in f.battery_settings. The settings include minSoc, minSocOnGrid, enable charge from grid and the charge times.
 
@@ -363,44 +368,41 @@ Given the data available, the modelling works as follows:
 
 The following parameters and default values are used to configure charge_needed and may be updated if required using name=value:
 ```
-contingency: [20,10,5,10]     # % of consumption. Single or [winter, spring, summer, autumn] values
-capacity: None                # Battery capacity in kWh (over-rides generated value if set)
-charge_current: None          # max battery charge current setting in A. None uses a value derrived from the inverter model
-discharge_current: None       # max battery discharge current setting in A. None uses a value derrived from the inverter model
-export_limit: None            # maximum export power in kW. None uses the inverter power rating
-discharge_loss: 0.98          # loss converting battery discharge power to grid power
-pv_loss: 0.95                 # loss converting PV power to battery charge power
-grid_loss: 0.975              # loss converting grid power to battery charge power
-inverter_power: None          # inverter power consumption in W (dynamically set)
-bms_power: 50                 # BMS power consumption in W
-force_charge_power: 5.00      # power used when Force Charge is scheduled
-allowed_drain: 4,             # % tolerance below min_soc before float charge starts
-float_current: 4,             # BMS float charge current in A
-bat_resistance: 0.070         # internal resistance of a battery in ohms
-volt_curve: lifepo4_curve     # battery OCV from 0% to 100% SoC
-nominal_soc: 55               # SoC for nominal open circuit voltage
-generation_days: 3            # number of days to use for average generation (1-7)
-consumption_days: 3           # number of days to use for average consumption (1-7)
-consumption_span: 'week'      # 'week' = last 7 days or 'weekday' = last 7 weekdays e.g. Saturdays
-use_today: 21.0               # hour when today's generation and consumption data will be used
-min_hours: 0.25               # minimum charge time to set (in decimal hours)
-min_kwh: 0.5                  # minimum charge to add in kwh
-solcast_adjust: 100           # % adjustment to make to Solcast forecast
-solar_adjust:  100            # % adjustment to make to Solar forecast
-forecast_selection: 1         # 1 = only update charge times if forecast is available, 0 = use best available data. Default is 1.
-annual_consumption: None      # optional annual consumption in kWh. If set, this replaces consumption history
-timed_mode: 0                 # 0 = None, 1 = use timed work mode, 2 = strategy mode
-special_contingency: 30       # contingency for special days when consumption might be higher
+contingency: [20,10,5,10]      # % of consumption. Single or [winter, spring, summer, autumn] values
+capacity: None                 # Battery capacity in kWh (over-rides generated value if set)
+charge_current: None           # max battery charge current setting in A. None uses a value derrived from the inverter model
+discharge_current: None        # max battery discharge current setting in A. None uses a value derrived from the inverter model
+export_limit: None             # maximum export power in kW. None uses the inverter power rating
+dc_ac_loss: 0.970              # loss converting battery DC power to AC grid power
+pv_loss: 0.950                 # loss converting PV power to DC battery charge power
+ac_dc_loss: 0.960              # loss converting AC grid power to DC battery charge power
+inverter_power: None           # inverter power consumption in W (dynamically set)
+bms_power: 50                  # BMS power consumption in W
+force_charge_power: 5.00       # power used when Force Charge is scheduled
+allowed_drain: 4,              # % tolerance below min_soc before float charge starts
+float_current: 4,              # BMS float charge current in A
+bat_resistance: 0.070          # internal resistance of a battery in ohms
+volt_curve: lifepo4_curve      # battery OCV from 0% to 100% SoC
+nominal_soc: 55                # SoC for nominal open circuit voltage
+generation_days: 3             # number of days to use for average generation (1-7)
+consumption_days: 3            # number of days to use for average consumption (1-7)
+consumption_span: 'week'       # 'week' = last 7 days or 'weekday' = last 7 weekdays e.g. Saturdays
+use_today: 21.0                # hour when today's generation and consumption data will be used
+min_hours: 0.25                # minimum charge time to set (in decimal hours)
+min_kwh: 0.5                   # minimum charge to add in kwh
+solcast_adjust: 100            # % adjustment to make to Solcast forecast
+solar_adjust:  100             # % adjustment to make to Solar forecast
+forecast_selection: 1          # 1 = only update charge times if forecast is available, 0 = use best available data. Default is 1.
+annual_consumption: None       # optional annual consumption in kWh. If set, this replaces consumption history
+timed_mode: 0                  # 0 = None, 1 = use timed work mode, 2 = strategy mode
+special_contingency: 30        # contingency for special days when consumption might be higher
 special_days: ['12-25', '12-26', '01-01']
-full_charge: None             # day of month (1-28) to do full charge or 'daily' or day of week: 'Mon', 'Tue' etc
-derate_temp: 21               # battery temperature in C when derating charge current is applied
-derate_step: 5                # step size for derating e.g. 21, 16, 11
-derating: [24, 15, 10, 2]     # derated charge current for each temperature step e.g. 21C, 16C, 11C, 6C
-force: 1                      # 1 = disable strategy periods when setting charge. 0 = fail if strategy period has been set.
-data_wrap: 6                  # data items to show per line
-target_soc: None              # target soc for charging
-shading: {}                   # effect of shading on Solcast / Solar (see below)
-save: 'charge_needed.txt'     # where to save calculation data for charge_compare(). '###' gets replaced with todays date.
+full_charge: None              # day of month (1-28) to do full charge or 'daily' or day of week: 'Mon', 'Tue' etc
+force: 1                       # 1 = disable strategy periods when setting charge. 0 = fail if strategy period has been set.
+data_wrap: 6                   # data items to show per line
+target_soc: None               # target soc for charging (over-rides calculated value)
+shading: {}                    # effect of shading on Solcast / Solar (see below)
+save: 'charge_needed.txt'      # where to save calculation data for charge_compare(). '###' gets replaced with todays date.
 ```
 
 These values are stored / available in f.charge_config.
@@ -784,6 +786,14 @@ This setting can be:
 
 # Version Info
 
+2.6.0<br>
+Rework charge de-rating with temperature, losses and other info provided by get_battery() to take new BMS behaviour into account.
+
+2.5.9<br>
+Change loss parameters to separate AC/DC, DC/AC conversion losses and battery charge / discharge losses.
+Update charge calibration for new BMS firmware.
+Increase de-rating temperature from 21C to 28C for new BMS firmware.
+
 2.5.8<br>
 Fix incorrect charging setup when force_charge=1.
 Rework charge_periods() to consolidate charge periods to reduce number of time segments when timed_mode=2.

foxesscloud-2.6.0.dist-info/RECORD

@@ -0,0 +1,7 @@
+foxesscloud/foxesscloud.py,sha256=2Nsr3hmtq6jcFLslcOrkBEvFKgjPiZdpsQO0O85S928,213371
+foxesscloud/openapi.py,sha256=sJopUqWbk8OrO9MXFAnEQvs14G6t4GWtX05lL-tXymI,206640
+foxesscloud-2.6.0.dist-info/LICENCE,sha256=-3xv8CElCJV8Bc8PbAsg3iyxMpAK8MoJneM3rXigxqI,1074
+foxesscloud-2.6.0.dist-info/METADATA,sha256=GUbw1ZyKzfPaxtmBxGnuM__RKDZqAJpeejiGb4KZfSM,56993
+foxesscloud-2.6.0.dist-info/WHEEL,sha256=oiQVh_5PnQM0E3gPdiz09WCNmwiHDMaGer_elqB3coM,92
+foxesscloud-2.6.0.dist-info/top_level.txt,sha256=IWOrKSNZCLU6IDXSX_b4_bqCfbZoWAT4CC0w0Lg7PuU,12
+foxesscloud-2.6.0.dist-info/RECORD,,

foxesscloud-2.5.8.dist-info/RECORD

@@ -1,7 +0,0 @@
-foxesscloud/foxesscloud.py,sha256=9WivKa8ysTZ52aTgPmTxJYDAnAz5TobibgvpJkrK_KM,211855
-foxesscloud/openapi.py,sha256=SCr8VOz1aP0Nk5vlUi3-kJ-AMNpW0STlVWQRbbABjBE,205502
-foxesscloud-2.5.8.dist-info/LICENCE,sha256=-3xv8CElCJV8Bc8PbAsg3iyxMpAK8MoJneM3rXigxqI,1074
-foxesscloud-2.5.8.dist-info/METADATA,sha256=Aon8oxksMZP4Z9NgWKOwsz7Dpo8raj8A2cIjakNL_xM,56304
-foxesscloud-2.5.8.dist-info/WHEEL,sha256=oiQVh_5PnQM0E3gPdiz09WCNmwiHDMaGer_elqB3coM,92
-foxesscloud-2.5.8.dist-info/top_level.txt,sha256=IWOrKSNZCLU6IDXSX_b4_bqCfbZoWAT4CC0w0Lg7PuU,12
-foxesscloud-2.5.8.dist-info/RECORD,,