earth套件讀取CWB_WRF數據

背景

動機與目的

• 中央氣象局雖然每天進行WRF數值預報，網路上卻沒有開發者使用來做氣流線等加值應用。雖然過去發展了高解析度軌跡近5日預報(see also [[2022-11-04-daily_traj]])、calpuff每日的預報等等，卻還沒有嘗試以earth.nullschool套件來加以延伸、展示。
• 除了每日、即時的展示之外，還有儲存、檢視的內部需求。相較windy，earth.nullschool可以按照使用者的需要，在url地址列簡單的指令，即可顯示過去特定時間的氣象與(或)濃度場。這個方案目前似乎是網路存取、檢視龐雜的氣象-空品模擬系統成果的最佳平台。

earth的發展與應用

• earth套件是Cameron Beccario(cambecc)早期為東京都環境局環境改善部撰寫的套件，後來發展成全球的服務網站，同時也被大陸地方選擇做為氣象預報可視化系統的設計項目。
• ESRI 2017年也將其繼續發展成類似油畫質感的動畫版本。
• 持續的商業版本除了windy之外，ventusky也有類似的應用。
• WebGL控制方案與3D繪圖效果(webgl-wind by XXHolic, 2022)。

整體工作流程

graph TD
    A(CWB_WRF3Km U10V10_d03) --> B((uv10_json.py))
    C(json template) --> B
    B --> D(current-wind-surface-level-cwb-3K.json)
    D --> E((products.js))
    E --> F((web rendering))

• uv10_json.py
  • source code
  • IO and program descriptions
• earth套件讀取gfs檔案的核心程式在products.js，必須讓它能夠讀取cwb轉成的json檔。
  • 找到有關gfs的程式碼：grep gfs $(find . -name "*.js")|more
• web rendering
  • intializing by node dev-server.js 80 or
  • reload browser
  • d3：http://125.229.149.182/#current/wind/surface/level/orthographic=-237.53,23.30,2000
  • d1：http://125.229.149.182/#current/wind/surface/level/orthographic=-236.33,23.30,800

diff of first paramter between gfs and cwbwrf_15Km files

• earth套件與gfs檔案的連結靠的是grib2json這支程式，cambecc也將其公開在github上。其下載、編譯、與應用的細節歷程可以參考FAQ->json。
• 經由下列簡單的python指令，可以列出gfs與cwb二者grib檔案內容的差異，作為修改(uv10_json.py)的重要指引。

fname='current-wind-surface-level-gfs-1.0.json'
with open(fname,'r') as fn:
  crnt=json.load(fn)
fname='M-A0061-006.grb2.json'
with open(fname,'r') as fn:
  cwb=json.load(fn)
c=set(cwb[0]['header'])
g=set(crnt[0]['header'])
a=[]
for i in c&g:a.append([i, crnt[0]['header'][i],cwb[0]['header'][i]])

grib file

grid definition

• 所有的實數必須是64位元(np.float64(…))，否則在json.dump時會報錯。
• 緯度的方向是由高到低。雖然dy是給正值，但程式內部gfs[:][‘data’]的排序在y方向是由北到南，需對其軸進行np.flip(…)。可以由風速較低值的範圍與海陸界限加以檢視，形狀會南北相反。

gfs[ir]['data']=list(np.flip(np.where(zz!=zz,0,zz),axis=0).flatten())

surface

set(surface1Value)
{0.0, 
1e-128, 1.5e-128, 2e-128, 2.5e-128, 3.e-128, 4e-128, 
5e-128, 7.e-128, 8.5e-128, 9.25e-128, 1e-127, 
2.0, 10.0}

time frame

installation of earth

download and compile

• git clone https://github.com/cambecc/earth.git
• cd earth
• npm install

initialize the website

• node dev-server.js 8080 (or 80 for usual port)

modifications of js

zoom limits

• from 3000 to 60,000
• function scaleExtent in file ./public/libs/earth/1.0.0/globes.js

FilePath function

