Domestic-load-modelling.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import numpy as np\n",
    "import math\n",
    "import random\n",
    "import array\n",
    "from numpy import zeros"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "house_inc = 0\n",
    "app_inc = 0\n",
    "inc = 0\n",
    "s=25\n",
    "m=0\n",
    "l=0\n",
    "healthcare=0\n",
    "school=2\n",
    "load_tot = []"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "hi=array.array('d',[])\n",
    "hour0_4=array.array('d',[])\n",
    "f = open(\"D:/HI.txt\")\n",
    "for i in f:\n",
    "    hi.append(float(i[0:6]))\n",
    "k=0\n",
    "hiavg=array.array('d',[])\n",
    "while k<8760:\n",
    "    dummy=sum(hi[k:k+4])\n",
    "    hiavg.append(round(dummy/4))\n",
    "    hiavg.append(round(dummy/4))\n",
    "    hiavg.append(round(dummy/4))\n",
    "    hiavg.append(round(dummy/4))\n",
    "    k+=4"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "for year in range (0,20):\n",
    "    l_fan = []\n",
    "    l_light = []\n",
    "    l_bulb = []\n",
    "    l_tv = []\n",
    "    l_socket = []\n",
    "    l_healthfan = []\n",
    "    l_healthlight = []\n",
    "    l_schoolfan = []\n",
    "    l_schoollight = []\n",
    "    load = []\n",
    "    healthfan=array.array('d',[])\n",
    "    healthlight=array.array('d',[])\n",
    "    schoolfan=array.array('d',[])\n",
    "    schoollight=array.array('d',[])\n",
    "    print ('year = ',year)\n",
    "    print ('small house = ',s)\n",
    "    print ('medium house = ',m)\n",
    "    print ('large house = ',l)\n",
    "    for iteration in range (int(s+m+2*l)):\n",
    "        bulb = [] \n",
    "        i=0\n",
    "        while i<8760:\n",
    "            if i%168<24 or i%168>71:\n",
    "                while i%24<4: #0-3                      \n",
    "                    if hiavg[i]<6:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    for p in range(4):\n",
    "                        bulb.append(hour0_4[p]) \n",
    "                    i+=1\n",
    "                while i%24<8 and i%24>3:     #4-7\n",
    "                    if hiavg[i]<3:\n",
    "                        hour0=math.floor(random.uniform(0,2))\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        for p in range(4):\n",
    "                            bulb.append(hour0_4[p])\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=math.floor(random.uniform(0,2))\n",
    "                        hour1=1\n",
    "                        bulb.append(hour0)\n",
    "                        bulb.append(hour1)\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                        for p in range(2):\n",
    "                            bulb.append(hour0_4[p])\n",
    "                    i+=1  \n",
    "                while i%24<12 and i%24>7:       #8-11\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=math.floor(random.uniform(0,2))\n",
    "                        hour1=0\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==3:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==4:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==5:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=1\n",
    "                    for p in range(4):\n",
    "                        bulb.append(hour0_4[p])  \n",
    "                while i%24<16 and i%24>11:            #12-15\n",
    "                    if hiavg[i]<4:\n",
    "                        hour0=0\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>3:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=1             \n",
    "                    for p in range(4):\n",
    "                        bulb.append(hour0_4[p])  \n",
    "                while i%24<20 and i%24>15:                 #16-19\n",
    "                    if hiavg[i]<3:\n",
    "                        hour0=0\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=1\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    i+=1      \n",
    "                    for p in range(4):\n",
    "                        bulb.append(hour0_4[p])  \n",
    "                while i%24<24 and i%24>19:                #20-23\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=1\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    if hiavg[i]>1:\n",
    "                        hour0=math.floor(random.uniform(0,2))\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    i+=1\n",
    "                    for p in range(4):\n",
    "                        bulb.append(hour0_4[p]) \n",
    "            else:\n",
    "                while i%24<4:                  #0-3\n",
    "                    if hiavg[i]<6:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    for p in range(4):\n",
    "                        bulb.append(hour0_4[p]) \n",
    "                    i+=1\n",
    "                while i%24<8 and i%24>3:            #4-7\n",
    "                    if hiavg[i]<3:\n",
    "                        hour0=math.floor(random.uniform(0,2))\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        for p in range(4):\n",
    "                            bulb.append(hour0_4[p])  \n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=math.floor(random.uniform(0,2))\n",
    "                        hour1=1\n",
    "                        bulb.append(hour0)\n",
    "                        bulb.append(hour1)\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                        for p in range(2):\n",
    "                            bulb.append(hour0_4[p])  \n",
    "                    i+=1\n",
    "                while i%24<12 and i%24>7:             #8-11\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=math.floor(random.uniform(0,2))\n",
    "                        hour1=0\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==3:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==4:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==5:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=1\n",
    "                    for p in range(4):\n",
    "                        bulb.append(hour0_4[p])  \n",
    "                while i%24<16 and i%24>11:                      #12-15\n",
    "                    if hiavg[i]<4:\n",
    "                        hour0=0\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>3:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=1             \n",
    "                    for p in range(4):\n",
    "                        bulb.append(hour0_4[p])  \n",
    "                while i%24<20 and i%24>15:                    #16-19\n",
    "                    if hiavg[i]<3:\n",
    "                        hour0=0\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=1\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    i+=1      \n",
    "                    for p in range(4):\n",
    "                        bulb.append(hour0_4[p])  \n",
    "                while i%24<24 and i%24>19:                  #20-23\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=1\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    if hiavg[i]>1:\n",
    "                        hour0=math.floor(random.uniform(0,2))\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    i+=1\n",
    "                    for p in range(4):\n",
    "                        bulb.append(hour0_4[p])\n",
    "        if iteration ==0:\n",
    "            for k in range (0,8760):\n",
    "                l_bulb.append(bulb[k])\n",
    "        else:\n",
    "            for k in range (0,8760):\n",
    "                l_bulb[k]=l_bulb[k]+bulb[k]\n",
    "\n",
    "    l_bulb = [i * 0.015 for i in l_bulb]\n",
    "    \n",
    "    for iteration in range (int(s+2*m+3*l)):\n",
    "        fan = []  \n",
    "        i=0\n",
    "        while i<8760:\n",
    "            if i%168<24 or i%168>71:\n",
    "                while i%24<4:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=0\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    for p in range(4):\n",
    "                        fan.append(hour0_4[p]) \n",
    "                    i+=4\n",
    "                while i%24<8 and i%24>3:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=0\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==3:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==4:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==5:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        fan.append(hour0_4[p])  \n",
    "                while i%24<12 and i%24>7:\n",
    "                    if hiavg[i]==1 or hiavg[i]==2:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==3:\n",
    "                        hour0=1\n",
    "                        hour1=0\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==4:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==5:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        fan.append(hour0_4[p])  \n",
    "                while i%24<16 and i%24>11:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2 or hiavg[i]==3:\n",
    "                        hour0=0\n",
    "                        hour1=1\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==4:\n",
    "                        hour0=0\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==5:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=4             \n",
    "                    for p in range(4):\n",
    "                        fan.append(hour0_4[p])  \n",
    "                while i%24<20 and i%24>15:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=0\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=1\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==3:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==4:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==5:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=4      \n",
    "                    for p in range(4):\n",
    "                        fan.append(hour0_4[p])  \n",
    "                while i%24<24 and i%24>19:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=0\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        fan.append(hour0_4[p]) \n",
    "            else:\n",
    "                while i%24<4:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=0\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    for p in range(4):\n",
    "                        fan.append(hour0_4[p]) \n",
    "                    i+=4\n",
    "                while i%24<8 and i%24>3:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=0\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==3:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==4:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==5:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        fan.append(hour0_4[p])  \n",
    "                while i%24<12 and i%24>7:\n",
    "                    if hiavg[i]==1 or hiavg[i]==2:\n",
    "                        hour0=0\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==3:\n",
    "                        hour0=1\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==4:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==5:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        fan.append(hour0_4[p])  \n",
    "                while i%24<16 and i%24>11:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=1\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2 or hiavg[i]==3:\n",
    "                        hour0=0\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==4:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==5:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=4             \n",
    "                    for p in range(4):\n",
    "                        fan.append(hour0_4[p])  \n",
    "                while i%24<20 and i%24>15:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=0\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=1\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==3:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==4:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==5:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=4      \n",
    "                    for p in range(4):\n",
    "                        fan.append(hour0_4[p])  \n",
    "                while i%24<24 and i%24>19:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=0\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        fan.append(hour0_4[p]) \n",
    "        if iteration ==0:\n",
    "            for k in range (0,8760):\n",
    "                l_fan.append(fan[k])\n",
    "        else:\n",
    "            for k in range (0,8760):\n",
    "                l_fan[k]=l_fan[k]+fan[k]\n",
    "\n",
    "    l_fan = [i * 0.075 for i in l_fan]\n",
    "    \n",
    "    for iteration in range (int(s+2*m+2*l)):\n",
    "        light = []\n",
    "        i=0\n",
    "        while i<8760:\n",
    "            if i%168<24 or i%168>71:\n",
    "                while i%24<4: #0-3\n",
    "                    if hiavg[i]<6:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    for p in range(4):\n",
    "                        light.append(hour0_4[p]) \n",
    "                    i+=1\n",
    "                while i%24<8 and i%24>3:     #4-7\n",
    "                    if hiavg[i]<3:\n",
    "                        hour0=0\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=0\n",
    "                        hour1=1\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    i+=1\n",
    "                    for p in range(4):\n",
    "                        light.append(hour0_4[p])  \n",
    "                while i%24<12 and i%24>7:       #8-11\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=math.floor(random.uniform(0,2))\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=math.floor(random.uniform(0,2))\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==3:\n",
    "                        hour0=math.floor(random.uniform(0,2))\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>3:\n",
    "                        hour0=math.floor(random.uniform(0,2))\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=1\n",
    "                    for p in range(4):\n",
    "                        light.append(hour0_4[p])  \n",
    "                while i%24<16 and i%24>11:            #12-15\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=0\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==3:\n",
    "                        hour0=0\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==4:\n",
    "                        hour0=0\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==5:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=1             \n",
    "                    for p in range(4):\n",
    "                        light.append(hour0_4[p])  \n",
    "                while i%24<20 and i%24>15:                 #16-19\n",
    "                    if hiavg[i]<3:\n",
    "                        hour0=0\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=1\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    i+=1      \n",
    "                    for p in range(4):\n",
    "                        light.append(hour0_4[p])  \n",
    "                while i%24<24 and i%24>19:                #20-23\n",
    "                    if hiavg[i]<3:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=1\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    i+=1\n",
    "                    for p in range(4):\n",
    "                        light.append(hour0_4[p]) \n",
    "            else:\n",
    "                while i%24<4:                  #0-3\n",
    "                    if hiavg[i]<6:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    for p in range(4):\n",
    "                        light.append(hour0_4[p]) \n",
    "                    i+=1\n",
    "                while i%24<8 and i%24>3:            #4-7\n",
    "                    if hiavg[i]<3:\n",
    "                        hour0=0\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=0\n",
    "                        hour1=1\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    i+=1\n",
    "                    for p in range(4):\n",
    "                        light.append(hour0_4[p])  \n",
    "                while i%24<12 and i%24>7:             #8-11\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=math.floor(random.uniform(0,2))\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=math.floor(random.uniform(0,2))\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==3:\n",
    "                        hour0=math.floor(random.uniform(0,2))\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==4:\n",
    "                        hour0=math.floor(random.uniform(0,2))\n",
    "                        hour1=0\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==5:\n",
    "                        hour0=math.floor(random.uniform(0,2))\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=1\n",
    "                    for p in range(4):\n",
    "                        light.append(hour0_4[p])  \n",
    "                while i%24<16 and i%24>11:                      #12-15\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=math.floor(random.uniform(0,2))\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==3:\n",
    "                        hour0=0\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>3:\n",
    "                        hour0=0\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=1             \n",
    "                    for p in range(4):\n",
    "                        light.append(hour0_4[p])  \n",
    "                while i%24<20 and i%24>15:                    #16-19\n",
    "                    if hiavg[i]<3:\n",
    "                        hour0=0\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=1\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    i+=1      \n",
    "                    for p in range(4):\n",
    "                        light.append(hour0_4[p])  \n",
    "                while i%24<24 and i%24>19:                  #20-23\n",
    "                    if hiavg[i]<3:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=1\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    i+=1\n",
    "                    for p in range(4):\n",
    "                        light.append(hour0_4[p])\n",
    "        if iteration ==0:\n",
    "            for k in range (0,8760):\n",
    "                l_light.append(light[k])\n",
    "        else:\n",
    "            for k in range (0,8760):\n",
    "                l_light[k]=l_light[k]+light[k]\n",
    "                            \n",
    "    l_light = [i * 0.040 for i in l_light]\n",
    "                            \n",
    "    for iteration in range (int(m+l)):\n",
    "        tv = []\n",
    "        i=0\n",
    "        while i<8760:\n",
    "            if i%168<24 or i%168>71:\n",
    "                while i%24<7:\n",
    "                    hour0=0\n",
    "                    hour1=0\n",
    "                    hour2=0\n",
    "                    hour3=0\n",
    "                    hour4=0\n",
    "                    hour5=0\n",
    "                    hour6=0\n",
    "                    hour0_4=[hour0,hour1,hour2,hour3,hour4,hour5,hour6]\n",
    "                    random.shuffle(hour0_4)\n",
    "                    for p in range(7):\n",
    "                        tv.append(hour0_4[p])\n",
    "                    i+=7\n",
    "                while i%24<13 and i%24>6:\n",
    "                    hour0=math.floor(random.uniform(0,2))\n",
    "                    hour1=0\n",
    "                    hour2=0\n",
    "                    hour3=0\n",
    "                    hour4=0\n",
    "                    hour5=0\n",
    "                    hour0_4=[hour0,hour1,hour2,hour3,hour4,hour5]\n",
    "                    i+=6\n",
    "                    for p in range(6):\n",
    "                        tv.append(hour0_4[p])\n",
    "                while i%24<18 and i%24>12:\n",
    "                    hour0=0\n",
    "                    hour1=1\n",
    "                    hour2=0\n",
    "                    hour3=math.floor(random.uniform(0,2))\n",
    "                    hour4=0\n",
    "                    hour0_4=[hour0,hour1,hour2,hour3,hour4]\n",
    "                    random.shuffle(hour0_4)\n",
    "                    i+=5\n",
    "                    for p in range(5):\n",
    "                        tv.append(hour0_4[p])\n",
    "                while i%24<22 and i%24>17:\n",
    "                    hour0=1\n",
    "                    hour1=1\n",
    "                    hour2=1\n",
    "                    hour3=math.floor(random.uniform(0,2))\n",
    "                    hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    random.shuffle(hour0_4)\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        tv.append(hour0_4[p])\n",
    "                while i%24<24 and i%24>21:\n",
    "                    hour0=0\n",
    "                    hour1=0\n",
    "                    hour0_4=[hour0,hour1]\n",
    "                    i+=2\n",
    "                    for p in range(2):\n",
    "                        tv.append(hour0_4[p])\n",
    "            else:\n",
    "                while i%24<7:\n",
    "                    hour0=0\n",
    "                    hour1=0\n",
    "                    hour2=0\n",
    "                    hour3=0\n",
    "                    hour4=0\n",
    "                    hour5=0\n",
    "                    hour6=0\n",
    "                    hour0_4=[hour0,hour1,hour2,hour3,hour4,hour5,hour6]\n",
    "                    random.shuffle(hour0_4)\n",
    "                    for p in range(7):\n",
    "                        tv.append(hour0_4[p])\n",
    "                    i+=7\n",
    "                while i%24<13 and i%24>6:\n",
    "                    hour0=math.floor(random.uniform(0,2))\n",
    "                    hour1=0\n",
    "                    hour2=1\n",
    "                    hour3=0\n",
    "                    hour4=0\n",
    "                    hour5=math.floor(random.uniform(0,2))\n",
    "                    hour0_4=[hour0,hour1,hour2,hour3,hour4,hour5]\n",
    "                    i+=6\n",
    "                    for p in range(6):\n",
    "                        tv.append(hour0_4[p])\n",
    "                while i%24<18 and i%24>12:\n",
    "                    hour0=0\n",
    "                    hour1=1\n",
    "                    hour2=0\n",
    "                    hour3=math.floor(random.uniform(0,2))\n",
    "                    hour4=1\n",
    "                    hour0_4=[hour0,hour1,hour2,hour3,hour4]\n",
    "                    random.shuffle(hour0_4)\n",
    "                    i+=5\n",
    "                    for p in range(5):\n",
    "                        tv.append(hour0_4[p])\n",
    "                while i%24<22 and i%24>17:\n",
    "                    hour0=1\n",
    "                    hour1=1\n",
    "                    hour2=1\n",
    "                    hour3=math.floor(random.uniform(0,2))\n",
    "                    hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    random.shuffle(hour0_4)\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        tv.append(hour0_4[p])\n",
    "                while i%24<24 and i%24>21:\n",
    "                    hour0=0\n",
    "                    hour1=0\n",
    "                    hour0_4=[hour0,hour1]\n",
    "                    i+=2\n",
    "                    for p in range(2):\n",
    "                        tv.append(hour0_4[p])\n",
    "                            \n",
    "        if iteration ==0:\n",
    "            for k in range (0,8760):\n",
    "                l_tv.append(tv[k])\n",
    "        else:\n",
    "            for k in range (0,8760):\n",
    "                l_tv[k]=l_tv[k]+tv[k]\n",
    "                            \n",
    "    l_tv = [i * 0.05 for i in l_tv]\n",
    "                            \n",
    "    for iteration in range (int(s+m+2*l)):\n",
    "        socket = []\n",
    "        i=0\n",
    "        while i<8760:\n",
    "            if i%168<24 or i%168>71:\n",
    "                while i%24<7:\n",
    "                    hour0=0\n",
    "                    hour1=0\n",
    "                    hour2=0\n",
    "                    hour3=0\n",
    "                    hour4=0\n",
    "                    hour5=0\n",
    "                    hour6=0\n",
    "                    hour0_4=[hour0,hour1,hour2,hour3,hour4,hour5,hour6]\n",
    "                    random.shuffle(hour0_4)\n",
    "                    for p in range(7):\n",
    "                        socket.append(hour0_4[p])\n",
    "                    i+=7\n",
    "                while i%24<13 and i%24>6:\n",
    "                    hour0=0\n",
    "                    hour1=0.007\n",
    "                    hour2=0.007\n",
    "                    hour3=0\n",
    "                    hour4=0\n",
    "                    hour5=0\n",
    "                    hour0_4=[hour0,hour1,hour2,hour3,hour4,hour5]\n",
    "                    i+=6\n",
    "                    for p in range(6):\n",
    "                        socket.append(hour0_4[p])\n",
    "                while i%24<18 and i%24>12:\n",
    "                    hour0=0\n",
    "                    hour1=0\n",
    "                    hour2=0\n",
    "                    hour3=0\n",
    "                    hour4=0\n",
    "                    hour0_4=[hour0,hour1,hour2,hour3,hour4]\n",
    "                    random.shuffle(hour0_4)\n",
    "                    i+=5\n",
    "                    for p in range(5):\n",
    "                        socket.append(hour0_4[p])\n",
    "                while i%24<22 and i%24>17:\n",
    "                    hour0=0\n",
    "                    hour1=0.007\n",
    "                    hour2=0.007\n",
    "                    hour3=math.floor(random.uniform(0,2))*0.007\n",
    "                    hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    random.shuffle(hour0_4)\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        socket.append(hour0_4[p])\n",
    "                while i%24<24 and i%24>21:\n",
    "                    hour0=0\n",
    "                    hour1=0\n",
    "                    hour0_4=[hour0,hour1]\n",
    "                    i+=2\n",
    "                    for p in range(2):\n",
    "                        socket.append(hour0_4[p])\n",
    "            else:\n",
    "                while i%24<7:\n",
    "                    hour0=0\n",
    "                    hour1=0\n",
    "                    hour2=0\n",
    "                    hour3=0\n",
    "                    hour4=0\n",
    "                    hour5=0\n",
    "                    hour6=0\n",
    "                    hour0_4=[hour0,hour1,hour2,hour3,hour4,hour5,hour6]\n",
    "                    random.shuffle(hour0_4)\n",
    "                    for p in range(7):\n",
    "                        socket.append(hour0_4[p])\n",
    "                    i+=7\n",
    "                while i%24<13 and i%24>6:\n",
    "                    hour0=math.floor(random.uniform(0,2))*0.1\n",
    "                    hour1=0\n",
    "                    hour2=1*0.007\n",
    "                    hour3=1*0.007\n",
    "                    hour4=0\n",
    "                    hour5=math.floor(random.uniform(0,2))*0.007\n",
    "                    hour0_4=[hour0,hour1,hour2,hour3,hour4,hour5]\n",
    "                    i+=6\n",
    "                    for p in range(6):\n",
    "                        socket.append(hour0_4[p])\n",
    "                while i%24<18 and i%24>12:\n",
    "                    hour0=0\n",
    "                    hour1=0\n",
    "                    hour2=0\n",
    "                    hour3=0\n",
    "                    hour4=0\n",
    "                    hour0_4=[hour0,hour1,hour2,hour3,hour4]\n",
    "                    random.shuffle(hour0_4)\n",
    "                    i+=5\n",
    "                    for p in range(5):\n",
    "                        socket.append(hour0_4[p])\n",
    "                while i%24<22 and i%24>17:\n",
    "                    hour0=1*0.007\n",
    "                    hour1=1*0.007\n",
    "                    hour2=0\n",
    "                    hour3=math.floor(random.uniform(0,2))*0.007\n",
    "                    hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    random.shuffle(hour0_4)\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        socket.append(hour0_4[p])\n",
    "                while i%24<24 and i%24>21:\n",
    "                    hour0=0\n",
    "                    hour1=0\n",
    "                    hour0_4=[hour0,hour1]\n",
    "                    i+=2\n",
    "                    for p in range(2):\n",
    "                        socket.append(hour0_4[p])\n",
    "                            \n",
    "        if iteration ==0:\n",
    "            for k in range (0,8760):\n",
    "                l_socket.append(socket[k])\n",
    "        else:\n",
    "            for k in range (0,8760):\n",
    "                l_socket[k]=l_socket[k]+socket[k]\n",
    "        \n",
    "    fridge=array.array('d',[])\n",
    "    i=0\n",
    "    while i<8760:\n",
    "        if i%2==0:\n",
    "            fridge.append(0.1)\n",
    "        else:\n",
    "            fridge.append(random.uniform(0.008,0.012))\n",
    "        i+=1                            \n",
    "    fridge = [i * l for i in fridge]\n",
    "    \n",
    "    for iteration in range (healthcare*2):\n",
    "        healthfan = []\n",
    "        i=0\n",
    "        while i<8760:\n",
    "            if i%168<24 or i%168>71:\n",
    "                while i%24<9:\n",
    "                    healthfan.append(0)\n",
    "                    i+=1\n",
    "                while i%24>8 and i%24<13:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    for p in range(4):\n",
    "                        healthfan.append(hour0_4[p]) \n",
    "                    i+=4\n",
    "                while i%24<17 and i%24>12:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        healthfan.append(hour0_4[p])  \n",
    "                while i%24<21 and i%24>16:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        healthfan.append(hour0_4[p])   \n",
    "                while i%24<24 and i%24>20 :\n",
    "                    healthfan.append(0)\n",
    "                    i+=1\n",
    "            else:\n",
    "                while i%24<9:\n",
    "                    healthfan.append(0)\n",
    "                    i+=1\n",
    "                while i%24>8 and i%24<13:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    for p in range(4):\n",
    "                        healthfan.append(hour0_4[p]) \n",
    "                    i+=4\n",
    "                while i%24<17 and i%24>12:\n",
    "                    hour0=0\n",
    "                    hour1=0\n",
    "                    hour2=0\n",
    "                    hour3=0\n",
    "                    hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    random.shuffle(hour0_4)        \n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        healthfan.append(hour0_4[p])  \n",
    "                while i%24<21 and i%24>16:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        healthfan.append(hour0_4[p])  \n",
    "                while i%24<24 and i%24>20:\n",
    "                    healthfan.append(0)\n",
    "                    i+=1\n",
    "        if iteration ==0:\n",
    "            for k in range (0,8760):\n",
    "                l_healthfan.append(healthfan[k])\n",
    "        else:\n",
    "            for k in range (0,8760):\n",
    "                l_healthfan[k]=l_healthfan[k]+healthfan[k]\n",
    "            \n",
    "    l_healthfan = [i * 0.075 for i in l_healthfan]\n",
    "    \n",
    "    for iteration in range (healthcare*2):\n",
    "        healthlight = []\n",
    "        i=0\n",
    "        while i<8760:\n",
    "            if i%168<24 or i%168>71:\n",
    "                while i%24<9:\n",
    "                    healthlight.append(0)\n",
    "                    i+=1\n",
    "                while i%24>8 and i%24<13:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    for p in range(4):\n",
    "                        healthlight.append(hour0_4[p]) \n",
    "                    i+=4\n",
    "                while i%24<17 and i%24>12:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        healthlight.append(hour0_4[p])  \n",
    "                while i%24<21 and i%24>16:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        healthlight.append(hour0_4[p])   \n",
    "                while i%24<24 and i%24>20 :\n",
    "                    healthlight.append(0)\n",
    "                    i+=1\n",
    "            else:\n",
    "                while i%24<9:\n",
    "                    healthlight.append(0)\n",
    "                    i+=1\n",
    "                while i%24>8 and i%24<13:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    for p in range(4):\n",
    "                        healthlight.append(hour0_4[p]) \n",
    "                    i+=4\n",
    "                while i%24<17 and i%24>12:\n",
    "                    hour0=0\n",
    "                    hour1=0\n",
    "                    hour2=0\n",
    "                    hour3=0\n",
    "                    hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    random.shuffle(hour0_4)        \n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        healthlight.append(hour0_4[p])  \n",
    "                while i%24<21 and i%24>16:\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        healthlight.append(hour0_4[p])  \n",
    "                while i%24<24 and i%24>20:\n",
    "                    healthlight.append(0)\n",
    "                    i+=1\n",
    "        if iteration ==0:\n",
    "            for k in range (0,8760):\n",
    "                l_healthlight.append(healthlight[k])\n",
    "        else:\n",
    "            for k in range (0,8760):\n",
    "                l_healthlight[k]=l_healthlight[k]+healthlight[k]\n",
    "            \n",
    "    l_healthlight = [i * 0.040 for i in l_healthlight]\n",
    "    \n",
    "    for iteration in range (school):\n",
    "        schoolfan = []\n",
    "        i=0\n",
    "        while i<8760:\n",
    "            if i%168<24 or i%168>71:\n",
    "                while i%24<4: #0-3\n",
    "                    if hiavg[i]<6:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    for p in range(4):\n",
    "                        schoolfan.append(hour0_4[p]) \n",
    "                    i+=4\n",
    "                while i%24<8 and i%24>3:     #4-7\n",
    "                    if hiavg[i]<6:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    for p in range(4):\n",
    "                        schoolfan.append(hour0_4[p])\n",
    "                    i+=4  \n",
    "                while i%24<12 and i%24>7:       #8-11\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        schoolfan.append(hour0_4[p])  \n",
    "                while i%24<16 and i%24>11:            #12-15\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=1\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=1\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    i+=4             \n",
    "                    for p in range(4):\n",
    "                        schoolfan.append(hour0_4[p])  \n",
    "                while i%24<20 and i%24>15:                 #16-19\n",
    "                    if hiavg[i]<6:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    i+=4      \n",
    "                    for p in range(4):\n",
    "                        schoolfan.append(hour0_4[p])  \n",
    "                while i%24<24 and i%24>19:                #20-23\n",
    "                    if hiavg[i]<6:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        schoolfan.append(hour0_4[p]) \n",
    "            else:  \n",
    "                schoolfan.append(0) \n",
    "                i+=1\n",
    "                \n",
    "        if iteration ==0:\n",
    "            for k in range (0,8760):\n",
    "                l_schoolfan.append(schoolfan[k])\n",
    "        else:\n",
    "            for k in range (0,8760):\n",
    "                l_schoolfan[k]=l_schoolfan[k]+schoolfan[k]\n",
    "            \n",
    "    l_schoolfan = [i * 0.075 for i in l_schoolfan]\n",
    "        \n",
    "    for iteration in range (school):\n",
    "        schoollight = []\n",
    "        i=0\n",
    "        while i<8760:\n",
    "            if i%168<24 or i%168>71:\n",
    "                while i%24<4:                  #0-3\n",
    "                    if hiavg[i]<6:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    for p in range(4):\n",
    "                        schoollight.append(hour0_4[p]) \n",
    "                    i+=4\n",
    "                while i%24<8 and i%24>3:     #4-7\n",
    "                    if hiavg[i]<6:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    for p in range(4):\n",
    "                        schoollight.append(hour0_4[p])\n",
    "                    i+=4  \n",
    "                while i%24<12 and i%24>7:       #8-11\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=math.floor(random.uniform(0,2))\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=1\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=math.floor(random.uniform(0,2))\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        schoollight.append(hour0_4[p])  \n",
    "                while i%24<16 and i%24>11:            #12-15\n",
    "                    if hiavg[i]==1:\n",
    "                        hour0=1\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=1\n",
    "                        hour3=1\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]==2:\n",
    "                        hour0=1\n",
    "                        hour1=1\n",
    "                        hour2=0\n",
    "                        hour3=math.floor(random.uniform(0,2))\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    if hiavg[i]>2:\n",
    "                        hour0=0\n",
    "                        hour1=math.floor(random.uniform(0,2))\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                        random.shuffle(hour0_4)\n",
    "                    i+=4            \n",
    "                    for p in range(4):\n",
    "                        schoollight.append(hour0_4[p])  \n",
    "                while i%24<20 and i%24>15:                 #16-19\n",
    "                    if hiavg[i]<6:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    i+=4    \n",
    "                    for p in range(4):\n",
    "                        schoollight.append(hour0_4[p])  \n",
    "                while i%24<24 and i%24>19:                #20-23\n",
    "                    if hiavg[i]<6:\n",
    "                        hour0=0\n",
    "                        hour1=0\n",
    "                        hour2=0\n",
    "                        hour3=0\n",
    "                        hour0_4=[hour0,hour1,hour2,hour3]\n",
    "                    i+=4\n",
    "                    for p in range(4):\n",
    "                        schoollight.append(hour0_4[p]) \n",
    "            else:  \n",
    "                schoollight.append(0) \n",
    "                i+=1 \n",
    "        if iteration ==0:\n",
    "            for k in range (0,8760):\n",
    "                l_schoollight.append(schoollight[k])\n",
    "        else:\n",
    "            for k in range (0,8760):\n",
    "                l_schoollight[k]=l_schoollight[k]+schoollight[k]\n",
    "            \n",
    "    l_schoollight = [i * 0.040 for i in l_schoollight]\n",
    "        \n",
    "    \n",
    "    if (healthcare == 0):\n",
    "        l_healthlight = [0]*8760\n",
    "        l_healthfan = [0]*8760\n",
    "        \n",
    "    if (school == 0):\n",
    "        l_schoollight = [0]*8760\n",
    "        l_schoolfan = [0]*8760\n",
    "    \n",
    "    if (l == 0 and m == 0):\n",
    "        l_tv = [0]*8760\n",
    "\n",
    "    load = [sum(x) for x in zip(l_fan,l_light,l_bulb,l_tv,l_socket,fridge,l_healthlight,l_healthfan,l_schoolfan,l_schoollight)]\n",
    "    \n",
    "    if year>0: #al_fan,l_light,l_bulb,l_tv,l_socket,fridge,l_healthlight,l_healthfan,l_schoolfan,l_schoollightnd year<16:\n",
    "        load = [x+z for x,z in zip(load_inc,load)]\n",
    "    \n",
    "    load_tot = load_tot + load\n",
    "    \n",
    "    data = [(l_fan[i],l_light[i],l_bulb[i],l_tv[i],l_socket[i],fridge[i],l_healthlight[i],l_healthfan[i],l_schoolfan[i],l_schoollight[i],load[i]) for i in range (0,8760)]\n",
    "    df = pd.DataFrame(data = data,columns =['fan','light','bulb','tv','socket','fridge','healthcare fan','healthcare light','school fan','school light','total load'])\n",
    "    \n",
    "    df.to_csv('load_madanpur_'+ str(year) +'.csv',index=False,encoding='utf-8')\n",
    "    \n",
    "    if year<16:\n",
    "        app_inc = int(random.uniform(4,7))\n",
    "        house_inc = 10 - app_inc\n",
    "\n",
    "        if house_inc==4:\n",
    "            h4=int(random.uniform(1,5))\n",
    "            if h4==1 or h4==2 or h4==3:\n",
    "                s=s+1\n",
    "            elif h4==4:\n",
    "                s=s+2\n",
    "        if house_inc==5:\n",
    "            h5=int(random.uniform(1,5))\n",
    "            if h5==1 or h5==2 or h5==3:\n",
    "                s=s+2\n",
    "            elif h5==4:\n",
    "                m=m+1\n",
    "        if house_inc==6:\n",
    "            h5=int(random.uniform(1,5))\n",
    "            if h5==1 or h5==2 or h5==3:\n",
    "                m=m+1\n",
    "                s=s+1\n",
    "            elif h5==4:\n",
    "                l=l+1\n",
    "\n",
    "        if app_inc==4:\n",
    "            inc = random.uniform(3,5)\n",
    "        if app_inc==5:\n",
    "            inc = random.uniform(4,6)\n",
    "        if app_inc==6:\n",
    "            inc = random.uniform(5,7)\n",
    "\n",
    "        print ('appliance increase = ',inc)\n",
    "         \n",
    "    load_inc = []\n",
    "    load_inc = [i * (inc/100) for i in load] "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "data = [(load_tot[i]) for i in range (0,175200)]\n",
    "df = pd.DataFrame(data = data,columns =['load'])\n",
    "df.to_csv('load_total_madanpur.csv',index=False,encoding='utf-8')"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}