| 
 
	积分2666贡献 精华在线时间 小时注册时间2014-12-1最后登录1970-1-1 
 | 
 
 发表于 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)
 
 | 
 |