求助，使用描述文件批量处理NC文件后出图和原文件不一致

zzza

批量研究出来了，但是合并的后一个文件画出图来和原始文件直接出图数值有差距，大概差了15左右，请问啥原因啊出图命令都是
set cint 5 d skip(msl/100,1,1)

river

这看不出来啥问题，你还叠加一个矢量图更看不出来了。你要先对比一下形势对不对，如果高压区都变成了低压区那就根本不对。

zzza

river 发表于 2024-3-20 12:47
这看不出来啥问题，你还叠加一个矢量图更看不出来了。你要先对比一下形势对不对，如果高压区都变成了低压区 ...

谢谢river回复，我又试了一下，第一个文件的形势场都是符合的，但是第二个文件就差了15左右，先读取第一个文件，第二个文件出图会少15左右，先读取第二个文件数据的话，第一个文件出图会高15左右，不知道是哪里出问题了，麻烦river再帮忙看看呀。
批量ctl
dset D:\fogeradata\sum\%ch.nc
CHSUB 1 48 201308
CHSUB 49 96 201403
options yrev
options template
tdef time  96 linear 00Z22AUG2013 60mn
vars 3
u10=>u10  0  t,y,x  10 metre U wind component
v10=>v10  0  t,y,x  10 metre V wind component
msl=>msl  0  t,y,x  Mean sea level pressure
endvars

测试批量描述的gs
'reinit'
'xdfopen D:\fogeradata\sum\mslpl.ctl'
'set t 13'
'set cint 5'
'd skip(msl/100,1,1)'
'gxprint  D:\fogeradata\sum\test13082220.png white'

'c'
'set t 73'
'set cint 5'
'd skip(msl/100,1,1)'
'gxprint  D:\fogeradata\sum\test14033008.png white'
'reinit'
;


测试原始文件的gs
'reinit'
'sdfopen D:\fogeradata\sum\201308.nc'
'set t 13'
'set cint 5'
'd skip(msl/100,1,1)'
'gxprint  D:\fogeradata\sum\ystest13082220.png white'
'c'
'close 1'
'sdfopen D:\fogeradata\sum\201403.nc'
'set t 25'
'set cint 5'
'd skip(msl/100,1,1)'
'gxprint  D:\fogeradata\sum\ystest14033008.png white'
'reinit'
;

zzza

图片是这样的

river

zzza 发表于 2024-3-25 11:23
谢谢river回复，我又试了一下，第一个文件的形势场都是符合的，但是第二个文件就差了15左右，先读取第一 ...

好像没啥问题啊，你要不把x y z都写全了试试

zzza

river 发表于 2024-3-26 09:15
好像没啥问题啊，你要不把x y z都写全了试试

                               
登录/注册后可看大图

写了之后变成了这样==
dset D:\fogeradata\sum\%ch.nc
CHSUB 1 48 201308
CHSUB 49 96 201403
CHSUB 97 240 201407
CHSUB 241 288 201408
CHSUB 289 336 201506
CHSUB 337 432 201507
CHSUB 433 528 201508
CHSUB 529 576 201606
CHSUB 577 648 201707
CHSUB 649 672 201708
CHSUB 673 768 201808
CHSUB 769 816 202204
options yrev
options template
xdef  lon 241 linear 70 0.25
ydef lat 201 linear 10 0.25
tdef time  816 linear 00Z22AUG2013 60mn
vars 3
u10=>u10  0  t,y,x  10 metre U wind component
v10=>v10  0  t,y,x  10 metre V wind component
msl=>msl  0  t,y,x  Mean sea level pressure
endvars

river

zzza 发表于 2024-4-1 10:48
写了之后变成了这样==
dset D:\fogeradata\sum\%ch.nc
CHSUB 1 48 201308

就是你这个修改完，资料画出来对了吗？

zzza

river 发表于 2024-4-1 12:49
就是你这个修改完，资料画出来对了吗？

不对，像上面那样修改加上经纬度之后打不开了

                               
登录/注册后可看大图

river

zzza 发表于 2024-4-2 16:14
不对，像上面那样修改加上经纬度之后打不开了

建议先写单个文件的描述文件，验证正确以后，再改写成批量的

zzza

river 发表于 2024-4-3 12:19
建议先写单个文件的描述文件，验证正确以后，再改写成批量的

我把描述文件中的经纬度写成了全称，没有之前的报错信息了，单个文件测试正常了，批量描述也可以，但是还是跟之前一样有出图15的差值
修改后的单个描述文件：
dset D:\fogeradata\msp\201308msp.nc
title single daily  Reanalysis
options yrev
dtype netcdf
xdef longitude 241 linear 70 0.25
ydef latitude 201 linear 10 0.25
tdef time 48 linear 00Z22AUG2013 60mn
vars 3
u10=>u10  0  t,y,x  10 metre U wind component
v10=>v10  0  t,y,x  10 metre V wind component
msl=>msl  0  t,y,x  Mean sea level pressure
endvars

但是上面描述文件没有垂直维度的描述，我想是否要写完整才行，再加上zdef level 1 levels 1000（或者zdef level 1 linear 0 1），但都会出现类似的错误提示信息
Scanning Descriptor File:  D:\fogeradata\sum\mslpltest.ctl
gadsdf: Lev dimension pressure is not an SDF dimension.
SDF Descriptor file D:\fogeradata\sum\mslpltest.ctl was not successfully opened & parsed.

这个zdef应该怎么写呢
不知道是不是因为没有zdef导致出图差值吗，还是没有这个也不影响，是其他地方错了
始末是这样的：我一开始是看了您批量描述的帖子跟着做，做的过程中发现会出现8楼回复的那个错误，我就在论坛搜到了妞胖胖这位贴主的帖子，https://bbs.06climate.com/forum. ... FA%C1%BF&page=1，我发现我的错误跟她的很相似，她解决了，所以我也用她的解决方法试了一下，批量描述成功了，但是测试发现出图不一样，所以实在不晓得是哪里错了====