clear all
close all
clc

% lonlims = [-79.52 -69.48]; % MABSB (degree)
% latlims = [35.37 42.52]; % MABSB (degree)
% depthlims = [0 500]; % in meters
% X1d=linspace(lonlims(1),lonlims(2),5); % (1, nx) row: 1 col: nx
% Y1d=linspace(latlims(1),latlims(2),10); % (1, ny)
% Z1d=linspace(depthlims(1),depthlims(2),15); % (1, nz)
% T1d=linspace(0,10,20); % (1, nt)
% Data=rand(5,10,15,20); % (nx, ny, nz, nt)
% pcolor(X1d,Y1d,Data(:,:,1,1))

run = 'run7';
[X1d,Y1d,Z1d,T1d,Data]=prep_lcs_data(run);
lonlims = [min(X1d) max(X1d)];
latlims = [min(Y1d) max(Y1d)];
statebdryfile    ='/Users/seroka/Research/asset_map/WindPower_ShapeFiles/FederalStateBoundary.dat';
boemfile         ='/Users/seroka/Research/asset_map/WindPower_ShapeFiles/2013_boem_osw_call_area.csv';
transectfile     ='/Users/seroka/Research/code/fy20/lcs/hovmoller/db_transect.csv';
longitude_coastal_roi_vertices = []; % longitude-coordinates of polygon vertices
latitude_coastal_roi_vertices = []; % latitude-coordinates of polygon vertices
wrfnc = fullfile(pwd,'wrfout_d02_2017-04-05_17:00:00_r.nc');
[transect]=hovmoller(X1d,Y1d,Z1d,T1d,Data,'lonlims',lonlims,'latlims',latlims,'transectfile',transectfile,'wrfout',wrfnc);

disp('that''s all!!!')
