python 劃的GFS圖及GEFS圖存檔（2015年至今重要過程）

45005TY

教學部分：最好使用anaconda navigator 的jupyter notebook
在cmd 打conda install numpy，matplotlib，pandas，pint,pooch,pyproj,scipy,traitlets,xarray,
還有最重要conda install metpy
開jupyter notebook copy and paste以下的code
便會出鎮樓圖

1. from datetime import datetime, timedelta
   
2. 
   
3. import matplotlib.pyplot as plt
   
4. 
   
5. import xarray as xr
   
6. from metpy.calc import wind_components
   
7. import metpy.calc as mpcalc
   
8. from metpy.cbook import get_test_data
   
9. from metpy.io import metar
   
10. from metpy.plots.declarative import (BarbPlot, ContourPlot, FilledContourPlot, MapPanel,
    
11.                                      PanelContainer, PlotObs)
    
12. from metpy.units import units
    
13. import numpy as np
    

复制代码

1. run_date="20180907"
   
2. init="18"
   
3. forecast="174"
   
4. url = f"https://thredds.rda.ucar.edu/thredds/dodsC/files/g/ds084.1/2018/{run_date}/gfs.0p25.{run_date}{init}.f{forecast}.grib2"
   
5. data = xr.open_dataset(url)

复制代码

1. ds = data.metpy.sel(lat=slice(25,19), lon=slice(122，128))

复制代码

cfill 即是底色

1. cfill = FilledContourPlot()
   
2. cfill.data = ds
   
3. cfill.field = 'Temperature_height_above_ground'
   
4. cfill.level = 0*units.m
   
5. cfill.colormap = 'gist_ncar'
   
6. cfill.time=None
   
7. cfill.colorbar = 'horizontal'
   
8. cfill.plot_units='K'
   
9. cfill.contours = list((range(233,313,2)))

复制代码

cntr2即是plot isobar

1. cntr2 = ContourPlot()
   
2. cntr2.data = ds
   
3. cntr2.field = 'MSLP_Eta_model_reduction_msl'
   
4. cntr2.linecolor = 'black'
   
5. cntr2.linestyle = 'solid'
   
6. cntr2.time=None
   
7. cntr2.clabels = True
   
8. cntr2.plot_units='hPa'
   
9. cntr2.contours = list(range(870,1014,4))

复制代码

barbs即是barbs

1. barbs = BarbPlot()
   
2. barbs.data = ds
   
3. barbs.time = None
   
4. barbs.field = ['u-component_of_wind_isobaric', 'v-component_of_wind_isobaric']
   
5. barbs.level = 1000* units.hPa
   
6. barbs.skip = (5, 5)
   
7. barbs.plot_units = 'knots'
   
8. barbs.color='green'

复制代码

1. panel = MapPanel()
   
2. panel.area = [122,128,19,25]
   
3. panel.projection = 'mer'
   
4. panel.layers = ['states', 'coastline', 'borders']
   
5. panel.layers_edgecolor='black'
   
6. panel.title = f'MSLP isobar at 2018091500Z ,{run_date}{init}Z run, Minimum pressure:{float(np.min(cfill.data.MSLP_Eta_model_reduction_msl))/100}hPa,Maximum Temperature:{int(np.max(cfill.data.Temperature_height_above_ground))}K'
   
7. panel.plots = [cfill,cntr2,barbs]

复制代码

1. pc = PanelContainer()
   
2. pc.size = (15,15)
   
3. pc.panels=[panel]

复制代码

最後就是出圖

1. pc.show()

复制代码

--------------------------------------------------------------------
2015年至今GFS及GEFS圖存檔
以下所有圖都是用python plot 的
鎮樓圖：山竹







（底色是溫度，以K為單位）
2016年1月極端寒潮


之後是重要過程前每日00Z的GFS圖存檔

45005TY

2015年蘇迪羅
7月30日00Z GFS
巔峰913hPa 趨向台灣


7月31日 00Z GFS 巔峰915hPa 趨向福建，浙江一帶

8月1日 00Z 巔峰886hPa 趨向台灣

8月2日00Z 巔峰887hPa北掠台灣


8月3日 00Z ： 888hPa北掠台灣

8月4日 00Z：

ForzaFerrari

具体是怎么画的呀？想请教一下

9916-Bart

我们建议标出中心最低气压
光等压线带来的冲击力还不够

Samuelle

https://thredds.rda.ucar.edu/thredds/dodsC/files/g/ds084.1/2018/{run_date}/gfs.0p25.{run_date}{init}.f{forecast}.grib2
这个网站打开后显示错误

45005TY

Samuelle 发表于 2024-6-13 21:34
https://thredds.rda.ucar.edu/thredds/dodsC/files/g/ds084.1/2018/{run_date}/gfs.0p25.{run_date}{init} ...

如果想檢視dataset, 可以加上.html例如：
https://thredds.rda.ucar.edu/thredds/dodsC/files/g/ds084.1/2018/20180907/gfs.0p25.2018090700.f144.grib2.html

但code不用加上.html

Samuelle

如果预报时间在100小时内，要变成0xx的形式

Samuelle

我们画的

Samuelle

雷伊

eoj329

當年18Z的預測,天鴿大約963hpa,並預測其將於正午時分會在珠江口以西登陸。