Skip to content

Commit

Permalink
Misc: pandas future warning
Browse files Browse the repository at this point in the history
  • Loading branch information
@ebranlard
ebranlard committed Aug 8, 2024
1 parent 7d39f05 commit e4a2413
Show file tree
Hide file tree
Showing 4 changed files with 56 additions and 9 deletions.
51 changes: 49 additions & 2 deletions pydatview/fast/postpro.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -689,6 +689,52 @@ def spanwiseColEDTwr(Cols):
def spanwiseColAD(Cols):
""" Return column info, available columns and indices that contain AD spanwise data"""
ADSpanMap=dict()
# From AeroDyn_AllBldNd: TODO Use it directly..
# "ALPHA ","AXIND ","AXIND_QS ","BEM_CT_QS","BEM_F_QS ","BEM_KP_QS","BEM_K_QS ","CD ", &
# "CD_QS ","CHI ","CL ","CLRNC ","CL_QS ","CM ","CMA ","CM_QS ", &
# "CN ","CPMIN ","CT ","CURVE ","CX ","CXA ","CY ","DEBUG1 ", &
# "DEBUG2 ","DEBUG3 ","DYNP ","FBN ","FBS ","FBT ","FBXA ","FBXI ", &
# "FBXL ","FBXP ","FBYA ","FBYI ","FBYL ","FBYP ","FBZA ","FBZI ", &
# "FBZL ","FBZP ","FD ","FL ","FN ","FT ","FX ","FXA ", &
# "FXI ","FXL ","FXP ","FY ","FYI ","FYL ","FYP ","FZI ", &
# "FZL ","FZP ","GAM ","GEOMPHI ","M ","MBN ","MBS ","MBT ", &
# "MBXA ","MBXI ","MBXL ","MBXP ","MBYA ","MBYI ","MBYL ","MBYP ", &
# "MBZA ","MBZI ","MBZL ","MBZP ","MM ","MXI ","MXL ","MXP ", &
# "MYI ","MYL ","MYP ","MZA ","MZI ","MZL ","MZP ","PHI ", &
# "RE ","SGCAV ","SIGCR ","STVX ","STVXA ","STVXI ","STVXL ","STVXP ", &
# "STVY ","STVYA ","STVYI ","STVYL ","STVYP ","STVZ ","STVZA ","STVZI ", &
# "STVZL ","STVZP ","THETA ","TNIND ","TNIND_QS ","TOE ","UA_FLAG ","UA_X1 ", &
# "UA_X2 ","UA_X3 ","UA_X4 ","UA_X5 ","UIN ","UIR ","UIT ","VDISX ", &
# "VDISXA ","VDISXI ","VDISXL ","VDISXP ","VDISY ","VDISYA ","VDISYI ","VDISYL ", &
# "VDISYP ","VDISZ ","VDISZA ","VDISZI ","VDISZL ","VDISZP ","VINDX ","VINDXA ", &
# "VINDXI ","VINDXL ","VINDXP ","VINDY ","VINDYA ","VINDYI ","VINDYL ","VINDYP ", &
# "VINDZA ","VINDZI ","VINDZL ","VINDZP ","VREL ","VUNDX ","VUNDXA ","VUNDXI ", &
# "VUNDXL ","VUNDXP ","VUNDY ","VUNDYA ","VUNDYI ","VUNDYL ","VUNDYP ","VUNDZ ", &
# "VUNDZA ","VUNDZI ","VUNDZL ","VUNDZP ","VX ","VY "/)
# CHARACTER(ChanLen), PARAMETER :: ParamUnitsAry(166) = (/ character(ChanLen) :: & ! This lists the units corresponding to the allowed parameters
# "(deg) ","(-) ","(-) ","(-) ","(-) ","(-) ","(-) ","(-) ", &
# "(-) ","(deg) ","(-) ","(m) ","(-) ","(-) ","(-) ","(-) ", &
# "(-) ","(-) ","(-) ","(deg) ","(-) ","(-) ","(-) ","(-) ", &
# "(-) ","(-) ","(Pa) ","(N/m) ","(N/m) ","(N/m) ","(N/m) ","(N/m) ", &
# "(N/m) ","(N/m) ","(N/m) ","(N/m) ","(N/m) ","(N/m) ","(N/m) ","(N/m) ", &
# "(N/m) ","(N/m) ","(N/m) ","(N/m) ","(N/m) ","(N/m) ","(N/m) ","(N/m) ", &
# "(N/m) ","(N/m) ","(N/m) ","(N/m) ","(N/m) ","(N/m) ","(N/m) ","(N/m) ", &
# "(N/m) ","(N/m) ","(m^2/s)","(1/0) ","(-) ","(N-m/m)","(N-m/m)","(N-m/m)", &
# "(N-m/m)","(N-m/m)","(N-m/m)","(N-m/m)","(N-m/m)","(N-m/m)","(N-m/m)","(N-m/m)", &
# "(N-m/m)","(N-m/m)","(N-m/m)","(N-m/m)","(N-m/m)","(N-m/m)","(N-m/m)","(N-m/m)", &
# "(N-m/m)","(N-m/m)","(N-m/m)","(N-m/m)","(N-m/m)","(N-m/m)","(N-m/m)","(deg) ", &
# "(-) ","(-) ","(-) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ", &
# "(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ", &
# "(m/s) ","(m/s) ","(deg) ","(-) ","(-) ","(deg) ","(-) ","(rad) ", &
# "(rad) ","(-) ","(-) ","(-) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ", &
# "(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ", &
# "(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ", &
# "(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ", &
# "(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ", &
# "(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ", &
# "(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) ","(m/s) "/)


