| 
 
	积分1285贡献 精华在线时间 小时注册时间2016-4-14最后登录1970-1-1 
 | 
 
| 
本帖最后由 主啊 于 2018-4-8 14:36 编辑
x
登录后查看更多精彩内容~您需要 登录 才可以下载或查看,没有帐号?立即注册 
  
 我改了很多次,无奈,各位前辈可以帮我看看是什么问题么????
 原始程序:---------------------------------------------------------------------------
 
 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
 
 
 改了变量类型还是不行
 
 
 
 
 
 
 
 
 
 | 
 |