帮忙看看chatgpt给的读取nc文件的fortran脚本

啥也没研究出来

如图需要读取一个名为msak.nc的9个变量，编译成功后，print每一个变量都为0，但是nc文件的值大都为1，极少数的值为0.3，格式是float32
找不到在哪里出错了，脚本部分也贴一下

include 'netcdf.inc'
integer :: status
integer :: ncid
integer :: varid
real,dimension(ny,nx,9) :: nc_data      ! 存储所有变量数据的数组

print*,"get mask nc"
print*,nx,ny
status = nf_open(trim('../input/mask.nc'), NF_NOWRITE, ncid)
status = nf_inq_varid(ncid, "so2", varid)
status = nf_get_var_real(ncid, varid, nc_data(:,:,1))
status = nf_inq_varid(ncid, "nox", varid)
status = nf_get_var_real(ncid, varid, nc_data(:,:,2))
status = nf_inq_varid(ncid, "co", varid)
status = nf_get_var_real(ncid, varid, nc_data(:,:,3))
status = nf_inq_varid(ncid, "voc", varid)
status = nf_get_var_real(ncid, varid, nc_data(:,:,4))
status = nf_inq_varid(ncid, "bc", varid)
status = nf_get_var_real(ncid, varid, nc_data(:,:,5))
status = nf_inq_varid(ncid, "oc", varid)
status = nf_get_var_real(ncid, varid, nc_data(:,:,6))
status = nf_inq_varid(ncid, "pm10", varid)
status = nf_get_var_real(ncid, varid, nc_data(:,:,7))
status = nf_inq_varid(ncid, "pm25", varid)
status = nf_get_var_real(ncid, varid, nc_data(:,:,8))
status = nf_inq_varid(ncid, "hn3", varid)
status = nf_get_var_real(ncid, varid, nc_data(:,:,9))
status = nf_close(ncid)

print*,nc_data(:,:,1)