wrfout文件进行轨迹模拟提示出错

guo871015

本人是wrf的初学者，用wrfout文件进行轨迹模拟，总是出错，状况如图，请教高手如何解决。

运行命令如下：
cd /home/zhangguo/wrf3/traj/

ifort get-namelist.f90 -o get-namelist.exe

./get-namelist.exe

ifort main.f90 interp.f90 ll2xy.f90 xy2ll.f90 hgt2z.f90 tt2time.f90 -o trajectory.exe

./trajectory.exe
在运行最后的命令时出现如图提示的错误。

本人分别将为wrf模拟后的wrfout文件分别提取height、tk、U、V、W五个变量，进行轨迹编译回推轨迹，出现如图问题。
get-namelist.f90文件是编译的轨迹程序，程序如下：
PROGRAM NAMELIST
REAL LON,LAT,H(7)
INTEGER STT,ENT,DI
INTEGER FNUM
CHARACTER ST*2,Fname*19
INTEGER MDH(3),DTIME(2), DREC(3)
REAL LCC(8)


OPEN(20,ACTION='READ',FILE='Lonlat', FORM='FORMATTED')

FNUM=2
MDH=(/5,5,0/)
DTIME=(/73,49/)
DREC=(/130,130,32/)
LCC=(/27.000,115.000,65.500,65.500,60.000,30.000,115.000,20.000/)
Tint=3600
H=(/2000,1750,1500,1250,1000,750,500/)
DI=-1

DO
READ(20,*,END=333) ST,LON,LAT, DI

Fname="namelist.input_B_"//ST
IF (DI .EQ. 1)Fname="namelist.input_F_"//ST
OPEN(21,ACTION='WRITE',FILE=Fname,FORM='FORMATTED')
WRITE(21,101) ST
WRITE(21,102) FNUM,(MDH(i),i=1,3)
WRITE(21,103) (DTIME(i),i=1,2)
WRITE(21,104) (DREC(i),i=1,3)
WRITE(21,105) (LCC(i),i=1,8)
WRITE(21,106) Tint
101 FORMAT(A2)
102 FORMAT(4I6.2)
103 FORMAT(2I6.2)
104 FORMAT(3I6)
105 FORMAT(8F7.2)
106 FORMAT(F6.0)
j=1
DO WHILE(j<=7)
i=37
DO WHILE(i<=119)
  STT=i; ENT=i-36
IF (DI .EQ. -1) THEN
  WRITE(21,100) LON,LAT,H(j),STT,ENT,DI
ELSE
  WRITE(21,100) LON,LAT,H(j),ENT,STT,DI
ENDIF
100 FORMAT (F8.4,1X,F7.4,1X,F6.0,1X,3(I3,1X))
  i=i+1
ENDDO
j=j+1
ENDDO
CLOSE(21)
ENDDO
333 CONTINUE
CLOSE(20)
END
请教各位高手是编译程序的错误还是其他问题，再次拜谢各位高手大侠。