Re_lookup.m

function data = Re_lookup(Re,Cl,ask,data_mat)
%-------------------------------------------------
%Cd = f(Re,Cl), alpha = f(Re,Cl)&#12434;&#27714;&#12417;&#12427;&#12383;&#12417;&#12398;&#38306;&#25968;
%ask = 1&#12398;&#26178;&#12399;Cd, ask = 2&#12398;&#26178;&#12399;alpha
%
%ZI = interp2(X,Y,Z,XI,YI,method)
%(XI,YI,ZI)&#35036;&#38291;&#12373;&#12428;&#12427;,&#27714;&#12417;&#12427;&#22580;&#25152;
%X=1:n, Y=1:m, size(Z)=[m,n]
%X = Relist, Y = Cllist, Z = data_mat&#12398;ask&#34892;&#21015;
%XI = Re, YI = Cl, ZI = data
%interp1&#12398;default&#12391;&#12399;&#22806;&#25407;&#12364;&#12391;&#12365;&#12378;&#12395;NA&#12364;&#20986;&#12427;&#12371;&#12392;&#12364;&#12354;&#12427;&#12398;&#12391;'cubic'&#12395;&#12377;&#12427;&#65288;&#35336;&#31639;&#26178;&#38291;&#12364;&#12363;&#12363;&#12427;)
%-------------------------------------------------

Relist = [10000 20000 30000 40000 50000 100000 	150000 200000 250000 300000 350000 400000];
Cllist = linspace(0,1.5,100);
if ask == 1
	for i = 1:12
		% Cd_mat_new(i,:) = interp1(data_mat(:,2,i),data_mat(:,3,i),Cllist);
		Cd_mat_new(i,:) = interp1(data_mat(:,2,i),data_mat(:,3,i),Cllist,'cubic',0);		
	end
	data = interp2(Relist,Cllist,Cd_mat_new',Re,Cl);
elseif ask == 2
	for i = 1:12
		% alpha_mat_new(i,:) = interp1(data_mat(:,2,i),data_mat(:,1,i),Cllist);
		alpha_mat_new(i,:) = interp1(data_mat(:,2,i),data_mat(:,1,i),Cllist,'cubic',0);
	end
	data = interp2(Relist,Cllist,alpha_mat_new',Re,Cl);
end