Thread Subject: a simple control flow to find the index and replace

Dear all,

Here is my code to supplement missing data in the specific time based on the index.

Can you advice me or provide me some hints how to write the part of control flow? I already try several way, but I need to pick up some clues.

Thanks so much.
Michael

load ('data_A_2007');
load ('data_B_2007');
load ('data_C_2007');

year = 2007;
doy=yeardays(year);
[L W] = size(year);

if L < W;
year = year';
[L W] = size(year);
end

ind = [1:1:L]';

start_time = datenum(year(ind(1)), 1, 1) + step/1440;
end_time = datenum(year(ind(L))+1, 1, 1);

time_vector = [start_time : step/1440 : end_time]';
base_vector = floor(time_vector(1,:));
day_vector = time_vector - base_vector - doy;

% summer time : doy >=91 and <=275; night time hour: 6 - 19
% winter time : doy <90 and >275; night time hour: 7 - 18

for day_vector(1:90*24*2 & 275*24*2:end)
night_index = find(floor(day_vector)< 7/24 | floor(day_vector)> 18/24);
else
night_index = find(floor(day_vector)< 6/24 | floor(day_vector)> 19/24);
end

night_missingdata_index = find(data_A_2007(night_index) & isnan(data_B_2007));
data_A_2007(night_missingdata_index) = data_C_2007(night_missingdata_index);

Hi all, are there any suggestions and tips for this post? Thank so much. Michael

"Kuo-Hsien"
> Hi all, are there any suggestions and tips for this post...

you need to curtail your (tedious) idiosyncratic snippet (including files that are not accessible to CSSMers) to the very (bare bone) problem - in easily copy/pasteable code - you want this NG to look at...

us

Thanks for your helps. I should make it more clearer and readable.

Here is my two problems:
(1) how to do some judgement based on a specific range of index?
(2) how to approach my purpose by utilizing if ... end or while ... end?

For example, I create a series of index<8760x1 double>, I like to sort out the index based on the conditions: (a) location from 1 to 1000, 8000 to the end and the rest; (b) the decimal value from 0 to 0.3, 0.7 to 1 and the rest.

index = repmat((0:1/24:0.9584)',[365 1]);
value = rand(1,8760)';

if index(1:1000 & 8000:end)
index_tail=find(floor(index)<0.3 | floor(index)>0.7)
else
index_mid=find(floor(index)>=0.3 | floor(index)<=0.7)
if index(1001:7999)
index_tail=find(floor(index)<0.3 | floor(index)>0.7)
else
index_mid=find(floor(index)>=0.3 | floor(index)<=0.7)
end
index_final = find((index_tail) & value>0.8);

Please advice me how to correct my code. If I did not make my question clear or misunderstand the logics, please feel free to let me know. Thanks for your advices.

Michael

Thanks for your helps. I should make it more clearer and readable.

Here is my two problems:
(1) how to do some judgement based on a specific range of index?
(2) how to approach my purpose by utilizing if ... end or while ... end?

For example, I create a series of index<8760x1 double>, I like to sort out the index based on the conditions: (a) location from 1 to 1000, 8000 to the end and the rest; (b) the decimal value from 0 to 0.3, 0.7 to 1 and the rest.

index = repmat((0:1/24:0.9584)',[365 1]);
value = rand(1,8760)';

if index(1:1000 & 8000:end)
index_tail=find(floor(index)<0.3 | floor(index)>0.7)
else
index_mid=find(floor(index)>=0.3 | floor(index)<=0.7)
if index(1001:7999)
index_tail=find(floor(index)<0.3 | floor(index)>0.7)
else
index_mid=find(floor(index)>=0.3 | floor(index)<=0.7)
end
index_final = find((index_tail) & value>0.8);

Please advice me how to correct my code. If I did not make my question clear or misunderstand the logics, please feel free to let me know. Thanks for your advices.

Michael