源代码——Cressman逐步订正法

倚天照海

不错，很好

涛声依旧否

顶起原创，顶。。。。

二斤幸福

下下来尝试一下

novelbean

好的，不错

Jack_TEA

好动西，先看看

hanqianli

利用版主提供的程序，带入以下数据进行插值，插值结果误差极大，请版主赐教：
51111,  115.35,   36.02,  274.42
51112,  116.59,   36.41,  284.35
51113,  117.06,   36.15,  184.61
51114,  118.09,   36.11,  263.00
52843,  119.07,   36.37,  271.77
52852,  120.42,   36.56,  267.32
53871,  122.22,   36.57,  252.61
53916,  116.51,   35.34,  278.68
55936,  118.50,   35.35,  267.55
55945,  119.32,   35.23,  259.71
56299,   92.16,   32.06,   92.77
51578,   88.55,   29.13,  149.68
52680,   89.36,   28.55,  134.74
53773,   89.05,   27.44,   83.58
54018,   94.43,   33.09,  113.03
55029,   96.45,   33.06,  131.42
56033,   98.08,   34.57,   81.87

部分插值结果：
     819.168800      810.563200      336.694100      344.277300
      116.099400      260.446100      249.674800      780.884000
      775.408000      780.460000      773.279100      783.133700
      262.149500      264.578700    0.000000E+00    0.000000E+00
    0.000000E+00    0.000000E+00    0.000000E+00    0.000000E+00
    0.000000E+00    0.000000E+00    0.000000E+00      116.795700
      337.594500      330.164000      814.167800      840.429400
     1759.299000     1747.698000      804.626400      804.944200
      344.863200      116.099400      263.569500      780.543300
      776.904200     1824.288000     1811.670000     1826.071000
      776.251600      782.601100      784.981900      260.438300
    0.000000E+00    0.000000E+00    0.000000E+00    0.000000E+00
    0.000000E+00    0.000000E+00    0.000000E+00    0.000000E+00
    0.000000E+00      113.259200      343.363700      359.300400
      832.994100      819.628700      816.086400      797.195400
      336.751000      344.702000      116.099400      246.216800
      783.753700     1826.057000     1821.608000     1825.256000
     3901.254000     1825.454000     1825.722000      783.951400
      266.081100    0.000000E+00    0.000000E+00    0.000000E+00
    0.000000E+00    0.000000E+00    0.000000E+00    0.000000E+00
    0.000000E+00    0.000000E+00    0.000000E+00      122.923200
      126.408700      362.742600      362.585500      342.292300
      341.962200      338.306600      111.362500      116.099400
      246.180500      784.084900     1826.899000     1822.166000
     1826.877000     1818.907000     3904.172000     1826.289000
      784.084100      265.766700    0.000000E+00    0.000000E+00
    0.000000E+00    0.000000E+00    0.000000E+00    0.000000E+00
    0.000000E+00    0.000000E+00    0.000000E+00    0.000000E+00
    0.000000E+00      123.400700      123.141800      116.976800
      346.213500      107.817600      113.242100      116.099400
    0.000000E+00      263.888300      781.675400      777.577500
附件grid.txt为插值结果原文件

hanqianli

           if(kk>0)         ge(i,j)=ge(i,j)+sum/kk              ! 格点初值
这句是不是应该改为           if(kk>0)         ge(i,j)=sum/kk              ! 格点初值

和清谈

hanqianli 发表于 2014-1-14 16:39

                               
登录/注册后可看大图

利用版主提供的程序，带入以下数据进行插值，插值结果误差极大，请版主赐教：
51111,  115.35,   36.02,   ...

你不妨修改一下精度、半径这些可调指数，还有要考虑下站点所处的位置，像中国西部的站点本来比较少，插值到比较精细的点时，误差肯定是比较大的

和清谈

hanqianli 发表于 2014-1-15 14:32

                               
登录/注册后可看大图

if(kk>0)         ge(i,j)=ge(i,j)+sum/kk              ! 格点初值
这句是不是应该改为           if(kk>0)         ge(i, ...

这是为了累加ge(i,j)

zengxiaoqing

必须要顶一下