2630 IF(ABS(H1-H(1) ) .LT. .000000001) KD=XKDNBR (NN) 



2640 IF(H1.LT.H(1)) KD=XKDBR(NN) 



2650 W=ftH0»(Hl*«3)/(KD»C0T»( ( (RHO/RHQW) -1 . )*»3) ) 



2660 AT=2.*((W/RHO)*«0.33) 



2670 IF(NN.EQ.3) AT=2 . 3* ( (W/RHO) *#0. 33) 



2680 IF(NN.EQ.4) AT= 1 . 03» ( ( W/ (0. 78«RH0 ) ) #*0. 33 ) 



2690 IF(NN.EQ.5) AT= 1 . 361 ♦ ( ( W/ (0. 2B»RH0) ) ♦♦O. 33) 



2700 IF(NN.EQ.6) AT=1 . 502* ( (W/ (0. 495#RH0) ) ♦»0. 33) 



2710 IF(NN.EQ.7) AT=1 . 29* ( ( W/ (0. 176*RH0) ) ♦♦0. 33) 



2720 IF(NN,EQ.8) AT=3. 68» ( ( W/ (6. 48»RH0) ) »«0. 33) 



2730 IF(NN.EQ.9) AT=0. 889* ( ( W/ (0. 083*RHO) ) ♦»0. 33) 



2740 IF(NN.EQ.10) AT=1 . 02» ( (W/ ( . 16*RH0) ) »»0. 33) 



2750 TOPW=3.0»XKDELTA(NN)#( (W/RHO) *«0. 33) 



2760 IFdUNITS.EQ.B) C=l. 



2770 IFdUNITS.EQ. 1) C=.3 



27B0C 



2790C ♦ ♦ INCREASE CREST ELEVATION (C) TO SATISFY HTMAX ♦ < 



2800C 



2810 996 IF(ICHECK.EQ.0) GO TO 500 



2820 IFdUNITS.EQ. 0) C=C+1. 



2830 IFdUNITS.EQ. 1) C=C + .3 



2840 500 HS=C+D 



2850 TANB=1./C0T 



2860 IFdUNITS.EQ. 0) SURF = TANB/SQRT (H2P/ (5. 1 2»T (2) #T (2) ) ) 



2870 IFdUNITS.EQ. 1) SURF = TANB/SQRT (H2P/ ( 1 . 56*T (2) ♦T (2) ) ) 



2880 550 RH = RA(NN)«SURF/d.+RB(NN)«SURF) 



2890 R=H2P»RH ■ 



2900 IF(ANSWIND.EQ."N") GO TO 551 '■ 



2910 IFdUNITS.EQ. 1) W IND=WIND/ , 447 



2920 WF=(WIND*«2)/1800. 



2930 IF(WIND.6T.60. ) WF=2. 



2940 WK=1.+(WF»(C/R+. 1)*SIN(THETA)) 



2950 R=R»WK 



2960 551 FR=C/R 



2970 CR = a.51 -0. lUTOPW/HS 



2980 KT = CR*d.-FR) 



2990 HT=H2P*KT 



3000C 



3010C * ♦ CHECK IF HTMAX CRITERIA MET * * 



3020C 



3030 604 IF(HT.LE. HTMAX) GO TO 605 



3040 ICHECK=l 



3050 GO TO 996 



3060C 



3070C ♦ ♦ VOLUME AND COST COMPUTATIONS * ♦ 



3080C 



3090 605 CONTINUE 



3100 TV0L=(T0PW*HS)+(CHS«*2)*C0T/2.)+( (HS*»2 ) * ( COTP/2 . ) ) 



3110 AH=C+(l.5»Hli 



3120 IF(AH.GT.HS) AH=HS 



3130 SAH = C+(0.5»Hn 



3140 AVOL=(TOPW«AT)+(AT*(AH/SIN(THETA) ))+(AT»(SAH/SIN(THETAP) ) ) 



3150 CVOL=TVOL-AVOL 



B17 



