- 积分
- 13650
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-9-17
- 最后登录
- 1970-1-1
|
发表于 2015-10-8 23:22:30
|
显示全部楼层
楼主您好!谢谢您的分享。我一直都用C语言编程,fortran程序能看懂,grads能看懂脚本,但自己没写过。
刚刚试了一下你的程序,想按照你的数据跟程序走一遍,学习一下。
发现open(12, file='sta.grd',status='replace',form='binary')有问题,编译通过,但是运行时出现如下错误信息:
[zhufx@M001-HPCC station_plots]$ gfortran sta2grd.f90 -lm -o sta2grd
[zhufx@M001-HPCC station_plots]$ ./sta2grd
At line 13 of file sta2grd.f90
Fortran runtime error: Bad FORM parameter in OPEN statement
我换成open(12, file='sta.grd',form='unformatted')后,编译、运行均正常,也得到sta.grd文件,
但是下一步stnmap时出了问题,如下:
ga-> !stnmap -i /home/zhufx/gradsfolder/station_plots/sta.ctl
Name of binary data set: /home/zhufx/gradsfolder/station_plots/sta.grd
Number of times in the data set: 1
Number of surface variables: 1
Number of level dependent variables: 0
Starting scan of station data binary file.
Binary data file open: /home/zhufx/gradsfolder/station_plots/sta.grd
Processing time = 1
Invalid station hdr found in station binary file
Possible causes: Invalid level count in hdr
Descriptor file mismatch
File not station data
Invalid relative time
levs = 1123862446 flag = 0 time = 51.43
恳请您帮我解决我的问题,非常感谢!ctl文件跟您的例子基本一样,如下:
DSET /home/zhufx/gradsfolder/station_plots/sta.grd
DTYPE station
STNMAP /home/zhufx/gradsfolder/station_plots/sta.map
UNDEF -9999
TITLE Rain Data Sample
TDEF 1 linear 01may2011 1mo
VARS 1
rain 0 99 Rain Data
ENDVARS
谢谢您的帮助! |
|