for sB in ['B1','B2','B3']:
ADSpanMap['^[A]*'+sB+r'N(\d*)Alpha_\[deg\]'] =sB+'Alpha_[deg]'
ADSpanMap['^[A]*'+sB+r'N(\d*)AxInd_\[-\]' ] =sB+'AxInd_[-]'
Expand All @@ -698,6 +744,7 @@ def spanwiseColAD(Cols):
ADSpanMap['^[A]*'+sB+r'N(\d*)BEM_kp_qs\[-\]' ]=sB+'BEM_kp_qs_[-]'
ADSpanMap['^[A]*'+sB+r'N(\d*)BEM_F_qs\[-\]' ]=sB+'BEM_F_qs_[-]'
ADSpanMap['^[A]*'+sB+r'N(\d*)BEM_CT_qs_\[-\]' ]=sB+'BEM_CT_qs_[-]'
ADSpanMap['^[A]*'+sB+r'N(\d*)Chi_\[deg\]' ]=sB+'Chi_[deg]'
ADSpanMap['^[A]*'+sB+r'N(\d*)Cl_\[-\]' ] =sB+'Cl_[-]'
ADSpanMap['^[A]*'+sB+r'N(\d*)Cd_\[-\]' ] =sB+'Cd_[-]'
ADSpanMap['^[A]*'+sB+r'N(\d*)Cm_\[-\]' ] =sB+'Cm_[-]'
Expand Down Expand Up @@ -1571,10 +1618,10 @@ def bin_mean_DF(df, xbins, colBin ):
raise Exception('The column `{}` does not appear to be in the dataframe'.format(colBin))
xmid = (xbins[:-1]+xbins[1:])/2
df['Bin'] = pd.cut(df[colBin], bins=xbins, labels=xmid ) # Adding a column that has bin attribute
df2 = df.groupby('Bin').mean() # Average by bin
df2 = df.groupby('Bin', observed=False).mean() # Average by bin
# also counting
df['Counts'] = 1
dfCount=df[['Counts','Bin']].groupby('Bin').sum()
dfCount=df[['Counts','Bin']].groupby('Bin', observed=False).sum()
df2['Counts'] = dfCount['Counts']
# Just in case some bins are missing (will be nan)
df2 = df2.reindex(xmid)
Expand Down
2 changes: 1 addition & 1 deletion pydatview/io/gnuplot_file.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -33,7 +33,7 @@ class GNUPlotFile(File):
@staticmethod
def defaultExtensions():
""" List of file extensions expected for this fileformat"""
return ['.dat']
return ['.dat','.raw']

@staticmethod
def formatName():
Expand Down
4 changes: 2 additions & 2 deletions pydatview/tools/fatigue.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -213,9 +213,9 @@ def bin_count(x, bins, meanBin=True):
df = pd.DataFrame(data=x, columns=['x'])
xmid = (bins[:-1]+bins[1:])/2
df['x_mid']= pd.cut(df['x'], bins= bins, labels = xmid ) # Adding a column that has bin attribute
df2 = df.groupby('x_mid').mean() # Average by bin
df2 = df.groupby('x_mid', observed=False).mean() # Average by bin
df['N'] = 1
dfCount = df[['N','x_mid']].groupby('x_mid').sum()
dfCount = df[['N','x_mid']].groupby('x_mid', observed=False).sum()
df2['N'] = dfCount['N']
# Just in case some bins are missing (will be nan)
df2 = df2.reindex(xmid)
Expand Down
8 changes: 4 additions & 4 deletions pydatview/tools/signal_analysis.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -260,15 +260,15 @@ def applySampler(x_old, y_old, sampDict, df_old=None):
if sample_time <= 0:
raise Exception('Error: sample time must be positive')

time_index = pd.TimedeltaIndex(x_old, unit="S")
x_new = pd.Series(x_old, index=time_index).resample("{:f}S".format(sample_time)).mean().interpolate().values
time_index = pd.TimedeltaIndex(x_old, unit="s")
x_new = pd.Series(x_old, index=time_index).resample("{:f}s".format(sample_time)).mean().interpolate().values

if df_old is not None:
df_new = df_old.set_index(time_index, inplace=False).resample("{:f}S".format(sample_time)).mean()
df_new = df_old.set_index(time_index, inplace=False).resample("{:f}s".format(sample_time)).mean()
df_new = df_new.interpolate().reset_index(drop=True)
return x_new, df_new
if y_old is not None:
y_new = pd.Series(y_old, index=time_index).resample("{:f}S".format(sample_time)).mean()
y_new = pd.Series(y_old, index=time_index).resample("{:f}s".format(sample_time)).mean()
y_new = y_new.interpolate().values
return x_new, y_new

Expand Down

0 comments on commit e4a2413

Please sign in to comment.