run_lct.m

一个研究声多普勒计程仪很好的工具箱

function run_lct

%run_lct gets the date, finds the folder for that date,
% and runs the lct functions to get the data and convert it into binary 
% At the end of run_lct batch processing is performed.
%get the date in the form yyyy/mm/dd for geobrowser
dat = floor(now);
[yesday,yyyy,mm,dd] = datmil(dat-1)

%create the new directory 
%(might need to add path for data)
%cd /usr2/users/jcote/NOPP
thePath = ['n' yyyy(3:4) mm dd]
while 0 %don't run path check now
if isequal(exist(thePath),7)
   break
else
   %check to make sure we haven't skipped a day
	[yesday,yyys,ms,ds] = datmil(dat-1);
	oldPath = [yyys(3:4) ms ds];
	if isequal(exist(oldPath),7)
   	break
	else
      thePath = oldPath;
      today = yesday;
	end

	makedir(thePath)
end
end
theFullPath = ['/usr2/users/jcote/NOPP/' thePath];
eval(['addpath ' theFullPath])
eval(['cd ' theFullPath])
disp(pwd)

%get the data and convert to binary
lct_get(yesday)
lct_conv

%these are the default (built - in) names 
%that are used when the data is downloaded
theCommandsFile = 'GB_commands.html';
thePayloadFile = 'GB_payload.html';
theTimeFile = 'GB_EIC_time.html';
theQualityFile = 'GB_rxqual.html';
theRDIFile = 'GB_payload.rdi';

%copy the file to a better name that can be carried on
theName = [thePath '.rdi'];
theFileName = fullfile(theFullPath,theName)
eval(['!cp -p ' theRDIFile ' ' theFileName])

%now run the processing in 2 batches
%need to change or regenerate the batch files for the corresponding date
bfile = [thePath '.bch']
bgfile = [thePath 'G.bch']
eval(['batch ' bfile ])
eval(['batch ' bgfile ])

quit