|  |  |  | @ -29,7 +29,9 @@ def forecasted_impacts(df_profile_features, df_forecasted_twl): | 
		
	
		
			
				|  |  |  |  | 
 | 
		
	
		
			
				|  |  |  |  |     # Join with df_profile features to find dune toe and crest elevations | 
		
	
		
			
				|  |  |  |  |     df_forecasted_impacts = df_forecasted_impacts.merge( | 
		
	
		
			
				|  |  |  |  |         df_profile_features.query("profile_type=='prestorm'")[["dune_toe_z", "dune_crest_z"]], how="left", left_index=True, right_index=True | 
		
	
		
			
				|  |  |  |  |         df_profile_features.query("profile_type=='prestorm'").reset_index('profile_type')[["dune_toe_z", | 
		
	
		
			
				|  |  |  |  |                                                                                          "dune_crest_z"]], how="left", | 
		
	
		
			
				|  |  |  |  |         left_on=['site_id'], right_on=['site_id'] | 
		
	
		
			
				|  |  |  |  |     ) | 
		
	
		
			
				|  |  |  |  | 
 | 
		
	
		
			
				|  |  |  |  |     # Compare R_high and R_low wirth dune crest and toe elevations | 
		
	
	
		
			
				
					|  |  |  | 
 |