diff --git a/docs/tutorials/13 - US Spatial Analysis for US Si & CdTe/1 - PCA baselines for Si & CdTe.ipynb b/docs/tutorials/13 - US Spatial Analysis for US Si & CdTe/1 - PCA baselines for Si & CdTe.ipynb
new file mode 100644
index 0000000..e70a4b7
--- /dev/null
+++ b/docs/tutorials/13 - US Spatial Analysis for US Si & CdTe/1 - PCA baselines for Si & CdTe.ipynb
@@ -0,0 +1,700 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# REEDs PCA Baselines for Si & CdTe"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "ReEDS installation projections from the Solar Futures Study, published by DOE on 2021.\n",
+ "\n",
+ "Scenario Interest: \n",
+ "o\t95-by-35+Elec.Adv+DR , a.k.a. \"Solar Futures Decarbonization + Electrification scenario\"\n",
+ "\n",
+ "This code performs three Methods to output files for the Mass Flows:\n",
+ "\n",
+ " - PCA original data by ReEDS
\n",
+ " - PCA data reordered based on ascending production up to 2035, and then descending,
\n",
+ " - PCA data calculated with an exponencial that mathces the cumulative 2035 and 2050 targets
\n",
+ "
"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "import PV_ICE\n",
+ "import numpy as np\n",
+ "import pandas as pd\n",
+ "import os,sys\n",
+ "import matplotlib.pyplot as plt\n",
+ "\n",
+ "plt.rcParams.update({'font.size': 22})\n",
+ "plt.rcParams['figure.figsize'] = (12, 8)"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 2,
+ "metadata": {},
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Your simulation will be stored in C:\\Users\\sayala\\Documents\\GitHub\\PV_ICE\\docs\\tutorials\\13 - US Spatial Analysis for US Si & CdTe\\TEMP\n"
+ ]
+ }
+ ],
+ "source": [
+ "import os\n",
+ "from pathlib import Path\n",
+ "\n",
+ "cwd = os.getcwd()\n",
+ "testfolder = os.path.join(cwd, 'TEMP')\n",
+ "\n",
+ "if not os.path.exists(testfolder):\n",
+ " os.makedirs(testfolder)\n",
+ "print (\"Your simulation will be stored in %s\" % testfolder)\n"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 3,
+ "metadata": {},
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Working on a Windows 10\n",
+ "Python version 3.11.7 | packaged by Anaconda, Inc. | (main, Dec 15 2023, 18:05:47) [MSC v.1916 64 bit (AMD64)]\n",
+ "Pandas version 2.1.0\n",
+ "PV_ICE version 0.3.4+7.g5a5d2dc.dirty\n"
+ ]
+ }
+ ],
+ "source": [
+ "# This information helps with debugging and getting support :)\n",
+ "import sys, platform\n",
+ "print(\"Working on a \", platform.system(), platform.release())\n",
+ "print(\"Python version \", sys.version)\n",
+ "print(\"Pandas version \", pd.__version__)\n",
+ "print(\"PV_ICE version \", PV_ICE.__version__)"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 4,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "WindowsPath('C:/Users/sayala/Documents/GitHub/PV_ICE/PV_ICE/baselines')"
+ ]
+ },
+ "execution_count": 4,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "baselinesFolder = Path().resolve().parent.parent.parent / 'PV_ICE' / 'baselines'\n",
+ "baselinesFolder"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## 1. Create Scenarios"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### A. Reading a standard ReEDS output file"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 6,
+ "metadata": {},
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Input file is stored in C:\\Users\\sayala\\Documents\\GitHub\\PV_ICE\\PV_ICE\\baselines\\SupportingMaterial\\December Core Scenarios ReEDS Outputs Solar Futures v3a.xlsx\n"
+ ]
+ }
+ ],
+ "source": [
+ "reedsFile = os.path.join(baselinesFolder, 'SupportingMaterial', 'December Core Scenarios ReEDS Outputs Solar Futures v3a.xlsx')\n",
+ "print (\"Input file is stored in %s\" % reedsFile)\n"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 7,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "REEDSInput = pd.read_excel(reedsFile,\n",
+ "# sheet_name=\"new installs PV (2)\")\n",
+ " sheet_name=\"new installs PV\")\n",
+ "\n",
+ "#index_col=[0,2,3]) #this casts scenario, PCA and State as levels\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### B. Save Input Files by PCA"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "#### Create a copy of the REEDS Input and modify structure for PCA focus"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 8,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/html": [
+ "\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " | \n",
+ " | \n",
+ " | \n",
+ " Capacity (GW) | \n",
+ "
\n",
+ " \n",
+ " | Scenario | \n",
+ " Year | \n",
+ " PCA | \n",
+ " | \n",
+ "
\n",
+ " \n",
+ " \n",
+ " \n",
+ " | 95-by-35.Adv | \n",
+ " 2010 | \n",
+ " p1 | \n",
+ " 0.000017 | \n",
+ "
\n",
+ " \n",
+ " | 2011 | \n",
+ " p1 | \n",
+ " 0.005158 | \n",
+ "
\n",
+ " \n",
+ " | 2012 | \n",
+ " p1 | \n",
+ " 0.005158 | \n",
+ "
\n",
+ " \n",
+ " | 2013 | \n",
+ " p1 | \n",
+ " 0.007146 | \n",
+ "
\n",
+ " \n",
+ " | 2014 | \n",
+ " p1 | \n",
+ " 0.007146 | \n",
+ "
\n",
+ " \n",
+ " | 2015 | \n",
+ " p1 | \n",
+ " 0.018253 | \n",
+ "
\n",
+ " \n",
+ " | 2016 | \n",
+ " p1 | \n",
+ " 0.018253 | \n",
+ "
\n",
+ " \n",
+ " | 2017 | \n",
+ " p1 | \n",
+ " 0.019460 | \n",
+ "
\n",
+ " \n",
+ " | 2018 | \n",
+ " p1 | \n",
+ " 0.019460 | \n",
+ "
\n",
+ " \n",
+ " | 2019 | \n",
+ " p1 | \n",
+ " 0.036320 | \n",
+ "
\n",
+ " \n",
+ " | 2020 | \n",
+ " p1 | \n",
+ " 0.036320 | \n",
+ "
\n",
+ " \n",
+ " | 2021 | \n",
+ " p1 | \n",
+ " 0.002434 | \n",
+ "
\n",
+ " \n",
+ " | 2022 | \n",
+ " p1 | \n",
+ " 0.002434 | \n",
+ "
\n",
+ " \n",
+ " | 2023 | \n",
+ " p1 | \n",
+ " 0.010066 | \n",
+ "
\n",
+ " \n",
+ " | 2024 | \n",
+ " p1 | \n",
+ " 0.010066 | \n",
+ "
\n",
+ " \n",
+ " | 2025 | \n",
+ " p1 | \n",
+ " 0.049445 | \n",
+ "
\n",
+ " \n",
+ " | 2026 | \n",
+ " p1 | \n",
+ " 0.049445 | \n",
+ "
\n",
+ " \n",
+ " | 2027 | \n",
+ " p1 | \n",
+ " 0.051342 | \n",
+ "
\n",
+ " \n",
+ " | 2028 | \n",
+ " p1 | \n",
+ " 0.051342 | \n",
+ "
\n",
+ " \n",
+ " | 2029 | \n",
+ " p1 | \n",
+ " 0.105931 | \n",
+ "
\n",
+ " \n",
+ " | 2030 | \n",
+ " p1 | \n",
+ " 0.105931 | \n",
+ "
\n",
+ " \n",
+ "
\n",
+ "