TITLE: 
"Python Models underlying the master thesis:
Managing Water and Energy on Small Touristic Islands: 
study case Caye Chapel"
-------------------------------------------------------------------------------------
DESCRIPTION:
The collection contains twelve different models. There 
is one model for each of the twelve (12) alternatives for 
the water-energy system described in the thesis "Managing 
Water and Energy on Small Touristic Islands: study case Caye 
Chapel". The input data for each model is contained in the 
.txt files. 
-------------------------------------------------------------------------------------
FILES (MODEL):
WEM_<A#>.ipynb = Water Energy Model (WEM) for alterntative <A#>.
-------------------------------------------------------------------------------------
FILES (INPUT DATA):
Precipitation_1981_2015.txt = Hourly precipitation data (mm/h) from 1981 to 2015.
Precipitation_1992.txt = Hourly precipitation data (mm/h) from 1992.
Precipitation_2015.txt = Hourly precipitation data (mm/h) from 2015.

Evaportation_1981_2015.txt = Hourly evaportation data (mm/h) from 1981 to 2015.
Evaportation_1992.txt = Hourly evaportation data (mm/h) from 1992.
Evaportation_2015.txt = Hourly evaportation data (mm/h) from 2015.

Wind_1981_2015.txt = Hourly wind speed data (m/s) from 1981 to 2015.
Wind_1992.txt = Hourly wind speed (m/s) from 1992.
Wind_2015.txt = Hourly wind speed data (m/s) from 2015.

Solar_1981_2015.txt = Hourly solar irradiance data (W/m2) from 1981 to 2015.
Solar_1992.txt = Hourly solar irradiance data (W/m2) from 1992.
Solar_2015.txt = Hourly solar irradiance data (W/m2) from 2015.

Water_Hourly_Patterns.txt = Hourly water demand patterns.
Water_Monthly_Patterns.txt = Monthly water demand patterns.
-------------------------------------------------------------------------------------
METEOROLOGICAL DATA SOURCES:
The precipitation data was retreived from an hourly time-averaged 2-dimensional data 
collection in Modern-Era Retrospective analysis of Research and Applications version 
2(MERRA-2). The variable extracted from the data collection was the bias-corrected 
total precipitation (kg/m2/s). 
Source: Global Modeling and Assimilation Office (GMAO). (2015a). MERRA-2 
tavg1_2d_flx_Nx: 2d,1-Hourly,Time-Averaged,Single-Level,Assimilation,Surface Flux 
Diagnostics V5.12.4, Greenbelt, MD, USA, Goddard Earth Sciences Data and Information 
Services Center (GES DISC). 10.5067/7MCPBJ41Y0K6

The precipitation data was retreived from an hourly time-averaged 2-dimensional data 
collection in Modern-Era Retrospective analysis of Research and Applications version 
2(MERRA-2). The variable extracted from the data collection was surface evaporation 
(kg/m2/s). 
Source: Global Modeling and Assimilation Office (GMAO). (2008). tavg1_2d_flx_Nx: 
MERRA 2D IAU Diagnostic, Surface Fluxes, Time Average 1-hourly V5.2.0, Greenbelt, 
MD, USA, Goddard Earth Sciences Data and Information Services Center (GES DISC). 
10.5067/4EQ54AKI405R

The wind speed data was retreived from an hourly time-averaged 2-dimensional data 
collection in Modern-Era Retrospective analysis of Research and Applications version 
2(MERRA-2). The variable extracted from the data collection was the surface wind speed 
(m/s). The data represents the wind speed at 10 m above the surface. 
Source: Global Modeling and Assimilation Office (GMAO). (2015). MERRA-2 tavg1_2d_flx_Nx: 
2d,1-Hourly,Time-Averaged,Single-Level,Assimilation,Surface Flux Diagnostics V5.12.4, 
Greenbelt, MD, USA, Goddard Earth Sciences Data and Information Services Center 
(GES DISC). 10.5067/7MCPBJ41Y0K6


The solar irradiance data was retreived from an hourly time-averaged 2-dimensional 
data collection in Modern-Era Retrospective analysis of Research and Applications 
version 2(MERRA-2). The variable extracted from the data collection was the incident 
short wave land (W/m2). 
Source: Global Modeling and Assimilation Office (GMAO), “MERRA-2 tavg1_2d_rad_Nx: 
2d,1-Hourly,Time-Averaged,Single-Level,Assimilation,Radiation Diagnostics V5.12.4, 
Greenbelt, MD, USA, Goddard Earth Sciences Data and Information Services Center 
(GES DISC),” 2015. Accessed: Feb. 11, 2022. [Online]. Available: 10.5067/Q9QMY5PBNV1T
----------------------------------------------------------------------------------------------------------------------
MODELS' NOMENCLATURE:

