fnl资料提取的dat数据用fortran无法处理

river

虹hhh 发表于 2023-4-10 00:03
找到问题了 ，用grads将nc转dat，时间这一维为1时（set t 1）正常,但是set t 2以后就提示Low Level I/O E ...

那就是只提取出来一层数据。你提取转换数据这一步就有问题，就不用说下一步了。

虹hhh

river 发表于 2023-4-10 12:08
那就是只提取出来一层数据。你提取转换数据这一步就有问题，就不用说下一步了。

这是我的grads程序，因为没有学习过grads，自学的可能程序上还有问题，麻烦老师看看错误在哪
'reinit'
'sdfopen C:\Users\ayqz\Desktop\energycalculate.nc'
'set gxout fwrite'
'set fwrite C:\Users\ayqz\Desktop\energy.dat'
mt=1
mz=1
'set undef -9.99e+33'
while (mt<=32)
'set t 'mt''
while(mz<=21)
'set z 'mz''
'set lon 90 150'
'set lat 20 60'
'd z '
mz=mz+1
endwhile
mt=mt+1
endwhile'disable fwrite'
'reinit'
时间 set t 2以后就没有数据了，lon 、lat、高度层设置其他值都有数据
ctl
dset C:\Users\ayqz\Desktop\energy.dat
undef -9.99e+33
title ERA5 z data
xdef  241 linear 90 0.25
ydef  161 linear 20 0.25
zdef 21 levels 1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100
tdef 32 linear 00Z28DEC2020 3hr
vars   1
z 21 99 geopotential
Endvars

river

虹hhh 发表于 2023-4-10 15:04
这是我的grads程序，因为没有学习过grads，自学的可能程序上还有问题，麻烦老师看看错误在哪
'reinit'
...

脚本看上去没啥大问题，我怀疑是脚本最后的 reinit 有问题，你换成 ; 试试，注意英文和中文标点的区别

river

虹hhh 发表于 2023-4-10 15:04
这是我的grads程序，因为没有学习过grads，自学的可能程序上还有问题，麻烦老师看看错误在哪
'reinit'
...

哦，不对不对，上面那个回复可能不是根本原因。你的循环有问题，mz是内循环啊，每次都要从1开始循环到21 ，你把它放在最前面，mz=mz+1是不是就没作用了。mz=1要写在while(mz<=21)前面。
我也有点儿懵，你自己试试吧

虹hhh

river 发表于 2023-4-11 13:26
哦，不对不对，上面那个回复可能不是根本原因。你的循环有问题，mz是内循环啊，每次都要从1开始循环到21  ...

问题解决了，谢谢老师