|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
想把KIVA3VR2移植到Window 的Visual Fortran 6.6下,有几个问题向想各位高人请教:
1)是否可把fdata.f中的dateh函数和timeh函数直接改为date函数和time函数?
2)是否可以直接屏蔽K3prep.f中gmova2、glina2、gplot、gwind2、gview2、gdone、gpage、gtext等unix作图函数?
3)K3PREP运算时"WARNING: cell 842 in region 1 is not convex:"是什么意思?我是用其提供的例子算的。convex cell check有关公式在哪可以查到?
v4138=(x8-x4)*((y1-y4)*(z3-z4)-(z1-z4)*(y3-y4))
& +(y8-y4)*((z1-z4)*(x3-x4)-(x1-x4)*(z3-z4))
& +(z8-z4)*((x1-x4)*(y3-y4)-(y1-y4)*(x3-x4))
4)程序中tinop变量有无用处,可否删掉?子程序call setctlc(tinop,msg,80)是否可以删掉?
5)向下列情况,变量赋值语句是否要改:
定义: common /kiva1/ aaa1(1),x(nv),y(nv),z(nv),u(nv),v(nv),w(nv),
1 ro(nv),vol(nv),p(nv),amu(nv),f(0:nv),fv(0:nv),temp(nv),sie(nv),
2 bcl(nv),bcf(nv),bcb(nv),alx(nv),aly(nv),alz(nv),afx(nv),
3 afy(nv),afz(nv),abx(nv),aby(nv),abz(nv),tke(nv),pit(nv),pit1(nv),
4 suvw(nv),cli(nv),clj(nv),clk(nv),cfi(nv),cfj(nv),cfk(nv),
5 cbi(nv),cbj(nv),cbk(nv),fsum14(nv),fsum34(nv),fsum84(nv),
6 fsv(nv),fvsv(nv),spd(nv,lnsp),drn(lsoot),dnuc(lsoot),con0(lsoot),
7 e1(nv),e2(nv),e3(nv),e4(nv),e5(nv),e6(nv),e7(nv),e8(nv),e9(nv),
8 e10(nv),e11(nv),e12(nv),e13(nv),e14(nv),e15(nv),e16(nv),e17(nv),
9 e18(nv),e19(nv),e20(nv),e21(nv),e22(nv),e23(nv),e24(nv),e25(nv),
x e26(nv),e27(nv),e28(nv),e29(nv),e30(nv),e31(nv),e32(nv),e33(nv),
1 e34(nv),e35(nv),e36(nv),e37(nv),e38(nv),e39(nv),e40(nv),e41(nv),
2 e42(nv),e43(nv),e44(nv),e45(nv),e46(nv),e47(nv),e48(nv),e49(nv),
3 e50(nv),e51(nv),e52(nv),e53(nv),e54(nv),e55(nv),e56(nv),e57(nv),
4 e58(nv),e59(nv),e60(nv),e61(nv),zzz1
parameter (ncom1 = 44*nv + lnsp*nv + 3*lsoot + 61*nv + 2)
赋值: do 10 n=1,ncom1
10 aaa1(n)=0.0
程序似乎通过一句语句给所有变量赋了值,在Visual Fortran 6.6中似乎要改?
6)程序中经常出现数组下界出界的情况,也既定义了x(nv),程运行时会出现x(0)的情况,而且不少变量都有这种情况出现,把数组定义改为x(0:nv),可以解决数组下界出界的情况,但会不会由此引起新的问题?
7)KIVA3VR2提供的例子在Red Hat Linux和红旗Linux下应用是否没有问题?Linux下什么Fortran编译器好使?哪有下载?
8)其他我还没有发现的问题,望各位高手指点!!!
不胜感谢!!!Thanks very much in deed!!! |
|