matrix points from a plot

Hello!

I got this code below thtr plot some profiles and her maximun value.

I want to have a matrix with the location of the max values for plot it with other varibles.

I mean something like:

Profile max_value depth_location

1 x.xx 15cm

2 x.xxx 10cm

3 x.xxx 20cm

....

[depht_maxvalue] = 15 10 20 ...

load 'biom.txt' % (Fortran file output)
tl=201; % Number of depth levels
[m,n]=size(biom); % derermining the size of the data
mm=m/tl; % total number of profiles
% sm=mm; % tottal profiles (for last profile); sm=mm-19; $for all profiles
% maxDepth = -400;

x= 1e6*biom; % <-- here variable concentration    microMole => 1e9 miliMol=>1e6
% nexttile %biom
% for i=sm:mm
for i = 1:mm
    
%     d=maxDepth ;
    hold on 
    
    plot(( x  ((i-1)*tl+1:(i-1)*tl+tl,1)), (-biom((i-1)*tl+1:(i-1)*tl+tl,2)), 'k', 'LineWidth',1.0)
    
    [x_max,idx] = max(x((i-1)*tl+1:(i-1)*tl+tl,1));
    plot(x_max,-biom((i-1)*tl+idx,2),'ro')
    
end
ylim([-50 -10])
xlabel ({'B' ; '(10^6 cells cm^-^2)'}) %label for variable and units
ylabel ({'Depth' ;'(cm)'}) 
set(gca,'XAxisLocation','top')

Answer 1

[Sulaymon Eshkabilov]

What you are trying to attain is to display the found max value and depth location to be displayed in the simulation order:

...
fprintf('Profile   max_value  depth-location \n')
for i = 1:mm    
    plot(x((i-1)*tl+1:(i-1)*tl+tl,1), -biom((i-1)*tl+1:(i-1)*tl+tl,2), 'k', 'LineWidth',1.0),  hold on 
    [x_max,idx] = max(x((i-1)*tl+1:(i-1)*tl+tl,1));
    plot(x_max,-biom((i-1)*tl+idx,2),'ro')
    fprintf('%d  %10.3f %10.0f cm  \n',[i, x_max, biom((i-1)*tl+idx,2)])  
end
...

Answer 2

[Voss]

The below includes code for making such a matrix and also a table of the same information (take whichever is more useful):

load 'biom.txt' % (Fortran file output)
tl=201; % Number of depth levels
[m,n]=size(biom); % derermining the size of the data
mm=m/tl; % total number of profiles

summary_matrix = zeros(mm,3);

summary_table = table( ...
    zeros(mm,1),zeros(mm,1),zeros(mm,1), ...
    'VariableNames',{'Profile' 'max_value' 'depth_location'});

x = 1e6*biom; % <-- here variable concentration    microMole => 1e9 miliMol=>1e6

hold on 

for i = 1:mm
    
    plot((x((i-1)*tl+1:(i-1)*tl+tl,1)),(-biom((i-1)*tl+1:(i-1)*tl+tl,2)),'k','LineWidth',1.0)
    [x_max,idx] = max(x((i-1)*tl+1:(i-1)*tl+tl,1));
    depth_of_max = biom((i-1)*tl+idx,2); % depth is positive in the summary matrix/table
    plot(x_max,-depth_of_max,'ro')       % and plotted negative
    
    summary_matrix(i,:) = [i x_max depth_of_max];
    summary_table(i,:) = {i x_max depth_of_max};
    
end

ylim([-50 -10])
xlabel ({'B' ; '(10^6 cells cm^-^2)'}) %label for variable and units
ylabel ({'Depth' ;'(cm)'}) 
set(gca,'XAxisLocation','top')

disp(summary_matrix);

   1.0e+03 *

    0.0010    0.0001         0
    0.0020    0.0001    0.0007
    0.0030    0.0007    0.1254
    0.0040    0.0023    0.1204
    0.0050    0.0042    0.1129
    0.0060    0.0060    0.1029
    0.0070    0.0076    0.0929
    0.0080    0.0092    0.0804
    0.0090    0.0110    0.0629
    0.0100    0.0159    0.0454
    0.0110    0.0319    0.0304
    0.0120    0.0824    0.0204
    0.0130    0.1975    0.0150
    0.0140    0.6344    0.0132
    0.0150    1.0471    0.0144
    0.0160    1.4309    0.0150
    0.0170    1.8585    0.0156
    0.0180    2.0602    0.0162
    0.0190    2.0714    0.0168
    0.0200    2.0328    0.0174

disp(summary_table);

    Profile    max_value    depth_location
    _______    _________    ______________

       1            0.1             0     
       2        0.10042           0.7     
       3         0.7237         125.4     
       4         2.2508         120.4     
       5         4.1584         112.9     
       6          5.988         102.9     
       7          7.619          92.9     
       8         9.1548          80.4     
       9         11.024          62.9     
      10         15.928          45.4     
      11         31.875          30.4     
      12         82.383          20.4     
      13         197.49            15     
      14         634.44          13.2     
      15         1047.1          14.4     
      16         1430.9            15     
      17         1858.5          15.6     
      18         2060.2          16.2     
      19         2071.4          16.8     
      20         2032.8          17.4