moitoi版主可否将你的ground贴出来一个学学

advard

包括，q1文件和ground文件。我用的不成功。谢谢！

moitoi

抱歉,我作的ground比较复杂不易贴出来.等我有时间，我给举几个examples,好吗？

betterman

水水水水水水水水水水水水水水水水水水水水水水水灌灌灌灌灌灌灌灌灌水水水水
水水水水水水水水水水水水水水水水水水灌灌灌灌灌灌灌灌灌灌灌灌灌灌灌水水水水
水水水水水水水水水水水水水水水灌灌灌灌灌灌灌灌灌灌灌灌灌灌灌灌灌灌水水水水
水水水水水水水水水水水灌灌灌灌水灌灌灌灌灌灌灌灌灌灌灌水水水水水水水水水水
水水水水水水水水灌灌灌灌灌灌灌水灌灌灌水水灌灌灌灌灌水水水水水水水水水水水
水水水灌灌灌灌灌灌灌灌灌灌灌灌水水水水水水灌灌灌灌水水水水水水水水水水水水
水灌灌灌灌灌灌灌灌灌灌灌灌灌灌水水水水水水灌灌灌灌水水水水水水水水水水水水
水灌灌灌灌灌灌灌灌灌灌灌灌水水水水水水水灌灌灌灌灌灌灌灌灌灌灌水水水水水水
水灌灌灌灌灌灌灌灌灌灌灌灌水水水水水水灌灌灌灌灌灌灌灌灌灌灌灌灌灌水水水水
水水灌灌灌灌灌灌灌灌灌灌水水水水水灌灌灌灌灌灌水水水灌灌灌灌灌灌灌水水水水
水水水水水水水水灌灌灌灌水水水水水灌灌灌灌水水水水水水灌灌灌灌灌水水水水水
水水水水水水水水灌灌灌灌水水水水灌灌灌灌水水灌灌水水水灌灌灌灌灌水水水水水
水水水水水水水水灌灌灌灌水水水水灌灌灌灌水水灌灌灌灌水灌灌灌灌灌水水水水水
水水水水水水水水灌灌灌灌水水水水灌灌灌灌水水灌灌灌灌水灌灌灌灌灌水水水水水
水水水水水水水水灌灌灌灌水水水水灌灌灌灌水水灌灌灌水水灌灌灌灌灌水水水水水
水水水水水水水水灌灌灌灌水水水水灌灌灌灌水水灌灌灌水水灌灌灌灌灌水水水水水
水水水水水水水水灌灌灌灌水水水水灌灌灌灌水灌灌灌灌水水灌灌灌灌灌水水水水水
水水水水水水水水灌灌灌灌水水水水灌灌灌灌水灌灌灌灌水水灌灌灌灌灌水水水水水
水水水水水水水水灌灌灌灌水水水水灌灌灌灌水灌灌灌灌水水灌灌灌灌灌水水水水水
水水水水水水水水灌灌灌灌水水水水灌灌灌灌水灌灌灌灌水水灌灌灌灌灌水水水水水
水水水水水水水水灌灌灌灌水水水水灌灌灌水水灌灌灌灌水水灌灌灌灌灌水水水水水
水水灌灌水水水灌灌灌灌灌水水水水灌灌灌水水灌灌灌水水水灌灌灌灌灌水水水水水
水水灌灌灌灌灌灌灌灌灌灌水水水水水灌灌水水灌灌水水水水灌灌灌灌灌水水水水水
水水水灌灌灌灌灌灌灌灌灌水水水水水水水水灌灌灌水水水水水灌灌灌灌水水水水水
水水水水水灌灌灌灌灌灌灌水水水水水水水水灌灌灌水灌灌灌灌水水水水水水水水水
水水水水水水灌灌灌灌灌灌水水水水水水水灌灌灌灌水水灌灌灌灌灌水水水水水水水
水水水水水水水水水灌灌灌水水水水水水灌灌灌灌灌水水水灌灌灌灌灌灌灌水水水水
水水水水水水水水水水水水水水水水灌灌灌灌灌灌水水水水水灌灌灌灌灌灌水水水水
水水水水水水水水水水水水水水水灌灌灌灌灌灌水水水水水水灌灌灌灌灌灌灌水水水
水水水水水水水水水水水水水水灌灌灌灌灌水水水水水水水水水灌灌灌灌灌灌水水水
水水水水水水水水水水水水水灌灌灌灌灌水水水水水水水水水水水灌灌灌灌水水水水
水水水水水水水水水水水水灌灌灌水水水水水水水水水水水水水水水灌灌灌水水水水
水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水

