gs中用fwrite写dat文件失败 出来的是类似二进制的文件

碧落

gs的目标是从grd文件中提取想要的数据，并生成dat文件，准备用fortran进一步计算，参照书本以及网上大家的脚本，尝试写了一个很简单的gs试一下效果：'reinit'
'open ./abc.ctl'

'set fwrite F:/work/dat/u3.dat'
'set gxout fwrite'
'set t 30'
'set lev 1000'
'd u3'
'disable fwrite'
'reinit'


运行中没有报错，然而生成的dat文件打开显示的是类似二进制文件那样的乱码：

生成的dat文件是这样的。。。

查了家园里其他人的类似脚本，貌似基本也是这个语法去写的gs，求助！希望写成正常十进制的dat。谢谢！
ps：数据文件是自己跑程序生成的，ctl文件也顺利画图无异样。

碧落

有一说一 我知道我的问题不入大神的眼 但是我也是认真看了书查了资料尝试无数次无果后过来提问 论坛的活跃我觉得毛毛雨小雨中雨的积极参与也很重要 感谢之前留言我的各位提供了各种思路 最后是在和我师兄的在线讨论中找到答案 感谢师兄提供参考文献和建议 解决方式不只有一条 我在这之后 也尝试了直接用fortran读 当然也成功了 小雨们加油 大胆提问 多多尝试 总有办法的！说了一堆废话 这小帖子里也未必有人看到 就当自勉吧

xfz147256

提取出二进制的，再放到Fortran转成十进制的应该就行了吧，grads可以直接写成十进制到文件吗？

霹雳电光鼠

这样提取出来的数据还是二进制的，直接拿fortran去读二进制数据就可以

river

写出来的本来就是二进制的，直接打开不乱吗才怪呢

伽蓝鸟

GrADS写得是single类型的，不是ASCII码。。

碧落

都说只提问不自己圆场的都是不负责的帖子，那我自问自答一下吧，希望帮得到其它童鞋。
上述大家都说的对的，我发问时的写法，写出来的依然是二进制。
然后，grads是可以直接写出二进制的

1. *******文件输出******
   
2. 'set t 'tii
   
3. file='./txt/'t2'_td3.txt'
   
4. 'set gxout print'
   
5. zi=zk
   
6. while ( zi <= zkk )
   
7. 'set z 'zi
   
8. 'set prnopts %g 402 1'
   
9. 'd td3'
   
10. rc=write(file,result)
    
11. zi=zi+1
    
12. endwhile
    
13. rc=close(file)

复制代码

set gxout时采用print即为输出ASCII码值，参考文献是：
https://max.book118.com/html/2017/0518/107570153.shtm
输出格式：
set prnopts format(%g) line_n(402) blank_n（1）'
其中%g为采用小数或标准指数格式中宽度较小的一种；line_n好像是列数，实操的时候确认一下输出即可，因为我需要输出数据矩阵，所以根据需要设置的。

碧落

上述代码是我gs中的一部分，直接截取出来的，大家主要看输出及其设置部分即可。

碧落

xfz147256 发表于 2021-4-19 14:15
提取出二进制的，再放到Fortran转成十进制的应该就行了吧，grads可以直接写成十进制到文件吗？

可以的 用set gxout print

碧落

霹雳电光鼠 发表于 2021-4-19 14:23
这样提取出来的数据还是二进制的，直接拿fortran去读二进制数据就可以

是的 后来经过反复确认 确实这样写出的依旧是二进制 改用print后实现ASCII输出了 后来也fortran编程实现了 谢谢！