Getting incorrect plots for cosmic-2 bending profiles

3 次查看（过去 30 天）

Imtiaz nabi 2023-12-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2057939-getting-incorrect-plots-for-cosmic-2-bending-profiles

评论： Mathieu NOE 2023-12-19

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)")

8 个评论
显示 6更早的评论隐藏 6更早的评论

Mathieu NOE 2023-12-8

在 MATLAB Online 中打开

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)")