- 积分
- 322
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-1-25
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
最近在用ivf编译一个pom程序,用的网上提供的2k版本。输入文件fort.40已经在前期做好,带入主程序后报错,具体报错如下:
源代码如下:
subroutine file2ic
C **********************************************************************
C * *
C * FUNCTION : Sets up my own problem. *
C * *
C * This example read IC from IC.dat file, generated by GRID.f in *
C * GRID-DATA directory. Only minimal number of fields are read, *
C * while others are calculated here. *
C * *
C **********************************************************************
C
implicit none
C
include 'pom2k.c'
C
real rad,re,dlat,dlon,cff
integer i,j,k,m
character*5 field
rad=0.01745329
re=6371.E3
C
write(6,'(/,'' Read grid and initial conditions '',/)')
C
C--- 1D ---
! open(40,file='fort.40',status='old')
read(40,'(a5)') field
write(6,'(a5)') field
read(40,'(8E12.5)') z
read(40,'(a5)') field
write(6,'(a5)') field
read(40,'(8E12.5)') zz
read(40,'(a5)') field
write(6,'(a5)') field
read(40,'(8E12.5)') dz
read(40,'(a5)') field
write(6,'(a5)') field
read(40,'(8E12.5)') dzz
C--- 2D ---
read(40,'(a5)') field
write(6,'(a5)') field
read(40,'(8E12.5)') east_e
read(40,'(a5)') field
write(6,'(a5)') field
read(40,'(8E12.5)') north_e
read(40,'(a5)') field
write(6,'(a5)') field
read(40,'(8E12.5)') h
C--- 3D ---
read(40,'(a5)') field
write(6,'(a5)') field
read(40,'(8E12.5)') t
read(40,'(a5)') field
write(6,'(a5)') field
read(40,'(8E12.5)') s
read(40,'(a5)') field
write(6,'(a5)') field
read(40,'(8E12.5)') rmean
C--- Constant wind stress read here
C (for time dep. read in loop 9000 & interpolate in time)
read(40,'(a5)') field
write(6,'(a5)') field
read(40,'(8E12.5)') wusurf
read(40,'(a5)') field
write(6,'(a5)') field
read(40,'(8E12.5)') wvsurf
C
其中输入文件fort.40具体内容:
Z
0.00000E+00-0.35714E-01-0.71429E-01-0.14286E+00-0.21429E+00-0.28571E+00-0.35714E+00-0.42857E+00
-0.50000E+00-0.57143E+00-0.64286E+00-0.71429E+00-0.78571E+00-0.85714E+00-0.92857E+00-0.10000E+01
等等,就不全贴上了。
小弟已经改了两天了,毫无头绪,请各位论坛的高手们帮忙看看是哪里出了问题吧,感谢!
|
-
|