A = load('my_dates.txt'); for iLength = 1:56 if iLength < 10 label = sprintf('0%s',num2str(iLength)); else label=num2str(iLength); end {label} fn =sprintf('%s.txt',label); B =load(fn); [m n] = size(B); for ib = 1 : m n_date = B(ib); start_date = num2str(n_date); im = find(A == n_date); event_date = A(m-3:m+iLength+7); nm = length(event_date); for ie = 1:nm if ie < 10 event_label = sprintf('0%s',num2str(ie)); else event_label = num2str(ie); end n_event_date = event_date(ie); s_event_date = num2str(n_event_date); s_year = s_event_date(1:4); s_date = s_event_date(1:8); s_hour = s_event_date(9:10);
fn = sprintf('mslp_data/%s/mslp_%s_%s',s_year,s_date,s_hour); MSLP = load(fn); gn = sprintf('events1/%s/mslp_%s_%s',label,start_date,event_label); dlmwrite(gn,MSLP) end endend
Best Answer