求用fortran或者matlab计算地闪密度的程序！急1不甚感激！

dl凋零

我有算法，但不会编程，我的思路是：把一地方分成0.1*0.1的网格，然后以网格中心为评估点。我好像没讲清楚，我的思路和我在一篇文献里看到的大致一样。如下：格子的中心位置为评估点所在经纬度。每个格子边长为3ｋｍ，每边共有９个格子。算法具体步骤如下：
评估点附近网格分布图：网格是9*9的网格
（１）在程序中声明４个二维数组：①ＡｒｒａｙＮｇ，②ＡｒｒａｙＡｖｇＬｉｇＣｕｒｔ，③ＡｒｒａｙＢｏｘｓＬｎｇ，④Ａｒｒａｙ－ＢｏｘｓＬａｔ。分别用来存放图１中每个格子的①Ｎｇ值，②雷电流平均强度，③格子中心点经度，④格子中心点纬度。
（２）已知用户输入的评估点经纬度，而且根据地球上任意两点的经纬度又可以得到两点间距离。计算公式如下：
ｄ（ｘ１，ｙ１，ｘ２，ｙ２）＝ｒａｒｃｃｏｓ［ｓｉｎ（ｘ１）ｓｉｎ（ｘ２）＋ｃｏｓ（ｘ１）ｃｏｓ（ｘ２）ｃｏｓ（ｙ１－ｙ２）］
其中：ｘ１，ｙ１是纬度、经度的弧度单位，ｒ 为地球半径。
据此，可以经过计算得出每个格子的经纬度范围，ＡｒｒａｙＢｏｘｓＬｎｇ，ＡｒｒａｙＢｏｘｓＬａｔ能够通过循环成功赋值。
（３）逐行读取经过处理的闪电定位系统（ＬＬＳ）数据，取闪电的经纬度存入局部变量ｌｎｇＴｍｐ、ｌａｔＴｍｐ。使用ｉｆ语句判断该行闪电定位数据是否属于图１中所示９×９＝８１个格子。如果属于，则ＡｒｒａｙＮｇ数组中第ｉ行ｊ列格子数据执行自加运算（ＡｒｒａｙＮｇ［ｉ，ｊ］＋＋；）；同时，地闪强度也累加（ＡｒｒａｙＡｖｇＬｉｇＣｕｒｔ［ｉ，ｊ］＋＝ ＬｉｇｈｔｎｉｎｇＣｕｒｒｅｎｔ；）。

（４）对数据源循环一遍后，使用ｆｏｒ循环语句重新对８１个格子进行循环赋值。使用公式为：平均电流强度＝ 总电流强度／闪击次数
                     Ｎｇ ＝ 总地闪次数／资料年数／格子面积
经过以上４步，即可得出每个格子的Ｎｇ值和雷

lqouc

同志，给你点提示好了。看你连发两个确实也挺着急的。
1.直接搜索论坛，有现成的话你就用，或者改改再用。
2.自己先把程序编出来，如果有问题，发上来大家给你改。
3.你直接开价吧，出钱的话就会有人写了。