爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 16036|回复: 3

[源代码] 提取wrfout输出的nc文件的格点数据

[复制链接]

新浪微博达人勋

发表于 2021-7-9 10:41:08 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
最近完成一个小脚本分享一下:
完成了输入一经纬度点,匹配最近网格点经纬度,输出气象要素,计算空气密度等输出文本,代码如下。(输入输出路径读取了配置文件)
另外想请教有没有好一点方法,把10m风速差值到70m,现在我是参考最简单风廓线公式来的。


config.ini

112 Bytes, 下载次数: 31, 下载积分: 金钱 -5

配置文件

getnc.py

4.15 KB, 下载次数: 43, 下载积分: 金钱 -5

源代码

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2021-8-3 17:46:37 | 显示全部楼层
大佬,
有没有方法提取wrfout里面的XLAT、XLONG Times PH PHB U V W T P PB PBLH UST等参数呢
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-8-10 18:26:21 | 显示全部楼层
小吴吴 发表于 2021-8-3 17:46
大佬,
有没有方法提取wrfout里面的XLAT、XLONG Times PH PHB U V W T P PB PBLH UST等参数呢

我这个脚本里就已经提取了XLAT、XLONG Times  其余参数对应修改就行。
   
    time = data.variables['Times'][:]
    lat = data.variables['XLAT'][:]
    lon = data.variables['XLONG'][:]
    temp = data.variables['T2'][:] - 273.15  #℃
    swdown = data.variables['SWDOWN'][:]  #w/m2
    press = data.variables['PSFC'][:]/100   #hpa
    wind_u=data.variables['U10'][:]
    wind_v=data.variables['V10'][:]
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-8-10 18:42:14 | 显示全部楼层
YanPenG 发表于 2021-8-10 18:26
我这个脚本里就已经提取了XLAT、XLONG Times  其余参数对应修改就行。
   
    time = data.variables ...

好的谢谢大哥
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表