fortran调用GMT API

GMT（Generic Mapping Tools）是一款开源的地图制图软件，提供了丰富的地图制图功能和API接口。Fortran可以通过调用GMT的API接口来实现地图制图。

以下是一个简单的Fortran程序，通过调用GMT的API接口绘制了一张简单的地图：

program gmt_example
  implicit none
  integer :: i, j
  integer :: iret
  real :: lon(3), lat(3), z(3)
  character(len=80) :: command

  ! 创建GMT工作区
  command = "gmt begin gmt_example"
  call system_trim(command, iret)

  ! 设置地图区域和投影方式
  command = "gmt basemap -R-180/180/-90/90 -Jm1"
  call system_trim(command, iret)

  ! 绘制三角形
  lon = (/ -100.0, -80.0, -90.0 /)
  lat = (/ 30.0, 20.0, 10.0 /)
  z = (/ 1.0, 2.0, 3.0 /)
  command = "gmt psxy -R -J -W1p -Gred -K -O <<EOF"
  do i = 1, 3
     write(command, "(3f10.5)") lon(i), lat(i), z(i)
  end do
  command = trim(command) // "EOF"
  call system_trim(command, iret)

  ! 添加标题和注释
  command = "gmt pstext -R -J -N -F+f16p+jTL -O <<EOF"
  command = trim(command) // "-100.0 30.0 Example"
  call system_trim(command, iret)

  ! 结束GMT工作区
  command = "gmt end"
  call system_trim(command, iret)
end program gmt_example

在程序中，我们首先使用`gmt begin`命令创建GMT工作区，然后使用`gmt basemap`命令设置地图区域和投影方式。接着，我们使用`gmt psxy`命令绘制一个三角形，并使用`gmt pstext`命令添加标题和注释。最后，我们使用`gmt end`命令结束GMT工作区。

需要注意的是，在Fortran中调用GMT API需要使用`call system_trim(command, iret)`来执行GMT命令，并检查返回值`iret`是否为0，以确保命令执行成功。

编译程序时需要链接GMT库，可以使用以下命令：

$ gfortran gmt_example.f90 -o gmt_example -lgmt

运行程序，将会生成一个名为`gmt_example.ps`的PostScript格式地图文件。