lyn 发表于 2004-6-17 08:30:57

【求助】e元件批处理命令

我做的非线性抛物方程,这几天的程序老是要吗没初值,要吗迭代一次就退出循环,今天看了看e元件批处理命令,
1。觉得不是书上教的那样?有点不理解?
tarta disp0 id0 nv disp0 bfd time0 time coor0 unod npa.io elem0 sys diag
if exist stop del stop
:1
bft disp0 time coor0 nv bfdt stop
if exist end del end
:2
enpa disp0 sys time nv coor0 bfd diag elem0 f displace end
if not exist end goto 2
call post.bat
if not exist stop goto 1
2。我的目录下没有生成end文件,为什么迭代后就退出,开始执行call post.bat?
3。我读了读enpa.for,发现程序中并没有用到enpa批处理命令最后两个文件displace ,end。为什么?
附件是我的enpa.for
请斑竹指点!!学生我非常感激!!

lyn 发表于 2004-6-17 11:52:15

回复: 【求助】e元件批处理命令

疑点可能找到,
1。
if (err .lt. 1.0e-8 .or. err .lt. 1.0e-8*ul) then
      open(11,file='end',status='unknown')
      close(11)
我觉得能不能写成
if (err .lt. 1.0e-8)then。。。。
2。
每次生成的end文件在下一次时步开始就被delete
if exist end del exist
所以我找不到这个文件

请教斑竹上面对不对?谢谢

lyn 发表于 2004-6-17 12:14:58

回复: 【求助】e元件批处理命令

但是即使我修改了程序,执行时还是按照原来的,怎么样才能执行修改后的程序?谢谢

lyn 发表于 2004-6-17 12:18:40

回复: 【求助】e元件批处理命令

我重新按照
fl32 /c enpa.for
编译链接时,出现许多warning,主要是参数类型不匹配,
是不是这个引起的?

ubah 发表于 2004-6-17 12:36:07

回复: 【求助】e元件批处理命令


if (err .lt. 1.0e-8 .or. err .lt. 1.0e-8*ul) then
      open(11,file='end',status='unknown')
      close(11)
是收敛准则,可以修改

算完后应该有end文件
页: [1]
查看完整版本: 【求助】e元件批处理命令