grads输出的txt格式文件可以控制输出格式吗？

冬日情愫+

我想将grads的输出文件格式变成txt格式的，但是文件中数据的输出格式是8*34+1（每行8个数据供34行，外加最后一行是1个数据的），但是我想文件中数据输出格式是21*13（每行21个数据13行），可以实现吗？
下面的数据是2.5*2.5的。
'reinit'
'sdfopen  e:\data\air.1998.nc'
'sdfopen  e:\data\rhum.1998.nc'
'set lon 80 130'
'set lat 15 45'
'set lev 850'
'set mpdset cnworld'
'enable print E:\data\sitase.gmf'
i=131
while(i<=134)
'set t 'i
'set gxout contour'
'define t0=air.1'
'define rh=rhum.2'
'define prs=lev'
'define es=(6.112*exp(17.67*(t0-273.15)/(t0-29.65)))'
'define q=rh*(0.62197*es/(prs-es))/100.'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(t0)-log(e)-4.805)'
'define theta=t0*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define esl=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'
'set cint 4'
'set gxout contour'
'set cthick 5'
'd esl'
'draw title sitase-'i' '
' print'
'set gxout print'
say i
file="E:\data\sitase.txt"
'd esl'
rc=write(file,result,append)
'c'
i=i+1
endwhile
'disable print'
;

mofangbao

这个可能有些困难，grads在d的时候直接把该维度上的资料都显示出来了，暂时还没想到好的方法，还是用fortran格式化一下吧

言深深

看来我的grads得从这边学起了

冬日情愫+

mofangbao 发表于 2011-11-18 08:37

                               
登录/注册后可看大图

这个可能有些困难，grads在d的时候直接把该维度上的资料都显示出来了，暂时还没想到好的方法，还是用fortra ...

目前我也只能用这种方法。

兰溪之水

O(∩_∩)O哈哈~见【兰溪金盆洗手帖】~

冬日情愫+

对呀，所以在此多谢兰溪帮我解决了问题。受益匪浅呀！

麻花

{:eb502:}