Add ezuetell preliminary analysis of August forecast

.gitignore

@@ -219,4 +219,8 @@ __marimo__/
 # Streamlit
 .streamlit/secrets.toml
 
+#Local data
 data/
+
+#Images
+.png

analysis/a08_qdm_may.ipynb

@@ -14,7 +14,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": 11,
    "id": "01ed0601",
    "metadata": {},
    "outputs": [
@@ -24,7 +24,7 @@
        "<module 'isku_utils' from '/home/emily_zuetell/projects/poreallas/analysis/isku_utils.py'>"
       ]
      },
-     "execution_count": 2,
+     "execution_count": 11,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -51,12 +51,12 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": null,
    "id": "06d22f1d",
    "metadata": {},
    "outputs": [],
    "source": [
-    "load_dotenv(dotenv_path=\"/home/emily_zuetell/projects/poreallas/local.env\")\n",
+    "load_dotenv()\n",
     "DATA_DIR = os.environ[\"DATA_DIR\"]\n",
     "EFFECTS_URI = os.environ[\"POREALLAS_EFFECTS_URI\"]\n",
     "IMPACT_REGION_POLYGONS = os.environ[\"POREALLAS_REGIONS_POLYGONS_URI\"]\n",
@@ -69,22 +69,14 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 4,
+   "execution_count": null,
    "id": "d23e9c64",
    "metadata": {},
-   "outputs": [
-    {
-     "name": "stderr",
-     "output_type": "stream",
-     "text": [
-      "/home/emily_zuetell/projects/poreallas/.venv/lib/python3.14/site-packages/zarr/core/group.py:3289: ZarrUserWarning: Object at zarr.json:Zone.Identifier is not recognized as a component of a Zarr hierarchy.\n",
-      "  warnings.warn(\n"
-     ]
-    }
-   ],
+   "outputs": [],
    "source": [
     "#Projection Effects\n",
-    "effect = xr.open_datatree(os.path.join(DATA_DIR, EFFECTS_URI), consolidated = False)\n",
+    "effect = xr.open_datatree(EFFECTS_URI, consolidated = False)\n",
+    "effect_qdm = xr.open_datatree(\"gs://poreallas-public-20260605/v20260731/parsed/effects.zarr\", backend_kwargs={\"storage_options\": {\"token\": \"anon\"}})\n",
     "baseline_period = analysis_utils.get_baseline_period(effect, years = 30)\n",
     "#Impact Regions\n",
     "_polygons = (\n",
@@ -101,23 +93,19 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 5,
-   "id": "e8fe1e54",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "effect_qdm = xr.open_datatree(\"gs://poreallas-public-20260605/v20260731/parsed/effects.zarr\", backend_kwargs={\"storage_options\": {\"token\": \"anon\"}})"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": 7,
+   "execution_count": null,
    "id": "cccf726d",
    "metadata": {},
    "outputs": [],
    "source": [
-    "impact_delta = analysis_utils.compute_impact(effect, socioeconomics, ensemble = True, baseline_period=baseline_period, hotonly = False, rate = True, age_weight = True)\n",
-    "impact_qdm = analysis_utils.compute_impact(effect_qdm, socioeconomics, ensemble = True, baseline_period=baseline_period, hotonly = False, rate = True, age_weight = True)"
+    "impact_delta = analysis_utils.compute_impact(effect, socioeconomics, \n",
+    "                                             chunks = {'number':-1, 'region': 'auto'},\n",
+    "                                             ensemble = True, baseline_period=baseline_period, hotonly = \"coldonly\", \n",
+    "                                             rate = False, age_weight = True)\n",
+    "impact_qdm = analysis_utils.compute_impact(effect_qdm, socioeconomics, \n",
+    "                                             chunks = {'number':-1, 'region': 'auto'},\n",
+    "                                             ensemble = True, baseline_period=baseline_period, hotonly = \"coldonly\", \n",
+    "                                             rate = False, age_weight = True)"
    ]
   },
   {