nc转grd出现错误

麦田_smile

我用的是2.5*2.5月平均风场资料。想把nc文件变为grd文件以便fortran的处理。但是出现错误（两者画出的图不一致），贴上来请教一下诸位
gs文件：
'reinit'
'sdfopen e:\uv\uwnd.mon.mean.nc'   
'set gxout fwrite'
'set fwrite e:\uv\u_850.grd'        
'set lon 0 357.5'
'set lat -90 90'
tt=121                     
while(tt<=781)                    
'set t ' tt                          
'set lev 1000'                     
'd uwnd'                     
'set lev 925'                       
'd uwnd'                     
'set lev 850'                       
'd uwnd'                     
'set lev 700'                     
'd uwnd'                        
'set lev 600'                       
'd uwnd'                     
'set lev 500'                       
'd uwnd'                        
'set lev 400'
'd uwnd'   
'set lev 300'
'd uwnd'   
'set lev 250'
'd uwnd'               
'set lev 200'
'd uwnd'
'set lev 150'
'd uwnd'
'set lev 100'
'd uwnd'
'set lev 70'
'd uwnd'
'set lev 50'
'd uwnd'
'set lev 30'
'd uwnd'
'set lev 20'
'd uwnd'
'set lev 10'
'd uwnd'
tt=tt+1
endwhile

'disable fwrite'
;

lqouc

'reinit'
'sdfopen e:\uv\uwnd.mon.mean.nc'
'set fwrite e:\uv\_u_850.grd'
'set gxout fwrite'
'set x 1 144
'set y 1 73'        
tt=121                     
while(tt<=781)                    
z1=1
while(z1<=17)
'set t 'tt''
'set z 'z1''
'd uwnd'                  
z1=z1+1
endwhile
tt=tt+1
endwhile
'disable fwrite'
'reinit'   
这样写肯定没问题，屡试不爽

麦田_smile

ctl文件
dset e:\uv\u_850.grd
undef 32766
title monthly mean U wind from the NCEP Reanalysis
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef    17 levels
    1000 925     850     700     600     500     400     300     250     200
     150 100  70 50 30 20 10
tdef 661 linear jan1958 1mon
vars 1
uwnd 1 -999 monthly mean U wind
endvars
;

麦田_smile

第一张图是nc画出的 第二是grd画出的

麦田_smile

这是第一张图 楼上上传失败了

lqouc

很奇怪啊，既然你会用循环，为什么不把z也循环，非要一个个写lev？
至于ctl感觉不大对，在对变量的描述中应该改成uwnd 17 99吧，你的风有17层呢。

麦田_smile

lqouc 发表于 2013-3-11 19:33

                               
登录/注册后可看大图

很奇怪啊，既然你会用循环，为什么不把z也循环，非要一个个写lev？
至于ctl感觉不大对，在对变量的描述中应 ...

哦 谢谢你的回复 我之前用过Z循环 不过的得到的文件竟然有1.2g 本来nc文件只有200m额。

麦田_smile

lqouc 发表于 2013-3-11 19:33

                               
登录/注册后可看大图

很奇怪啊，既然你会用循环，为什么不把z也循环，非要一个个写lev？
至于ctl感觉不大对，在对变量的描述中应 ...

这是把lev也加进循环的gs
'reinit'
'sdfopen e:\uv\uwnd.mon.mean.nc'
'set fwrite e:\uv\_u_850.grd'
'set gxout fwrite'
'set x 0 357.5'
'set y -90 90'        
tt=121                     
while(tt<=781)                    
'set t ' tt
z1=1
while(z1<=17)
'set z ' z1
'd uwnd'                  
z1=z1+1
endwhile
tt=tt+1
endwhile
'disable fwrite'
'reinit'                        
;

麦田_smile

lqouc 发表于 2013-3-11 19:53

                               
登录/注册后可看大图

'reinit'
'sdfopen e:\uv\uwnd.mon.mean.nc'
'set fwrite e:\uv\_u_850.grd'

谢谢您 请问ctl文件除了你指出的外，还有别的问题吗

lqouc

麦田_smile 发表于 2013-3-11 20:00

                               
登录/注册后可看大图

谢谢您 请问ctl文件除了你指出的外，还有别的问题吗

别的我也看不出来啊，莫非还是不对？