1 view (last 30 days)

format long

station = input('Enter the station you want data for with out the K identifier (BOS not KBOS) ','s');

Day1 = input('Enter first day ','s'); %enter as without leading zero ex: 1 not 01

Month1 = input('Enter first month ','s'); %enter as without leading zero ex: 1 not 01

Year1 = input('Enter first year ','s'); %enter as without leading zero ex: 1 not 01

Day2 = input('Enter second day ','s'); %enter as without leading zero ex: 1 not 01

Month2 = input('Enter second month ','s'); %enter as without leading zero ex: 1 not 01

Year2 = input('Enter second year ','s'); %enter as without leading zero ex: 1 not 01

Date1 = strcat(Month1, '-',Day1,'-',Year1);

Date2 = strcat(Month2, '-',Day2,'-',Year2);

DateNum1 = datenum(Date1);

DateNum2 = datenum(Date2);

t1 = datevec(DateNum1,'mmmm dd, yyyy HH:MM:SS');

t2 = datevec(DateNum2,'mmmm dd, yyyy HH:MM:SS');

MinBetDat = ((etime(t2,t1))/60) +1440; %+1440 to count end date

Dataurl = strcat('https://mesonet.agron.iastate.edu/cgi-bin/request/asos.py?station=',station,'&data=skyc1&data=skyc2&data=skyc3&year1=',Year1,'&month1=',Month1,'&day1=',Day1,'&year2=',Year2,'&month2=',Month2,'&day2=',Day2,'&tz=Etc%2FUTC&format=onlycomma&latlon=no&missing=M&trace=T&direct=no&report_type=1&report_type=2');

% Reference URL: 'https://mesonet.agron.iastate.edu/cgi-bin/request/asos.py?station=BOS&data=skyc1&data=skyc2&data=skyc3&year1=2019&month1=1&day1=1&year2=2019&month2=4&day2=8&tz=Etc%2FUTC&format=onlycomma&latlon=no&missing=M&trace=T&direct=no&report_type=1&report_type=2'

str = urlread(Dataurl); %pulls data from url as a single string

DATES1 = []; %blank date array

CLOUDS = []; %Blank array to separate the sky condition into, NOT FUNCTIONAL YET

Data = strsplit(str,{'\n',','}); %parses sinlge url string first separating by each new line and then by commas such that each element has it's own indexed number

L = length(Data); % probably not necessary, but I'm not going to change it

Data = Data(6:end); %removes column headers

Data(1:5:end) = []; %removes BOS identifier in each line

k=1; %leftover from previous versions

DATES1 = Data(1:4:end); %separates dates from data into their own array, makes referencing the dates easier in the following lines

k = length(DATES1);

Days = MinBetDat/1440;

DATES1 = datenum(DATES1);

DATES1 = datevec(DATES1);

D = ones(1,length(DATES1));

DATES2 = mat2cell(DATES1,D);

k = 6; %index of second line so that it actually has something to compare against.

t = 1;

tO = 1;

tB = 1;

tS = 1;

tF = 1;

tV = 1;

tC = 1;

time = zeros(1,length(DATES2));

timeOVC = zeros(1,length(DATES2));

timeBKN = zeros(1,length(DATES2));

timeSCT = zeros(1,length(DATES2));

timeFEW = zeros(1,length(DATES2));

timeVV = zeros(1,length(DATES2));

timeCLR = zeros(1,length(DATES2));%Faster to preallocate the times array

while k<=length(Data)

if any(strcmp(Data(k:k+2), 'OVC')) || any(strcmp(Data(k:k+2), 'BKN')) || any(strcmp(Data(k:k+2), 'SCT')) || any(strcmp(Data(k:k+2), 'FEW')) || any(strcmp(Data(k:k+2), 'VV')) || any(strcmp(Data(k:k+2), 'CLR'))%checks if the current line contains any cloud conditions

if any(strcmp(Data(k-4:k-2), 'OVC')) %checks if the line above had any cloud conditions, if it does then it checks which sky condition was in the lines above and adds the time difference to that respective sky con matrix

timeOVC(tO)= ((etime(DATES2{t+1},DATES2{t})/60));

t = t+1;

tO = tO+1;

elseif any(strcmp(Data(k-4:k-2), 'BKN'))

timeBKN(tB)= ((etime(DATES2{t+1},DATES2{t})/60));

t = t+1;

tB = tB + 1;

