立即注册 登录
气象家园 返回首页

Kimmy-24686的个人空间 http://bbs.06climate.com/?442 [收藏] [复制] [分享] [RSS]

日志

哭死,折腾了两个小时,这个弱智问题,求大家帮帮忙,看看哪里出问题了,

热度 1已有 170 次阅读2017-3-17 14:52 | 哭死, 折腾了两个小时, 这个弱智问题, 求大家帮帮忙, 看看哪里出问题了

fortran十进制直接成二进制,配相应的ctl,在grads中显示出来的值是错的,求求大侠帮帮忙吧,
parameter (N=46)
integer i
real yr,x(N),y(N),z(N)
open(10,file='D:\2017\ice_jx_year_jibing_1970-2015-dong.TXT')
    do i=1,N 
      read(10,*) yr,x(i),y(i) 
    enddo
   close(10)
   open(100,file='D:\2017\ice_jx_year_jibing_1970-2015-dong.grd',form='binary')
   do i=1,N
   write(100)x(i), y(i)
   ENDDO
 end
下面是相应的ctl,用grads是可以打开的,可是画出来的图却显示全是0,
dset d:\2017\ice_jx_year_jibing_1970-2015-dong.grd
undef -999.0
title SVD time coefficient of rainfall and NP MLD
xdef  1 linear  119.5    1.0
ydef   1 linear  19.739  0.333
zdef   1 levels  5    
tdef  46 linear   jan1969  1yr
vars 2
tcf 1 999 mean Jul rainfall (mm)
tcg 1 999 mean Jan mld
endvars
 

发表评论 评论 (1 个评论)

回复 Soaring 2017-3-18 08:31
read和write那边语句有问题,不能这么读的,还有你的yr是为了忽略连续3个数据中的第一个?

facelist doodle 涂鸦板

您需要登录后才可以评论 登录 | 立即注册

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

返回顶部