NO.	NOMENCLATURE		DESCRIPTION									UNITS
GENERAL PARAMETER			
1	days			Number of days that will be modeled.						days
2	hours			Number of hours per day.							hours
INPUT			
3	water_demand_pot_daily	Average inhabitants' daily water demand.					m3/d
4	PF_hourly		Hourly peak factor.								adim
5	PF_monthly		Monthly peak factor.								adim
6	water_demand_pot_hourly	Hourly inhabitants' water demand.						m3/h
7	water_demand_irrig_daily	Average irrigation daily water demand.					m3/d
8	irrig_t_hours		Number of hours that the irrigation system operates per day.			hours
9	irrig_t_starting	Daily time at which the irrigation system starts to operate.			hours
10	water_demand_irrig_hourly	Hourly irrigation water demand.						m3/h
11	P			Precipitation hourly data.							m/h
12	Ev			Evaporation hourly data.							m/h
13	Ws			Wind speed hourly data @10 m relative to the surface.				m/s
14	Sr			Solar irradiance.								W/m2
(DWTP) DRINKING WATER TREATMENT PLANT			
15	seawater_t_production_starting	Daily time at which the desalination facility starts to operate.	hours
16	seawater_t_production_hours	Number of hours that the desalination facility operates per day.	hours
17	dwtp_Q_capacity		Maximum hourly capacity.							m3/h
(DWT) DRINKING WATER STORAGE TANK			
18	dwt_V_max		Maximum storage capacity.							m3
19	dwt_V_per_low		Indicator for low percentage of stored water.					adim
20	dwt_V_full		Available water in the water tank.						m3
21	dwt_V_per		Percentage of available water in the water tank.				%
22	dwt_V_empty		Available storage capacity in the water tank.					m3
23	dwt_Q_out_pot		Water flow going to the potable services.					m3/h
24	dwt_Q_out_irrig		Water flow going to the IWT.							m3/h
25	dwt_Q_in		Water flow coming into the tank. 						m3/h
(WWTP) WASTEWATER TREATMENT PLANT			
26	ww_factor		Percentage of drinking water that becomes wastewater.				%
27	wwtp_factor		Percentage of treated wastewater produced per cubic 
				meter of wastewater entering the WWTP.						%
28	delay_dwt_wwtp		Amount of timespets required for one cubic meter of 
				water to go from the DWT to the WWTP.						adim
29	delay_wwtp_wwt		Amount of timespets required for one cubic meter of 
				wastewater to be converted into treated wastewater.				adim
30	wwtp_Q_capacity		Maximum hourly flow capacity.							m3/h
31	wwtp_Q_in		Water flow coming in. 								m3/h
32	wwtp_Q_out		Water flow going out.								m3/h
33	wwtp_Q_rejected		Water flow rejected (do nt go inside the WWTP).					m3/h
(WWT) TREATED WASTEWATER STORAGE TANK			
34	wwt_V_max		Maximum storage capacity.							m3
35	wwt_V_min		Minimum storage capacity.							m3
36	wwt_Q_in		Water flow coming in. 								m3/h
37	wwt_Q_out_spill		Water flow spilled.								m3/h
38	wwt_Q_out_irrig		Water flow going to the IWT.							m3/h
39	wwt_V_full		Available water in the water tank.						m3
(WB) GOLF COURSE'S WATER BODIES			
40	wb_hh  			"The set level is defined as: wb_hh - wb_h_range.				m
41	wb_h_range		"The overflow level is defined as: wb_hh + wb_h_range"				m
42	wb_area			Surface area.									m
43	wb_h_min		Minimum water level.								m
44	wb_h			Water level.									m
45	wb_Q_in			Water flow coming in. 								m3/h
46	wb_Q_out		Water flow overflowed.								m3/h
47	wb_Q_in_demand		Water flow that is required to keep the set water level.			m3/h
(RWB) RAINWATER BUFFER TANK			
48	delay_wb_rwb		Amount of timespets required for one cubic meter of rainwater to 
				flow from the WB to the RWB.							adim
49	rwb_V_max		Maximum storage capacity.							m3
50	rwb_V_min		Minimum storage capacity.							m3
51	rwb_V_full		Available water in the water tank.						m3
52	rwb_V_empty		Available storage capacity in the water tank.					m3
53	rwb_Q_in		Water flow coming in. 								m3/h
54	rwb_Q_out_spill		Water flow spilled.								m3/h
55	rwb_Q_out_irrig		Water flow going the RWT.							m3/h
(RWTP) RAINWATER TREATMENT PLANT			
56	rwtp_factor		Percentage of treated rainwater produced per cubic meter 
				of rainwater entering the RWTP.							%
57	delay_rwtp_rwt		Amount of timespets required for one cubic meter of rainwater 
				to be converted into treated rainwater.						adim
