Getting incorrect plots for cosmic-2 bending profiles

Hey there, I am trying to reproduce results of a paper (see image below) (https://www.mdpi.com/2073-4433/14/1/121). I am trying to plot the bending profile of the files used in this paper (file attached) but for some reason my results are far off than the ones in the paper. I've searched everywhere but seems like I am missing something.

The thing that I want to get (at least one line from the following)

Please check the following code

File = "atmPrf_C2E3.2022.016.01.45.R04_0001.0001_nc";
info = ncinfo(File);
impactHeigth = ncread(File,"Impact_height");
Temp = ncread(File,'Ref');
tongLat = ncread(File, "Lat");
TongLong = ncread(File, "Lon");
pres = ncread(File,"Pres");
bendAngle = ncread(File, "Bend_ang");
bendAngle_sigma = ncread(File, 'Bend_ang_stdv');
bendAngle_ci = ncread(File, 'Bend_ang_conf');
MSL = ncread(File, "MSL_alt");
maxHeight = 40;
minHeight = 5;
alt = MSL((MSL >= minHeight) & (MSL <= maxHeight));
temp = Temp(1:length(alt));

figure 
plot(bendAngle, impactHeigth)
title("COSMIC-2 RO Bending Angle W.R.T Mean Sea Altitude")
xlabel("Bending Angle")
ylabel("Mean Sea Altitude (km)")

[Mathieu NOE]

I see that your reference plot x axis is bending angle change and not just bending angle itself

so I wonder if you should not take the diff or gradient of this angle

I get a new plot that shows better the bump around altitude = 33 km but at lower altitudes I see still a large discrepancy between your data and the reference plot - maybe it's not related to how we process it in matlab but the nc file itself ?

File = "atmPrf_C2E3.2022.016.01.45.R04_0001.0001_nc";
info = ncinfo(File)
impactHeigth = ncread(File,"Impact_height");
Temp = ncread(File,'Ref');
tongLat = ncread(File, "Lat");
TongLong = ncread(File, "Lon");
pres = ncread(File,"Pres");
bendAngle = ncread(File, "Bend_ang");
bendAngle_sigma = ncread(File, 'Bend_ang_stdv');
bendAngle_ci = ncread(File, 'Bend_ang_conf');
MSL = ncread(File, "MSL_alt");
maxHeight = 40;
minHeight = 5;
alt = MSL((MSL >= minHeight) & (MSL <= maxHeight));
temp = Temp(1:length(alt));
figure 
plot(bendAngle, impactHeigth)
ylim([10 max(impactHeigth)]);
title("COSMIC-2 RO Bending Angle W.R.T Mean Sea Altitude")
xlabel("Bending Angle")
ylabel("Mean Sea Altitude (km)")

figure 
plot(gradient(bendAngle), impactHeigth)
ylim([10 max(impactHeigth)]);
title("COSMIC-2 RO Bending Angle W.R.T Mean Sea Altitude")
xlabel("Gradient of Bending Angle")
ylabel("Mean Sea Altitude (km)")

[Imtiaz nabi]

perhaps it was the gradient!!!! Thank you please take a look at the following figure now. One more thing, Can we smoothen the graph by filtering it somehow? umm like by applying a butterworth filter if possible.

[Mathieu NOE]

try smoothdata

[Imtiaz nabi]

Hey, I tried that but it didn't work!

[Mathieu NOE]

did you get an error message or is the result not what you expected ?

[Imtiaz nabi]

The result was simply not what I expected

[Imtiaz nabi]

One more question, how would one calculate the % change in bending angle?

[Mathieu NOE]

I would say :

percentage = 100 * (angle - mean(angle)) / mean(angle)