关于清风大佬的双线性插值程序读取负值的疑惑

979162804

最近刚好用到清风大佬的这个fortran自编双线性插值fortran程序，我自己加入了读取数据等步骤，然后发现当格点数据为正时，插值到站点以后的结果都是正确的（虽然有一点点误差），但是当我用到格点数据为负时，出来的数据误差特别大，仔细看了下好像是负值不支持？我自己改了一天还是没改出个结果，特来请教

mofangbao

说句鸡汤：在快要坚持不下去的时候再坚持一下，很大可能就成功了。

979162804

汗，我重新把权重部分理了一遍，竟然通了，只需要一些地方加abs进行绝对值就可以了！！，还是感谢清风大佬的程序，真的很牛

979162804

mofangbao 发表于 2020-2-28 12:21
说句鸡汤：在快要坚持不下去的时候再坚持一下，很大可能就成功了。

大佬我想问下，我现在已经把micaps格点数据插值到站点上（用fortran处理），现在需要进行读取时间段来做时间平均。想到了用ncl可能比较方便点毕竟有函数。
现在想了3个办法，但前面两个都遇到了麻烦
1.将插值完的站点数据批量转为netcdf格式，用ncl来处理（站点数据转为nc格式比较困难，不太懂这方面）
2.继续用fortran来处理，并参考大佬您的那个fortran批量读取的教程。（fortran好像对这种超多文件的时间处理比较繁琐？）
3.全部推倒重来，直接用ncl读入micaps数据进行插值到站点，并直接输出为nc文件，方便后续进行时间平均等步骤

979162804

mofangbao 发表于 2020-2-28 12:21
说句鸡汤：在快要坚持不下去的时候再坚持一下，很大可能就成功了。

另外我现在插值完的站点数据格式是
站点号stid   lat   lon   数据

mofangbao

979162804 发表于 2020-3-6 15:32
大佬我想问下，我现在已经把micaps格点数据插值到站点上（用fortran处理），现在需要进行读取时间段来做 ...

第二种是比较常规的处理方式，遇到问题应该是你对Fortran应用还不熟练，这种场景应该是很基础的了，多思考尝试吧，应该能搞出来的。