-
Notifications
You must be signed in to change notification settings - Fork 0
/
Main.m
38 lines (24 loc) · 1.91 KB
/
Main.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
% -------------------------------------------------------------------------
% Main file for processing data and running simulation
%
% Instructions:
% Run the following functions sequentially. They will load, process, and
% plot the data. The final function will use the information to run a
% simulation over the time horizon 2022-2080 for new wind capacity for
% the United States.
% -------------------------------------------------------------------------
clc
clearvars
latlim = [ 23, 50]; % Latitude limits for bounding box
lonlim = [-127,-65]; % Longitude limits for bounding box
[grid,gridR] = FullArea(150, latlim, lonlim); % Process reference background data for full US area
[wind,windR] = WindMap(150, latlim, lonlim, grid); % Process wind speed measurements at 110-m above ground level
[crops,cropsR] = CropArea(150, latlim, lonlim); % Process agricultural development data
[water,waterR] = WaterArea(150, latlim, lonlim); % Process water area data
[city,cityR] = CityArea(150, latlim, lonlim); % Process urban development data
[forest,forestR] = ForestArea(150, latlim, lonlim, grid); % Process forest area data
[elev,elevR] = ElevationMap(150, latlim, lonlim, grid); % Process geospatial elevation data
[comb,combW] = CombineAreas(wind,crops,water,city,forest,elev); % Combine map layers
[~,~,Lats,Lons] = SumArea(grid, gridR, combW); % Find remaining areas available for wind farm development
[Outputs] = PowerCalcs(Lons, Lats, grid, latlim, lonlim); % Simulate new wind farm development
MakePlots(grid,wind,crops,water,city,forest,elev,comb,combW,gridR,latlim, lonlim)