关于循环数值传递

fcode

修改见下面代码以 !// 为注释的内容

[Fortran] 纯文本查看 复制代码

Program hello
  Implicit None
  Integer , parameter :: NN = 96232 , MM = 201801
  Real :: distance, row_3, k, dd(MM)!//dd用一维数组即可
  Integer :: row_1, row_2, i, ii
  Real, Allocatable :: lonx(:), laty(:), ilon(:), ilat(:)
!--------------------------------------------------------------------------
!a文件坐标点
!--------------------------------------------------------------------------
  Allocate (lonx(nn), laty(nn))

  Open (Unit=11, File='D:\气象数据\tmin\pointmatch_all.txt')
  Do i = 1, nn
    Read (11, *) row_1, row_2, lonx(i), laty(i), row_3
  End Do
  Close (Unit=11)
!----------------------------------------------------------------------------
!b文件坐标点
!----------------------------------------------------------------------------
  Allocate (ilon(mm), ilat(mm))
  Open (Unit=12, File='D:\气象数据\tmin\test.txt')
  Do ii = 1, mm
    Read (12, *) ilon(ii), ilat(ii), row_3
  End Do
  Close (Unit=12)
!----------------------------------------------------------------------------
!两个循环用来计算距离
!----------------------------------------------------------------------------
 !Allocate (dd(mm))!放入坐标点距离
  Open (Unit=13, File='D:\气象数据\tmin\test5.txt')

  Do i = 1, nn
    Do ii = 1, mm !//不减1
      dd(ii) = distance(lonx(i), laty(i), ilon(ii), ilat(ii))
    End Do
    Write (*, *) dd(mm)
!--------------------------------------------------------------------------
!输出B文件中最近点坐标
!--------------------------------------------------------------------------
    k = minloc(dd(:), 1) !//去掉1:ii
    !// 去掉If (ii/=k) Cycle
    Write (13, *) lonx(i), laty(i), ilon(k), ilat(k), dd(k)
  End Do
  Deallocate (lonx, laty, ilon, ilat)
End Program hello