【Matlab】HWM14：繪製一天全球水平風場 (南北、東西) 的變化，300公里高

tags: Matlab 高層大氣科學

Common 目標

• 利用 Matlab 內建的 HWM14 套件
• 繪製任一天全球水平風場(南北或東西)的變化情形
• 高度：300 公里

Introduction 簡介

• Horizontal Wind Model (HWM) 水平風場模式

安裝 Aerospace 工具套件

• HWM14 在 Aerospace 此工具套件中
• 直接在 Matlab 操作

步驟教學

1. 上方工具列 >> HOME >> Add-Ons

   

2. 搜尋 Horizontal Wind Model，因為 Aerospace 資料筆數過多

3. 選擇 Aerospace Toolbox，並下載安裝

   

HWM14

atmoshwm() 方法

• 高度：公尺
• 時間：秒

Code 程式碼

• 全球海岸線：coast.mat

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 hr=1
    for i=1:length(gim_lat)
        for j=1:length(gim_lon)            
            temp = atmoshwm(gim_lat(i),gim_lon(j),300000,'day',1,'seconds',hr*3600,'apindex',80,'model','total','version','14');
            % 1是南北向，2是東西向
            %data(i,j) = temp(1);
            data(i,j) = temp(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)')
    
    if hr<10
        hour=['0' num2str(hr)];
    else
        hour=num2str(hr);
    end
    
    title(['01/01 ' hour '00UT'])
    %----- Plot the Global coast map -----
    load coast
    hold on
    plot(long,lat,'k')
    
    %print('-dpng', ['HWM14_SN_' hour '.png'])
    print('-dpng', ['HWM14_EW_' hour '.png'])
 end

Result 結果

方向	圖
南北
東西

Reference 參考資料

1. 高層大氣科學 上課講義：
   http://class-qry.acad.ncku.edu.tw/syllabus/online_display.php?syear=0106&sem=2&co_no=C44A900