剖面诊断分析绘图diag_crossection
<p>[TOC]</p>
<h1>相当位温, 沿剖面风和风速</h1>
<p>相当位温, 沿剖面风和风速产品定制代码示例如下:</p>
<pre><code class="language-python">draw_wind_theta_wsp(cross_wsp, cross_theta, cross_u, cross_v, cross_terrain, hgt,
st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65),
h_pos=[0.125, 0.665, 0.25, 0.2],
wsp_contourf_kwargs={}, theta_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={},
**pallete_kwargs)</code></pre>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">类型</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">cross_wsp</td>
<td style="text-align: left;"><font color="red">是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">风速</td>
</tr>
<tr>
<td style="text-align: left;">cross_theta</td>
<td style="text-align: left;"><font color="red">是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">相当位温</td>
</tr>
<tr>
<td style="text-align: left;">cross_u</td>
<td style="text-align: left;"><font color="red">是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">U风</td>
</tr>
<tr>
<td style="text-align: left;">cross_v</td>
<td style="text-align: left;"><font color="red">是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">cross_terrain</td>
<td style="text-align: left;"><font color="red">是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">hgt</td>
<td style="text-align: left;"><font color="red">是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">高度数据</td>
</tr>
<tr>
<td style="text-align: left;">st_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">起始点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">ed_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">结束点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">lon_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">经度剖面</td>
</tr>
<tr>
<td style="text-align: left;">lat_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">纬度剖面</td>
</tr>
<tr>
<td style="text-align: left;">map_extent</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">(50, 150, 0, 65)</td>
<td style="text-align: left;">绘图经纬度</td>
</tr>
<tr>
<td style="text-align: left;">h_pos</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td>
<td style="text-align: left;">小图区域</td>
</tr>
<tr>
<td style="text-align: left;">wsp_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">theta_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中barbs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">terrain_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv
from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d
from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d
from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D
from metdig.onestep.complexgrid_var.theta import read_theta3d
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
rh = get_model_3D_grid(...)
u = get_model_3D_grid(...)
v = get_model_3D_grid(...)
tmp = get_model_3D_grid(...)
hgt = get_model_3D_grid(...)
psfc = get_model_3D_grid(...)
wsp = mdgcal.other.wind_speed(u, v)
cross_rh = mdgcal.cross_section(rh, st_point, ed_point)
cross_u = mdgcal.cross_section(u, st_point, ed_point)
cross_v = mdgcal.cross_section(v, st_point, ed_point)
cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point)
cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point)
cross_wsp = mdgcal.cross_section(wsp, st_point, ed_point)
cross_td = mdgcal.dewpoint_from_relative_humidity(cross_tmp, cross_rh)
pressure = mdgstda.gridstda_full_like_by_levels(cross_rh, cross_tmp['level'].values)
cross_theta = mdgcal.equivalent_potential_temperature(pressure, cross_tmp, cross_td)
cross_terrain = pressure - cross_psfc
cross_terrain.attrs['var_units'] = ''
# plot
draw_cross.draw_time_rh_uv_theta(rh, u, v, theta, ...)</code></pre>
<h1>相当位温, 锋生函数, 沿剖面风</h1>
<p>相当位温, 锋生函数, 沿剖面风产品定制代码示例如下:</p>
<pre><code class="language-python">draw_wind_theta_fg( cross_u, cross_v, cross_theta, cross_fg, cross_terrain, hgt,
st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65),
h_pos=[0.125, 0.665, 0.25, 0.2],
fg_contourf_kwargs={}, theta_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={},
**pallete_kwargs)</code></pre>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">类型</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">cross_u</td>
<td style="text-align: left;"><font color="red">是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">U风</td>
</tr>
<tr>
<td style="text-align: left;">cross_v</td>
<td style="text-align: left;"><font color="red">是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">cross_theta</td>
<td style="text-align: left;"><font color="red">是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">相当位温</td>
</tr>
<tr>
<td style="text-align: left;">cross_fg</td>
<td style="text-align: left;"><font color="red">是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
</tr>
<tr>
<td style="text-align: left;">cross_terrain</td>
<td style="text-align: left;"><font color="red">是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">hgt</td>
<td style="text-align: left;"><font color="red">是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">高度数据</td>
</tr>
<tr>
<td style="text-align: left;">st_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">起始点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">ed_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">结束点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">lon_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">经度剖面</td>
</tr>
<tr>
<td style="text-align: left;">lat_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">纬度剖面</td>
</tr>
<tr>
<td style="text-align: left;">map_extent</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">(50, 150, 0, 65)</td>
<td style="text-align: left;">绘图经纬度</td>
</tr>
<tr>
<td style="text-align: left;">h_pos</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td>
<td style="text-align: left;">小图区域</td>
</tr>
<tr>
<td style="text-align: left;">fg_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">theta_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中barbs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">terrain_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python"></code></pre>
<h1>相当位温、垂直运动速度、沿剖面风</h1>
<p>相当位温、垂直运动速度、沿剖面风产品定制代码示例如下:</p>
<pre><code class="language-python">draw_wind_theta_w(cross_u, cross_v, cross_theta, cross_w, cross_terrain, hgt,
st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65),
h_pos=[0.125, 0.665, 0.25, 0.2],
w_contourf_kwargs={}, theta_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">cross_w</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">垂直运动速度</td>
</tr>
<tr>
<td style="text-align: left;">cross_theta</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">相当位温</td>
</tr>
<tr>
<td style="text-align: left;">cross_u</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">U风</td>
</tr>
<tr>
<td style="text-align: left;">cross_v</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">cross_terrain</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">hgt</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">高度数据</td>
</tr>
<tr>
<td style="text-align: left;">st_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">起始点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">ed_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">结束点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">lon_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">经度剖面</td>
</tr>
<tr>
<td style="text-align: left;">lat_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">纬度剖面</td>
</tr>
<tr>
<td style="text-align: left;">h_pos</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td>
<td style="text-align: left;">小图区域</td>
</tr>
<tr>
<td style="text-align: left;">map_extent</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">(50, 150, 0, 65)</td>
<td style="text-align: left;">绘图经纬度</td>
</tr>
<tr>
<td style="text-align: left;">w_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">theta_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中barbs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">terrain_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
rh = get_model_3D_grid(...)
u = get_model_3D_grid(...)
v = get_model_3D_grid(...)
tmp = get_model_3D_grid(...)
hgt = get_model_grid(...)
psfc = get_model_grid(...)
w=read_w3d(...)
cross_rh = mdgcal.cross_section(rh, st_point, ed_point)
cross_u = mdgcal.cross_section(u, st_point, ed_point)
cross_v = mdgcal.cross_section(v, st_point, ed_point)
cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point)
cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point)
cross_w = mdgcal.cross_section(w, st_point, ed_point)
cross_td = mdgcal.dewpoint_from_relative_humidity(cross_tmp, cross_rh)
pressure = mdgstda.gridstda_full_like_by_levels(cross_rh, cross_tmp['level'].values)
cross_theta = mdgcal.equivalent_potential_temperature(pressure, cross_tmp, cross_td)
cross_terrain = pressure - cross_psfc
cross_terrain.attrs['var_units'] = ''
# plot
draw_cross.draw_wind_theta_w(cross_u, cross_v, cross_theta, cross_w, cross_terrain, hgt, ...)</code></pre>
<h1>相当位温、水平散度、沿剖面风</h1>
<p>相当位温、水平散度、沿剖面风产品定制代码示例如下:</p>
<pre><code class="language-python">draw_wind_theta_div(cross_div, cross_theta, cross_u, cross_v, cross_terrain, hgt,
st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65),
h_pos=[0.125, 0.665, 0.25, 0.2],
div_contourf_kwargs={}, theta_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">cross_div</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">水平散度</td>
</tr>
<tr>
<td style="text-align: left;">cross_theta</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">相当位温</td>
</tr>
<tr>
<td style="text-align: left;">cross_u</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">U风</td>
</tr>
<tr>
<td style="text-align: left;">cross_v</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">cross_terrain</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">hgt</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">高度数据</td>
</tr>
<tr>
<td style="text-align: left;">st_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">起始点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">ed_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">结束点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">lon_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">经度剖面</td>
</tr>
<tr>
<td style="text-align: left;">lat_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">纬度剖面</td>
</tr>
<tr>
<td style="text-align: left;">h_pos</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td>
<td style="text-align: left;">小图区域</td>
</tr>
<tr>
<td style="text-align: left;">map_extent</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">(50, 150, 0, 65)</td>
<td style="text-align: left;">绘图经纬度</td>
</tr>
<tr>
<td style="text-align: left;">div_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">theta_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中barbs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">terrain_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv
from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d
from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d
from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D
from metdig.onestep.complexgrid_var.theta import read_theta3d
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
rh = get_model_3D_grid(...)
div,u,v = read_div_uv_3d(...)
tmp = get_model_3D_grid(...)
hgt = get_model_grid(...)
psfc = get_model_grid(...)
# +form 3D psfc
_, psfc_bdcst = xr.broadcast(tmp, psfc.squeeze())
psfc_bdcst = psfc_bdcst.where(psfc_bdcst &gt; -10000, drop=True) # 去除小于-10000
cross_rh = mdgcal.cross_section(rh, st_point, ed_point)
cross_u = mdgcal.cross_section(u, st_point, ed_point)
cross_v = mdgcal.cross_section(v, st_point, ed_point)
cross_div = mdgcal.cross_section(div, st_point, ed_point)
cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point)
cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point)
cross_td = mdgcal.dewpoint_from_relative_humidity(cross_tmp, cross_rh)
pressure = mdgstda.gridstda_full_like_by_levels(cross_rh, cross_tmp['level'].values)
cross_theta = mdgcal.equivalent_potential_temperature(pressure, cross_tmp, cross_td)
cross_terrain = pressure - cross_psfc
# plot
draw_cross.draw_wind_theta_div(cross_div, cross_theta, cross_u, cross_v, cross_terrain, hgt, ...)</code></pre>
<h1>云水比、雪水比、水平风、温度时间剖面</h1>
<p>云水比、雪水比、水平风、温度时间剖面产品定制代码示例如下:</p>
<pre><code class="language-python">draw_time_wind_qcld_qsn_tmp(qcld, qsn, tmp, u, v, terrain, mean_area=None,
qcld_contour_kwargs={},qice_contour_kwargs={}, tmp_contourf_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">qcld</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">云水混合比</td>
</tr>
<tr>
<td style="text-align: left;">qsn</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">雨水混合比</td>
</tr>
<tr>
<td style="text-align: left;">tmp</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">温度</td>
</tr>
<tr>
<td style="text-align: left;">u</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">v</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">mean_area</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;"></td>
</tr>
<tr>
<td style="text-align: left;">qcld_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">qice_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">tmp_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中brabs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">terrain_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv
from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d
from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d
from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D
from metdig.onestep.complexgrid_var.theta import read_theta3d
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
u = get_model_3D_grids(...)
v = get_model_3D_grids(...)
qcld = get_model_3D_grids(...)
qsn = get_model_3D_grids(...)
psfc = get_model_3D_grids(...)
tmp = get_model_3D_grids(...)
_, pressure = xr.broadcast(v, v['level'])
terrain = pressure - psfc.values
# plot
draw_cross.draw_time_wind_qcld_qsn_tmp(qcld, qsn, tmp, u, v, terrain, ...)</code></pre>
<h1>云水比、冰水比、水平风、温度时间剖面</h1>
<p>云水比、冰水比、水平风、温度时间剖面产品定制代码示例如下:</p>
<pre><code class="language-python">draw_time_wind_qcld_qice_tmp(qcld, qice, tmp, u, v, terrain,mean_area=None,
qcld_contour_kwargs={},qice_contour_kwargs={}, tmp_contourf_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">qcld</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">云水混合比</td>
</tr>
<tr>
<td style="text-align: left;">qice</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">云冰混合比</td>
</tr>
<tr>
<td style="text-align: left;">tmp</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">温度</td>
</tr>
<tr>
<td style="text-align: left;">u</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">v</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">mean_area</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;"></td>
</tr>
<tr>
<td style="text-align: left;">qcld_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">qice_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">tmp_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中brabs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">terrain_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv
from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d
from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d
from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D
from metdig.onestep.complexgrid_var.theta import read_theta3d
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
u = get_model_3D_grids(...)
v = get_model_3D_grids(...)
qcld = get_model_3D_grids(...)
qice = get_model_3D_grids(...)
psfc = get_model_3D_grids(...)
tmp = get_model_3D_grids(...)
_, pressure = xr.broadcast(v, v['level'])
terrain = pressure - psfc.values
# plot
draw_cross.draw_time_wind_qcld_qice_tmp(qcld, qice, tmp, u, v, terrain, ...)</code></pre>
<h1>温度、温度平流、沿剖垂直环流</h1>
<p>温度、温度平流、沿剖垂直环流产品定制代码示例如下:</p>
<pre><code class="language-python">draw_wind_w_tmpadv_tmp(cross_tmpadv, cross_tmp, cross_t, cross_w, cross_terrain, hgt,
st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65),
h_pos=[0.125, 0.665, 0.25, 0.2],
tmpadv_contourf_kwargs={}, tmp_contour_kwargs={}, wind_quiver_kwargs={},terrain_contourf_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">cross_tmpadv</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">温度平流</td>
</tr>
<tr>
<td style="text-align: left;">cross_tmp</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">温度</td>
</tr>
<tr>
<td style="text-align: left;">cross_t</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">温度</td>
</tr>
<tr>
<td style="text-align: left;">cross_w</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">风速</td>
</tr>
<tr>
<td style="text-align: left;">cross_terrain</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">hgt</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">高度</td>
</tr>
<tr>
<td style="text-align: left;">st_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">起始点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">ed_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">结束点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">lon_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">经度剖面</td>
</tr>
<tr>
<td style="text-align: left;">lat_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">纬度剖面</td>
</tr>
<tr>
<td style="text-align: left;">h_pos</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td>
<td style="text-align: left;">小图区域</td>
</tr>
<tr>
<td style="text-align: left;">map_extent</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">(50, 150, 0, 65)</td>
<td style="text-align: left;">绘图经纬度</td>
</tr>
<tr>
<td style="text-align: left;">tmpadv_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">tmp_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">wind_quiver_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中quiver函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">terrain_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv
from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d
from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d
from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D
from metdig.onestep.complexgrid_var.theta import read_theta3d
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
u = get_model_3D_grid(...)
v = get_model_3D_grid(...)
w=read_w3d(...)
tmp = get_model_3D_grid(...)
hgt = get_model_grid(...)
psfc = get_model_grid(...)
tmpadv = mdgcal.var_advect(tmp, u, v)
cross_u = mdgcal.cross_section(u, st_point, ed_point)
cross_v = mdgcal.cross_section(v, st_point, ed_point)
cross_w = mdgcal.cross_section(w, st_point, ed_point)
cross_t, cross_n = mdgcal.cross_section_components(cross_u, cross_v)
cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point)
cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point)
cross_tmpadv = mdgcal.cross_section(tmpadv, st_point, ed_point)
pressure = mdgstda.gridstda_full_like_by_levels(cross_tmp, cross_tmp['level'].values)
cross_terrain = pressure - cross_psfc
cross_terrain.attrs['var_units'] = ''
# plot
draw_cross.draw_wind_w_tmpadv_tmp(cross_tmpadv, cross_tmp, cross_t, cross_w*ratio, cross_terrain, hgt, ...)</code></pre>
<h1>温度、温度平流、沿剖面风</h1>
<p>温度、温度平流、沿剖面风产品定制代码示例如下:</p>
<pre><code class="language-python">draw_wind_tmpadv_tmp(cross_tmpadv, cross_tmp, cross_u, cross_v, cross_terrain, hgt,
st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65),
h_pos=[0.125, 0.665, 0.25, 0.2],
tmpadv_contourf_kwargs={}, tmp_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">cross_tmpadv</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">温度平流</td>
</tr>
<tr>
<td style="text-align: left;">cross_tmp</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">温度</td>
</tr>
<tr>
<td style="text-align: left;">cross_u</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">U风</td>
</tr>
<tr>
<td style="text-align: left;">cross_v</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">cross_terrain</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">hgt</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">高度</td>
</tr>
<tr>
<td style="text-align: left;">st_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">起始点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">ed_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">结束点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">lon_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">经度剖面</td>
</tr>
<tr>
<td style="text-align: left;">lat_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">纬度剖面</td>
</tr>
<tr>
<td style="text-align: left;">h_pos</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td>
<td style="text-align: left;">小图区域</td>
</tr>
<tr>
<td style="text-align: left;">map_extent</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">(50, 150, 0, 65)</td>
<td style="text-align: left;">绘图经纬度</td>
</tr>
<tr>
<td style="text-align: left;">tmpadv_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">tmp_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中barbs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">terrain_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv
from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d
from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d
from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D
from metdig.onestep.complexgrid_var.theta import read_theta3d
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
u = get_model_3D_grid(...)
v = get_model_3D_grid(...)
tmp = get_model_3D_grid(...)
hgt = get_model_grid(...)
psfc = get_model_grid(...)
tmpadv = mdgcal.var_advect(tmp, u, v)
cross_u = mdgcal.cross_section(u, st_point, ed_point)
cross_v = mdgcal.cross_section(v, st_point, ed_point)
cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point)
cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point)
cross_tmpadv = mdgcal.cross_section(tmpadv, st_point, ed_point)
pressure = mdgstda.gridstda_full_like_by_levels(cross_tmp, cross_tmp['level'].values)
cross_terrain = pressure - cross_psfc
cross_terrain.attrs['var_units'] = ''
# plot
draw_cross.draw_wind_tmpadv_tmp(cross_tmpadv, cross_tmp, cross_u, cross_v, cross_terrain, hgt, ...)</code></pre>
<h1>温度、垂直涡度平流、沿剖面风</h1>
<p>温度、垂直涡度平流、沿剖面风产品定制代码示例如下:</p>
<pre><code class="language-python">draw_wind_vortadv_tmp(cross_vortadv, cross_tmp, cross_u, cross_v, cross_terrain, hgt,
st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65),
h_pos=[0.125, 0.665, 0.25, 0.2],
vortadv_contourf_kwargs={}, tmp_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">cross_vortadv</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">垂直涡度平流</td>
</tr>
<tr>
<td style="text-align: left;">cross_tmp</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">温度</td>
</tr>
<tr>
<td style="text-align: left;">cross_u</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">U风</td>
</tr>
<tr>
<td style="text-align: left;">cross_v</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">cross_terrain</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">hgt</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">高度</td>
</tr>
<tr>
<td style="text-align: left;">st_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">起始点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">ed_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">结束点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">lon_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">经度剖面</td>
</tr>
<tr>
<td style="text-align: left;">lat_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">纬度剖面</td>
</tr>
<tr>
<td style="text-align: left;">h_pos</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td>
<td style="text-align: left;">小图区域</td>
</tr>
<tr>
<td style="text-align: left;">map_extent</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">(50, 150, 0, 65)</td>
<td style="text-align: left;">绘图经纬度</td>
</tr>
<tr>
<td style="text-align: left;">vortadv_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">tmp_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中barbs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">terrain_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv
from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d
from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d
from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D
from metdig.onestep.complexgrid_var.theta import read_theta3d
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
u = get_model_3D_grid(...)
v = get_model_3D_grid(...)
tmp = get_model_3D_grid(...)
hgt = get_model_grid(...)
psfc = get_model_grid(...)
vort = mdgcal.vorticity(u, v)
vortadv = mdgcal.var_advect(vort, u, v)
cross_u = mdgcal.cross_section(u, st_point, ed_point)
cross_v = mdgcal.cross_section(v, st_point, ed_point)
cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point)
cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point)
cross_vortadv = mdgcal.cross_section(vortadv, st_point, ed_point)
pressure = mdgstda.gridstda_full_like_by_levels(cross_tmp, cross_tmp['level'].values)
cross_terrain = pressure - cross_psfc
cross_terrain.attrs['var_units'] = ''
# plot
draw_cross.draw_wind_vortadv_tmp(cross_vortadv, cross_tmp, cross_u, cross_v, cross_terrain, hgt, ...)</code></pre>
<h1>饱和相当位温, 准地转湿位涡, 沿剖面准地转风</h1>
<p>饱和相当位温, 准地转湿位涡, 沿剖面准地转风产品定制代码示例如下:</p>
<pre><code class="language-python">draw_wind_thetaes_mpvg(cross_mpvg, cross_theta, cross_u, cross_v, cross_terrain, hgt,
st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65),
h_pos=[0.125, 0.665, 0.25, 0.2],
mpv_contourf_kwargs={}, theta_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">cross_mpvg</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
</tr>
<tr>
<td style="text-align: left;">cross_theta</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">相当位温</td>
</tr>
<tr>
<td style="text-align: left;">cross_u</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">U风</td>
</tr>
<tr>
<td style="text-align: left;">cross_v</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">cross_terrain</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">hgt</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">高度</td>
</tr>
<tr>
<td style="text-align: left;">st_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">起始点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">ed_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">结束点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">lon_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">经度剖面</td>
</tr>
<tr>
<td style="text-align: left;">lat_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">纬度剖面</td>
</tr>
<tr>
<td style="text-align: left;">h_pos</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td>
<td style="text-align: left;">小图区域</td>
</tr>
<tr>
<td style="text-align: left;">map_extent</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">(50, 150, 0, 65)</td>
<td style="text-align: left;">绘图经纬度</td>
</tr>
<tr>
<td style="text-align: left;">mpv_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">theta_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中barbs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">terrain_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python"></code></pre>
<h1>相当位温、湿位涡、沿剖面风</h1>
<p>相当位温、湿位涡、沿剖面风产品定制代码示例如下:</p>
<pre><code class="language-python">draw_wind_theta_mpv(cross_mpv, cross_theta, cross_u, cross_v, cross_terrain, hgt,
st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65),
h_pos=[0.125, 0.665, 0.25, 0.2],
mpv_contourf_kwargs={}, theta_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">cross_mpv</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">湿位涡</td>
</tr>
<tr>
<td style="text-align: left;">cross_theta</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">相当位温</td>
</tr>
<tr>
<td style="text-align: left;">cross_u</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">U风</td>
</tr>
<tr>
<td style="text-align: left;">cross_v</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">cross_terrain</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">hgt</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">高度</td>
</tr>
<tr>
<td style="text-align: left;">st_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">起始点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">ed_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">结束点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">lon_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">经度剖面</td>
</tr>
<tr>
<td style="text-align: left;">lat_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">纬度剖面</td>
</tr>
<tr>
<td style="text-align: left;">h_pos</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td>
<td style="text-align: left;">小图区域</td>
</tr>
<tr>
<td style="text-align: left;">map_extent</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">(50, 150, 0, 65)</td>
<td style="text-align: left;">绘图经纬度</td>
</tr>
<tr>
<td style="text-align: left;">mpv_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">theta_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中barbs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">terrain_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv
from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d
from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d
from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D
from metdig.onestep.complexgrid_var.theta import read_theta3d
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
theta = read_theta3d(...)
mpv, _div, u, v = read_pv_div_uv(...)
hgt = get_model_grid(...)
psfc = get_model_grid(...)
cross_theta = mdgcal.cross_section(theta, st_point, ed_point)
cross_u = mdgcal.cross_section(u, st_point, ed_point)
cross_v = mdgcal.cross_section(v, st_point, ed_point)
cross_mpv = mdgcal.cross_section(mpv, st_point, ed_point)
cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point)
pressure = mdgstda.gridstda_full_like_by_levels(cross_theta, cross_theta['level'].values)
cross_terrain = pressure - cross_psfc
cross_terrain.attrs['var_units'] = ''
# plot
draw_cross.draw_wind_theta_mpv(cross_mpv, cross_theta, cross_u, cross_v, cross_terrain, hgt, ...)</code></pre>
<h1>相当位温、绝对涡度、沿剖面风</h1>
<p>相当位温、绝对涡度、沿剖面风产品定制代码示例如下:</p>
<pre><code class="language-python">draw_wind_theta_absv(cross_absv, cross_theta, cross_u, cross_v, cross_terrain, hgt,
st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65),
h_pos=[0.125, 0.665, 0.25, 0.2],
absv_contourf_kwargs={}, theta_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">cross_absv</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">绝对涡度</td>
</tr>
<tr>
<td style="text-align: left;">cross_theta</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">相当位温</td>
</tr>
<tr>
<td style="text-align: left;">cross_u</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">U风</td>
</tr>
<tr>
<td style="text-align: left;">cross_v</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">cross_terrain</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">hgt</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">高度</td>
</tr>
<tr>
<td style="text-align: left;">st_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">起始点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">ed_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">结束点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">lon_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">经度剖面</td>
</tr>
<tr>
<td style="text-align: left;">lat_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">纬度剖面</td>
</tr>
<tr>
<td style="text-align: left;">h_pos</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td>
<td style="text-align: left;">小图区域</td>
</tr>
<tr>
<td style="text-align: left;">map_extent</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">(50, 150, 0, 65)</td>
<td style="text-align: left;">绘图经纬度</td>
</tr>
<tr>
<td style="text-align: left;">absv_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">theta_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中barbs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">terrain_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv
from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d
from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d
from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D
from metdig.onestep.complexgrid_var.theta import read_theta3d
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
rh = get_model_3D_grid(...)
u = get_model_3D_grid(...)
v = get_model_3D_grid(...)
tmp = get_model_3D_grid(...)
hgt = get_model_grid(...)
psfc = get_model_grid(...)
absv = mdgcal.absolute_vorticity(u, v)
cross_rh = mdgcal.cross_section(rh, st_point, ed_point)
cross_u = mdgcal.cross_section(u, st_point, ed_point)
cross_v = mdgcal.cross_section(v, st_point, ed_point)
cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point)
cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point)
cross_absv = mdgcal.cross_section(absv, st_point, ed_point)
cross_td = mdgcal.dewpoint_from_relative_humidity(cross_tmp, cross_rh)
pressure = mdgstda.gridstda_full_like_by_levels(cross_rh, cross_tmp['level'].values)
cross_theta = mdgcal.equivalent_potential_temperature(pressure, cross_tmp, cross_td)
cross_terrain = pressure - cross_psfc
cross_terrain.attrs['var_units'] = ''
# plot
draw_cross.draw_wind_theta_absv(cross_absv, cross_theta, cross_u, cross_v, cross_terrain, hgt, ...)</code></pre>
<h1>相当位温、相对湿度、沿剖面风</h1>
<p>相当位温、相对湿度、沿剖面风产品定制代码示例如下:</p>
<pre><code class="language-python">draw_wind_theta_rh(cross_rh, cross_theta, cross_u, cross_v, cross_terrain, hgt,
st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65),
h_pos=[0.125, 0.665, 0.25, 0.2],
rh_contourf_kwargs={}, theta_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">cross_rh</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">相对湿度</td>
</tr>
<tr>
<td style="text-align: left;">cross_theta</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">相当位温</td>
</tr>
<tr>
<td style="text-align: left;">cross_u</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">U风</td>
</tr>
<tr>
<td style="text-align: left;">cross_v</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">cross_terrain</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">hgt</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">高度</td>
</tr>
<tr>
<td style="text-align: left;">st_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">起始点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">ed_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">结束点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">lon_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">经度剖面</td>
</tr>
<tr>
<td style="text-align: left;">lat_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">纬度剖面</td>
</tr>
<tr>
<td style="text-align: left;">h_pos</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td>
<td style="text-align: left;">小图区域</td>
</tr>
<tr>
<td style="text-align: left;">map_extent</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">(50, 150, 0, 65)</td>
<td style="text-align: left;">绘图经纬度</td>
</tr>
<tr>
<td style="text-align: left;">rh_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">theta_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中barbs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">terrain_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv
from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d
from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d
from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D
from metdig.onestep.complexgrid_var.theta import read_theta3d
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
rh = get_model_3D_grid(...)
u = get_model_3D_grid(...)
v = get_model_3D_grid(...)
tmp = get_model_3D_grid(...)
hgt = get_model_grid(...)
psfc = get_model_grid(...)
cross_rh = mdgcal.cross_section(rh, st_point, ed_point)
cross_u = mdgcal.cross_section(u, st_point, ed_point)
cross_v = mdgcal.cross_section(v, st_point, ed_point)
cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point)
cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point)
cross_td = mdgcal.dewpoint_from_relative_humidity(cross_tmp, cross_rh)
pressure = mdgstda.gridstda_full_like_by_levels(cross_rh, cross_tmp['level'].values)
cross_theta = mdgcal.equivalent_potential_temperature(pressure, cross_tmp, cross_td)
cross_terrain = pressure - cross_psfc
cross_terrain.attrs['var_units'] = ''
# plot
draw_cross.draw_wind_theta_rh(cross_rh, cross_theta, cross_u, cross_v, cross_terrain, hgt, ...)</code></pre>
<h1>相当位温、绝对湿度、沿剖面风</h1>
<p>相当位温、绝对湿度、沿剖面风产品定制代码示例如下:</p>
<pre><code class="language-python">draw_wind_w_theta_spfh(cross_spfh, cross_theta, cross_t, cross_w, cross_terrain, hgt,
st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65),
h_pos=[0.125, 0.665, 0.25, 0.2],
spfh_contourf_kwargs={}, theta_contour_kwargs={}, wind_quiver_kwargs={},terrain_contourf_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">cross_spfh</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">绝对湿度</td>
</tr>
<tr>
<td style="text-align: left;">cross_theta</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">相当位温</td>
</tr>
<tr>
<td style="text-align: left;">cross_t</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">温度</td>
</tr>
<tr>
<td style="text-align: left;">cross_w</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">风</td>
</tr>
<tr>
<td style="text-align: left;">cross_terrain</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">hgt</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">高度</td>
</tr>
<tr>
<td style="text-align: left;">st_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">起始点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">ed_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">结束点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">lon_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">经度剖面</td>
</tr>
<tr>
<td style="text-align: left;">lat_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">纬度剖面</td>
</tr>
<tr>
<td style="text-align: left;">h_pos</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td>
<td style="text-align: left;">小图区域</td>
</tr>
<tr>
<td style="text-align: left;">map_extent</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">(50, 150, 0, 65)</td>
<td style="text-align: left;">绘图经纬度</td>
</tr>
<tr>
<td style="text-align: left;">spfh_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">theta_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">wind_quiver_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中quiver函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">terrain_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv
from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d
from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d
from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D
from metdig.onestep.complexgrid_var.theta import read_theta3d
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
rh = get_model_3D_grid(...)
u = get_model_3D_grid(...)
v = get_model_3D_grid(...)
tmp = get_model_3D_grid(...)
hgt = get_model_grid(...)
vvel = get_model_3D_grid(...)
spfh = get_model_3D_grid(...)
psfc = get_model_grid(...)
#calcu
w = mdgcal.vertical_velocity(vvel, tmp, spfh)
cross_rh = mdgcal.cross_section(rh, st_point, ed_point)
cross_u = mdgcal.cross_section(u, st_point, ed_point)
cross_v = mdgcal.cross_section(v, st_point, ed_point)
cross_w = mdgcal.cross_section(w, st_point, ed_point)
cross_t, cross_n = mdgcal.cross_section_components(cross_u, cross_v)
cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point)
cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point)
cross_td = mdgcal.dewpoint_from_relative_humidity(cross_tmp, cross_rh)
pressure = mdgstda.gridstda_full_like_by_levels(cross_rh, cross_tmp['level'].values)
cross_spfh = mdgcal.specific_humidity_from_dewpoint(pressure, cross_td)
cross_theta = mdgcal.equivalent_potential_temperature(pressure, cross_tmp, cross_td)
cross_terrain = pressure - cross_psfc
cross_terrain.attrs['var_units'] = ''
# plot
draw_cross.draw_wind_w_theta_spfh(cross_spfh, cross_theta, cross_t, cross_w*ratio, cross_terrain, hgt, ...)</code></pre>
<h1>相当位温、绝对湿度、沿剖面风`</h1>
<p>相当位温、绝对湿度、沿剖面风产品定制代码示例如下:</p>
<pre><code class="language-python">draw_wind_theta_spfh(cross_spfh, cross_theta, cross_u, cross_v, cross_terrain, hgt,
st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65),
h_pos=[0.125, 0.665, 0.25, 0.2],
spfh_contourf_kwargs={}, theta_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">cross_spfh</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">绝对湿度</td>
</tr>
<tr>
<td style="text-align: left;">cross_theta</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">相当位温</td>
</tr>
<tr>
<td style="text-align: left;">cross_u</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">U风</td>
</tr>
<tr>
<td style="text-align: left;">cross_v</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">cross_terrain</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">hgt</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">高度</td>
</tr>
<tr>
<td style="text-align: left;">st_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">起始点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">ed_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">结束点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">lon_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">经度剖面</td>
</tr>
<tr>
<td style="text-align: left;">lat_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">纬度剖面</td>
</tr>
<tr>
<td style="text-align: left;">h_pos</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td>
<td style="text-align: left;">小图区域</td>
</tr>
<tr>
<td style="text-align: left;">map_extent</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">(50, 150, 0, 65)</td>
<td style="text-align: left;">绘图经纬度</td>
</tr>
<tr>
<td style="text-align: left;">spfh_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">theta_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中barbs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">terrain_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv
from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d
from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d
from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D
from metdig.onestep.complexgrid_var.theta import read_theta3d
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
rh = get_model_3D_grid(...)
u = get_model_3D_grid(...)
v = get_model_3D_grid(...)
tmp = get_model_3D_grid(...)
hgt = get_model_grid(...)
psfc = get_model_grid(...)
cross_rh = mdgcal.cross_section(rh, st_point, ed_point)
cross_u = mdgcal.cross_section(u, st_point, ed_point)
cross_v = mdgcal.cross_section(v, st_point, ed_point)
cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point)
cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point)
cross_td = mdgcal.dewpoint_from_relative_humidity(cross_tmp, cross_rh)
pressure = mdgstda.gridstda_full_like_by_levels(cross_rh, cross_tmp['level'].values)
cross_spfh = mdgcal.specific_humidity_from_dewpoint(pressure, cross_td)
cross_theta = mdgcal.equivalent_potential_temperature(pressure, cross_tmp, cross_td)
cross_terrain = pressure - cross_psfc
cross_terrain.attrs['var_units'] = ''
# plot
draw_cross.draw_wind_theta_spfh(cross_spfh, cross_theta, cross_u, cross_v, cross_terrain, hgt, ...)</code></pre>
<h1>温度、相对湿度、水平风场</h1>
<p>温度、相对湿度、水平风场产品定制代码示例如下:</p>
<pre><code class="language-python">draw_wind_tmp_rh_vvel(cross_rh, cross_tmp, cross_u, cross_v, cross_vvel, cross_terrain, hgt,
st_point=None, ed_point=None, lon_cross=None, lat_cross=None, map_extent=(50, 150, 0, 65),
h_pos=[0.125, 0.665, 0.25, 0.2],
rh_contourf_kwargs={}, tmp_contour_kwargs={}, uv_barbs_kwargs={},terrain_contourf_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">cross_rh</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">水平散度</td>
</tr>
<tr>
<td style="text-align: left;">cross_tmp</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">相当位温</td>
</tr>
<tr>
<td style="text-align: left;">cross_u</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">U风</td>
</tr>
<tr>
<td style="text-align: left;">cross_v</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">cross_vvel</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
</tr>
<tr>
<td style="text-align: left;">cross_terrain</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">hgt</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">高度</td>
</tr>
<tr>
<td style="text-align: left;">st_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">起始点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">ed_point</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">结束点经纬度</td>
</tr>
<tr>
<td style="text-align: left;">lon_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">经度剖面</td>
</tr>
<tr>
<td style="text-align: left;">lat_cross</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">纬度剖面</td>
</tr>
<tr>
<td style="text-align: left;">h_pos</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">[0.125, 0.665, 0.25, 0.2]</td>
<td style="text-align: left;">小图区域</td>
</tr>
<tr>
<td style="text-align: left;">map_extent</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">(50, 150, 0, 65)</td>
<td style="text-align: left;">绘图经纬度</td>
</tr>
<tr>
<td style="text-align: left;">rh_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">tmp_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中barbs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">terrain_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv
from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d
from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d
from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D
from metdig.onestep.complexgrid_var.theta import read_theta3d
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
rh = get_model_3D_grid(...)
u = get_model_3D_grid(...)
v = get_model_3D_grid(...)
tmp = get_model_3D_grid(...)
hgt = get_model_grid(...)
psfc = get_model_grid(...)
cross_rh = mdgcal.cross_section(rh, st_point, ed_point)
cross_u = mdgcal.cross_section(u, st_point, ed_point)
cross_v = mdgcal.cross_section(v, st_point, ed_point)
cross_tmp = mdgcal.cross_section(tmp, st_point, ed_point)
cross_psfc = mdgcal.cross_section(psfc_bdcst, st_point, ed_point)
cross_u_t, cross_v_n = mdgcal.cross_section_components(cross_u, cross_v)
_, pressure = xr.broadcast(cross_rh, cross_tmp['level'])
cross_terrain = pressure - cross_psfc
cross_terrain.attrs['var_units'] = ''
# plot
draw_cross.draw_wind_tmp_rh(cross_rh, cross_tmp, cross_u, cross_v, cross_u_t, cross_v_n, cross_terrain, hgt, ...)</code></pre>
<h1>相当位温、相对湿度、水平风</h1>
<p>相当位温、相对湿度、水平风产品定制代码示例如下:</p>
<pre><code class="language-python">draw_time_rh_uv_theta(rh, u, v, theta,terrain,
rh_contourf_kwargs={}, uv_barbs_kwargs={}, theta_contour_kwargs={}, terrain_contourf_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">rh</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">相对湿度</td>
</tr>
<tr>
<td style="text-align: left;">theta</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">相当位温</td>
</tr>
<tr>
<td style="text-align: left;">u</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">U风</td>
</tr>
<tr>
<td style="text-align: left;">v</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">terrain</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
</tr>
<tr>
<td style="text-align: left;">rh_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">theta_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中barbs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">terrain_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv
from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d
from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d
from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D
from metdig.onestep.complexgrid_var.theta import read_theta3d
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
tmp = get_model_3D_grids(...)
u = get_model_3D_grids(...)
v = get_model_3D_grids(...)
rh = get_model_3D_grids(...)
td = mdgcal.dewpoint_from_relative_humidity(tmp, rh)
pressure = mdgstda.gridstda_full_like_by_levels(rh, rh['level'].values)
theta = mdgcal.equivalent_potential_temperature(pressure, tmp, td)
# plot
draw_cross.draw_time_rh_uv_theta(rh, u, v, theta, ...)</code></pre>
<h1>散度、垂直涡度、绝对湿度、水平风时间剖面</h1>
<p>散度、垂直涡度、绝对湿度、水平风时间剖面产品定制代码示例如下:</p>
<pre><code class="language-python">draw_time_div_vort_spfh_uv(div, vort, spfh, u, v, terrain,
spfh_contourf_kwargs={}, uv_barbs_kwargs={}, div_contour_kwargs={}, vort_contourf_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">div</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">水平散度</td>
</tr>
<tr>
<td style="text-align: left;">vort</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">垂直涡度</td>
</tr>
<tr>
<td style="text-align: left;">spfh</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">绝对湿度</td>
</tr>
<tr>
<td style="text-align: left;">u</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">U风</td>
</tr>
<tr>
<td style="text-align: left;">v</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">terrain</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">spfh_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中brabs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">div_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">vort_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv
from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d
from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d
from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D
from metdig.onestep.complexgrid_var.theta import read_theta3d
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
div, u, v = read_div_uv_4d(...)
vort, u, v = read_vort_uv_4d(...)
spfh = read_spfh_4D(...)
psfc = get_model_3D_grids(...)
_, pressure = xr.broadcast(v, v['level'])
terrain = pressure - psfc.values.repeat(pressure['level'].size, axis=1)
# plot
draw_cross.draw_time_div_vort_spfh_uv(div, vort, spfh, u, v, terrain, ...)</code></pre>
<h1>温度、温度平流、水平风时间剖面</h1>
<p>温度、温度平流、水平风时间剖面产品定制代码示例如下:</p>
<pre><code class="language-python">draw_time_wind_tmpadv_tmp(tmpadv, tmp, u, v, terrain,mean_area=None,
tmpadv_contour_kwargs={'levels': np.arange(-15, 15, 1)}, tmp_contourf_kwargs={}, uv_barbs_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">tmpadv</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">温度平流</td>
</tr>
<tr>
<td style="text-align: left;">tmp</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">温度</td>
</tr>
<tr>
<td style="text-align: left;">u</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">U风</td>
</tr>
<tr>
<td style="text-align: left;">v</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">terrain</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">mean_area</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">tmpadv_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">tmp_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中brabs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv
from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d
from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d
from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D
from metdig.onestep.complexgrid_var.theta import read_theta3d
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
u = get_model_3D_grids(...)
v = get_model_3D_grids(...)
psfc = get_model_3D_grids(...)
tmp = get_model_3D_grids(...)
tmpadv = mdgcal.var_advect(tmp, u, v)
_, pressure = xr.broadcast(v, v['level'])
terrain = pressure - psfc.values
terrain.attrs['var_units'] = ''
# plot
draw_cross.draw_time_wind_tmpadv_tmp(tmpadv, tmp, u, v, terrain, ...)</code></pre>
<h1>温度、垂直涡度平流、水平风时间剖面</h1>
<p>温度、垂直涡度平流、水平风时间剖面产品定制代码示例如下:</p>
<pre><code class="language-python">draw_time_wind_vortadv_tmp(vortadv, tmp, u, v, terrain,mean_area=None,
vortadv_contour_kwargs={}, tmp_contourf_kwargs={}, uv_barbs_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">vortadv</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">水平散度</td>
</tr>
<tr>
<td style="text-align: left;">tmp</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">相当位温</td>
</tr>
<tr>
<td style="text-align: left;">u</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">U风</td>
</tr>
<tr>
<td style="text-align: left;">v</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">terrain</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">mean_area</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">None</td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">vortadv_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">tmp_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中brabs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv
from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d
from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d
from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D
from metdig.onestep.complexgrid_var.theta import read_theta3d
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
vort, u, v = read_vort_uv_4d(...)
psfc = get_model_3D_grids(...)
tmp = get_model_3D_grids(...)
vortadv = mdgcal.var_advect(vort, u, v)
_, pressure = xr.broadcast(v, v['level'])
terrain = pressure - psfc.values
terrain.attrs['var_units'] = ''
# plot
draw_cross.draw_time_wind_vortadv_tmp(vortadv, tmp, u, v, terrain, ...)</code></pre>
<h1>散度、垂直涡度、相对湿度、水平风时间剖面</h1>
<p>散度、垂直涡度、相对湿度、水平风时间剖面产品定制代码示例如下:</p>
<pre><code class="language-python">draw_time_div_vort_rh_uv(div, vort, rh, u, v, terrain,
rh_contourf_kwargs={}, uv_barbs_kwargs={}, div_contour_kwargs={}, vort_contourf_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">div</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">水平散度</td>
</tr>
<tr>
<td style="text-align: left;">vort</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">相当位温</td>
</tr>
<tr>
<td style="text-align: left;">rh</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">相对湿度</td>
</tr>
<tr>
<td style="text-align: left;">u</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">U风</td>
</tr>
<tr>
<td style="text-align: left;">v</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">terrain</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">rh_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中brabs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">div_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">vort_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv
from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d
from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d
from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D
from metdig.onestep.complexgrid_var.theta import read_theta3d
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
div, u, v = read_div_uv_4d(...)
vort, u, v = read_vort_uv_4d(...)
rh = get_model_3D_grids(...)
psfc = get_model_3D_grids(...)
_, pressure = xr.broadcast(v, v['level'])
terrain = pressure - psfc.values.repeat(pressure['level'].size, axis=1)
# plot
draw_cross.draw_time_div_vort_rh_uv(div, vort, rh, u, v, terrain, ...)</code></pre>
<h1>温度 相对湿度 水平风 气压垂直速度</h1>
<p>温度 相对湿度 水平风 气压垂直速度产品定制代码示例如下:</p>
<pre><code class="language-python">draw_time_rh_uv_tmp_vvel(rh, u, v, tmp, vvel, terrain,
rh_contourf_kwargs={}, uv_barbs_kwargs={}, tmp_contour_kwargs={},vvel_contour_kwargs={},
**pallete_kwargs)</code></pre>
<p>根据输入的数据,生成叠加产品并输出到指定目录。</p>
<p><font size="4"><b>参数说明</b></font></p>
<table>
<thead>
<tr>
<th style="text-align: left;">参数</th>
<th style="text-align: left;">必选</th>
<th style="text-align: left;">默认</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">rh</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">相对湿度</td>
</tr>
<tr>
<td style="text-align: left;">u</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">U风</td>
</tr>
<tr>
<td style="text-align: left;">v</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">V风</td>
</tr>
<tr>
<td style="text-align: left;">vvel</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
</tr>
<tr>
<td style="text-align: left;">tmp</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">温度</td>
</tr>
<tr>
<td style="text-align: left;">terrain</td>
<td style="text-align: left;"><font color='red'>是</font></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">地形</td>
</tr>
<tr>
<td style="text-align: left;">rh_contourf_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contourf函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">tmp_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">uv_barbs_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中barbs函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;">vvel_contour_kwargs</td>
<td style="text-align: left;">否</td>
<td style="text-align: left;">{}</td>
<td style="text-align: left;">matplotlib中contour函数的可变参数</td>
</tr>
<tr>
<td style="text-align: left;"><code>**pallete_kwargs</code></td>
<td style="text-align: left;"></td>
<td style="text-align: left;"></td>
<td style="text-align: left;">详见[graphics层画板设置draw_compose](<a href="https://www.showdoc.com.cn/metdig/7407616276019947">https://www.showdoc.com.cn/metdig/7407616276019947</a> "graphics层画板设置draw_compose")</td>
</tr>
</tbody>
</table>
<table>
<thead>
<tr>
<th style="text-align: left;">返回值</th>
<th style="text-align: left;">说明</th>
</tr>
</thead>
<tbody>
<tr>
<td style="text-align: left;">ret</td>
<td style="text-align: left;">图像字典</td>
</tr>
</tbody>
</table>
<p><font size="4"><b>调用示例</b></font></p>
<pre><code class="language-python">from metdig.io import get_model_grid
from metdig.io import get_model_3D_grid
from metdig.io import get_model_3D_grids
from metdig.onestep.complexgrid_var.pv_div_uv import read_pv_div_uv
from metdig.onestep.complexgrid_var.div_uv import read_div_uv_4d,read_div_uv_3d
from metdig.onestep.complexgrid_var.vort_uv import read_vort_uv_4d
from metdig.onestep.complexgrid_var.spfh import read_spfh_4D,read_spfh_3D
from metdig.onestep.complexgrid_var.theta import read_theta3d
from metdig.onestep.complexgrid_var.w import read_w3d
from metdig.products import diag_crossection as draw_cross
import metdig.cal as mdgcal
import metdig.utl as mdgstda
# get data
tmp = get_model_3D_grids(...)
u = get_model_3D_grids(...)
v = get_model_3D_grids(...)
rh = get_model_3D_grids(...)
psfc = get_model_3D_grids(...)
pressure = xr.broadcast(v, v['level'])
terrain = pressure - psfc.values.repeat(pressure['level'].size, axis=1)
terrain.attrs['var_units'] = ''
# plot
draw_cross.draw_time_rh_uv_tmp(rh, u, v, tmp, terrain, ...)</code></pre>