- 积分
- 226
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2020-2-26
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
我画了一张K指数的填色图,但是最后出来的图基本上就一个颜色,想问一下这个要怎么修改。
我的脚本是这样
f1= addfile('E:/ncdata/85008t.nc')
f2=addfile('E:/ncdata/85008q1.nc')
T850 = f1['t'][4,:,:] - 273.16
rh = f2['r'][4,:,:]
Td850 = T850-((14.55+0.114*T850)*(1-0.01*rh) + pow((2.5+0.007*T850)*(1-0.01*rh),3) + (15.9+0.117*T850)*pow((1-0.01*rh),14))
f3 =addfile('E:/ncdata/70008t.nc')
f4=addfile('E:/ncdata/70008rh.nc')
T700 = f3['t'][4,:,:] - 273.16
rh1 = f4['r'][4,:,:]
Td700 = T700-((14.55+0.114*T700)*(1-0.01*rh1) + pow((2.5+0.007*T700)*(1-0.01*rh1),3) + (15.9+0.117*T700)*pow((1-0.01*rh1),14))
f5=addfile('E:/ncdata/50008t.nc')
T500 = f5['t'][4,:,:] - 273.16
K = T850-T500+Td850-(T700-Td700)
axesm()
mlayer = shaperead('F:\Program Files (x86)\MeteoInfo\Map\country1.shp')
geoshow(mlayer, edgecolor='black')
layer = contourfm(K, cmap='grads_rainbow')
colorbar(layer)
title('K index 2019.8.5.08')
|
|