• 原本(gfs1p0degPath)指定檔名的函數，寫死一定要開current-wind-surface-level-gfs-1.0.json
  • 引數只提供attr, type, surface, level
• 增加src, res 2個引數，來源（gfs或cwb）與解析度('1p90', '15K', '3K')
• in file ./public/libs/earth/1.0.0/products.js

default settings

• in file: ./public/libs/earth/1.0.0/micro.js:14:
• 可設定：
  1. 時間
  2. 貼上特定屬性的shade等值圖
  3. 高度(地面或定壓層)
  4. 座標系統投影方式
  5. 初始中心點經緯度與zoom倍率

var DEFAULT_CONFIG = "current/wind/surface/level/orthographic";

• 修正為

var DEFAULT_CONFIG = "current/wind/surface/level/orthographic=-238.80,23.73,7500";

縣市界

要在earth套件中更改海岸線、湖泊線，或是加上縣市界，可以參考詳見[shp檔轉json][shp_json]說明

• topojson檔案放在./public/data/topo目錄下
• 有關的js檔包括
  1. ./public/libs/earth/1.0.0/earth.js (開啟topojson檔)
  2. ./public/libs/earth/1.0.0/globes.js (變數定義)

氣流線之調整

earth套件的動態氣流線在其程式中稱為particle，可能因為看起來像是顆粒的軌跡動態。particle相關設定有：

1. ./public/libs/earth/1.0.0/[products.js][js1]： particles: {velocityScale: 1/3000000, maxIntensity: 1 }，這2個設定前者為風速的尺度(1/3M ~ 1/60K)，其值越小，流線越短。後者為密集度(1 ~ 17)，其值越小會有越多處有particle起點，即使風速很低，適用在小比例尺圖面。
2. ./public/libs/earth/1.0.0/earth.js

    var MAX_PARTICLE_AGE = 100;               // max number of frames a particle is drawn before regeneration
    var PARTICLE_LINE_WIDTH = 1.0;            // line width of a drawn particle
    var PARTICLE_MULTIPLIER = 7;              // particle count scalar (completely arbitrary--this values looks nice)
    var PARTICLE_REDUCTION = 0.75;            // reduce particle count to this much of normal for mobile devices

下載與執行

• CWB WRF程式結果每6小時更新，分別為每天的2/8/14/20時等4次。每次預報84-6=78小時(0~6小時warm up)
  • 由於CWB是陸續更新，而其檔名系統又以預報起始時間為0，因此有可能發生錯亂覆蓋的情形。
  • 解決方式：解析結果按照檔案內的時間(grbs[1].validDate)來命名
  • 並將結果按照年月日另建目錄擺放(參earth/public/test/products-test.html內容)
  • eg: "/data/weather/2013/11/20/0800-a-b-c-gfs-1.0.json"
  • 下載、解讀之自動化腳本(earth_cwbwrf.cs)乃在node03執行，以節省頻寬。
• 網頁程式安放在iMacKuang，以避免防火牆干擾。
  • current目錄下之檔案：每小時執行連結。

node03作業

腳本

• grbuv10_json.py的概略說明可以參考地面風wrfout檔轉json->grib2 to json directly。
  • node03版本加入了後處理項目：
    1. 在mac上創新目錄
    2. 將結果傳遞到mac的指定目錄

#kuang@node03 /nas1/Data/javascripts/D3js/earth/public/data/weather/current
#$ cat earth_cwbwrf.cs
weather=/nas1/Data/javascripts/D3js/earth/public/data/weather
cd $weather/current
for dom in 4;do
for ((i=6; i<=84; i=i+6)); do
  ii=`printf "%02d" $i`
  fn=M-A006${dom}-0$ii.grb2
  rm -f $fn
  /usr/bin/wget https://cwbopendata.s3.ap-northeast-1.amazonaws.com/MIC/$fn
  ./grbuv10_json.py $fn
done
done

crontab內容

crontab -l|grep earth
# earth CWB_WRF
55 2,8,14,20 * * * /nas1/Data/javascripts/D3js/earth/public/data/weather/current/earth_cwbwrf.cs

