From owner-freebsd-questions@FreeBSD.ORG Mon Jun 17 17:16:53 2013 Return-Path: Delivered-To: freebsd-questions@freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [8.8.178.115]) by hub.freebsd.org (Postfix) with ESMTP id 80A02170 for ; Mon, 17 Jun 2013 17:16:53 +0000 (UTC) (envelope-from christopher.maness@gmail.com) Received: from mail-vb0-x22b.google.com (mail-vb0-x22b.google.com [IPv6:2607:f8b0:400c:c02::22b]) by mx1.freebsd.org (Postfix) with ESMTP id 4481D128C for ; Mon, 17 Jun 2013 17:16:53 +0000 (UTC) Received: by mail-vb0-f43.google.com with SMTP id e12so2154985vbg.2 for ; Mon, 17 Jun 2013 10:16:52 -0700 (PDT) DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20120113; h=mime-version:sender:date:x-google-sender-auth:message-id:subject :from:to:content-type; bh=JbIGRhyAWyvQgfdVhoZYDbh7BHOuVpPeBVXuEbFtfnM=; b=mi0PYHdDxEtrrAZFE9VoOkWTYp/oy1ypWKKhQYpVO79WOoWD6LXGwbnvdu0ScwROgU zdFKMhlUbTGBE+Ib7Fxaot4eDnPpyGewS1jDleb0B+gIuupUAwwTM1EgLf8JuvEqLXUy jTxd0W6jr1pGn5cxVSQ1CKt2ARnBiV9Qgv4IpGS4xpVZ04qaRp2opfN2XVm/RIoDYJa6 FuPm3lNZFWdkZNY7lofQjNFNCgzM1OCMu2ZYCqh3qD5yzCuvJm0VEkRbcEnI/7tjp3Er K1lm2HY1ORvHX0YO+HaIBmH0VkSWIGrvnf3tmIIElhr1ga2XxrD8lKIvd2mshf1ixXKk DfRg== MIME-Version: 1.0 X-Received: by 10.220.43.14 with SMTP id u14mr583574vce.35.1371489412667; Mon, 17 Jun 2013 10:16:52 -0700 (PDT) Sender: christopher.maness@gmail.com Received: by 10.58.220.105 with HTTP; Mon, 17 Jun 2013 10:16:52 -0700 (PDT) Date: Mon, 17 Jun 2013 10:16:52 -0700 X-Google-Sender-Auth: v7opXDQoW9xLEXJClJIBr8sg1Cc Message-ID: Subject: Any BASIC Gurus around? From: Chris Maness To: "freebsd-questions@freebsd.org" Content-Type: text/plain; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable X-Content-Filtered-By: Mailman/MimeDel 2.1.14 X-BeenThere: freebsd-questions@freebsd.org X-Mailman-Version: 2.1.14 Precedence: list List-Id: User questions List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Mon, 17 Jun 2013 17:16:53 -0000 I am having trouble getting this old USNO basic program running in bwBASIC. The error output is not clear to me where the problem is. Here is the code, if someone wouldn't mind running it and suggesting edits: 10 DEF FNARCOS(ARG)=3D1.570796-ATN(ARG/SQR(1.-ARG*ARG)) 20 DEF FNARCSIN(ARG)=3DATN(ARG/SQR(1.-ARG*ARG)) 30 DEF FNDEG(ARG)=3DINT(ARG)+((ARG-INT(ARG))*1O.)/6. 40 DEF FNDMS(ARG)=3DINT(ARG)+6. * (ARG - INT(ARG)) / 10! 50 RD=3D57.29578 60 DR=3D1./RD 70 DIM A(4) 80 DIM B(2) 90 A(1)=3D-.01454 100 A(2)=3D-.10453 110 A(3)=3D-.20791 120 A(4)=3D.00233 130 CE=3D.91775 140 SE=3D.39715 150 INPUT "LONGITUDE IN DEG."; LO 160 LI=3DABS(LO) 170 INPUT "LATITUDE IN DEG."; F 180 F=3DF*DR 190 S1=3DSIN(F) 200 CI=3DCOS(F) 210 INPUT "YEAR (4 DIGITS)"; IY 220 INPUT "MONTH (NUMERAL)"; IM 230 INPUT "DAY (NUMERAL)"; ID 240 C=3D360. 250 J=3D367*IY-INT(7*(IY+INT((IM+9)/12))/4)+INT(275*IM/9)+ID-730531. 260 INPUT "UNIVERSAL TIME =3D 0, ZONE TIME =3D 1, LOCAL MEAN TIME =3D 2"; Z 270 DT=3D0. 280 IF Z=3D0. THEN LET DT=3D-LO/C 290 IF Z=3D1. THEN LET DT=3D-(LI-15*INT((LI+7.5)/15))/C*SGN(LO) 300 INPUT "HOUR (4 DIGIT NUMERAL ON 24 HOUR CLOCK)"; H 310 ZO#=3DJ-.5 320 IF H>0 THEN GOTO 870 330 PRINT "DATA FOR ";IY; ", MONTH ";IM; ", DAY";ID 340 FOR L=3D1 TO 4 350 ON L GOTO 370, 650, 650, 360 360 C=3D347.81 370 M=3D.5+DT 380 K=3D1 390 M=3DM-DT 400 E=3DM-LO/360. 410 GOSUB 430 420 GOTO 530 430 D#-ZO#+E 440 IF ABS(E)>=3D1 THEN LET E=3DE-SGN(E) 450 GOSUB 1220 460 IF L=3D4 THEN GOSUB 1720 470 T=3DT+LO+360.*E 480 T=3DT-INT(T/360.)*360. 490 U=3DT-AS 500 IF ABS(U) > 180! THEN LET U=3DU-360.*SGN(U) 510 U=3DU/C 520 RETURN 530 M=3DM-U+DT 540 IF L<4 THEN LET K=3DK+1 550 ON K GOTO 600,560,600,580,600,620 560 IF M>=3DO. AND M<1. THEN GOTO 620 570 GOTO 590 580 IF M>=3DO. THEN GOTO 620 590 M=3DM-SGN(M) 600 K=3DK+1 610 GOTO 390 620 H=3DFNARCSIN(COS(F-DS))*RD 630 IF L=3D4 THEN LET H=3DH-.95*COS(H) 640 GOSUB 2160 650 GOSUB 2000 660 B(1)=3DM-H 670 B(2)=3DM+H 680 FOR I=3D1 TO 2 690 K=3D2*I-3 700 FOR N=3D1 TO 6 710 B(I)=3DB(I)-DT 720 E=3DB(I)-LO/360. 730 GOSUB 430 740 GOSUB 2000 750 B(I)=3DB(I)+K*H-U+DT 760 IF L<4 THEN LET N=3DN+1 770 ON N GOTO 820,780,820,800,820,830 780 IF B(I)>=3DO. AND B(I)<1. THEN GOTO 830 790 GOTO 810 800 IF B(I)>=3DO. THEN GOTO 830 810 B(I)=3DB(I)-SGN(B(I)) 820 NEXT N 830 NEXT I 840 ON L GOSUB 1350,1400,1400,1610 850 NEXT L 860 GOTO 150 870 INPUT"SKY CONDITION =3D 1,2,3,10,";SK 880 PRINT"DATA FOR ";IY;", MONTH ";IM;", DAY ";ID;", AT ";H;" HOURS" 890 E=3DFNDEG(H/100.)/24.-DT-LO/360. 900 D#=3DZO#+E 910 N=3D1 920 GOSUB 1220 930 T=3DT+360.*E+LO 940 IF N=3D2 THEN GOSUB 1720 950 H=3DT-AS 960 GOSUB 2060 970 Z=3DH*DR 980 H=3DH-.95*(N-1)*COS(H*DR) 990 GOSUB 2160 1000 GOSUB 2200 1010 HA=3DINT(ABS(HA)+.5)*SGN(HA) 1020 ON N GOTO 1030, 1090 1030 IS=3D133775.*M/SK 1040 PRINT"SUN AZIMUTH (DEG.) ";AZ 1050 PRINT"SUN ALTITUDE (DOG.) ";HA 1060 PRINT"SUN ILLUMINANCE (LUX) ";IS 1070 N=3D2 1080 GOTO 940 1090 E=3DFNARCOS(COS(V-LS)*CB) 1100 P=3D.892*EXP(-3.343/((TAN(E/2.))^.632))+.0344*(SIN(E)-E*COS(E)) 1110 P=3D.418*P/(1.-.005*COS(E)-.03*SIN(Z)) 1120 IL=3DP*M/SK 1130 IS=3DIS+IL+.0005/SK 1140 PRINT"MOON AZIMUTH (DEG.) ";AZ 1150 PRINT"MOON ALTITUDE (DEG.) ";HA 1160 PRINT"MOON ILLUMINANCE (LUX) ";IL 1170 IL=3DINT(50.*(1.-COS(E=BB+.5) 1180 PRINT" (";IL;"% OF MOON ILLUMINATED)" 1190 PRINT"TOTAL ILLUMINANCE (LUX) ";IS 1200 GOTO 300 1210 END 1220 TD#=3D280.46#+.98565#*D# 1230 T=3DTD#-INT(TD#/360#)*360# 1240 IF T=3D4800. OR R1. THEN GOTO 2040 2020 H=3DFNARCOS(H)*RD/C 2030 RETURN 2040 H=3D1.5 2050 RETURN 2060 CD=3DCOS(PS) 2070 CS=3DCOS(H*DR) 2080 Q=3DSD*CI-CD*SI*CS 2090 P=3D-CD*SIN(H*DR) 2100 AZ=3DATN(P/Q)*RD 2110 IF Q