请用fortran的编写一个100*100*100的格点，服从高斯分布的代码。

以下是使用 Fortran 编写的生成 100*100*100 格点，服从高斯分布的代码示例：

program generate_gaussian_grid
    implicit none
    
    integer, parameter :: N = 100  ! 格点数量
    
    real(kind=8), parameter :: mu = 0.0  ! 均值
    real(kind=8), parameter :: sigma = 1.0  ! 标准差
    
    real(kind=8), dimension(N, N, N) :: grid  ! 存储生成的格点
    integer :: i, j, k
    
    call random_seed()  ! 初始化随机数种子
    
    ! 生成格点
    call generate_gaussian_grid(grid, mu, sigma)
    
    ! 打印部分格点值
    do i = 1, 5
        do j = 1, 5
            do k = 1, 5
                write(*, '(F8.4)') grid(i, j, k)
            end do
        end do
        write(*, *)
    end do
    
end program generate_gaussian_grid


subroutine generate_gaussian_grid(grid, mu, sigma)
    implicit none
    real(kind=8), dimension(:,:,:), intent(out) :: grid  ! 存储生成的格点
    real(kind=8), intent(in) :: mu, sigma  ! 正态分布的参数
    integer :: i, j, k
    
    real(kind=8) :: u1, u2, z
    
    do i = 1, size(grid, 1)
        do j = 1, size(grid, 2)
            do k = 1, size(grid, 3)
                call random_number(u1)
                call random_number(u2)
                
                z = sqrt(-2.0 * log(u1)) * cos(2.0 * acos(-1.0) * u2)
                
                grid(i, j, k) = mu + sigma * z
            end do
        end do
    end do
    
end subroutine generate_gaussian_grid

在这个代码中，`generate_gaussian_grid` 子程序使用直接生成高斯分布随机数的方法生成了一个 100*100*100 的格点。在每次循环中，通过调用 `random_number` 生成两个均匀分布的随机数 u1 和 u2。然后，根据直接生成高斯分布的公式，计算出服从标准正态分布的随机数 z。最后，根据给定的均值 mu 和标准差 sigma，计算出最终的服从高斯分布的随机数，并存储到 grid 数组中。

在主程序中，你可以根据需要修改打印部分格点值的代码，以输出你感兴趣的部分或整个格点的值。