如何在运行时分配Netcdf fortran数组的大小？

program dynamic_array_allocation
  use netcdf
  implicit none
  
  integer :: ncid, varid, ndims, dimids(3), dimlen(3)
  real, allocatable :: data(:,:,:)
  
  ! 打开Netcdf文件
  call nf_open("example.nc", NF_NOWRITE, ncid)
  
  ! 获取变量ID
  call nf_inq_varid(ncid, "temperature", varid)
  
  ! 获取维度信息
  call nf_inq_varndims(ncid, varid, ndims)
  call nf_inq_vardimid(ncid, varid, dimids)
  
  ! 获取维度大小
  do i = 1, ndims
    call nf_inq_dimlen(ncid, dimids(i), dimlen(i))
  end do
  
  ! 分配数组内存
  allocate(data(dimlen(1), dimlen(2), dimlen(3)))
  
  ! 读取数据
  call nf_get_var(ncid, varid, data)
  
  ! 关闭Netcdf文件
  call nf_close(ncid)
  
  ! 使用数组进行后续操作
  ! ...
  
  ! 释放数组内存
  deallocate(data)
  
end program dynamic_array_allocation