Esope source files

   1 : C YKFPT     SOURCE    MAGN      06/03/01    21:15:10     5317           
   2 :       SUBROUTINE YKFPT(XRO,IK1,VMU,IK2,CP,IK3,XLB,IK4,UET,YP,H,NEL)
   3 : C
   4 : C
   5 : C
   6 : C
   7 : C
   8 :       IMPLICIT INTEGER(I-N)
   9 :       IMPLICIT REAL*8 (A-H,O-Z)
  10 :       DIMENSION H(*),UET(*),XRO(*),VMU(*),CP(*),XLB(*)
  11 :       DIMENSION TPLUS(101)
  12 : C
  13 :       DO 10 K=1,NEL
  14 :       PRANDT = VMU(k*(1-IK2)+IK2)*CP(k*(1-IK3)+IK3)/XLB(k*(1-IK4)+IK4)
  14 :       CALL KALTET (TPLUS,PRANDT)
  15 : C
  16 :       YPLUS = XRO(k*(1-IK1)+IK1)*YP*UET(K)/VMU(k*(1-IK2)+IK2)
  17 :       YPLUS = ABS(YPLUS) + 1.E-5
  18 :       IF (YPLUS.GT.100.0) THEN
  19 :         ATHE = (TPLUS(101)-TPLUS(91))/(LOG(1.D2)-LOG(9.D1))
  20 :         BTHE = TPLUS(101)-(ATHE*LOG(1.D2))
  21 :         TPIU = ATHE*(LOG(YPLUS))+BTHE
  22 :       ELSE
  23 :       IY = INT(YPLUS)
  24 :       RY = YPLUS-IY
  25 :       TPIU = TPLUS(IY+1)+RY*(TPLUS(IY+2)-TPLUS(IY+1))
  26 :       ENDIF
  27 :       H(K)=UET(K)/ABS(TPIU + 1.E-5)*XRO(k*(1-IK1)+IK1)*CP(k*(1-IK3)+IK3)
  28 :  10   CONTINUE
  29 :       RETURN
  30 :       END
  31 : 
  32 : 
  33 : 
  34 : 
  35 : 
  36 :  
  37 :  
  38 :