【发布时间】：2016-03-11 05:39:11
【问题描述】：
我尝试运行一个 fortran 程序，以从 MODIS 8 每日叶面积指数 LAI 和光合有效辐射 FPAR 分数中提取数据。首先，我使用 *.sh 将 MODIS 原始 *.hdf 文件投影到 *.dat 文件。然后使用 g++ 将 *.dat 重新格式化为 *.txt，这个阶段是我数据准备的最后阶段，将所有 *.txt 数据组合成特定格式 *.direct 和 *.8daily 从头到尾全年通过使用 fortran 代码。
Program read_8daily
Implicit None
integer ir,ic
integer nnr, nnc
character filename1*80,filename2*80, string*20, year*4, day*3
character hour*2, outfile*200, ch2*2, ch5*5, chim*2
integer rec,count,i,j,sta,im, iday, ihour,iyear,d
integer, dimension (300,300):: fpar0, lai0
real, dimension (300,300):: fpar, lai
integer, dimension (300,300)::modis1, modis2
integer startday,endday,starthour,endhour
integer rec_no
nnc = 328 ! For Red River Basin
nnr = 395
do iyear = 2010, 2010 ! Change year
write(year, '(i4.4)') iyear
startday = 1 ! change start date (if)
endday = 361 ! change end date (if)
print *, 'EVALUATING from:', iyear, startday
print *, '.............To:', iyear, endday
print *, 'Reading modis data...'
OPEN (15, file='FPAR.'//year//'.8daily',
: status='unknown')
OPEN (16, file='LAI.'//year//'.8daily',
: status='unknown')
OPEN (25, file='FPAR.'//year//'.direct',
: form = 'unformatted', access = 'direct',
: recl = 4*nnc*nnr, status='unknown')
OPEN (26, file='LAI.'//year//'.direct',
: form = 'unformatted', access = 'direct',
: recl = 4*nnc*nnr, status='unknown')
print *,'Writing', '..modis.'//year//'.daily'
rec_no = 0 !record number for every monthly data file
do iday = startday, endday, 8 !8-day snow cover product
write(day, '(i3.3)') iday
filename1= year//day//'.fpar.txt'
c 文件名1='/home/maheswor/koshi/LAIFPAR/input/koshi_' c://year//'_'//day//'_fpar.txt'
filename2= year//day//'.lai.txt'
c filename2='/home/maheswor/koshi/LAIFPAR/input/koshi_' c ://year//'_'//day//'_lai.txt'
OPEN(7, file=filename1)
OPEN(8, file=filename2)
read(7,*) ((modis1(ir, ic),ic=1,nnc), ir=1,nnr)
read(8,*) ((modis2(ir, ic),ic=1,nnc), ir=1,nnr)
rec_no = rec_no + 1
print *, 'record number:', rec_no
do ir = 1, nnr
do ic=1, nnc
fpar0(ir,ic) = modis1(ir,ic)
lai0(ir,ic) = modis2(ir,ic)
end do
end do
print *, 'begin writing'
50 格式(1x, 300i10) 做 ir = 1 ,nnr 写(15,50) (fpar0(ir,ic),ic=1,nnc) 写(16,50) (lai0(ir,ic),ic=1,nnc) 结束做
write(25,rec = rec_no) ((fpar0(ir,ic),ic=1,nnc),ir=1,nnr)
write(26,rec = rec_no) ((lai0(ir,ic), ic=1,nnc),ir=1,nnr)
end do !do iday = startday, enday
close(7)
close(8)
close(25)
close(26)
close(15)
close(16)
end do !do iyear
end program
我在 cygwin 中运行了该脚本。我在照片中遇到了以下问题。请检查并帮助。 gfortran run and error
【问题讨论】：
-
提示：nnc 和 nnr 大于 300。
-
对所有 Fortran 问题使用 fortran 标签以获得更多关注。