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

river

zzza 发表于 2024-4-3 14:17
我把描述文件中的经纬度写成了全称，没有之前的报错信息了，单个文件测试正常了，批量描述也可以，但是还 ...

你用你这个描述文件单独描述你测试用的那两个资料，出来的图都是正常的了，是这个意思吧？还有就是直接 d msl/100 ，不使用skip函数再看看
实在不行，你这个资料不太大的话，传你测试的那两个资料上来我看看

zzza

river 发表于 2024-4-4 12:26
你用你这个描述文件单独描述你测试用的那两个资料，出来的图都是正常的了，是这个意思吧？还有就是直接 d ...

是的，单独描述测试正常，批量也看起来没有问题，我用msl/100也试了，但是出图还是有差值，文件上传了，麻烦您方便的话帮我看看，十分感谢
单独的ctl
dset D:\fogeradata\sum\201308.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
批量描述的ctl
dset D:\fogeradata\sum\%ch.nc
CHSUB 1 48 201308
CHSUB 49 96 201403
options template
options yrev
dtype netcdf
xdef longitude 241 linear 70 0.25
ydef latitude 201 linear 10 0.25
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

river

zzza 发表于 2024-4-4 21:19
是的，单独描述测试正常，批量也看起来没有问题，我用msl/100也试了，但是出图还是有差值，文件上传了， ...

找到问题了，先上图






你对比一下橙红色的部分，数据存储过程中使用的影响因子和偏移量是不同的，所以打开的时次不同，默认的影响因子和偏移量是不同的。先打开t=13的话，就是默认使用了第一个资料的影响因子和偏移量，导致第二个资料里t=73时次出现偏差。反之，同理。
我暂时没有找到通过修改描述文件解决问题的方法。但有一个笨办法，比较麻烦。需要你先写一个完全版的描述文件，让他的影响因子和偏移量默认为1.0和0.0，然后需要你把每个资料的影响因子和偏移量都找出来，然后display的时候再计算进去。
dset D:\fogeradata\sum\%ch.nc
CHSUB 1 48 201308
CHSUB 49 96 201403
title single daily  Reanalysis
undef -32767
dtype netcdf
*UNPACK scale_factor add_offset 0.000378815 -0.0524355
*UNPACK scale_factor add_offset 0.000430352 -0.424341
*UNPACK scale_factor add_offset 0.0751844 102148
options template
options yrev
xdef  241 linear 70 0.25
ydef  201 linear 10 0.25
tdef  96 linear 00Z22AUG2013 60mn
zdef 1 linear 0 1
vars 3
u10  1  t,y,x  10 metre U wind component
v10  1  t,y,x  10 metre V wind component
msl  1  t,y,x  Mean sea level pressure
endvars


'reinit'
'open D:\fogeradata\sum\mslpl.ctl'
'set t 13'
'set cint 5'
'd skip( (msl*0.0640412+100364)/100,1,1)'
'gxprint  D:\fogeradata\sum\test13082220.png white'


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

zzza

river 发表于 2024-4-5 18:26
找到问题了，先上图

原来是这样，谢谢您啦，找到问题了我再研究一下，十分十分感谢

river

zzza 发表于 2024-4-6 10:50
原来是这样，谢谢您啦，找到问题了我再研究一下，十分十分感谢

研究出来了记得会来分享经验

zzza

river 发表于 2024-4-6 14:30
研究出来了记得会来分享经验

没有研究出来什么办法，直接把文件合并了

river

zzza 发表于 2024-4-20 20:38
没有研究出来什么办法，直接把文件合并了

合并以后，画图就都正常了？
那也挺好，能解决问题就是好办法

zzza

river 发表于 2024-4-20 22:07
合并以后，画图就都正常了？
那也挺好，能解决问题就是好办法

是的，正常了，谢谢river耐心解答呀