大家帮我看看这段subroutine？

cnzjuyilin

一段管道，内部布满细小的铜丝，重点模拟其中温度场的变化情况。
采取的方法是：将整个流动区域设为porous，以此来模拟流场内的压力损失。然后在同样的空间位置再建立一套相同大小的网格，设为固体，两套网格之间会发生热量传递。能量方程源项的子程序如下：
      IF(INTFLG(61).NE.1) THEN
        OPEN(87,FILE=';CHKSOLI.txt';,STATUS=';UNKNOWN';)
        OPEN(85,FILE=';CHKFLUI.txt';,STATUS=';UNKNOWN';)
        INTFLG(61)=1
      ENDIF

      IF(INTFLG(10).NE.ITER) THEN
        OPEN(UNIT=86,FILE=';T.rest';,STATUS=';UNKNOWN';)
  10    READ(86,*,END=20) I,v1,K
        TMON(I) = v1
        GOTO 10
  20    CONTINUE
        CLOSE(86)
        INTFLG(10)=ITER
      ENDIF
      ISOL=6
      IPORO=1
      IOFSET=5830
      H=200.
      AV=100.0/5099e-6
C--------SOLID-------------------C
      IF(ICTID.EQ.ISOL) THEN
       ISTG=IP-IOFSET
       TG=TMON(ISTG)+273.
       S1P=H*AV*TG
       S2P=H*AV
C-- CHECK WRITE
       IF(INTFLG(51).NE.ITER) THEN
         WRITE(6,*) ';SOLID DATA ';,IP,ISTG,T,TG,S1P,S2P
         INTFLG(51)=ITER
       ENDIF      
      ENDIF
C--------FLUID-------------------C
      IF(ICTID.EQ.IPORO) THEN
       IPSTS=IP+IOFSET
       TS=TMON(IPSTS)+273.
       S1P=H*AV*TS
       S2P=H*AV
C-- CHECK WRITE
       IF(INTFLG(52).NE.ITER) THEN
         WRITE(6,*) ';FLUID DATA ';,IP,IPSTS,T,TS,S1P,S2P
         INTFLG(52)=ITER
       ENDIF
      ENDIF      
      RETURN
      END
C
对starcd里的子程序不是很了解，哪位详细解释一下？
谢谢！