grads批量合并有规律的二进制文件（如NC逐日资料等）和批量提取特定时间的特定要素

river

海洋里的鱼 发表于 2019-4-14 15:25
楼主，请问这样的错误是怎么回事呢？
这是我的gs文件
'reinit'

说得哼明确了，if 的函数结构有问题，嵌套出问题了，最主要是放的位置不对。你自己看看吧

RainyYuan

楼主好，想请教一个问题，就是我现在想把多个nc文件整合成一个，nc文件是某个变量从2000-2017年的季度平均，是每年一个文件。然后grd文件已经出来了，编写好ctl后放入grads中，发现只有当t等于1时才能出结果，不知道是哪一步做错。
'reinit'
'set gxout fwrite'
'set fwrite D:\aa\MAM.grd'
tt=2000
while(tt<=2017)
'sdfopen

D:\aa\NDVI\g4.timeAvg.MOD13C2_006_CMG_0_05_Deg_Monthly_NDVI.'tt'.SEASON_MAM.100E_25N_110E_35N.nc'
'set x 1 200'
'set y 1 200'
'd mod13c2_006_cmg'
'close 1'
tt=tt+1
endwhile
'disable fwrite'
;

dset D:\aa\MAM.grd
title mean monthly NDVI （2000-2017）
undef -3000
xdef 200 linear 100.025 0.05
ydef 200 linear 25.025 0.05
zdef 1 linear 0 1
tdef 18 linear jan2000 1yr
vars 1
mod13c2_006_cmg 18 99 mean monthly NDVI
endvars

river

RainyYuan 发表于 2019-5-8 11:45
楼主好，想请教一个问题，就是我现在想把多个nc文件整合成一个，nc文件是某个变量从2000-2017年的季度平均 ...

你的资料没有介绍清楚，原始资料只有一层？每个资料只有一个时次？
如果是这样你的描述文件为什么写有18层  mod13c2_006_cmg 18 99 mean monthly NDVI

冬日情愫+

以前经常用楼主这种方法整合多个nc文件,由于不知道怎么描述多个nc文件一个ctl描述。学习一下nc文件的ctl描述文件如何写，自己尝试写了一下发现不对，在论坛上搜索一下，看到楼主的帖子，追根溯源看看相关经典好贴！

river

冬日情愫+ 发表于 2019-8-14 10:21
以前经常用楼主这种方法整合多个nc文件,由于不知道怎么描述多个nc文件一个ctl描述。学习一下nc文件的ctl描 ...

你们还在搞科研啊，羡慕。我都放弃了······

youhuiqi7455

你好，请教一个问题，假如我的文件名是198001，198002....198101 这样的37年逐月资料该怎么写文件名呢？

river

youhuiqi7455 发表于 2020-4-30 18:37
你好，请教一个问题，假如我的文件名是198001，198002....198101 这样的37年逐月资料该怎么写文件名呢？

看我另一个帖子，用批量描述的方法更简单更方便

一片酸甜叶

非常感谢！虽然一开始明明是照着写的，但是不知道为什么会被“unknown command 2”刷屏，但后面自己捣鼓捣鼓着又好了

Hooly！

你好，想请教一下啊，我用这个画程序出来的降水的图有些地方和原nc文件不一样，有缺测，数值也有些地方对不上，可以帮我看一下错在哪吗？谢谢啦

'reinit'
'sdfopen C:\grads\CNrain\CN05.prec.19790101-20200816.nc'
'set fwrite c:\grads\cnrain\jul.grd'
'set gxout fwrite'
'set lat 26 36'
'set lon 110 123'
yy=1979
while(yy<=2020)
if ((math_mod(yy,4)=0&math_mod(yy,100)!=0)|math_mod(yy,400)=0)
tt=182
while(tt<=213)
'set t 'tt''
'd p'
tt=tt+1
endwhile
else
tt=181
while(tt<=212)
'set t 'tt''
'd p'
tt=tt+1
endwhile
endif
yy=yy+1
endwhile
'disable fwrite'
;

这是ctl文件
dset c:\grads\cnrain\jul.grd
title daily prec in jul from 1979 to 2020
undef -999000000.000000
xdef 27 linear 110 0.5
ydef 21 linear 26 0.5
zdef 1 levels 0 1
tdef 1302 linear 00Z01jul1979 1440mn
vars 1
p 0 99  mean t,y,x,p
endvars

river

Hooly！ 发表于 2020-9-16 21:35
你好，想请教一下啊，我用这个画程序出来的降水的图有些地方和原nc文件不一样，有缺测，数值也有些地方对不 ...

我不懂你这个是什么逻辑？看着又不是合并文件，循环年份是为了什么？如果只是提取自己需要的区域的资料，判断是否闰年有什么作用？