3.2版本涡度散度计算不出来？

棉花猪仔

3.2版本涡度散度计算不出来？
用的是站点数据插值后hdivg    hcurl

1. fnu = 'D:/Project/smfkx/webapps/fkx/data/picture/test/data/925_u.dat'
   
2. fnv = 'D:/Project/smfkx/webapps/fkx/data/picture/test/data/925_v.dat'
   
3. ncol = numasciicol(fnu)
   
4. nrow = numasciirow(fnu)
   
5. ncol2 = numasciicol(fnv)
   
6. nrow2 = numasciirow(fnv)
   
7. a = asciiread(fnu,shape=(nrow,ncol))
   
8. a2 = asciiread(fnv,shape=(nrow2,ncol2))
   
9. lon = a[:,0]
   
10. lat = a[:,1]
    
11. u = a[:,2]
    
12. lon2 = a2[:,0]
    
13. lat2 = a2[:,1]
    
14. v = a2[:,2]
    
15. 
    
16. x = arange(109.1, 117.7,0.09)
    
17. y = arange(20, 26, 0.09)
    
18. prg = griddata((lon, lat), u, xi=(x, y), method='idw', radius=10)[0]
    
19. prg2 = griddata((lon2, lat2), v, xi=(x, y), method='idw', radius=10)[0]
    
20. 
    
21. divg = hdivg(prg, prg2)
    
22. vort = hcurl(prg, prg2)
    
23. 
    
24. print('divg: ' + str(divg))
    
25. print('vort: ' + str(vort))
    
26. 
    

复制代码

925_u.dat (767 Bytes, 下载次数: 0)

2021-10-29 12:43 上传

点击文件名下载附件

925_v.dat (788 Bytes, 下载次数: 1)

2021-10-29 12:43 上传

点击文件名下载附件

MeteoInfo

在 MeteoInfo -> pylib -> mipylib -> meteolib -> calc 目录中更新一下下面的文件：
kinematics.py (2.81 KB, 下载次数: 4)

2021-10-29 15:12 上传

点击文件名下载附件

脚本程序参考如下代码：

1. fnu = 'D:/Temp/ascii/925_u.dat'
   
2. fnv = 'D:/Temp/ascii/925_v.dat'
   
3. ncol = numasciicol(fnu)
   
4. nrow = numasciirow(fnu)
   
5. ncol2 = numasciicol(fnv)
   
6. nrow2 = numasciirow(fnv)
   
7. a = asciiread(fnu,shape=(nrow,ncol))
   
8. a2 = asciiread(fnv,shape=(nrow2,ncol2))
   
9. lon = a[:,0]
   
10. lat = a[:,1]
    
11. u = a[:,2]
    
12. lon2 = a2[:,0]
    
13. lat2 = a2[:,1]
    
14. v = a2[:,2]
    
15. 
    
16. x = arange(109.1, 117.7,0.09)
    
17. y = arange(20, 26, 0.09)
    
18. ug = griddata((lon, lat), u, xi=(x, y), method='idw', radius=10)[0]
    
19. vg = griddata((lon2, lat2), v, xi=(x, y), method='idw', radius=10)[0]
    
20. 
    
21. divg = meteo.divergence(ug, vg, x, y)
    
22. vort = meteo.vorticity(ug, vg, x, y)
    
23. 
    
24. print('divg: ' + str(divg))
    
25. print('vort: ' + str(vort))
    
26. 
    
27. subplot(2,1,1,axestype='map')
    
28. geoshow('cn_province')
    
29. contourf(x, y, vort, 20)
    
30. title('Vorticity')
    
31. colorbar(aspect=30)
    
32. subplot(2,1,2,axestype='map')
    
33. geoshow('cn_province')
    
34. contourf(x, y, divg, 20)
    
35. title('Divergence')
    
36. colorbar(aspect=30)

fnu = 'D:/Temp/ascii/925_u.dat'<br /> fnv = 'D:/Temp/ascii/925_v.dat'<br /> ncol = numasciicol(fnu)<br /> nrow = numasciirow(fnu)<br /> ncol2 = numasciicol(fnv)<br /> nrow2 = numasciirow(fnv)<br /> a = asciiread(fnu,shape=(nrow,ncol))<br /> a2 = asciiread(fnv,shape=(nrow2,ncol2))<br /> lon = a[:,0]<br /> lat = a[:,1]<br /> u = a[:,2]<br /> lon2 = a2[:,0]<br /> lat2 = a2[:,1]<br /> v = a2[:,2]<br /> <br /> x = arange(109.1, 117.7,0.09)<br /> y = arange(20, 26, 0.09)<br /> ug = griddata((lon, lat), u, xi=(x, y), method='idw', radius=10)[0]<br /> vg = griddata((lon2, lat2), v, xi=(x, y), method='idw', radius=10)[0]<br /> <br /> divg = meteo.divergence(ug, vg, x, y)<br /> vort = meteo.vorticity(ug, vg, x, y)<br /> <br /> print('divg: ' + str(divg))<br /> print('vort: ' + str(vort))<br /> <br /> subplot(2,1,1,axestype='map')<br /> geoshow('cn_province')<br /> contourf(x, y, vort, 20)<br /> title('Vorticity')<br /> colorbar(aspect=30)<br /> subplot(2,1,2,axestype='map')<br /> geoshow('cn_province')<br /> contourf(x, y, divg, 20)<br /> title('Divergence')<br /> colorbar(aspect=30)

棉花猪仔

MeteoInfo 发表于 2021-10-29 15:13
在 MeteoInfo -> pylib -> mipylib -> meteolib -> calc 目录中更新一下下面的文件：

可以了老师，非常感谢！