在wrfinput中增加新的数组

cici破晓碎寒 · 发表于 2021-7-12 09:24:35

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

您需要登录才可以下载或查看，没有帐号？立即注册

x

大家好，最近我在修正wrf中的辐射分量，需要对不同的辐分量乘除修正数组，这个数组已经用高程数据算出，是二维的，grd格式，并且插值到模式模拟的网格点了，现在需要读入这个数组在模式运行过程中加以纠正辐射通量，我的想法是在wrfinput.nc中加入这几个数组，但是我没有找到相关读取数据的.F程序，请问大家有什么好的建议吗？谢谢！

ljchen1989 · 发表于 2021-7-12 14:00:24

这样改动还挺大的吧
ncl或者nco可还行

cici破晓碎寒 · 发表于 2021-7-12 14:23:23

ljchen1989 发表于 2021-7-12 14:00
这样改动还挺大的吧
ncl或者nco可还行

是在NCL中把这个数组添到wrfinput.nc中吗？这样的话，后续模式读取的时候也要做一些处理吧？

ljchen1989 · 发表于 2021-7-12 16:59:18

cici破晓碎寒发表于 2021-7-12 14:23
是在NCL中把这个数组添到wrfinput.nc中吗？这样的话，后续模式读取的时候也要做一些处理吧？

所以你为啥要在wrfinput里面加变量？加变量还是改变量？

cici破晓碎寒 · 发表于 2021-7-12 17:03:30

加数组，是一个已经用高程数据算好的并插值到模式格点(nx,ny)上的数组，加了之后在运行模式的时候，在辐射模块用这个数组修正相应的通量

astiny · 发表于 2021-7-12 18:47:33

本帖最后由 astiny 于 2021-7-12 18:49 编辑

你最好还是不要直接改wrfinput吧，在别的地方开一个接口把数据接进去。如果硬要搞的话，可以试试用mozbc,这个是wrfchem里用来给初始场和边界条件加入气溶胶和化学物质的，感觉应该也可以用来做你想要的那种效果。
原版：https://www.acom.ucar.edu/wrf-chem/download.shtml
github上改版的支持HVC的版本：https://github.com/fengx7/mozbc_for_WRFv3.9

cici破晓碎寒 · 发表于 2021-7-13 09:12:22

astiny 发表于 2021-7-12 18:47
你最好还是不要直接改wrfinput吧，在别的地方开一个接口把数据接进去。如果硬要搞的话，可以试试用mozbc,这 ...

谢谢呀！！！我参考一下~~~

dqjmorpheus · 发表于 2021-7-13 14:46:14

有个简单粗暴的方法：转成文本格式，然后在所用的辐射方案里直接打开，就用fortran的open命令打开文件。

dqjmorpheus · 发表于 2021-7-13 14:48:46

dqjmorpheus 发表于 2021-7-13 14:46
有个简单粗暴的方法：转成文本格式，然后在所用的辐射方案里直接打开，就用fortran的open命令打开文件。

这样每个辐射方案的时步都会打开和关闭这个文件，如果数据量不是很大的话，不会影响I/O速度。

cici破晓碎寒 · 发表于 2021-7-13 20:42:02

dqjmorpheus 发表于 2021-7-13 14:48
这样每个辐射方案的时步都会打开和关闭这个文件，如果数据量不是很大的话，不会影响I/O速度。

主要是，模式在物理方案计算的时候，已经并行了，分成许多块，我的数据是整个domain的

		自动登录	找回密码
密码			立即注册

在wrfinput中增加新的数组

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

浏览过的版块