Update app.py

app.py

@@ -140,7 +140,7 @@ def parse_weather_data(data):
 def calculate_total_new_snow(df):
     """
     Calculate total new snow by:
-    1. Summing 3-hour snowfall amounts within each day period
+    1. Using ONLY the 3-hour snowfall amounts
     2. Using 9 AM as the daily reset point
     3. Filtering out obvious anomalies (>9 inches in 3 hours)
     
@@ -153,7 +153,7 @@ def calculate_total_new_snow(df):
     # Sort by datetime to ensure correct calculation
     df = df.sort_values('datetime')
     
-    # Create a copy of the dataframe with just the columns we need
+    # Create a copy of the dataframe with ONLY datetime and 3-hour snowfall
     snow_df = df[['datetime', 'snowfall_3hr']].copy()
     
     # Create a day group that starts at 9 AM instead of midnight
@@ -162,17 +162,21 @@ def calculate_total_new_snow(df):
     )
     
     def process_daily_snow(group):
-        """Sum up the 3-hour snowfall amounts for each day period"""
+        """Sum up ONLY the 3-hour snowfall amounts for each day period"""
         # Sort by time to ensure proper sequence
         group = group.sort_values('datetime')
         
-        # Filter out obvious anomalies (more than 9 inches in 3 hours)
-        MAX_THREE_HOUR_SNOW = 9.0  # Maximum reasonable snow in 3 hours
-        valid_snow = group['snowfall_3hr'].where(group['snowfall_3hr'] <= MAX_THREE_HOUR_SNOW, 0.0)
+        # Print debugging information
+        print(f"\nSnowfall amounts for {group['day_group'].iloc[0]}:")
+        for _, row in group.iterrows():
+            if pd.notna(row['snowfall_3hr']):
+                print(f"{row['datetime'].strftime('%Y-%m-%d %H:%M')}: {row['snowfall_3hr']} inches")
         
-        # Sum up all valid 3-hour amounts
-        daily_total = valid_snow.sum()
+        # Sum only the valid 3-hour amounts, treating NaN as 0
+        valid_amounts = group['snowfall_3hr'].fillna(0)
+        daily_total = valid_amounts.sum()
         
+        print(f"Daily total: {daily_total} inches")
         return daily_total
     
     # Calculate daily snow totals
@@ -211,6 +215,8 @@ def create_daily_snow_plot(df, ax):
                    ha='center', va='bottom')
 
 
+
+
     
 def create_wind_rose(df, ax):
     """Create a wind rose plot"""