求助！如何基于 ERA-Interim 再分析数据资料计算出水汽通量数据，后导入中Arcgis？

周一的07

一、目前的成果：已基于 ERA-Interim 再分析数据资料使用grads画出了夏季多年平均水汽输送通量场，如图所示。
二、想要达到的目的：想将计算出的水汽通量资料导入Arcgis中，最终得到多年平均水汽输送通量栅格数据。（用于提取指定点的水汽输送通量数据）
三、思考：
查阅到一篇相关的文献，如附件pdf所示：基于GIS的汉江流域水汽输送与降水时空模式研究_梁文谦。
该篇文献中提到，首先使用grads计算得到旬水汽通量矢量场，后将将汉江流域19个气象站点的旬水汽通量值、旬可降水量、旬实测降水量数据导入到ArcGIS中，并利用普通克里金插值法得出整个汉江流域的水汽输送通量、可降水量、实测降水量空间分布。疑问：我们已得到了大范围流域的水汽通量数据，作者是通过什么办法得到研究流域的19个气象站点具体的相关数据的？我的另一思考：我已得到整层的水汽输送通量，那我是否可以通过先得到水汽输送通量的栅格数据，再基于研究流域的矢量边界，剪裁或者掩膜提取出研究流域的水汽输送通量的空间分布呢？（但目前我不知道该如何将水汽通量数据在gis中转成栅格数据）
四、问题：
1、如何将水汽通量数据在gis中转成栅格数据？
2、是否是先在grads计算出水汽输送通量，后直接将数据导入gis中，那么水汽通量应该保存为什么格式的数据呢？
在写毕业论文时卡在这儿了，希望各位前辈解答疑惑，不胜感激！

周一的07

咦，这个附件要如何设置，可以不需要扣除金钱就能下载。

river

首先，很遗憾我不会用GIS，所以用GIS转换栅格数据啥的我不太懂，不过思路是一样的，都是导入站点的数据，然后通过插值来得到空间分布。
其次，根据我个人的理解，来捋一捋作者的思路。文章一共使用了两种资料，一是ncep再分析资料，这个资料应该是用来计算整层水汽通量、可降水量，然后通过插值到站点的函数（gr2stn）插值到汉江流域19个站点上，导出这个资料备用。二是汉江流域19个站点实测资料，这个资料就是用来计算汉江流域19个站点实测降水量的，计算完备用。  最后就是把上述两个备用的站点资料导入到arcis软件进行相关的计算以后，通过克里金插值得到空间分布。

周一的07

river 发表于 2021-10-9 18:24
首先，很遗憾我不会用GIS，所以用GIS转换栅格数据啥的我不太懂，不过思路是一样的，都是导入站点的数据，然 ...

感谢river回复！我会再顺着这个思路再思考一下，如果后面成功做出来了，会分享一下。

river

周一的07 发表于 2021-10-10 11:51
感谢river回复！我会再顺着这个思路再思考一下，如果后面成功做出来了，会分享一下。

加油！上面那些就是我个人的理解，我建议最好可以和做过相关工作的师兄师姐沟通沟通。

影伊

你好，请教一下。我在画水汽通量场的时候出错。不知道错在了哪里。用到了ERA-interim资料的surface pressure、q（相对湿度）、u、v这四个要素计算的，公式如下
'define uq=(1000*shum.3)*u.1'
'define vq=(1000*shum.3)*v.2'
'define uqs=vint(0.01*pres.4(lev=0),uq,300)'
'define vqs=vint(0.01*pres.4(lev=0),vq,300)'

周一的07

影伊 发表于 2022-1-21 10:09
你好，请教一下。我在画水汽通量场的时候出错。不知道错在了哪里。用到了ERA-interim资料的surface pressur ...

看看GrADS软件里是如何报错的

影伊

周一的07 发表于 2022-1-21 10:29
看看GrADS软件里是如何报错的

没有报错，但是画出来的箭头朝下

周一的07

影伊 发表于 2022-1-21 12:12
没有报错，但是画出来的箭头朝下

没报错的话证明编程应该没有问题