import pandas as pd

df = pd.read_csv('bcdata_wndf.csv')


def grp_yrs(df, period):
    print(period)
    grp = df.groupby(['lat', 'lon', 'tmax_window', 'wind_thld', 'tmax_thld', 'gcm']).mean()
    grp = grp.reset_index()
    grp['year'] = period
    return grp

#1981-1990s-2009
#2020-2030s-2049
#2040-2050s-2069
#2070-2080s-2099

df1990s = grp_yrs(df[(df.year >= 1981) & (df.year <= 2010)], '1990s')
df2030s = grp_yrs(df[(df.year >= 2020) & (df.year <= 2049)], '2030s')
df2050s = grp_yrs(df[(df.year >= 2040) & (df.year <= 2069)], '2050s')
df2080s = grp_yrs(df[(df.year >= 2070) & (df.year <= 2099)], '2080s')

df = pd.concat([df1990s, df2030s, df2050s, df2080s])
df.to_csv('boxdata_wndf.csv', index=False, float_format='%4f')