58	rwtp_Q_capacity		Maximum hourly flow capacity.							m3/h
59	rwtp_Q_in		Water flow coming in. 								m3/h
60	rwtp_Q_out		Water flow going out.								m3/h
(RWT) RAIN WATER STORAGE TANK			
61	rwt_V_max		Maximum storage capacity.							m3
62	rwt_V_min		Minimum storage capacity.							m3
63	rwt_Q_in 		Water flow coming in. 								m3/h
64	rwt_Q_out_spill		Water flow spilled.								m3/h
65	rwt_Q_out_irrig		Water flow going to the IWT.							m3/h
66	rwt_V_full		Available water in the water tank.						m3
(IWT) IRRIGATION AND RECOVERY WATER STORAGE TANK			
67	iwt_V_max		Maximum storage capacity.							m3
68	iwt_V_min		Minimum storage capacity.							m3
69	iwt_V_per_low		Indicator for low percentage of stored water.					m3
70	iwt_V_per_med		Indicator for medium percentage of stored water.				m3
71	iwt_V_per_high		Indicator for high percentage of stored water.					m3
72	iwt_P_min		Minimum precipitation's treshold to shut down the irrigation system.		m
73	iwt_V_full		Available water in the water tank.						m3
74	iwt_V_empty		Available storage capacity in the water tank.					m3
75	iwt_V_per		Percentage of available water in the water tank.				m3
76	iwt_Q_in_wwt		Water flow coming in from the WWT. 						m3/h
77	iwt_Q_in_rwt		Water flow coming in from the RWT. 						m3/h
78	iwt_Q_in_dwt		Water flow coming in from the DWT. 						m3/h
79	iwt_Q_out_wb		Water flow going to the WB.							m3/h
80	iwt_Q_out_irrig		Water flow going to the irrigation system.					m3/h
81	iwt_Q_out_spill		Water flow spilled.								m3/h
PUMPING STATIONS			
82	dwt_Q_out_irrig_pump	Flow capacity of pumping station no. 1 						m3/h
83	wwt_Q_out_irrig_pump	Flow capacity of pumping station no. 4						m3/h
84	rwb_Q_out_irrig_pump	Flow capacity of pumping station no. 5						m3/h
85	rwt_Q_out_irrig_pump	Flow capacity of pumping station no. 6						m3/h
86	rwt_Q_out_irrig_pump	Flow capacity of pumping station no. 7						m3/h
87	iwt_Q_out_irrig_pump	Flow capacity of pumping station no. 8						m3/h
88	iwt_Q_out_wb_pump	Flow capacity of pumping station no. 9						m3/h
URBAN WATER CYCLE'S ELECTRICITY DEMAND			
89	pump_eff		Pumping station's efficiency.							%
90	pump_Q			Pumping station's hourly flow.							m3/h
91	pump_set		Number of pumps in parallel for each pumping station.				adim
92	pump1_E			Electricity demand from the pumping station no.1				kW/h
93	pump2_E			Electricity demand from the pumping station no.2				kW/h
94	pump3_E			Electricity demand from the pumping station no.3				kW/h
95	pump4_E			Electricity demand from the pumping station no.4				kW/h
96	pump5_E			Electricity demand from the pumping station no.5				kW/h
97	pump6_E			Electricity demand from the pumping station no.6				kW/h
98	pump7_E			Electricity demand from the pumping station no.7				kW/h
99	pump8_E			Electricity demand from the pumping station no.8				kW/h
100	pump9_E			Electricity demand from the pumping station no.9				kW/h
101	dwtp_E			Electricity demand from the DWTP.						kW/h
102	wwtp_E			Electricity demand from the WWTP.						kW/h
103	rwtp_E			Electricity demand from the RWTP.						kW/h
104	water_E_demand		Urban water cycle's electricity demand.						kW/h
(WT) WIND TURBINES			
105	wt_z			Hub's height.									m
106	wt_capacity		Rated capacity.									kWh
107	wt_quantity		Number of wind turbines.							adim
108	wt_E_out		Electricity delivered by the wind turbine(s).					kWh
(PV) PV PANELS			
109	pv_A			Area of PV panels.								m2
110	pv_yield		Yield.										%
111	pv_tf			Transposition factor.								adim
112	pv_pf			Performance factor.								%
113	pv_E_out		Electricity delivered by the PV panels.						kWh
(IRE) INTERMITTENT RENEWABLE ENERGY			
114	ire_E			Intermittent renewable energy.							kWh
115	ire_E_grid		Intermittent renewable energy going into the urban water cycle's grid.		kWh
116	ire_E_stor		Intermittent renewable energy going into the battery.				kWh
117	ire_E_rejected		Intermittent renewable energy rejected.						kWh
(BT) BATTERY			
118	b_eff			Battery's efficiency.								%
119	b_stor_capacity		Maximum storage capacity.							kW
120	b_in_capacity		Maximum inlet capacity.								kWh
121	b_out_capacity		Maximum outlet capacity.							kWh
122	b_E_full		Available energy stored in the battery.						kW
123	b_E_empty		Available storage capacity in the battery.					kW
124	b_E_out			Electricity delivered by the battery.						kWh
125	b_E_in			Energy coming into the battery.							kWh

