【Matlab】nrlmsise-00：全球300公里高的氧原子 [O] 的分布情形

tags: Matlab 高層大氣科學

Common 簡介

• 中性大氣經驗模式 NRLMSISE-00

  

Data Source 資料來源

• NRLMSISE-00 Atmosphere Model：連結
• Matlab Add On (附加功能、插件) - NRLMSISE-00 Atmosphere Model：連結
• 改寫後的 NRLMSISE-00 Atmosphere Model：下載

Example：[O]氧原子分布情形

• Q：試著跑出2010年1月1日下午16點，全球300公里高之氧原子[O]的分布情形

步驟1：使用 NRLMSISE-00 Atmosphere Model (改寫)

• 改寫後的 NRLMSISE-00 Atmosphere Model：下載

• 輸入的資料

  

• 輸出的資料

  

步驟2：matlab code

clear all
close all

gim_lat=87.5:-2.5:-87.5;
gim_lon=-180:5:180;

data = zeros([length(gim_lat) length(gim_lon)]);

for i=1:length(gim_lat)
    for j=1:length(gim_lon)
        % 呼叫改寫後的 function
        temp = nrlmsise00_routine(1,57600,300,gim_lat(i),gim_lon(j),16,150,150,4);
        data(i,j) = temp.d(2);
    end
end

% fig
[a,b]=contourf(gim_lon,gim_lat,data,100);
set(b,'edgecolor','none')
colorbar
colormap('jet')
axis([-180 180 -90 90])
set(gca,'xtick',-180:30:180,'ytick',-90:30:90,'fontsize',12)
xlabel('Longitude (\circE)')
ylabel('Latitude (\circN)')

title('2010/01/01 16UT 300km [O] concentration')
%----- Plot the Global coast map -----
load coast
hold on
plot(long,lat,'k')

print -dpng 'O concentration.png'

Result 結果

Reference 參考資料

1. 高層大氣科學 上課講義：
   http://class-qry.acad.ncku.edu.tw/syllabus/online_display.php?syear=0106&sem=2&co_no=C44A900
2. NRLMSISE-00 Atmosphere Model：https://ccmc.gsfc.nasa.gov/modelweb/models/nrlmsise00.php
3. NRLMSISE-00 Atmosphere Model：https://www.mathworks.com/matlabcentral/fileexchange/56253-nrlmsise-00-atmosphere-model