- 积分
- 2666
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-12-1
- 最后登录
- 1970-1-1
![未绑定新浪微博用户 新浪微博达人勋](source/plugin/sina_login/img/gray.png)
|
发表于 2020-4-28 11:52:45
|
显示全部楼层
回帖奖励 +1 金钱
本帖最后由 radon 于 2020-4-28 11:56 编辑
简单写一段fortran程序,功能是将txt转为二进制。配合.ctl描述文件,grads就可以读取
注意grads的纬度排列默认是从南向北,如果txt文件中不是,要调顺序,或在ctl文件中标明
integer lon,lat
parameter(lon=72,lat=18)
character(len=200) tempc
real dat(lon,lat)
integer i,j
open(110,file='SURF_CLI_CHN_TEM_DAY_GRID_0.5-MIN-20040704.txt',form='formatted',status='old',action='read')
do i=1,6
read(110,*) tempc
enddo
do j=1,lat
read(110,*) (dat(i,j),i=1,lon)
enddo
close(110)
open(210,file='TEM_20040704.dat',form='unformatted',access='direct',status='unknown',action='write',recl=lat*lon)
write(210,rec=1) ((dat(i,j),i=1,lon),j=1,lat)
close(210)
|
|