grads fwrite写入出错

Rplay

各位大神，小弟最近对一次天气过程做了模拟，用arwpost处理后得到了相应的文件，小弟希望先计算假相当位温，然后将假相当位温的计算结果输入一个二进制文件以便于fortran进行处理。小弟的计算结果貌似是正确的，但是写入文件时一直有问题。屡次尝试均以失败告终。希望各位grads高手给小弟看一下，gs文件中哪里出现了问题。不胜感激。

'reinit'
'open fd01'
'set grads off'
'set grid off'
'set fwrite THSE.grd'
'set gxout fwrite'
'set t 1'
'set lev 1000 100'
'define es=6.11*exp(21.87*tc/(tc+265.49))'
'define qs=0.622*es/lev-0.378*es)'
'define qa=rh*qs/100.0'
'define thse=(tc+273.15)*exp(0.28586*log(1000.0/lev)+2500.0*qa/(338.52-0.24*(tc+273.15)+1.24*td))'
zz=1
while(zz<20)
'set z 'zz
'd thse'
zz=zz+1
endwhile
'disable fwrite'

风飘摇

'define qs=0.622*es/lev-0.378*es)'这句是不是少半个括号？

Rplay

风飘摇 发表于 2013-8-5 16:24

                               
登录/注册后可看大图

'define qs=0.622*es/lev-0.378*es)'这句是不是少半个括号？

是的。。。但是有这半个括号依然有问题

river

失败是怎么失败，有没有错误提示？楼主要提供一定的线索让其他人去帮你定位问题所在。这只是一个短短gs文件，别人还能从头看到尾。那要是一个几万行的甚至几十万行的程序，楼主也说遇到问题，屡屡尝试解决不得，然后把长长的程序贴上来，让别人从头看到末尾？那是不现实的。

才仔

超出能力范围了

Rplay

river 发表于 2013-8-5 19:00

                               
登录/注册后可看大图

失败是怎么失败，有没有错误提示？楼主要提供一定的线索让其他人去帮你定位问题所在。这只是一个短短gs文件 ...

谢谢提醒。。。失败之处就是写进文件去数据的不是我计算得到的数据，都是很大很大的，大小在1E35这个样子的

river

Rplay 发表于 2013-8-5 22:23

                               
登录/注册后可看大图

谢谢提醒。。。失败之处就是写进文件去数据的不是我计算得到的数据，都是很大很大的，大小在1E35这个样子 ...

你画图出来看见的，还是fortran读出来的？

Rplay

river 发表于 2013-8-5 22:43

                               
登录/注册后可看大图

你画图出来看见的，还是fortran读出来的？

图画出来和fortran读出来都很大

Rplay

Rplay 发表于 2013-8-5 22:51

                               
登录/注册后可看大图

图画出来和fortran读出来都很大

我甚至已经用了最简单的几行语句：
'reinit'
'open fd01'
'set fwrite V.grd'
'set gxout fwrite'
'set t 1'
zz=1
while(zz<20)
'd V'
zz=zz+1
endwhile
'disable fwrite'
依然数据很大

river

Rplay 发表于 2013-8-6 23:47

                               
登录/注册后可看大图

已经找到问题所在，谢谢各位

不带你这样的，找到问题就感谢一下没事儿了。是什么问题，怎么解决的，说出来。也可以让别人看见，以后作为参考