关于grads批量提取nc4文件求助

hjx125128

我的gs文件
'reinit'
dy=20030101
while(dy<=20031231)

yy=000
while(yy<=2100)

if(yy<=900)

'sdfopen G:\GLDAS_NOAH_025_3H\2003\GLDAS_NOAH025_3H.A'dy'.0'yy'.021.nc4'
'set gxout fwrite'
'set fwrite G:\nc4dat\2003\GLDAS_NOAH025_3H.A'dy'.0'yy'.021.dat'
'set x 1044 1112'
'set y 348 368'
'set z 1 1'
'set t 1 1'
'd tveg_tavg'
'd wind_f_inst'
'd rainf_f_tavg'
'd tair_f_inst'
'd qair_f_inst'
'd swdown_f_tavg'
'd lwdown_f_tavg'
else
'reinit'
'sdfopen G:\GLDAS_NOAH_025_3H\2003\GLDAS_NOAH025_3H.A'dy'.0'yy'.021.nc4'
'set gxout fwrite'
'set fwrite G:\nc4dat\2003\GLDAS_NOAH025_3H.A'dy'.0'yy'.021.dat'
'set x 1044 1112'
'set y 348 368'
'set z 1 1'
'set t 1 1'
'd tveg_tavg'
'd wind_f_inst'
'd rainf_f_tavg'
'd tair_f_inst'
'd qair_f_inst'
'd swdown_f_tavg'
'd lwdown_f_tavg'
endif

yy=yy+300
endwhile

dy=dy+1

grads里的错误提示

hjx125128

river

先去看看各个时间的替换格式，哪有你这样写的······就算是fortran之类的这样写应该也不太合适吧·······你自己想一想也不对，光dy的循环就有一千多个文件，你真有那么多资料？
建议你在论坛搜索一下批量描述看看，会有启发的。还有你先用一个文件确定grads可不可以直接打开你的nc4的资料啊

hjx125128

river 发表于 2017-10-29 13:44
先去看看各个时间的替换格式，哪有你这样写的······就算是fortran之类的这样写应该也不太合适吧·· ...

已经提取出来了，我的数据是GLADS数据每三小时，放到一个文件夹里。这样一年就有2920个文件夹{:lol:}，现在已经成功的把它提取成dat形式了，想通过fortran把它转换成十进制以便后续的操作，但是我是刚学fortran，所以写不出来相应的fortran程序，大神有什么建议吗

hjx125128

hjx125128 发表于 2017-10-30 11:59
已经提取出来了，我的数据是GLADS数据每三小时，放到一个文件夹里。这样一年就有2920个文件夹{:lo ...

还有grads是可以直接操作nc4数据的

river

hjx125128 发表于 2017-10-30 11:59
已经提取出来了，我的数据是GLADS数据每三小时，放到一个文件夹里。这样一年就有2920个文件夹{:lo ...

哦，看来你研究的是大项目啊。Fortran还是好入门的。你不要一直看书，你要多看看别人现成的程序，理解意思，然后自己实践。只有实践才能最快的掌握fortran

小懒你不能废呀

请问直接用sdfopen语句就能打开nc4数据吗？我的无法打开是什么原因呢？

Linux

hjx125128 发表于 2017-10-30 12:01
还有grads是可以直接操作nc4数据的

您好，请问您是如何用grads打开nc4文件的。可不可以说明一下