1 REM ===================================== 2 REM DRAWDOWN.BAS 3 REM 4 REM By John M. Reed 5 REM Commercial Rights Reserved 6 REM 7 REM Edited for COGS Disk by J Phillips 8 REM 26 May 1984 9 REM ===================================== 10 REM Last Modified: october 1985 11 REM 12 REM Program revised and corrected by 13 REM Milovan S. Beljin 14 REM International Ground Water Modeling Center 15 REM 4600 Sunset Ave. 16 REM Indianapolis, Indiana 46208 17 REM 20 BKG=15: REM Inverse background color for IBM-PC 40 REM ********************** 42 REM TITLE PAGE 44 REM ********************** 50 CLS: LOCATE 4,1,1,0,13 52 PRINT" ÍËÍÍ» É " 54 PRINT" º º ÌÍ ÍÍ» É É É ÉÍ͹ ÉÍÍ» É É É ÌÍÍ»" 56 PRINT" º º º ÉÍͺ º º º º º º º º º º º º" 58 PRINT" ÍÊÍͼ ¼ ÈÍͼ Èͼͼ ÈÍͼ ÈÍͼ Èͼͼ ¼ ¼" 60 PRINT"ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ" 62 PRINT" By John M. Reed ////////////³ ³////////////" 64 PRINT" Modified by Milovan S. Beljin ³ ³ " 66 PRINT" ³ ³ " 68 PRINT" ³ ³ " 70 PRINT"Calculates water well drawdown due | ³ ³ | " 72 PRINT"to pumping of a nearby well. V ³ ³ V " 74 PRINT" ------------³ ³------------" 76 PRINT"Uses the Theis non-equilibrium . . . . . . ³ ³ . . . . . ." 78 PRINT"well formula. . . . . . .³ ³. . . . . . " 80 PRINT" ............³ ³............" 82 PRINT" ............³ ³............" 84 LOCATE 23,1,1,0,13 86 PRINT " to continue; to exit; for help.........."; 88 Q$=INKEY$: IF Q$="" THEN 88 90 IF Q$=CHR$(27) THEN 20000 92 IF Q$="H" OR Q$="h" THEN 1500 94 IF Q$<>" " THEN 88 96 REM ********************** 97 REM DATA ENTRY 98 REM ********************** 100 CLS 102 PRINT TAB(24)"DRAWDOWN CALCULATION--DATA ENTRY" 104 PRINT: PRINT: PRINT 106 INPUT "Flow Rate (Gal/Minute).................= ",G 112 INPUT "Transmissivity (Gal/Day/Ft)............= ",T 114 INPUT "Radius of Effect (Feet)................= ",R 116 INPUT "Storage Cofficient (0.0-1.0)...........= ",S 118 INPUT "Time (Days)............................= ",P 120 U=(1.87*(R*R)*S)/(T*P) 125 PRINT:PRINT 170 REM ********************** 172 REM COMPUTE W(U) 174 REM ********************** 200 IF U > 10 THEN WU = 0 :GOTO 310 210 IF U > 7.45 THEN WU = 9.999999E-06 :GOTO 310 220 E0 = .5772156649# + LOG(U) 230 KN = 1: ESUM = 0 240 FOR NN = 1 TO 100 250 KN = KN * NN 260 E1 = (-U)^NN/(NN*KN) 270 IF ABS(E1) < 9.999999E-06 GOTO 300 280 ESUM = ESUM + E1 290 NEXT NN 300 WU = -E0 - ESUM 310 D=114.6*G*WU/T 790 REM ********************** 792 REM SCREEN OR PRINTER 794 REM ********************** 800 CLS: LOCATE 12,1,1,0,13 805 PRINT "Send results to screen or printer (S/P)........."; 810 P$=INKEY$: IF P$="" THEN 810 815 IF P$="s" THEN P$="S" 820 IF P$="p" THEN P$="P" 825 IF P$<>"S" AND P$<>"P" THEN 810 830 PRINT P$ 835 IF P$ = "S" THEN OPEN "SCRN: " FOR OUTPUT AS #1: GOTO 900 840 REM --- PRINTER PROMPT --------- 845 CLS: PRINT "READY TO PRINT.": PRINT " Adjust Paper." 850 PRINT " Turn on Printer.": PRINT " Press a Key to Continue....."; 855 Q$ = INKEY$: IF Q$ = "" THEN 855 860 PRINT Q$ 865 OPEN "LPT1:" FOR OUTPUT AS #1 870 PRINT #1, CHR$(27);"E";: REM EPSON command for emphasized type 890 REM ********************** 892 REM DISPLAY RESULTS 894 REM ********************** 900 IF P$ = "S" THEN CLS: LOCATE 3,1 910 PRINT #1, "=============================================" 915 PRINT #1, "D R A W D O W N C A L C U L A T I O N " 920 PRINT #1, "Using the Theis non-equilibrium well formula " 925 PRINT #1, "=============================================" 930 PRINT #1, " " 935 PRINT #1, "DATA ENTERED" 940 PRINT #1, " Flow Rate (Gal/Minute).................=";G 950 PRINT #1, " Transmisivity (Gal/Day/Ft).............=";T 955 PRINT #1, " Radius of Effect (Feet)................=";R 960 PRINT #1, " Storage Cofficient (0.0-1.0)...........=";S 965 PRINT #1, " Time (Days)............................=";P 990 PRINT #1, " " 1000 PRINT #1, "CALCULATED RESULTS" 1005 PRINT #1, " Value of U........ = ";U 1010 PRINT #1, " Well Function W(U) = ";WU 1030 PRINT #1, " DRAWDOWN in feet = "; 1040 PRINT #1, USING "####.##";D 1050 CLOSE #1 1060 IF P$="S" THEN GOSUB 2100 1070 GOTO 50 1490 REM ********************** 1492 REM HELP SCREENS 1494 REM ********************** 1496 REM --- Page 1 ----- 1500 CLS: PRINT TAB(31)"D R A W N D O W N" 1505 PRINT 1510 PRINT "This program calculates drawdown in a water well due to pumping of another" 1512 PRINT "nearby water well." 1515 PRINT 1520 PRINT "It uses the Theis non-equilibrium well formula which has the form:" 1525 PRINT " " 1530 PRINT " " 1535 PRINT " s = 114.6*Q*W(u)/T u=1.87*r*r*S/T*t 1540 PRINT " " 1545 PRINT " where 1550 PRINT " " 1555 PRINT " s = drawdown; 1560 PRINT " Q = discharge rate (gal/min); 1565 PRINT " T = transmissivity (gal/day/ft); 1570 PRINT " r = radius of effect (ft); 1575 PRINT " S = storage coefficient; 1580 PRINT " t = time since pumping started (days); 1585 PRINT " W(u)= well function. 1980 GOSUB 2100 1990 REM --- RETURN TO PROGRAM ----- 2000 GOTO 50 2090 REM --- PAUSE SUBROUTINE ----- 2100 LOCATE 24,1,1,0,13 2105 PRINT "Press a key to continue....."; 2110 Q$ = INKEY$: IF Q$ = "" THEN 2110 2120 RETURN 19990 REM ********************** 19992 REM EXIT ROUTINE 19994 REM ********************** 20000 CLS: PRINT TAB(28) "E X I T P R O G R A M" 20005 LOCATE 8,20: PRINT "EXIT TO:" 20010 M$(1) = " 1. MENU, to run another program. " 20015 M$(2) = " 2. BASIC " 20020 M$(3) = " 3. MS-DOS System level " 20025 M$(4) = " 4. CANCEL Exit Command " 20030 FOR Y = 1 TO 4: LOCATE 2*Y + 8, 20: PRINT M$(Y): NEXT Y 20035 LOCATE 23,1: PRINT "Enter number or use arrows to make selection; then press " 20037 LOCATE ,,0: REM Cursor OFF 20040 Y=1: MVE=NO 20045 WHILE Q$ <> CHR$(13) 20050 LOCATE 2*Y + 8, 20: COLOR 0,BKG: PRINT M$(Y): COLOR 7,0 20055 Q$=INKEY$: IF Q$="" THEN 20055 20060 IF LEN(Q$) = 2 AND ASC(RIGHT$(Q$,1)) = 72 THEN MVE=YES: X=Y-1: IF X=0 THEN X=4 20065 IF LEN(Q$) = 2 AND ASC(RIGHT$(Q$,1)) = 80 THEN MVE=YES:X=Y+1: IF X=5 THEN X=1 20067 IF Q$ >="1" AND Q$ <= "4" THEN X=VAL(Q$): MVE=YES 20070 IF MVE=YES THEN LOCATE 2*Y + 8, 20: PRINT M$(Y): Y=X: MVE=NO 20075 WEND 20080 LOCATE ,,1: REM Cursor ON 20085 CLS 20100 IF Y = 1 THEN LOCATE 12,1: PRINT "Loading.....MENU": RUN "MENU" 20110 IF Y = 2 THEN END 20120 IF Y = 3 THEN SYSTEM 20130 IF Y = 4 THEN 50: REM Title Page 30000 CLS: END