CERBmathloop

function COPSmathloop(filename)
%mathloop that only works for prr
%filename = 'bsakrmP101216B.lcd.1'


domath(filename,'r','LuZ320','EdZ320',0,'rrs320')

domath(strcat('m',filename),'r','LuZ340','EdZ340',0,'rrs340')

domath(strcat('m',filename),'r','LuZ380','EdZ380',0,'rrs380')

domath(strcat('m',filename),'r','LuZ395','EdZ395',0,'rrs395')

domath(strcat('m',filename),'r','LuZ412','EdZ412',0,'rrs412')

domath(strcat('m',filename),'r','LuZ443','EdZ443',0,'rrs443')
domath(strcat('m',filename),'r','LuZ465','EdZ465',0,'rrs465')
%movefile(strcat('mm',filename),strcat('m',filename))
domath(strcat('m',filename),'r','LuZ490','EdZ490',0,'rrs490')
%movefile(strcat('mm',filename),strcat('m',filename))
domath(strcat('m',filename),'r','LuZ510','EdZ510',0,'rrs510')
domath(strcat('m',filename),'r','LuZ532','EdZ532',0,'rrs532')
%movefile(strcat('mm',filename),strcat('m',filename))
domath(strcat('m',filename),'r','LuZ555','EdZ555',0,'rrs555')
domath(strcat('m',filename),'r','LuZ560','EdZ560',0,'rrs560')
domath(strcat('m',filename),'r','LuZ625','EdZ625',0,'rrs625')
%movefile(strcat('mm',filename),strcat('m',filename))
domath(strcat('m',filename),'r','LuZ665','EdZ665',0,'rrs665')
%movefile(strcat('mm',filename),strcat('m',filename))
domath(strcat('m',filename),'r','LuZ670','EdZ670',0,'rrs670')
domath(strcat('m',filename),'r','LuZ683','EdZ683',0,'rrs683')
domath(strcat('m',filename),'r','LuZ710','EdZ710',0,'rrs710')
domath(strcat('m',filename),'r','LuZ780','EdZ780',0,'rrs780')
%command = (strcat('cp ',strcat('mmmmm',filename)))
%command  = strcat(command,strcat('m',filename));

% system('cp mmmmmbsakrmp101216B.lcd.1 mbsakrmP10121gB.lcd.1')
%  system(command)
%system('rm mmmmmbsakrmp101216B.lcd.1')
% system('rm mmmmbsakrmp101216B.lcd.1')
% system('rm mmmbsakrmp101216B.lcd.1')
% system('rm mmbsakrmp101216B.lcd.1')