Esope source files

   1 : C KALTET    SOURCE    CHAT      05/01/13    00:51:01     5004           
   2 :       SUBROUTINE KALTET(TPLUS,PRANDT)                                           
   3 : C                                                                               
   4 : C     PROGRAMMA PER LA TABULAZIONE DI TPLUS                                     
   5 : C     PER UN ASSEGNATO VALORE DI PRANDTL                                        
   6 : C                                                                               
   7 :       IMPLICIT INTEGER(I-N)
   8 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
   9 :       DIMENSION TPLUS(101)                                                      
  10 :       APLUS=26.0                                                                
  11 :       BPLUS=34.96+28.79*(LOG10(PRANDT))+33.95*(LOG10(PRANDT))**2              
  12 :      *+6.33*(LOG10(PRANDT))**3-1.186*(LOG10(PRANDT))**4                       
  13 :       TPLUS(1)=0.0D0                                                            
  14 :       YPLUS=0.0D0                                                               
  15 :       ZL=ZLM(YPLUS,APLUS)                                                       
  16 :       UP=UPRIMO(ZL)                                                             
  17 :       ST=SIGMAT(PRANDT,YPLUS,APLUS,BPLUS)                                       
  18 :       ZETAV=PRANDT/(1.0+ZL*UP*PRANDT/ST)                                        
  19 :       DO 10 I = 2,101                                                           
  20 :       YPLUS=YPLUS+1.0D0                                                         
  21 :       ZL=ZLM(YPLUS,APLUS)                                                       
  22 :       UP=UPRIMO(ZL)                                                             
  23 :       ST=SIGMAT(PRANDT,YPLUS,APLUS,BPLUS)                                       
  24 :       ZETAN=PRANDT/(1.0+ZL*UP*PRANDT/ST)                                        
  25 :       TPLUS(I)=TPLUS(I-1)+(ZETAV+ZETAN)/2.0                                     
  26 :       ZETAV=ZETAN                                                               
  27 :    10 CONTINUE                                                                  
  28 :       RETURN                                                                    
  29 :       END
  30 :                                                                        
  31 :                                                                                 
  32 :