Copy source dataframe

probabilistic-analysis/probabilistic_assessment.py

@@ -471,7 +471,8 @@ def process(beach_name, beach_scenario, n_runs, start_year, end_year,
                 pbar_year.set_description('Year: {}'.format(year))
 
                 # Create output dataframe
-                df_out = df_in[['beach', 'block', 'profile']]
+                df_out = df_in[['beach', 'block', 'profile']].copy()
+                df_out['block'] = df_out['block'].astype(str)
 
                 # Add info on non-erodable sections
                 df_out = df_out.assign(min_chainage=np.nan)
@@ -486,7 +487,7 @@ def process(beach_name, beach_scenario, n_runs, start_year, end_year,
                 df_out = df_out.assign(segment_gaps=False)
                 for b, p, in zip(segment_gaps['block'],
                                  segment_gaps['profile']):
-                    idx = (str(df_out['block']) == str(b)) & (df_out['profile']
+                    idx = (df_out['block'] == str(b)) & (df_out['profile']
                                                               == p)
                     df_out.loc[idx, 'segment_gaps'] = True