polis中的范例不好使

advard

polis中的范例,C:\phoenics\D_POLIS\D_TUTS\GRID\GRTUT3.HTM的case5不能正确生成bfc网格,q1文件如下:请指教
************************************************************
   Q1 created by VDI menu, Version 3.3, Date 14/12/00
CPVNAM=VDI;SPPNAM=Core
************************************************************
IRUNN   =       1 ;LIBREF =       0
************************************************************
  Group 1. Run Title
TEXT(No title has been set for this run.     )
************************************************************
  Group 2. Transience
STEADY  =    T
************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
    * Set overall domain extent:
    *        xulast  yvlast  zwlast
            name
************************************************************
  Group 6. Body-Fitted coordinates
TALK=T;RUN(1,1);VDU=WINDOWED
BFC=T
GSET(P,A,0,-.1,0)
GSET(P,B,0,0,0)
GSET(P,C,0,0,0.25)
GSET(P,D,0,-.07,0.25)
GSET(P,E,0,-.1,0.3)
GSET(P,F,0,-.1,0.4)
GSET(P,G,0,-.15,0.4)
GSET(P,H,0,-.15,0.2)
GSET(L,L1,A,B,10)
GSET(L,L2,B,C,9)
GSET(L,L3,C,D,4)
GSET(L,L4,D,E,2)
GSET(L,L5,E,F,4)
GSET(L,L6,F,G,6)
GSET(L,L7,G,H,6,1.0,ARC,0,-.2,.3)
GSET(L,L8,D,H,6)
GSET(V,V1,S,A,SPLINE)
GSET(V,0,-.09,.08)
GSET(V,0,-.08,.14)
GSET(V,V1,E,H)
GSET(L,L9,A,H,9,1.0CRVV1)
GSET(F,F1,A,-,B,-,C,D,H,-)
GSET(F,F2,H,-,D,E,F,-,G,-)
GSET(D,20,20,15,.1,.2,.4)
GSET(M,F1,-J+K,1,11,1,TRANS)
GSET(M,F2,-J+K,1,11,10,TRANS)
GSET(I)
GSET(P,A,-0.707*0.05,0.6,0.35-0.707*0.05)
GSET(P,B, 0.707*0.05,0.6,0.35-0.707*0.05)
GSET(P,C, 0.707*0.05,0.6,0.35+0.707*0.05)
GSET(P,D,-0.707*0.05,0.6,0.35+0.707*0.05)
GSET(L,L1,A,B,6,1.0,ARC,0,0.6,0.35-0.05)
GSET(L,L2,B,C,6,1.0,ARC,0.05,0.6,0.35)
GSET(L,L3,C,D,6,1.0,ARC,0,0.6,0.35+0.05)
GSET(L,L4,D,A,6,1.0,ARC,-0.05,0.6,0.35)
GSET(F,F1,A,-,B,-,C,-,D,-)
GSET(M,F1,+I+K,8,21,10)
GSET(C,J14,F,J21,8,13,10,15,+,0.0,-0.3,0.0,INC,1.0)
GSET(C,I21,F,I1,RZ,-3.14159*2,0,0,INC,1)
GSET(T,J14,F,J11,8,13,10,15,1.0)
GSET(B,7,12,5,13,10,15,TRANS)
GSET(I)

    * Match a grid mesh
    **********
NONORT  =    T
************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS  =    T
    * Non-default variable names
NAME(149) =DEN1 ; NAME(150) =PRPS
STORE(PRPS,DEN1)

************************************************************
  Group 8. Terms & Devices
************************************************************
  Group 9. Properties
SETPRPS(1,  0)
RHO1    = 1.189000E+00
CP1     = 1.005000E+03
ENUL    = 1.544000E-05 ;ENUT   = 0.000000E+00
DVO1DT  = 3.410000E-03
************************************************************
  Group 10.Inter-Phase Transfer Processes
************************************************************
  Group 11.Initialise Var/Porosity Fields
FIINIT(DEN1) =  1.001000E-10 ;FIINIT(PRPS) =  0.000000E+00
   No PATCHes used for this Group


INIADD  =    T
************************************************************
  Group 12. Convection and diffusion adjustments
   No PATCHes used for this Group
************************************************************
  Group 13. Boundary & Special Sources
   No PATCHes used for this Group

EGWF    =    T
************************************************************
  Group 14. Downstream Pressure For PARAB
************************************************************
  Group 15. Terminate Sweeps
LSWEEP  =       1
RESFAC  = 1.000000E-03
************************************************************
  Group 16. Terminate Iterations
************************************************************
  Group 17. Relaxation
************************************************************
  Group 18. Limits
************************************************************
  Group 19. EARTH Calls To GROUND Station
USEGRD  =    T  ;USEGRX =    T
RADIA   = 0.000000E+00 ;RADIB  = 0.000000E+00
************************************************************
  Group 20. Preliminary Printout
ECHO    =    T
************************************************************
  Group 21. Print-out of Variables
************************************************************
  Group 22. Monitor Print-Out
NPRMON  =  100000
NPRMNT  =       1
TSTSWP  =      -1
************************************************************
  Group 23.Field Print-Out & Plot Control
NPRINT  =  100000
ISWPRF  =       1 ;ISWPRL =  100000
   No PATCHes used for this Group
************************************************************
  Group 24. Dumps For Restarts
NOWIPE  =    T

GVIEW(P,7.770572E-01,4.280243E-01,4.614949E-01)
GVIEW(UP,-3.143628E-01,8.991129E-01,-3.045850E-01)

> DOM,    SIZE,        1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    MONIT,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> DOM,    RELAX,       5.000000E-01
************************************************************
  

mix

有的范例只能在dos模式下用

advard

谢谢!非常感谢!

horse

polis中的例子基本上都可以在sattelite界面下运行，你可以试试看。

advard

faint, 不试我怎么知道不行.