Fortran读入nc数据出错

hsuchueng · 发表于 2021-7-23 15:07:21

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 hsuchueng 于 2021-7-23 15:09 编辑

大家好，最近用Fortran读nc文件时遇到了些问题。
我利用下面这串代码想读入海温和经纬度信息，经纬度变量能够成功读进来，并且和nc文件里面的值一样。但是海温SST却无法读进来，SST_RES里的数据全是0.00000000E+0000这样的，nc文件中的SST数据是float型请问这该如何解决？
real (r8), dimension(1:320,1:384,1:365) :: SST_RES
real (r8), dimension(1:384) :: lat_RES
real (r8), dimension(1:320) :: lon_RES
integer (int_kind) :: status_RES, ncid_RES, varid_RES
status_RES = nf90_open(".plus1.nc",nf90_nowrite,ncid_RES)
status_RES = nf90_inq_varid(ncid_RES,"sst",varid_RES)
status_RES = nf90_get_var(ncid_RES,varid_RES,SST_RES)
status_RES = nf90_inq_varid(ncid_RES,"lat",varid_RES)
status_RES = nf90_get_var(ncid_RES,varid_RES,lat_RES)
status_RES = nf90_inq_varid(ncid_RES,"lon",varid_RES)
status_RES = nf90_get_var(ncid_RES,varid_RES,lon_RES)
status_RES = nf90_close(ncid_RES)

蒙同学_/╱°. · 发表于 2022-5-11 17:25:11

我也想学会Fortran啊~~~~一切从0开始

hsuchueng · 发表于 2021-7-23 16:40:07

啊啊啊，我崩溃了，不小心一试竟然发现问题了，虽然还不知道为什么。
是这样的，我的SST数据是365x384x320的，但是用nf90_get_var读入SST数据，他的维度不知为何就变换成320x384x365了，所以我定义的数组SST_RES的维度就和读入的不一致，最终导致错误了。虽然我还是不知道为什么用nf90_get_var之后会把我数据的维度给转换了。

人可山 · 发表于 2021-7-24 10:16:04

进来学习的，没啥贡献

蓝色蜗牛角 · 发表于 2021-7-25 21:23:21

哈哈哈我看到金币进来白嫖了

kaiofchina · 发表于 2021-8-25 14:41:14

帮楼主顶一下

baishizhou · 发表于 2021-8-26 10:46:04

建议先用matlab看一下该数据头的信息，尤其是温度数据的数据类型，以及缺省值是什么，有没有offset 和scale

有消息 · 发表于 2021-9-5 09:08:11

帮楼主顶一下

顷刻花 · 发表于 2022-4-8 21:19:57

dddddddddddd

susus · 发表于 2022-4-30 17:59:18

帮楼主顶一下

		自动登录	找回密码
密码			立即注册

[求助] Fortran读入nc数据出错

登录后查看更多精彩内容~

回帖奖励 +3 金钱

回帖奖励 +3 金钱

回帖奖励 +3 金钱

回帖奖励 +3 金钱

回帖奖励 +3 金钱