advard

谢谢板主,举个简单的例子就可以,我用ground计算,运行一次earth,则ground.for文件
变成原来的样子.不知道为什么?

moitoi

在Q1文件中的第十三组中，用户可用PATCH和COVAL PHOENICS函数命令提供要加的非线性源项的区域和类型及边界条件，并指令EARTH访问GROUND中的第十三组，用户可在此组中加入计算中待求的非线性系数和值的程序。例如在X方向动量方程中有下例形式的源项：
      S=-au(-1)[0.0-u]
注：(-1)表示指数项（考过来的公式显示不出来，只能如此表示了）
式中系数为非线性，值为零。在Q1文件第十三组中输入
PATCH（NONL，PHASEM，1，NX，1，NY，1，NZ，1，1）
COVAL（NONL，U1，GRND，0.0）
PATCH中第一项为源项名，第二项为源类型，后面各项为坐标、时间区域。COVAL中第一项为源项名，第二项为变量名，第三项GRND是PHOENICS 指令，表示系数不是常数，通知EARTH 访问GROUND第十三组中的第一项计算系数项；第四项是值系数项。
在Q1文件中设置USEGRD=T
在GROUND第十三组中的第一项中用户要加入以下程序段
    C --------GROUP 13. Source and Boundary condition
    13  GO TO (131,132,133,134,135,136,137,138,139,1310),ISC
    131 CONTINUE
    C-------------------SECTION 1   Grnd -----
IF (NPATCH(1:4).EQ.’NONL’) THEN
  IF(INDVAR.EQ.U1) THEN
  L0U=L0F(U1)
  L0CO=L0F(CO)
  NXNY=NX*NY
  GB=R(2)-1
  DO 1300 I=1,NXNY
  XX=F(L0U+I)**GB
  1300  F(L0CO+I)=R(1)*XX
  ENDIF
ENDIF
RETURN
132CONTINUE
而后编译ground , 和其他项连接，生成earth.exe, 在RUN选项中选用private.运行生成的earth.exe.

advard

板主您说的是plant,而我的意思是直接到ground中编程.不用plant自动生成.因为有时候plant做不到.

moitoi

是ground!!!!

moitoi

plant 仅在q1文件写命领语句,不需要在ground中写程序.直接到ground中编程.需在用ground,q1文件中都要编程!!!!

mix

是否使用了plant，重新编译，会冲掉旧版本的ground文件？

moitoi

开始GROUND是空的只要你的q1不变就行了

mix

这么说来，只要通过plant修改q1文件
ground就会重新被改写
所以大家要做好备份

下面引用由moitoi在 2003/03/07 04:50pm 发表的内容：
开始GROUND是空的只要你的q1不变就行了

流体

下面引用由moitoi在 2003/02/24 06:30pm 发表的内容：
在Q1文件中的第十三组中，用户可用PATCH和COVAL PHOENICS函数命令提供要加的非线性源项的区域和类型及边界条件，并指令EARTH访问GROUND中的第十三组，用户可在此组中加入计算中待求的非线性系数和值的程序。例如 ...

斑竹，可否将你的这几行程序稍做一下注释说明！！（要求有点过分了。）
L0U=L0F(U1)
L0CO=L0F(CO)
。。。。。。

phoenics

不好意思在斑竹面前献丑了：
L0U=L0F(U1)
就是将U1在当前slab中存储的起始位置读出来。就象计算机原理的中的堆栈指针。

XX=F(L0U+I)**GB  LOU是个地址+I，即现在单元的地址，F取值即将现在单元的值取出来
1300  F(L0CO+I)=R(1)*XX 将LOCO+I位置的值改为R（1）*XX

总之：LOF取地址，F付值。
斑竹我的理解是否正确，望指教。

merryfish

^_^，真有耐心，看来是真该你才能做这个班长了。我是不行了。。。。

李莉你好

^_^