mac上的任務

lnk_curr.cs腳本

• mac的任務就是保持current檔案隨時都是當下的時間
• 執行方式是每小時進行UTC時間的計算，並將對的檔案連結到current
  • UTC = LST - 8H
  • 因為不是每小時都有模擬結果。需要進行計算，將LST逐時對照到UTC逐6小時。
  • 文字的數字改成10進位計算：#10$h，計算完後再轉成2碼數字

#kuang@114-32-164-198 /Users/Data/javascripts/D3js/earth/public/data/weather/current
#$ cat lnk_curr.cs
now=$(date -v-8H -j  +%Y%m%d%H)
y=$(date -j -f "%Y%m%d%H" "${now}" +%Y)
m=$(date -j -f "%Y%m%d%H" "${now}" +%m)
d=$(date -j -f "%Y%m%d%H" "${now}" +%d)
h=$(date -j -f "%Y%m%d%H" "${now}" +%H)
h=$(( 10#$h / 6 * 6|bc -l ))
h=$(printf "%02d" $h)
weather=/Users/Data/javascripts/D3js/earth/public/data/weather
fn=$weather/$y/$m/$d/${h}00-wind-surface-level-cwb-3K.json
if [ -e $fn ]; then
  ln -sf $fn $weather/current/current-wind-surface-level-cwb-3K.json
else
  echo $fn
fi

crontab

0 * * * * /Users/Data/javascripts/D3js/earth/public/data/weather/current/lnk_curr.cs

成果檢討

3公里解析度範圍與流線場

2021070206中央氣象局WRF預報之地面風氣流線
同時間的GFS地面風

• CWB預報有颱風、GFS並沒有發現
• 很多地型效應、台灣東側背風渦漩、小島背風面較低風速等等現象，在CWB風場中都有展現。
• 內陸風場細節：GFS無法模擬

TODO

• other timeframe, fields
• webGL controler
• fine resolution shape file
  • administrative boundies
  • isopleth for terrain elevation(geotiff to KML, KML to shape file)

resource

earth package and D3js implements

• cambecc(2016), earth building, launching and etc on GitHub.
• cambecc(2017), grib2json on GitHub.
• Roger Veciana i Rovira(2018), Drawing wind barbs with D3js from a GeoTIFF, on bl.ocks.org
• 陈晖, 范玉鑫, 陈杨,及 吴天亭(2016), 吉林省WRF模式数值预报可视化系统设计, 现代农业科技2016年4期.
• Kuan-Jung, Huang(2019), 透過 D3.js 調用外部資料集, mybaseball52.medium.com

about GFS downloading

• curl command in Github
  • curl "http://nomads.ncep.noaa.gov/cgi-bin/filter_gfs.pl?file=gfs.t00z.pgrb2.1p00.f000&lev_10_m_above_ground=on&var_UGRD=on&var_VGRD=on&dir=%2Fgfs.${YYYYMMDD}00" -o gfs.t00z.pgrb2.1p00.f000
• nomads website has been upgraded since cambecc released earth.
• curl command become:
  • curl "https://nomads.ncep.noaa.gov/cgi-bin/filter_gfs_1p00.pl?file=gfs.t00z.pgrb2.1p00.f000&lev_10_m_above_ground=on&var_UGRD=on&var_VGRD=on&leftlon=0&rightlon=360&toplat=90&bottomlat=-90&dir=%2Fgfs.${YYYYMMDD}%2F00%2Fatmos" -o gfs.t00z.pgrb2.1p00.f000
• GFS PARAMETERS & UNITS
• GRIB Filters and View the URL

about HTTPS

• 詮力科技(2019),為您的網站加上「-ssl憑證-」，成為https網頁,十二月 13,2019/技術探討
• Mangle Kuo(2021),設定macOS本地端HTTPs/SSL證書
• Bharath(2021) How to install the Securly SSL certificate on Mac OSX ?,support.securly.com

about json TypeError

• 咸魚(2020), TypeError: Object of type ‘float32’ is not JSON serializable解决方案, blog.csdn.net