elseif any(strcmp(Data(k-4:k-2), 'SCT'))

timeSCT(tS)= ((etime(DATES2{t+1},DATES2{t})/60));

t = t+1;

tS = tS + 1;

elseif any(strcmp(Data(k-4:k-2), 'FEW'))

timeFEW(tF)= ((etime(DATES2{t+1},DATES2{t})/60));

t = t+1;

tF = tF + 1;

elseif any(strcmp(Data(k-4:k-2), 'CLR'))

timeCLR(tC) = ((etime(DATES2{t+1},DATES2{t})/60));

t = t+1;

tC = tC + 1;

elseif any(strcmp(Data(k-4:k-2), "VV"))

timeVV(tV)= ((etime(DATES2{t+1},DATES2{t})/60));

t = t+1;

tV = tV + 1;

else

time(t) = 0;

t = t+1;

end

elseif any(strcmp(Data(k:k+2), 'OVC')) || any(strcmp(Data(k:k+2), 'BKN')) || any(strcmp(Data(k:k+2), 'SCT')) || any(strcmp(Data(k:k+2), 'FEW')) || any(strcmp(Data(k:k+2), 'VV')) || any(strcmp(Data(k:k+2), 'CLR'))%checks if the current line contains any cloud conditions

if any(strcmp(Data(k-4:k-2), 'OVC')) %checks if the line above had any cloud conditions, if it does then it checks which sky condition was in the lines above and adds the time difference to that respective sky con matrix

timeOVC(tO)= ((etime(DATES2{t+1},DATES2{t})/60));

t = t+1;

tO = tO+1;

elseif any(strcmp(Data(k-4:k-2), 'BKN'))

timeBKN(tB)= ((etime(DATES2{t+1},DATES2{t})/60));

t = t+1;

tB = tB + 1;

elseif any(strcmp(Data(k-4:k-2), 'SCT'))

timeSCT(tS)= ((etime(DATES2{t+1},DATES2{t})/60));

t = t+1;

tS = tS + 1;

elseif any(strcmp(Data(k-4:k-2), 'FEW'))

timeFEW(tF)= ((etime(DATES2{t+1},DATES2{t})/60));

t = t+1;

tF = tF + 1;

elseif any(strcmp(Data(k-4:k-2), 'CLR'))

timeCLR(tC) = ((etime(DATES2{t+1},DATES2{t})/60));

t = t+1;

tC = tC + 1;

elseif any(strcmp(Data(k-4:k-2), "VV"))

timeVV(tV) = ((etime(DATES2{t+1},DATES2{t})/60));

t = t+1;

tV = tV + 1;

else

time(t) = 0;

t = t+1;

end

end

progress_clouds = (k/length(Data))* 100

k = k+4;

end

OVC = sum(timeOVC)

BKN = sum(timeBKN)

SCT = sum(timeSCT)

FEW = sum(timeFEW)

VV = sum(timeVV)

CLR = sum(timeCLR)

totaltime = OVC + BKN + SCT + FEW + VV;

CloudyTime = (totaltime/MinBetDat)*100

This is my code, this program is designed to pull weather observations from a user generated url and parse the data down to only sky conditions and then calculate the total time spent at each condition when it is the worst one present. When I run it I can get values for all of the sky conditions except VV which is always returned as 0? I've checked it againts an excel sheet with the same data and there are several instances where VV is the worst but it doesn't perform the calculation. Does anyone have any insight as to why it won't pick it up?

Geoff Hayes
on 3 Jun 2019

Edited: Geoff Hayes
on 3 Jun 2019

Joshua - from this code

strcmp(Data(k-4:k-2), "VV"))

you are comparing the two character string "VV" with a three character string. If k is 5 then you are extracting Data(1,3) which is three characters long...and so will never be identical to "VV". Either you need to ensure that you only extract two characters from Data or you will could use strfind or contains to see if the shorter string is contained within the longer. (Or use any of the regular expressions to do the same)

EDIT

Data is a cell array of strings and so Data(k-4:k-2) is not a single string but three strings.

Geoff Hayes
on 3 Jun 2019

Looking at the data, the "VV" string seems to really be "VV ". Have you tried adding an extra space to your code to see if you can string compare on that?

any(strcmp(Data(k-4:k-2), "VV "))

Opportunities for recent engineering grads.

Apply TodayFind the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!
## 0 Comments

Sign in to comment.