小编Jes*_*seH的帖子

我是使用Fortran编写的程序的研究员。我具有非常基本的编码技能，因此需要一些帮助来使一些代码正确编译。

在显示代码之前，我将提供一些背景知识。我正在处理大量数据，这将需要64位编译和大于2 gb的内存。我在代码中注意到的第一件事是，许多变量被写为“实数”，但是在我的研究中，我发现“双精度”允许更大的变量，并且将是一个更灵活的选择，因此我更改了所有从“真实”变量到“双精度”变量。

fortran构建文件“ dist.f”的编译中还包含一个文件，称为“ geocoord.inc”。我发现变量已保存到一个公共块中，但是再一次，我需要一些可以容纳大量数据的东西。正如我一直相信的那样，模块将是一个更好的程序。在转换此包含文件以使其与模块程序正常工作时，我需要一些建议，我将在下面列出。

Dist.f：

c Convert latitude and longitude to kilometers relative
c to center of coordinates by short distance conversion.

subroutine dist(xlat, xlon, xkm, ykm)

implicit none

c   Parameters:
double precision    xlat, xlon  ! (input)
double precision        xkm, ykm    ! (output)

c   Local variables:
double precision lat1, lat2, lat3
double precision    q
double precision    xx
double precision    yp

include "geocoord.inc"

c Set up short distance conversion by subr. SETORG
  q=60*xlat-olat
  yp=q+olat
  lat1=datan(rlatc*dtan(RAD*yp/60.0))
  lat2=datan(rlatc*dtan(RAD*OLAT/60.0))
  LAT3=(LAT2+LAT1)/2.
  xx=60*xlon-olon  ! …