clear T=load('S3_T_region.txt'); row_T=size(T,1); column_T=size(T,2); %2 means the second dimension row_extract=43; colu_extract=51; T_region=zeros((row_extract-1),(colu_extract-1)); for i=1:(row_extract-1) for j=1:(colu_extract-1) T_region(i,j)=T(i,j); end end [X,Y]=meshgrid(linspace(0,300,(colu_extract-1)), linspace(250,0,(row_extract-1))); %surface(X,Y,T,'EdgeColor','none') surface(X,Y,T_region,'EdgeColor','none') ylabel('Height [km]') xlabel('Off-trench distance [km]') %contour(X,Y,T,6,'LineWidth',1,'LineColor','b','showtext','on') colormap jet %cool cbh=colorbar('southoutside'); cbh.Label.String = 'Temperature [degree]';