读取两类资料并从中选取经纬度近似的点输出变量，流程.....

麒麟

程序大致如此，全面的在附件中可以找到。想实现读取两类资料并从中选取经纬度近似的点输出变量。遇见问题在：
1。  if((abs(lonin(l,m)-line(2)).le.0.06).and.(abs(latin(l,m)-line(3)).le.0.06))then                                   !判断两类资料的经纬度间隔最小以输出变量        
                       write(6,'(3i5,3f12.3)')l,m,n,latin(l,m),lonin(l,m), pre_wrf(l,m,n)
这样输出受到外层循环i,j,k的影响，使得输出文件中不停的循环。
请教大家，如何在寻找经纬度相近两点的情况下，输出时间（120次）合适的变量？

百事可乐

楼主，这是你的代码逻辑。别人不清楚的话，很难帮你改正哟。

li913

你的文字叙述部分，我没看懂，建议整理一下你的问题，表述清楚。是否：从文件A读取坐标a，在文件B中搜索与a最近的点b并输出相关信息？
如果文件不大，建议先读取所有文件，再处理。

麒麟

li913 发表于 2017-2-19 16:55
你的文字叙述部分，我没看懂，建议整理一下你的问题，表述清楚。是否：从文件A读取坐标a，在文件B中搜索与a ...

恩，是的，我的问题没说清楚。你理解的意思对了，就是想从文件A中读取所要a点的坐标，再到文件B中搜索与a点坐标最经的坐标点b并输出相关信息。下面语句已经能够输出我想要的信息，但由于处在读取a点坐标的循环内，输出的信息会一致循环。其实这整个程序分别是两个独立程序合在一起的。
if((abs(lonin(l,m)-line(2)).le.0.06).and.(abs(latin(l,m)-line(3)).le.0.06))then
write(6,'(3i5,3f12.3)')l,m,n,latin(l,m),lonin(l,m), pre_wrf(l,m,n)

麒麟

百事可乐 发表于 2017-2-19 15:13
楼主，这是你的代码逻辑。别人不清楚的话，很难帮你改正哟。

怪我没说清楚，我想实现想从文件A中读取所要a点的坐标，再到文件B中搜索与a点坐标最经的坐标点b并输出相关信息。

li913

就问题本身，不难，可是你没有描述清楚，期望的结果是什么样子。