爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 4930|回复: 0

[脚本编辑] 紧急求助,grads生成的数据用工作站上的fortran读取显示一堆-1.0318071E-22

[复制链接]
回帖奖励 3 金钱 回复本帖可获得 3 金钱奖励! 每人限 1 次
发表于 2016-8-14 17:14:38 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
各位老师,我在工作上用grads生成的数据,用fortran读取时,显示一堆-1.0318071E-22(实际上正确值是-999.9),在open语句中加了little_endian也没用,有人知道是怎么回事吗?着急中,谢谢!
下面是grads生成的数据的ctl描述文件
dset ^obs_daily.bin
undef -999.9
title forcing
xdef 126 linear 72.75 0.5
ydef 74 linear 17.75 0.5
tdef 185  linear 01Mar1982 1dy
edef 45 names 1962 1963 1964 1965 1966 1967 1968 1969 1970 1971 1972 1973 1974 1975 1976 1977 1978 1979 1980 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006
zdef 1 linear 1 1
vars 1
prate  0 99 precipitation rate [mm/day]
endvars

下面是读取的fortran程序
implicit none
integer,parameter :: nx=126,ny=74
integer :: irec
integer :: iday,iyr,iyrr
integer :: imon
integer :: ix,iy,i,j
real :: obs(nx,ny,185,45)

write(*,*) 'read_obs and prior'
open(unit=1,file='obs_daily.bin', &
    form='unformatted',access='direct', &
  recl=nx*ny,status='old',convert='little_endian')
irec=1
iyrr=1
do iyr=1,45
   if(iyr.eq.(1982-1961))then
          irec=irec+185
   else
       do iday=1,185
           read(1,rec=irec)obs(:,:,iday,iyrr)
           print*,obs(:,:,iday,iyrr)
           pause
           irec=irec+1
       enddo
       iyrr=iyrr+1
   endif
enddo
close(1)
end
相同的fortran程序在微机上能正确读取,但是我需要在工作站上运行程序,求教各位知道是什么原因吗?
先谢谢了!

密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

快速回复 返回顶部 返回列表