1. 回复关键词获取相关的代码:
怀俄明探空站数据解算PWV和Tm:怀俄明探空站数据解算PWV和Tm
怀俄明多线程下载探空站数据(包括检查和下载遗漏数据的代码):怀俄明多线程下载
对IGRAv2进行质量控制得到PWV和Tm的 matlab 代码:IGRAv2进行质量控制得到PWV和Tm
算 IGRAv2 探空站的 Tm 和 PWV:IGRAv2计算Tm和PWV
提取探空站 IGRAv2 全部数据:
ERA5 解算合集(温度、气压、PWV、水汽压和 Tm)代码获取:ERA5合集
2. 代码仅仅适用于怀俄明下载链接的新网站
如下图,从新网站下载探空站数据存的txt文件:
![图片](https://i-blog.csdnimg.cn/blog_migrate/3ef126b7a9e67748202b1d989e6631dd.png)
我使用第一列(气压)、第二列(位势高)、第三列(温度)、第五列(相对湿度)计算得到PWV和Tm。计算方法跟前面博文IGRA计算PWV和Tm类似,这里不再赘述。
使用Matlab计算IGRAv2探空站的Tm和PWV
这两个代码无缝连接,直接使用。如果是别的代码下载的怀俄明数据,请自行修改。
3. 数据已进行质量控制:
1.顶层至少达到300 hPa
2.地面有数据可用
3.在低于(高于)1000hPa的台站,地面以上至少有五(四个)个标准压力水平
4.连续记录温度或湿度之间的压力间差值大于200 hPa的剔除。
注意怀俄明站点的变化
53772站点在不同时间他的经度竟然是不一样的,纬度没有任何变化,离谱。如果经纬度和位势高度都变可以理解为换位置了,可是纬度和位势高没有任何变化,由于不确定哪个经度正确,这个站点我弃用了。
![图片](https://i-blog.csdnimg.cn/blog_migrate/a66e1c729a5e74ea866326c5d4cd1834.jpeg)
4. Matlab代码压缩包内容如下:
![图片](https://i-blog.csdnimg.cn/blog_migrate/c0776525b2cae68b7f09216fad8a9c5a.png)
data:为数据集,下载好的探空站数据放里面。
![图片](https://i-blog.csdnimg.cn/blog_migrate/d5d2752ad83ae76a0aafcf34bc1b285d.png)
![图片](https://i-blog.csdnimg.cn/blog_migrate/c9718bc44359aa5ba26719f7d1a72b5e.png)
探空站经纬度.xlsx:这个是探空站的经纬度和ID。
![图片](https://i-blog.csdnimg.cn/blog_migrate/2930dac68cf1160598998cc0c7af4c35.png)
把数据集改成我的格式,直接运行readsounding.m即可。结果存为sounding_PWV_Tm.mat,每一行是一个站点(第一列为探空站id,第二列为时间和PWV),如下图所示。
![图片](https://i-blog.csdnimg.cn/blog_migrate/bdfa9309b151ce0c17e3e0bb7fa86140.png)
![图片](https://i-blog.csdnimg.cn/blog_migrate/e9240973d19e4f097e6c0362c89fed41.png)
前提是跟我数据格式一样.前提是跟我数据格式一样.前提是跟我数据格式一样.因为每个人存数据的方式不一致,使用的网址也不同,格式不一致请参考并自行修改。代码的获取关键词:怀俄明探空站数据解算PWV和Tm。