ALCALC.f90#40 ALPHPD = (1.0+R*R)/(1.0+3.0*R*R) BBCALC.f90#43 BSUB1 = R*ASUB1 #44 BSUB2 = -R*ASUB2 BECALC.f90#40 BETAPD = 1.0 + R*R CAVITY_HT.f90#2 SUBROUTINE CAVITY_HT(H,W,L,R,LC,LR,HR,X,ZCAV,YCAV) #36 REAL H , W , L , R , HR , LR , LC #52 IF ( X.GE.0. .AND. X.LE.R ) THEN ! 656374 #53 YCAV = (W/2.+R/3.) - (X-R)**2/(3.*R) ! 151928 #54 ELSEIF ( X.GT.R .AND. X.LE.(L+LR) ) THEN #55 YCAV = (W/2.+R/3.)*(1.-((X-R)/(L+LR-R))**2)**0.5 ! 504446 #69 IF ( X.GE.0. .AND. X.LE.0.5*R ) THEN ! 9660 #70 ZCAV = HR + 4.*(X-0.5*R)**2*(H-HR)/(R**2) ! 2564 #71 ELSEIF ( X.GT.0.5*R .AND. X.LE.(L+LR) ) THEN #72 ZCAV = HR*(1.-((X-0.5*R)/(L+LR-0.5*R))**2)**0.5 ! 7096 CRCALC.f90#41 R = 2.0 GCUBIC.f90#27 REAL :: A1 , A2 , A3 , THIRD , A1SQ , A1CUBE , A1BY3 , Q , R , & #38 R = (2.*A1CUBE-9.*A1*A2+27.*A3)/54. #41 RSQ = R*R #46 THETA = ACOS(R/SQRT(QCUBE)) #52 ARG = (SQRT(RSQ-QCUBE)+ABS(R))**THIRD ! 0 #53 ROOT(1) = -SIGN(1.0,R)*(ARG+Q/ARG) - A1BY3 LUMP.f90#20 COMMON /PLU / S , X , Y , Z , R , U , V , W , USC , PHI , DEN , & #23 F(1) = DEN*USC*R*R ! 6088 MAIN1.f90#670 REAL :: SKEW , R , ALPHPD , BETAPD , ASUB1 , ASUB2 , BSUB1 , & NUMRISE.f90#104 COMMON /PLU / S , X , Y , Z , R , U , V , W , USC , PHI , DEN , & #124 R = REFF #158 FDUM = W*R*R/TP #166 XMAXM = 8.*R*WBYU*(1.+3./WBYU)**2 #183 ZNF = 6.*R*W/UAM ! 0 #297 WRITE (IO6,114) NNP , X , Y , ZC , R , U , V , W , USC , PHI , & RATE.f90#40 COMMON /PLU / S , X , Y , Z , R , U , V , W , USC , PHI , DEN , & #67 RHS(1) = 2.0*R*ALPHA*RA*ABS(USC-UA*U/USC) & #68 & + 2.0*R*BETA*RA*ABS(UA*SIN(PHI)) #71 RHS1A = 2.0*R*RA*UA*DRDXA #74 RHS(2) = -R*R*DEN*W*DUDZ #75 RHS(3) = GRAVI*R*R*(RA-DEN) #76 RHS(4) = -DPDZ*DEN*W*R*R - RP*R*(TP**4-TA**4) UNLUMP.f90#27 COMMON /PLU / S , X , Y , Z , R , U , V , W , USC , PHI , DEN , & #38 R = SQRT(F(1)/USC/DEN) WAKE_DIM.f90#3 SUBROUTINE WAKE_DIM(X,H,W,R,HWAKE,WWAKE) #38 XBYR = X/R ! 24347K #42 XPBYR = -(H/(CWKHT*R))**3 #47 IF ( XBYR.GT.ZERO ) HWAKE = CWKHT*R*DXBYR**THIRD #51 IF ( XBYR.GT.ZERO ) WWAKE = HALF*W + (THIRD*R)*XBYR3RD WAKE_SIG.f90#3 SUBROUTINE WAKE_SIG(X,XD,XOLD,TURBZ,TURBY,SZOLD,SYOLD,H,W,R,ZK,YK,& #41 CALL WAKE_DIM(X,H,W,R,HTWAKE,HWWAKE) !7144148 ZSTREAM.f90#2 SUBROUTINE ZSTREAM(H,W,L,R,LR,HR,X,Y,Z,DZDX) #36 REAL H , W , L , R , HR , LR , EXPZ1 , EXPZ2 , EXPZG , EXPX #66 IF ( X.LT.-R ) THEN #78 ZSLOPE = 2.*(HR-H)*(X+R)*ONEBYR**2. ! 3398 #80 DZDX2 = ZSLOPE/((Z-H+R)*ONEBYR)**EXPZ1 ! 1812 #81 ELSEIF ( R.LE.H .AND. Z.LE.0.67*H ) THEN #83 ELSEIF ( R.LE.H .AND. Z.GT.0.67*H ) THEN #85 ELSEIF ( R.GT.H .AND. Z.LE.0.67*(2*H-R) ) THEN #87 ELSEIF ( R.GT.H .AND. Z.GT.0.67*(2*H-R) ) THEN #102 DZDX2 = ZSLOPE/((Z-H+R)*ONEBYR)**EXPZ1 ! 16930 #111 ZSLOPE = (HR-H)*(R-2.*X)/((L-RBY2+LR)**2) ! 532546 #113 DZDX2 = ZSLOPE/((Z-H+R)*ONEBYR)**EXPZ2 ! 407962 #126 ZSLOPE = ZSLOPELR/((X-(L+LR-R))*ONEBYR)**EXPX #128 DZDX2 = ZSLOPE/((Z-H+R)*ONEBYR)**EXPZ2 !2747174 #138 IF ( YPOS.GT.(WBY2+R/3.) ) THEN !7941052
HyperKWIC - Version 1.00DD executed at 20:00 on 1 Mar 2018 | Personal or Academic or Evaluation User | Free for Non-Commercial, Non-Government Use