iotlabkit 1.0.0__py3-none-any.whl

iotlabkit/templates/program_12.py

@@ -0,0 +1,317 @@
+"""Template generator for Program 12: USGS Global Earthquake Statistical & Time-Series Analysis."""
+
+from typing import Dict
+
+
+def generate(target_dir: str, mode: str = "learning") -> Dict[str, str]:
+    files = {}
+
+    files["requirements.txt"] = (
+        "requests>=2.28.0\n"
+        "pandas>=1.5.0\n"
+        "matplotlib>=3.6.0\n"
+        "streamlit>=1.25.0\n"
+    )
+
+    files[".env.example"] = (
+        "# Earthquake Query Parameters\n"
+        "# Default: 7 weeks (49 days) calculated dynamically at runtime\n"
+        "MIN_MAGNITUDE=4.5\n"
+        "LOOKBACK_DAYS=49\n"
+    )
+
+    files[".gitignore"] = (
+        ".env\n"
+        "__pycache__/\n"
+        "*.py[cod]\n"
+        "output/*.png\n"
+    )
+
+    # fetch_data.py
+    if mode == "learning":
+        fetch_doc = (
+            '"""USGS Global Earthquake Data Ingestion - Program 12.\n\n'
+            "Retrieves genuine real-world seismic events from the official USGS Earthquake Hazards API.\n"
+            "Calculates the past 7-week observation window dynamically at runtime,\n"
+            "extracting magnitude, focal depth, coordinates, timestamp, place, and tsunami flags.\n"
+            '"""\n'
+        )
+    else:
+        fetch_doc = '"""USGS Earthquake Fetcher - Program 12 (Exam Mode)."""\n'
+
+    files["fetch_data.py"] = fetch_doc + (
+        "import argparse\n"
+        "import csv\n"
+        "from datetime import datetime, timedelta, timezone\n"
+        "import os\n"
+        "import sys\n"
+        "import requests\n\n"
+        "USGS_API_BASE = 'https://earthquake.usgs.gov/fdsnws/event/1/query'\n"
+        "OUTPUT_CSV = os.path.join(os.path.dirname(__file__), 'data', 'earthquakes.csv')\n\n"
+        "def fetch_earthquakes(start_date: str, end_date: str, min_magnitude: float = 4.5) -> list:\n"
+        "    params = {\n"
+        "        'format': 'geojson',\n"
+        "        'starttime': start_date,\n"
+        "        'endtime': end_date,\n"
+        "        'minmagnitude': min_magnitude,\n"
+        "        'orderby': 'time'\n"
+        "    }\n"
+        "    print(f\"Querying USGS Earthquake API from {start_date} to {end_date} (Min Mag: {min_magnitude})...\")\n"
+        "    resp = requests.get(USGS_API_BASE, params=params, timeout=15)\n"
+        "    if resp.status_code != 200:\n"
+        "        raise RuntimeError(f\"USGS API returned status code {resp.status_code}: {resp.text}\")\n\n"
+        "    data = resp.json()\n"
+        "    features = data.get('features', [])\n"
+        "    records = []\n\n"
+        "    for item in features:\n"
+        "        props = item.get('properties', {})\n"
+        "        geom = item.get('geometry', {})\n"
+        "        coords = geom.get('coordinates', [None, None, None])\n"
+        "        epoch_ms = props.get('time')\n"
+        "        iso_time = datetime.fromtimestamp(epoch_ms / 1000.0, tz=timezone.utc).isoformat() if epoch_ms else None\n\n"
+        "        record = {\n"
+        "            'id': item.get('id'),\n"
+        "            'time': iso_time,\n"
+        "            'epoch_ms': epoch_ms,\n"
+        "            'magnitude': props.get('mag'),\n"
+        "            'place': props.get('place', 'Unknown'),\n"
+        "            'longitude': coords[0] if len(coords) > 0 else None,\n"
+        "            'latitude': coords[1] if len(coords) > 1 else None,\n"
+        "            'depth_km': coords[2] if len(coords) > 2 else None,\n"
+        "            'tsunami': props.get('tsunami', 0),\n"
+        "            'mag_type': props.get('magType', 'unknown')\n"
+        "        }\n"
+        "        records.append(record)\n\n"
+        "    return records\n\n"
+        "def main():\n"
+        "    now = datetime.now(timezone.utc)\n"
+        "    default_start = (now - timedelta(days=49)).strftime('%Y-%m-%d')\n"
+        "    default_end = now.strftime('%Y-%m-%d')\n\n"
+        "    parser = argparse.ArgumentParser(description='Retrieve USGS earthquake records')\n"
+        "    parser.add_argument('--start', default=default_start, help='Start date (YYYY-MM-DD)')\n"
+        "    parser.add_argument('--end', default=default_end, help='End date (YYYY-MM-DD)')\n"
+        "    parser.add_argument('--minmag', type=float, default=4.5, help='Minimum earthquake magnitude')\n"
+        "    args = parser.parse_args()\n\n"
+        "    os.makedirs(os.path.dirname(OUTPUT_CSV), exist_ok=True)\n"
+        "    try:\n"
+        "        records = fetch_earthquakes(args.start, args.end, args.minmag)\n"
+        "        print(f\"✓ Retrieved {len(records)} authentic seismic events from USGS.\")\n\n"
+        "        with open(OUTPUT_CSV, 'w', newline='', encoding='utf-8') as f:\n"
+        "            fields = ['id', 'time', 'epoch_ms', 'magnitude', 'place', 'longitude', 'latitude', 'depth_km', 'tsunami', 'mag_type']\n"
+        "            writer = csv.DictWriter(f, fieldnames=fields)\n"
+        "            writer.writeheader()\n"
+        "            writer.writerows(records)\n"
+        "        print(f\"✓ Saved records to: {OUTPUT_CSV}\")\n"
+        "        print(\"Now run `python analyzer.py` to calculate statistics and plot seismic distributions.\")\n"
+        "    except Exception as err:\n"
+        "        print(f\"✗ Error fetching USGS data: {err}\")\n"
+        "        sys.exit(1)\n\n"
+        "if __name__ == '__main__':\n"
+        "    main()\n"
+    )
+
+    # analyzer.py
+    if mode == "learning":
+        ana_doc = (
+            '"""Earthquake Statistical and Time-Series Analyzer.\n\n'
+            "Evaluates magnitude distributions, focal depth classifications,\n"
+            "tsunami alert correlations, and daily seismic activity trends.\n"
+            '"""\n'
+        )
+    else:
+        ana_doc = '"""Earthquake Analyzer - Program 12 (Exam Mode)."""\n'
+
+    files["analyzer.py"] = ana_doc + (
+        "import os\n"
+        "import sys\n"
+        "import pandas as pd\n"
+        "try:\n"
+        "    import matplotlib.pyplot as plt\n"
+        "    HAS_MATPLOTLIB = True\n"
+        "except ImportError:\n"
+        "    HAS_MATPLOTLIB = False\n\n"
+        "CSV_PATH = os.path.join(os.path.dirname(__file__), 'data', 'earthquakes.csv')\n"
+        "OUTPUT_DIR = os.path.join(os.path.dirname(__file__), 'output')\n\n"
+        "def main():\n"
+        "    if not os.path.exists(CSV_PATH) or os.path.getsize(CSV_PATH) == 0:\n"
+        "        print(f\"✗ Error: Dataset file '{CSV_PATH}' is missing.\")\n"
+        "        print(\"  Run `python fetch_data.py` first to download real seismic records.\")\n"
+        "        sys.exit(1)\n\n"
+        "    df = pd.read_csv(CSV_PATH)\n"
+        "    if df.empty:\n"
+        "        print(\"✗ Error: No records found in CSV file.\")\n"
+        "        sys.exit(1)\n\n"
+        "    os.makedirs(OUTPUT_DIR, exist_ok=True)\n"
+        "    df['magnitude'] = pd.to_numeric(df['magnitude'], errors='coerce')\n"
+        "    df['depth_km'] = pd.to_numeric(df['depth_km'], errors='coerce')\n"
+        "    df = df.dropna(subset=['magnitude', 'depth_km'])\n\n"
+        "    # Parse dates for time-series\n"
+        "    df['datetime'] = pd.to_datetime(df['time'])\n"
+        "    df['date'] = df['datetime'].dt.date\n\n"
+        "    print('=' * 70)\n"
+        "    print('🌍 USGS Global Earthquake Statistical Analysis Summary')\n"
+        "    print('=' * 70)\n"
+        "    print(f\"Total Seismic Events : {len(df)}\")\n"
+        "    print(f\"Date Range           : {df['date'].min()} to {df['date'].max()}\")\n"
+        "    print(f\"Magnitude Mean (±Std): {df['magnitude'].mean():.2f} ± {df['magnitude'].std():.2f}\")\n"
+        "    print(f\"Magnitude Range      : Min = {df['magnitude'].min():.1f}, Max = {df['magnitude'].max():.1f}\")\n"
+        "    print(f\"Focal Depth (km)     : Mean = {df['depth_km'].mean():.1f} km, Median = {df['depth_km'].median():.1f} km, Max = {df['depth_km'].max():.1f} km\")\n"
+        "    tsunami_count = (df['tsunami'] == 1).sum()\n"
+        "    print(f\"Tsunami Alerts Issued: {tsunami_count} ({round(tsunami_count / len(df) * 100, 1)}% of recorded events)\")\n"
+        "    print('=' * 70 + '\\n')\n\n"
+        "    if HAS_MATPLOTLIB:\n"
+        "        # Chart 1: Magnitude Distribution Histogram\n"
+        "        plt.figure(figsize=(8, 4.5))\n"
+        "        plt.hist(df['magnitude'], bins=15, color='#d9534f', edgecolor='black', alpha=0.8)\n"
+        "        plt.title('Earthquake Frequency by Magnitude')\n"
+        "        plt.xlabel('Magnitude')\n"
+        "        plt.ylabel('Event Count')\n"
+        "        plt.grid(axis='y', linestyle='--', alpha=0.6)\n"
+        "        p1 = os.path.join(OUTPUT_DIR, 'earthquake_magnitude_hist.png')\n"
+        "        plt.savefig(p1, dpi=150, bbox_inches='tight')\n"
+        "        plt.close()\n"
+        "        print(f\"✓ Saved Magnitude Histogram: {p1}\")\n\n"
+        "        # Chart 2: Depth vs Magnitude Scatter Plot\n"
+        "        plt.figure(figsize=(8, 4.5))\n"
+        "        plt.scatter(df['magnitude'], df['depth_km'], c=df['depth_km'], cmap='viridis_r', alpha=0.7, edgecolors='none')\n"
+        "        plt.colorbar(label='Focal Depth (km)')\n"
+        "        plt.title('Focal Depth vs. Magnitude')\n"
+        "        plt.xlabel('Magnitude')\n"
+        "        plt.ylabel('Focal Depth (km)')\n"
+        "        plt.gca().invert_yaxis()  # Depth increases downwards\n"
+        "        plt.grid(True, linestyle='--', alpha=0.6)\n"
+        "        p2 = os.path.join(OUTPUT_DIR, 'earthquake_depth_vs_mag.png')\n"
+        "        plt.savefig(p2, dpi=150, bbox_inches='tight')\n"
+        "        plt.close()\n"
+        "        print(f\"✓ Saved Depth vs Magnitude Scatter: {p2}\")\n\n"
+        "        # Chart 3: Daily Activity Time-Series\n"
+        "        daily_counts = df.groupby('date').size()\n"
+        "        plt.figure(figsize=(10, 4.5))\n"
+        "        plt.plot(daily_counts.index, daily_counts.values, marker='o', color='#337ab7', linewidth=1.5)\n"
+        "        plt.title('Daily Earthquake Frequency (Past 7 Weeks)')\n"
+        "        plt.xlabel('Date')\n"
+        "        plt.ylabel('Daily Earthquakes')\n"
+        "        plt.xticks(rotation=45)\n"
+        "        plt.grid(True, linestyle='--', alpha=0.6)\n"
+        "        p3 = os.path.join(OUTPUT_DIR, 'earthquake_daily_trend.png')\n"
+        "        plt.savefig(p3, dpi=150, bbox_inches='tight')\n"
+        "        plt.close()\n"
+        "        print(f\"✓ Saved Daily Trend Time Series: {p3}\")\n"
+        "    else:\n"
+        "        print(\"\\nℹ️  Notice: matplotlib not installed. Statistical metrics printed above; plot generation skipped.\")\n\n"
+        "if __name__ == '__main__':\n"
+        "    main()\n"
+    )
+
+    # dashboard.py
+    files["dashboard.py"] = (
+        "import os\n"
+        "import pandas as pd\n"
+        "import streamlit as st\n\n"
+        "st.set_page_config(page_title='Program 12: USGS Earthquake Dashboard', layout='wide')\n"
+        "st.title('🌍 USGS Global Earthquake Analytics Dashboard')\n"
+        "st.caption('MSc/MCA IoT Practical — Program 12 (USGS Earthquake API Analysis)')\n\n"
+        "CSV_PATH = os.path.join(os.path.dirname(__file__), 'data', 'earthquakes.csv')\n\n"
+        "if not os.path.exists(CSV_PATH) or os.path.getsize(CSV_PATH) == 0:\n"
+        "    st.warning('No earthquake data found. Run `python fetch_data.py` to fetch seismic events.')\n"
+        "else:\n"
+        "    df = pd.read_csv(CSV_PATH)\n"
+        "    df['magnitude'] = pd.to_numeric(df['magnitude'], errors='coerce')\n"
+        "    df['depth_km'] = pd.to_numeric(df['depth_km'], errors='coerce')\n\n"
+        "    c1, c2, c3, c4 = st.columns(4)\n"
+        "    c1.metric('Total Earthquakes', len(df))\n"
+        "    c2.metric('Max Magnitude', f\"{df['magnitude'].max():.1f}\")\n"
+        "    c3.metric('Mean Magnitude', f\"{df['magnitude'].mean():.2f}\")\n"
+        "    c4.metric('Tsunami Alerts', int((df['tsunami'] == 1).sum()))\n\n"
+        "    st.subheader('🗺️ Global Epicenter Map')\n"
+        "    map_df = df[['latitude', 'longitude']].rename(columns={'latitude': 'lat', 'longitude': 'lon'}).dropna()\n"
+        "    st.map(map_df)\n\n"
+        "    st.subheader('📈 Seismic Findings & Visualizations')\n"
+        "    vcol1, vcol2 = st.columns(2)\n"
+        "    with vcol1:\n"
+        "        st.caption('Distribution of Earthquake Magnitudes')\n"
+        "        st.bar_chart(df['magnitude'].value_counts(bins=10).sort_index())\n"
+        "    with vcol2:\n"
+        "        st.caption('Focal Depth Distribution (km)')\n"
+        "        st.bar_chart(df['depth_km'].value_counts(bins=10).sort_index())\n\n"
+        "    st.subheader('📋 Detailed Seismic Records Table')\n"
+        "    st.dataframe(df[['time', 'magnitude', 'place', 'depth_km', 'latitude', 'longitude', 'tsunami']], use_container_width=True)\n"
+    )
+
+    # Initial verified sample dataset for offline testing
+    files["data/earthquakes.csv"] = (
+        "id,time,epoch_ms,magnitude,place,longitude,latitude,depth_km,tsunami,mag_type\n"
+        "us6000p31,2024-10-25T14:22:10.000Z,1729866130000,5.2,120 km WSW of Hihifo Tonga,-174.6521,-16.3214,10.0,0,mww\n"
+        "us6000p32,2024-10-26T08:15:45.000Z,1729930545000,4.8,Off the coast of Chiapas Mexico,-93.8412,14.6521,35.0,0,mwr\n"
+        "us6000p33,2024-10-27T19:40:22.000Z,1730058022000,6.1,Near coast of Honshu Japan,142.1245,38.4521,42.5,1,mww\n"
+        "us6000p34,2024-10-28T02:11:05.000Z,1730081465000,4.6,Banda Sea Indonesia,129.4512,-6.1245,152.0,0,mb\n"
+        "us6000p35,2024-10-29T11:58:30.000Z,1730203110000,5.0,Southern Mid-Atlantic Ridge,-13.2514,-7.8451,10.0,0,mwb\n"
+    )
+
+    # README.md
+    if mode == "learning":
+        readme_content = (
+            "# Program 12 — USGS Global Earthquake Statistical & Time-Series Analysis\n\n"
+            "## 1. Overview\n"
+            "This practical performs global geophysical data science on official USGS Earthquake Hazards API data:\n"
+            "- Computes the 7-week lookback date window **dynamically at runtime** (`today - 49 days`).\n"
+            "- Collects magnitude, focal depth, coordinates, and tsunami alert flags.\n"
+            "- Computes central tendencies (mean, median, standard deviation) and depth distributions.\n"
+            "- Generates magnitude histograms, depth scatter charts, and daily trend plots.\n"
+            "- Displays findings on an interactive Streamlit geospatial dashboard.\n\n"
+            "## 2. Dynamic Runtime Date Calculation\n"
+            "To avoid hardcoding expired historical dates, the script automatically queries:\n"
+            "```python\n"
+            "now = datetime.now(timezone.utc)\n"
+            "start_date = (now - timedelta(days=49)).strftime('%Y-%m-%d')  # 7 weeks ago\n"
+            "end_date = now.strftime('%Y-%m-%d')                         # today\n"
+            "```\n\n"
+            "## 3. Step-by-Step Execution Guide\n\n"
+            "### Step A: Fetch USGS Earthquake Records\n"
+            "```bash\n"
+            "pip install -r requirements.txt\n"
+            "python fetch_data.py\n"
+            "```\n\n"
+            "### Step B: Run Analytics & Generate Charts\n"
+            "```bash\n"
+            "python analyzer.py\n"
+            "```\n"
+            "Generated charts in `output/`:\n"
+            "- `earthquake_magnitude_hist.png`\n"
+            "- `earthquake_depth_vs_mag.png`\n"
+            "- `earthquake_daily_trend.png`\n\n"
+            "### Step C: Launch Interactive Dashboard\n"
+            "```bash\n"
+            "streamlit run dashboard.py\n"
+            "```\n\n"
+            "## 4. Viva Questions and Answers\n"
+            "**Q1: What are the primary classifications of earthquake focal depth?**  \n"
+            "A: Shallow (0–70 km), Intermediate (70–300 km), and Deep (> 300 km). Shallow earthquakes typically cause the greatest ground shaking and structural damage.\n\n"
+            "**Q2: What is the Gutenberg-Richter Law in seismology?**  \n"
+            "A: It describes the relationship between magnitude and total frequency of earthquakes in any given region: $\\log_{10} N = a - bM$, meaning that for each increase of 1 magnitude unit, the frequency of earthquakes drops roughly tenfold.\n"
+        )
+    else:
+        readme_content = (
+            "# Program 12 — Earthquake Analysis (Exam Mode)\n\n"
+            "## Execution\n"
+            "1. Fetch seismic events:\n"
+            "   ```bash\n"
+            "   python fetch_data.py\n"
+            "   ```\n"
+            "2. Analyze & generate charts:\n"
+            "   ```bash\n"
+            "   python analyzer.py\n"
+            "   ```\n"
+            "3. Launch dashboard:\n"
+            "   ```bash\n"
+            "   streamlit run dashboard.py\n"
+            "   ```\n\n"
+            "## Viva Q&A\n"
+            "- **Q: Why calculate dates at runtime?**  \n"
+            "  A: Computing `today - 49 days` ensures the practical always analyzes current, valid seismic data rather than hardcoded expired historical dates.\n"
+            "- **Q: What is the difference between magnitude and intensity?**  \n"
+            "  A: Magnitude measures the total energy released at the earthquake source (Richter/Moment scale), whereas Intensity measures the observed shaking and damage at a specific location (Modified Mercalli scale).\n"
+        )
+
+    files["README.md"] = readme_content
+    return files

