- 积分
- 1285
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-4-14
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 主啊 于 2018-4-8 14:36 编辑
我改了很多次,无奈,各位前辈可以帮我看看是什么问题么????
原始程序:---------------------------------------------------------------------------
program binary
implicit none
integer,parameter::nnt=31,stii=2425
integer:: ii=1,jj=1,i,j
integer nlev,flag
real tim
character*7 sta1(stii)
real lon1(stii),lat1(stii),rain1(nnt,stii)
!===========================================
character*7 name !name字符型
real lat,lon,year,mon,day,water
!===========================================
Open(1,file='real16061108.txt',status='old') !status = old防止覆盖,有格式顺序
open(2,file='sta.grd',status='replace',form='binary') ! 二进制型
!=============================================================
do while(.not.eof(1))
if(ii.NE.32)then
READ(1,*) name,lon,lat,year,mon,day,water
rain1(ii,jj)=water
ii=ii+1
end if
if(ii.EQ.32)then
sta1(jj)=name
lon1(jj)=lon
lat1(jj)=lat
jj=jj+1
ii=1
end if
end do
!=================================================================
do i=1,nnt
do j=1,stii
tim=0.0
nlev=1
flag=1 !无高空数据tim=0 , nlev = 0,有地面数据所以flag = 1
Write(2) sta1(j),lon1(j),lat1(j),tim,nlev,flag,rain1(i,j)
print *,sta1(j),lon1(j),lat1(j),tim,nlev,flag,rain1(i,j)
pause
end do
nlev=0
Write(2) sta1(j-1),lon1(j-1),lat1(j-1),tim,nlev,flag
end do
close(1)
close(2)
end
CTL文件:------------------------------------------------------------------------------------
DSET F:/3-mon/sta.grd
DTYPE station
STNMAP F:/3-mon/sta.map
UNDEF -32766
TITLE Rain Data Sample
*这里的时间可以根据上面示例文件中的第一行的时间填写
TDEF 31 linear 01JUL2016 1dy
VARS 1
rain1 0 99 Rain Data
ENDVARS
改了变量类型还是不行
|
|