tut2.4运行问题

meconglin

　　在运行tut2.4的例子时，系统提示子程序droico.f不能通过。子程序严格按照例子中程序运行，但系统提示的语法错误不知如何修改，在这里要请教大家，下面是系统提示的截图
  
下面是droico.f的程序内容C***********************************************************************
      SUBROUTINE DROICO(IDRT,RNDP,XD,YD,ZD,UD,VD,WD,DRD,DENDR,TD,
     &                  ICSRGP,OMEGRP)
C     Droplet initial conditions
C***********************************************************************
C--------------------------------------------------------------------------*
C     STAR VERSION 3.20.000                                                *
C--------------------------------------------------------------------------*
      INCLUDE ';comdb.inc';

      GLOBAL/USR001/INTFLG(100)

      INCLUDE ';usrdat.inc';
      EQUIVALENCE( UDAT05(001), IDR )
C-----------------------------------------------------------------------
C
C    This subroutine enables the user to specify droplet initial
C    conditions.
C
C    ** Parameters to be returned to STAR: IDRT, RNDP, XD,YD,ZD,
C                                          UD,VD,WD, DRD, DENDR, TD
C                                          ICSRGP,OMEGRP
C
C    Droplet type IDRT is passed to this subroutine as a negative
C    number. If its value is not changed to valid value
C    (between 1 and 10) this subroutine will not be called again in
C    this time-step and injection of particles will be
C    held until the next time-step.
C
C-----------------------------------------------------------------------
C-----------------------------------------------------------------------
     INTEGER*4 NN,ISEED,TWO14
     SAVE NN
     data NN/0/
C-----------------------------------------------------------------------
     IF(.NOT.(TIME.LT.(0.02-0.5*DT).OR.(TIME.GT.(0.03+0.5*DT).AND.TIME.LT.(0.05-0.5*DT)))) RETURN
C----NUMBER OF PARCELS INJECTED PER TIME STEP
     MPAPTS=10
     IF(IDN.GT.NPAPTS) RETURN
     PI=3.1415927
     DEGRAD=PI/180.
C----DROPLET MASS INJECTION RATE (kg/s)
     FLOWDR=1.0E-04
C----DROPLET TYPE
     IDRT=1
C----DROPLET DIAMETER (m)
     DRD=2.5E-05
C----DROPLET DENSITY (kg/m3)
     DENDR=995.
C----DROPLET TEMPERATURE (K)
     TD=300.
C----MASS OF SINGLE DROPLET
     DRMAS=DENDR*PI/6.*DRD**3
C----NUMBER OF DROPLETS REPRESENTED BY DROPLET
     RNDP=FLOWDR/(NPAPTS*DRMAS)*DT
C----INITIAL POSITION OF DROPLET (m)
     XD=2.9E-02
     YD=1.0E-07
     ZD=1.4E-02
C----GET RANDOM NUMBER RN1 AND RN2
     ISEED=150623
     TWO14=16384
     IF(NN.LT.1500.AND.NN.GY.-1500) NN=ISEED
     NN=NN*ISEED
     RN1=MOD(NN,16384)/16384
     NN=NN*ISEED
     RN2=MOD(NN,16384)/16384
C----INITIAL VELOCITY OF DROPLET (m/s)
     CONANG=30.
     FI=(90.+CONANG*RN1)*DEGRAD
     TH=180.*CONANG*RN2*DEGRAD
     UMAG=10.0
     UD= UMAG*COS(FI)*SIN(TH)
     VD= UMAG*SIN(FI)
     WD= UMAG*COS(FI)*COS(TH)
C-------------------------------------------------------------------------C

      RETURN
      END
C

supergsx

這個截圖不太能看出error在哪 ?
只知道你的 droico.f set 有問題
你必須去 check error message 才能知道 或是 修改/check droico.f 檔案內容

meconglin

谢谢你的指点，的确是子程序出了问题，我在补充中将程序粘上了，还肯请帮忙修改一下。[br][br][以下内容由 meconglin 在 2007年04月10日 09:25am 时添加] [br]
我用的是star_cd的3.2版本

changhai

子程序是版本自带的，
C:\Program Files\CD adapco\STAR-CD\PROSTAR\3.22\tutorials\star\tut2\tut2.4
能找到的