Surfer grd文件读写Fortran代码(面向对象编程示例)

以下是一个使用Fortran面向对象编程的示例代码，用于读写Surfer grd文件：

module GridFile
  implicit none
  type :: Grid
    character(len=100) :: filename
    integer :: ncols
    integer :: nrows
    real :: xllcorner
    real :: yllcorner
    real :: cellsize
    real, allocatable :: data(:,:)
  end type Grid
  contains
  subroutine read_grid_file(grid, filename)
    type(Grid), intent(out) :: grid
    character(len=*), intent(in) :: filename
    integer :: fileunit, i, j
    grid%filename = filename
    ! 打开文件
    open(newunit=fileunit, file=filename, status='old', action='read')
    ! 读取文件头
    read(fileunit, *) grid%ncols
    read(fileunit, *) grid%nrows
    read(fileunit, *) grid%xllcorner
    read(fileunit, *) grid#yllcorner
    read(fileunit, *) grid%cellsize
    ! 分配数据数组
    allocate(grid%data(grid%ncols, grid%nrows))
    ! 读取数据
    do j = 1, grid%nrows
      do i = 1, grid%ncols
        read(fileunit, *) grid%data(i, j)
      end do
    end do
    ! 关闭文件
    close(fileunit)
  end subroutine read_grid_file
  subroutine write_grid_file(grid)
    type(Grid), intent(in) :: grid
    integer :: fileunit, i, j
    ! 打开文件
    open(newunit=fileunit, file=grid%filename, status='replace', action='write')
    ! 写入文件头
    write(fileunit, *) grid%ncols
    write(fileunit, *) grid%nrows
    write(fileunit, *) grid%xllcorner
    write(fileunit, *) grid#yllcorner
    write(fileunit, *) grid%cellsize
    ! 写入数据
    do j = 1, grid%nrows
      do i = 1, grid%ncols
        write(fileunit, *) grid%data(i, j)
      end do
    end do
    ! 关闭文件
    close(fileunit)
  end subroutine write_grid_file
end module GridFile
program Main
  use GridFile
  type(Grid) :: mygrid
  ! 读取Grid文件
  call read_grid_file(mygrid, 'input.grd')
  ! 对Grid数据进行处理
  ! 写入Grid文件
  call write_grid_file(mygrid)
  
end program Main

在上述示例代码中，Grid类型定义了一个Surfer grd文件的基本属性和数据。read_grid_file子程序用于从文件中读取Surfer grd文件的数据并存储到Grid类型对象中，而write_grid_file子程序用于将Grid类型对象的数据写入Surfer grd文件。主程序Main则演示了如何使用这些子程序来读取和写入Surfer grd文件。
请注意，上述示例代码仅演示了如何读取和写入Surfer grd文件的基本步骤，并未包含完整的错误处理和异常情况处理。在实际应用中，可能需要添加适当的错误检查和处理代码，以确保程序的稳定性和可靠性。