0

I hope you all are fine.

I got a function (Y axis) that gets values according to a datetime pandas array (X axis) as shown below.

Effective Irradiance PVLIB-Python

I need to get all the maximum values. Those n maximum values depends on the number of days that are between a start and end pd.timestamp. Then, I also need to save those values in an array; something like this:

max_values = [
    (A_value_1, B_value_1),
    (A_value_2, B_value_1),
        ...,
    (A_value_n, A_value_n)]

Note: The B_value_#is the same that I am asking for but with another function (A_values: Irradiance and B_values: Cell temperature, just a contextualization for the PVLIB-Python community.

For those who are familiar with PVLIB-Python, I am trying to get the I-V curve by using the pvsystem.calcparams_desoto function (this is because I am using a CEC module and it doesn't have the parameters required by the SAPM function. Part of this code is snipped below:

    IL, I0, Rs, Rsh, nNsVth = pvsystem.calcparams_desoto(
    effective_irrad_calc,
    pvtemps['temp_cell'],
    module['alpha_sc'],
    module['a_ref'],
    module['I_L_ref'],
    module['I_o_ref'],
    module['R_sh_ref'],
    module['R_s'],
    EgRef=1.121,
    dEgdT=-0.0002677
)

curve_info = pvsystem.singlediode(
    photocurrent=IL,
    saturation_current=I0,
    resistance_series=Rs,
    resistance_shunt=Rsh,
    nNsVth=nNsVth,
    ivcurve_pnts=100,
    method='lambertw'
)

# plot the calculated curves:
plt.figure()
for i, case in conditions.iterrows():
    label = (
        "$G_{eff}$ " + f"{case['Geff']} $W/m^2$\n"
        "$T_{cell}$ " + f"{case['Tcell']} $C$"
    )
    plt.plot(curve_info['v'][i], curve_info['i'][i], label=label)
    v_mp = curve_info['v_mp'][i]
    i_mp = curve_info['i_mp'][i]
    # mark the MPP
    plt.plot([v_mp], [i_mp], ls='', marker='o', c='k')

plt.legend(loc=(1.0, 0))
plt.xlabel('Module voltage [V]')
plt.ylabel('Module current [A]')
plt.title(parameters['Name'])
plt.show()
plt.gcf().set_tight_layout(True)

Note: This code snipped was taken at: https://github.com/pvlib/pvlib-python/blob/master/docs/examples/plot_singlediode.py

Thanks in advance. :)

Nelson Salazar
  • 21
  • 6

1 Answers1

1

First just making a small example:

df = pd.DataFrame(index=pd.date_range('2020-01-01', freq='1h', periods=24*5))
df['irradiance'] = np.clip(np.cos(np.linspace(np.pi, 2*np.pi*5+np.pi, 24*5)), 
                           a_min=0, a_max=None)

Then using pandas' resample function with a daily frequency and the max function, would give you the daily value for the peak:

max_value = df['irradiance'].resample('1d').max()

Similarly, the location of the max value can be found:

max_location = df['irradiance'].resample('1d').apply(lambda s: s.idxmax())

Now knowing the location/index of the max value you can apply this to another column of your DataFrame.

Adam R. Jensen
  • 617
  • 5
  • 12