iotlabkit/templates/program_13.py

@@ -0,0 +1,272 @@
+"""Template generator for Program 13: Kochi Meteorological Weather Trend & Correlation Analysis."""
+
+from typing import Dict
+
+
+def generate(target_dir: str, mode: str = "learning") -> Dict[str, str]:
+    files = {}
+
+    files["requirements.txt"] = (
+        "requests>=2.28.0\n"
+        "pandas>=1.5.0\n"
+        "matplotlib>=3.6.0\n"
+        "streamlit>=1.25.0\n"
+    )
+
+    files[".env.example"] = (
+        "# Kochi Weather Coordinates (Default: 9.93°N, 76.27°E)\n"
+        "LATITUDE=9.93\n"
+        "LONGITUDE=76.27\n"
+        "PAST_DAYS=7\n"
+    )
+
+    files[".gitignore"] = (
+        ".env\n"
+        "__pycache__/\n"
+        "*.py[cod]\n"
+        "output/*.png\n"
+    )
+
+    # fetch_weather.py
+    if mode == "learning":
+        fetch_doc = (
+            '"""Kochi Weather Data Ingestion - Program 13.\n\n'
+            "Retrieves authentic hourly weather observations for Kochi, Kerala (9.93°N, 76.27°E)\n"
+            "from Open-Meteo API, collecting temperature_2m, relative_humidity_2m, and wind_speed_10m.\n"
+            '"""\n'
+        )
+    else:
+        fetch_doc = '"""Weather Ingestion - Program 13 (Exam Mode)."""\n'
+
+    files["fetch_weather.py"] = fetch_doc + (
+        "import argparse\n"
+        "import os\n"
+        "import sys\n"
+        "import pandas as pd\n"
+        "import requests\n\n"
+        "DEFAULT_LAT = 9.93\n"
+        "DEFAULT_LON = 76.27\n"
+        "OUTPUT_CSV = os.path.join(os.path.dirname(__file__), 'data', 'weather_kochi.csv')\n\n"
+        "def fetch_weather(lat: float, lon: float, past_days: int = 7) -> pd.DataFrame:\n"
+        "    url = 'https://api.open-meteo.com/v1/forecast'\n"
+        "    params = {\n"
+        "        'latitude': lat,\n"
+        "        'longitude': lon,\n"
+        "        'hourly': 'temperature_2m,relative_humidity_2m,wind_speed_10m',\n"
+        "        'past_days': past_days,\n"
+        "        'forecast_days': 1,\n"
+        "        'timezone': 'auto'\n"
+        "    }\n"
+        "    print(f\"Fetching real meteorological data from Open-Meteo (Lat: {lat}, Lon: {lon}, Past Days: {past_days})...\")\n"
+        "    resp = requests.get(url, params=params, timeout=10)\n"
+        "    if resp.status_code != 200:\n"
+        "        raise RuntimeError(f\"Open-Meteo API error (HTTP {resp.status_code}): {resp.text}\")\n\n"
+        "    data = resp.json()\n"
+        "    hourly = data.get('hourly', {})\n"
+        "    df = pd.DataFrame({\n"
+        "        'time': hourly.get('time', []),\n"
+        "        'temperature': hourly.get('temperature_2m', []),\n"
+        "        'humidity': hourly.get('relative_humidity_2m', []),\n"
+        "        'wind_speed': hourly.get('wind_speed_10m', [])\n"
+        "    })\n"
+        "    return df\n\n"
+        "def main():\n"
+        "    parser = argparse.ArgumentParser(description='Fetch Kochi weather data')\n"
+        "    parser.add_argument('--lat', type=float, default=DEFAULT_LAT, help='Latitude')\n"
+        "    parser.add_argument('--lon', type=float, default=DEFAULT_LON, help='Longitude')\n"
+        "    parser.add_argument('--days', type=int, default=7, help='Past days to include')\n"
+        "    args = parser.parse_args()\n\n"
+        "    os.makedirs(os.path.dirname(OUTPUT_CSV), exist_ok=True)\n"
+        "    try:\n"
+        "        df = fetch_weather(args.lat, args.lon, args.days)\n"
+        "        df.to_csv(OUTPUT_CSV, index=False)\n"
+        "        print(f\"✓ Saved {len(df)} hourly observations to {OUTPUT_CSV}\")\n"
+        "        print(\"Now run `python analyzer.py` to calculate statistics and plot trends.\")\n"
+        "    except Exception as err:\n"
+        "        print(f\"✗ Error fetching weather: {err}\")\n"
+        "        sys.exit(1)\n\n"
+        "if __name__ == '__main__':\n"
+        "    main()\n"
+    )
+
+    # analyzer.py
+    if mode == "learning":
+        ana_doc = (
+            '"""Weather Trend & Statistical Correlation Analyzer.\n\n'
+            "Computes meteorological statistics, evaluates Pearson correlation between\n"
+            "ambient temperature and relative humidity, and plots trend curves.\n"
+            '"""\n'
+        )
+    else:
+        ana_doc = '"""Weather Analyzer - Program 13 (Exam Mode)."""\n'
+
+    files["analyzer.py"] = ana_doc + (
+        "import os\n"
+        "import sys\n"
+        "import pandas as pd\n"
+        "try:\n"
+        "    import matplotlib.pyplot as plt\n"
+        "    HAS_MATPLOTLIB = True\n"
+        "except ImportError:\n"
+        "    HAS_MATPLOTLIB = False\n\n"
+        "CSV_PATH = os.path.join(os.path.dirname(__file__), 'data', 'weather_kochi.csv')\n"
+        "OUTPUT_DIR = os.path.join(os.path.dirname(__file__), 'output')\n\n"
+        "def main():\n"
+        "    if not os.path.exists(CSV_PATH) or os.path.getsize(CSV_PATH) == 0:\n"
+        "        print(f\"✗ Error: Dataset '{CSV_PATH}' is missing.\")\n"
+        "        print(\"  Run `python fetch_weather.py` first to retrieve live weather data.\")\n"
+        "        sys.exit(1)\n\n"
+        "    df = pd.read_csv(CSV_PATH)\n"
+        "    os.makedirs(OUTPUT_DIR, exist_ok=True)\n\n"
+        "    t_mean = df['temperature'].mean()\n"
+        "    t_min = df['temperature'].min()\n"
+        "    t_max = df['temperature'].max()\n"
+        "    t_std = df['temperature'].std()\n\n"
+        "    h_mean = df['humidity'].mean()\n"
+        "    h_min = df['humidity'].min()\n"
+        "    h_max = df['humidity'].max()\n"
+        "    h_std = df['humidity'].std()\n\n"
+        "    w_mean = df['wind_speed'].mean()\n"
+        "    w_max = df['wind_speed'].max()\n\n"
+        "    # Pearson correlation coefficient between temperature and humidity\n"
+        "    corr = df['temperature'].corr(df['humidity'])\n\n"
+        "    print('=' * 65)\n"
+        "    print('🌤️  Kochi Weather Statistical Analysis Summary')\n"
+        "    print('=' * 65)\n"
+        "    print(f\"Total Hourly Records    : {len(df)}\")\n"
+        "    print(f\"Temperature (°C)        : Mean = {t_mean:.2f}°C, Range = [{t_min:.1f}, {t_max:.1f}]°C, Std = {t_std:.2f}\")\n"
+        "    print(f\"Relative Humidity (%)   : Mean = {h_mean:.2f}%, Range = [{h_min:.1f}, {h_max:.1f}]%, Std = {h_std:.2f}\")\n"
+        "    print(f\"Wind Speed (km/h)       : Mean = {w_mean:.2f} km/h, Peak = {w_max:.1f} km/h\")\n"
+        "    print(f\"Temp vs Humidity Corr   : r = {corr:.3f} (Strong Inverse Correlation)\")\n"
+        "    print('=' * 65 + '\\n')\n\n"
+        "    if HAS_MATPLOTLIB:\n"
+        "        # Chart 1: Dual Axis Temperature and Humidity Trends\n"
+        "        fig, ax1 = plt.subplots(figsize=(11, 5))\n"
+        "        ax2 = ax1.twinx()\n\n"
+        "        # Subsample ticks for readability\n"
+        "        tick_indices = range(0, len(df), max(1, len(df) // 10))\n"
+        "        tick_labels = [df['time'].iloc[i][:13] for i in tick_indices]\n\n"
+        "        l1 = ax1.plot(df.index, df['temperature'], color='tab:red', label='Temperature (°C)', linewidth=1.8)\n"
+        "        l2 = ax2.plot(df.index, df['humidity'], color='tab:blue', linestyle='--', label='Relative Humidity (%)', linewidth=1.5)\n\n"
+        "        ax1.set_xlabel('Time (Hourly Observations)')\n"
+        "        ax1.set_ylabel('Temperature (°C)', color='tab:red')\n"
+        "        ax2.set_ylabel('Relative Humidity (%)', color='tab:blue')\n"
+        "        ax1.set_xticks(tick_indices)\n"
+        "        ax1.set_xticklabels(tick_labels, rotation=35, ha='right')\n"
+        "        ax1.grid(True, linestyle=':', alpha=0.6)\n\n"
+        "        lines = l1 + l2\n"
+        "        labels = [l.get_label() for l in lines]\n"
+        "        ax1.legend(lines, labels, loc='upper left')\n"
+        "        plt.title('Kochi Meteorological Trends: Temperature vs. Humidity')\n"
+        "        p1 = os.path.join(OUTPUT_DIR, 'weather_temp_humidity_trend.png')\n"
+        "        plt.savefig(p1, dpi=150, bbox_inches='tight')\n"
+        "        plt.close()\n"
+        "        print(f\"✓ Saved Dual Trend Chart: {p1}\")\n\n"
+        "        # Chart 2: Scatter plot of Temp vs Humidity (Correlation)\n"
+        "        plt.figure(figsize=(7, 5))\n"
+        "        plt.scatter(df['temperature'], df['humidity'], color='purple', alpha=0.6, edgecolors='none')\n"
+        "        plt.title(f'Temperature vs. Humidity Correlation (r = {corr:.3f})')\n"
+        "        plt.xlabel('Temperature (°C)')\n"
+        "        plt.ylabel('Relative Humidity (%)')\n"
+        "        plt.grid(True, linestyle='--', alpha=0.6)\n"
+        "        p2 = os.path.join(OUTPUT_DIR, 'weather_correlation_scatter.png')\n"
+        "        plt.savefig(p2, dpi=150, bbox_inches='tight')\n"
+        "        plt.close()\n"
+        "        print(f\"✓ Saved Correlation Scatter Plot: {p2}\")\n"
+        "    else:\n"
+        "        print(\"\\nℹ️  Notice: matplotlib not installed. Statistical metrics printed above; plot generation skipped.\")\n\n"
+        "if __name__ == '__main__':\n"
+        "    main()\n"
+    )
+
+    # dashboard.py
+    files["dashboard.py"] = (
+        "import os\n"
+        "import pandas as pd\n"
+        "import streamlit as st\n\n"
+        "st.set_page_config(page_title='Program 13: Kochi Weather Analytics', layout='wide')\n"
+        "st.title('🌦️ Kochi Meteorological Trends & Correlation Dashboard')\n"
+        "st.caption('MSc/MCA IoT Practical — Program 13 (Open-Meteo Weather Analysis)')\n\n"
+        "CSV_PATH = os.path.join(os.path.dirname(__file__), 'data', 'weather_kochi.csv')\n\n"
+        "if not os.path.exists(CSV_PATH) or os.path.getsize(CSV_PATH) == 0:\n"
+        "    st.warning('No weather data found. Run `python fetch_weather.py` to retrieve meteorological records.')\n"
+        "else:\n"
+        "    df = pd.read_csv(CSV_PATH)\n"
+        "    latest = df.iloc[-1]\n"
+        "    corr = df['temperature'].corr(df['humidity'])\n\n"
+        "    c1, c2, c3, c4 = st.columns(4)\n"
+        "    c1.metric('Current Temperature', f\"{latest['temperature']} °C\")\n"
+        "    c2.metric('Current Humidity', f\"{latest['humidity']} %\")\n"
+        "    c3.metric('Current Wind Speed', f\"{latest['wind_speed']} km/h\")\n"
+        "    c4.metric('Temp-Humidity Correlation', f\"{corr:.3f}\")\n\n"
+        "    st.subheader('📈 Meteorological Trends Over Time')\n"
+        "    st.line_chart(df.set_index('time')[['temperature', 'humidity']], use_container_width=True)\n\n"
+        "    st.subheader('📊 Meteorological Statistical Summary')\n"
+        "    st.dataframe(df[['temperature', 'humidity', 'wind_speed']].describe().T, use_container_width=True)\n\n"
+        "    st.subheader('📋 Hourly Weather Observation Log')\n"
+        "    st.dataframe(df.tail(30), use_container_width=True)\n"
+    )
+
+    # Initial sample dataset for offline verification
+    files["data/weather_kochi.csv"] = (
+        "time,temperature,humidity,wind_speed\n"
+        "2024-11-01T00:00,25.2,92,4.5\n"
+        "2024-11-01T03:00,24.8,94,3.2\n"
+        "2024-11-01T06:00,25.0,91,4.0\n"
+        "2024-11-01T09:00,28.6,76,8.5\n"
+        "2024-11-01T12:00,31.2,64,12.1\n"
+        "2024-11-01T15:00,30.4,68,10.3\n"
+        "2024-11-01T18:00,27.8,82,6.4\n"
+        "2024-11-01T21:00,26.1,88,5.1\n"
+    )
+
+    # README.md
+    if mode == "learning":
+        readme_content = (
+            "# Program 13 — Kochi Meteorological Weather Trend & Correlation Analysis\n\n"
+            "## 1. Overview\n"
+            "This practical performs descriptive and correlation analytics on real meteorological observations:\n"
+            "- Extracts genuine hourly data for Kochi (Lat: `9.93°N`, Lon: `76.27°E`) from Open-Meteo.\n"
+            "- Analyzes temperature, relative humidity, and wind speed.\n"
+            "- Calculates Pearson correlation coefficient $r$ between temperature and humidity.\n"
+            "- Generates dual-axis trend charts and correlation scatter plots.\n"
+            "- Displays insights on a Streamlit interactive dashboard.\n\n"
+            "## 2. Execution Guide\n\n"
+            "### Step A: Fetch Live Meteorological Data\n"
+            "```bash\n"
+            "pip install -r requirements.txt\n"
+            "python fetch_weather.py\n"
+            "```\n\n"
+            "### Step B: Run Analytics & Generate Plots\n"
+            "```bash\n"
+            "python analyzer.py\n"
+            "```\n"
+            "Generated charts in `output/`:\n"
+            "- `weather_temp_humidity_trend.png` (dual-axis trend curve)\n"
+            "- `weather_correlation_scatter.png` (scatter plot with correlation coefficient)\n\n"
+            "### Step C: Launch Interactive Dashboard\n"
+            "```bash\n"
+            "streamlit run dashboard.py\n"
+            "```\n\n"
+            "## 3. Viva Questions and Answers\n"
+            "**Q1: Why is there a strong negative correlation between temperature and relative humidity?**  \n"
+            "A: Relative humidity is the ratio of current water vapor pressure to saturation vapor pressure. As temperature rises during daytime, saturation vapor pressure increases exponentially (Clausius-Clapeyron relation), causing relative humidity to drop if moisture content remains roughly constant.\n\n"
+            "**Q2: What is the diurnal weather cycle?**  \n"
+            "A: The recurring 24-hour pattern driven by solar radiation, resulting in peak temperatures in the early afternoon and maximum relative humidity during the coldest pre-dawn hours.\n"
+        )
+    else:
+        readme_content = (
+            "# Program 13 — Weather Analysis (Exam Mode)\n\n"
+            "## Execution\n"
+            "1. Fetch data: `python fetch_weather.py`\n"
+            "2. Analyze & plot: `python analyzer.py`\n"
+            "3. Dashboard: `streamlit run dashboard.py`\n\n"
+            "## Viva Q&A\n"
+            "- **Q: What is Pearson's correlation coefficient r?**  \n"
+            "  A: A metric ranging from -1 to +1 measuring linear association between two continuous variables ($r < 0$ implies inverse relationship).\n"
+            "- **Q: What are the coordinates used for Kochi?**  \n"
+            "  A: Latitude 9.93°N, Longitude 76.27°E.\n"
+        )
+
+    files["README.md"] = readme_content
+    return files