Т > NGRCalc 1.03 Т SТ 28-Feb-1997 v1.00 Program written by Len Killip, G0APZ, for his own interest. (Т 2OТ 15-Oct-1999 v1.01 JGH: Tidied output slightly, QUITs on exit if possible, <0Т windows not defined FNТ 18-Apr-2009 v1.02 JGH: UTM grid for Channel Islands, command line access PAТ 25-Jun-2009 v1.03 JGH: Works in interactive mode on Windows Z;Т Note: Some issues on SDL and Brandy dТ nТ Len's original comments xТ ----------------------- ┌DТ This program converts National Grid References to latitude and ▄IТ longitude, and calculates the deviation between true north and grid √Т north. ═Т ╙EТ Also the reverse, lat and long to easting and northing and NGR. Є1Т Also the bearing and distance between NGRs. ЎCТ The "mathematical engine" is derived from the Ordnance Survey х<Т "Geodetic Information Paper No 1: 2/1996 (Version 2.0) рCТ 10 and 8 figure as well as 6 figure NGRs may be entered. Also эAТ Eastings and Northings (to nearest millimetre if desired!). ФGТ Using OS worked examples, eastings and northings entered directly ПAТ to the nearest millimetre(!) result agrees to within 0.0004 ЗHТ seconds lat. and is exact to four decimal places of seconds long. BТ (Exact to four decimal places in both cases using BASIC VI). FТ Deviation checks correct at sheet corners, and exact at Caister. DТ Bearings checked on p19 Fram to Caister 7.488 seconds v. 7.487 "FТ Caister to Fram 10.833 seconds v. 10.832. I'll settle for these! ,:Т Scale factor for mid-point of that line also checks. 6: @2A$=єOS_GetEnv:ver$="1.03":Г os%=32 ▄ *FLOAT 64 JГ A$<>"":Рcmd(A$):Рexit(0) T8Г os%=32:х≥ "ShowWindow",@hwnd%,3:Р_setfocus(@hwnd%) ^: hК 27+128 rРborder |Н ┘ Рerror ├: ░Рconstants_OSGB : єrepeat=╧ ўУ ╦Рchoice1 бЩrepeat<>╧ лЮ ж: Ю щРchoice1 Й4ш:Я'┼15)"GRID REFERENCE CONVERTOR"'┼15)д24,"=")' Т!Я┼18)"CURRENT GRID: ";grid$'' Ч?Я┼5)"(A) CONVERT GRID REFERENCE TO LATITUDE AND LONGITUDE"' 6Я┼5)"(B) BEARING AND DISTANCE BETWEEN TWO POINTS"' ?Я┼5)"(C) CONVERT LATITUDE AND LONGITUDE TO GRID REFERENCE"' !Я┼5)"(D) SELECT A BASE GRID"' &Я┼5)"(E) EXIT"' 0Я┼5)"Press a key:"; :х▌ єanswer("ABCDE") й Dи 1 Nlatlongrepeat=╧ Xх∙ latlongrepeat=╧ b ш:Я'' lРlatlong vн ─и 2 ┼bdrepeat=╧ ■х∙ bdrepeat=╧ · ш:Я'' ╗Р_bd ╡н ╪и 3 ф lltongr=╧ пх∙ lltongr=╧ з ш:Я'' Д Р_lltongr Нн Ьи 4 Рchoosegrid и 5 Я:quit%=╧:Н ┘ quit%=ё Гquit%:х≤ *Ю 4к >А H: RщРchoosegrid \Я''┼5)"SELECT GRID:"' f0Я┼5)"(A) OSGB - Great Britain & Isle of Man" p$Я┼5)"(B) UTM - Channel Islands" z'Я┼5)"(C) ITM - Ireland (new grid)" └'Я┼5)"(D) OSI - Ireland (old grid)" ▌Я'┼5)"Press a key:"; ≤х▌ єanswer("ABCD") й ╒и 1:Рconstants_OSGB ╛и 2:Рconstants_UTM Іи 3:Рconstants_ITM юи 4:Рconstants_OSI йк тА ч: Х щРlatlong РРchoice2 Э ш:Я'' Я ┼10)"GRID REFERENCE: "G$' Р_Mvrho(E,N) Р_toll $ Р_C_EN .$Рoutput(lambda,phi,C):Рlldisplay 8<Я '┼10) "Do you want another Lat Long calculation? Y/N"; B&Г єanswer("YN")<>1 latlongrepeat=ё LА V: ` щР_bd jЯ ┼20)"FIRST POSITION"' tРchoice2 ~ E1=E:N1=N ┬Р_Mvrho(E1,N1):y1%=y:F1=F ▓ Р_toll ° Р_C_EN іЯ ┼20)"SECOND POSITION"' ╟Рchoice2 ╨ E2=E:N2=N дР_Mvrho(E2,N2):y2%=y:F2=F н Р_toll ьх▌ N2 й Би N1 Л alpha=╞/2 Ж alpha=≥((E2-E1)/(N2-N1)) к х▌ E2 й и E1 (alpha=0 2к <Г N2N1 ▄ Zalpha=alpha+2*╞ dм nм x%Nm=(N1+N2)/2:q=N1-N2:Em=(E1+E2)/2 ┌Г y1%<0 ┌ y2%<0 ▄ ▄Рalphasplit √л ═Р_converge ╙м ЄР_distance ЎР_bd_output х6Я ┼15)"Another bearing/dist calculation? Key Y/N"; р!Г єanswer("YN")<>1 bdrepeat=ё эА Ф: ПщР_lltongr З7Я ┼15)"Enter latitude north in form deg, min, sec"' 6Х ┼15)deg%, min%,sec':degN=deg%:minN=min%:secN=sec 7Я ┼15)"Enter longitude in form deg, min, sec, E/W"' ;Х ┼15)deg%, min%, sec,ew$':degE=deg%:minE=min%:secE=sec "4Р_lltoEN:ш:Я''┼15)"Easting= ";E;" Northing= ";N' ,=Я ┼15) "Latitude "; degN; " Degs "; minN; " Mins "; secN; 6Я " Secs North"' @>Я ┼15) "Longitude "; degE; " Degs "; minE; " Mins "; secE; JЯ " Secs "; ew$' TР_ENtongr(E,N) ^Р_ngrdisplay hА r: | щР_lltoEN ├7Т On entry, degN, minN, secN = latitude north ░<Т degE, minE, secE, ew$ = longitude west/east +phi=єconvert(degN,minN,secN):phi=╡(phi) є(v=єA(phi):rho=єB(phi):itasqd=v/rho-1 ў6lambda=єconvert(degE, minE, secE):lambda=╡(lambda) ╦Lew$=юew$,1):Г ew$="W" └ ew$="w" ▄ lambda=-lambda:ew$="West" ▀ ew$="East" бP=lambda-lambda0:phip=phi лM=єM(b,n,phip,phi0) ж I=M+N0 ЮII=(v/2)*╣(phi)*⌡(phi) Й8III=(v/24)*╣(phi)*(⌡(phi))^3*(5-(Ї(phi))^2+9*itasqd) Т@IIIA=(v/720)*╣(phi)*(⌡(phi)^5)*(61-58*(Ї(phi)^2)+(Ї(phi)^4)) ЧN=I+P^2*II+P^4*III+P^6*IIIA IV=v*⌡(phi) )V=(v/6)*(⌡(phi))^3*(v/rho-(Ї(phi)^2)) VI=(v/120) * (⌡(phi)^5) &EVI=VI*(5-18*(Ї(phi)^2)+(Ї(phi)^4)+14*itasqd-58*(Ї(phi)^2)*itasqd) 0E=E0+P*IV+P^3*V+P^5*VI :А D: NщР_ENtongr(E,N) XГgrid$="ITM":E=E+100000 bГ E<0:E=E+.1:Т IF E<0 E=E-1 lГ N<0:N=N+.1:Т IF N<0 N=N-1 vРen_string(E) ─ east$=en$ ┼Рen_string(N) ■north$=en$ ·: ╗Рgrid_letters ╡: ╪,easting$=бeast$,5): northing$=бnorth$,5) фГ E<0 ▄ п(easting$=бц(1000000-╩(аeast$,2))),5) зм ДГ N<0 ▄ Н)northing$=бц(500000-╩(аnorth$,2))),5) Ьм !NGR$=A$+B$+easting$+northing$ А : щР_ngrdisplay *6Я ┼15)"10-digit National Grid Reference is ";NGR$' 4Рcurtail(easting$) >easting$=en$ HРcurtail(northing$) Rnorthing$=en$ \!NGR$=A$+B$+easting$+northing$ f5Я ┼15)"6-digit National Grid Reference is ";NGR$' pDЯ ┼15)"4-digit National Grid Reference is ";юNGR$,4);аNGR$,6,2)' z: └:Я '┼15)"Another lat/long to NGR calculation? Key Y/N"; ▌: ≤ Г єanswer("YN")<>1 lltongr=ё ╒А ╛: ІщРcurtail(EN$) юen%=╗(╩(EN$)/100) йen$=б"0000"+ц(en%),3) тА ч: Х щРchoice2 Р$Я ┼15)"DO YOU WISH TO ENTER AS"' Э,Я ┼15)"(A) NATIONAL GRID REFERENCE, OR"' &Я ┼15)"(B) EASTING AND NORTHING?"' Я ┼15)"KEY A OR B"' х▌ єanswer("AB") й $и 1 .?Я ┼2)"Enter the NGR, with grid letters, in 10, 8, 6 or 4 "; 8Я "numeral form "' BХ ┼15) G$' LРngr Vи 2 `:Я ┼10) "ENTER FULL EASTING AND NORTHING, IN FORM E,N"' jFЯ ┼5 ) "(six or seven figures in each, plus decimals if need be)"' tХ ┼20) E,N' ~G$=цE+","+цN ┬к ▓А °: і-Т This puts NGR into easting and northing ╟ щРngr ╨Г (╘ G$ ─ 1)=1:G$="I"+G$ д$Г юG$,1)>"`":G$=Ґ(≈G$-32)+аG$,2) н1Г аG$,2,1)>"`":G$=юG$,1)+Ґ(≈аG$,2)-32)+аG$,3) ьРconstants_OSGB БГ юG$,1)="I":Рconstants_ITM ЛГ юG$,1)="W":Рconstants_UTM Ж: X% = (╘(G$) - 2) / 2 Y% = 10^(5-X%) !X$ = юG$, 1): Y$ = аG$, 2, 1) 1E = ╩(аG$, 3, X%)) * Y%: N = ╩(бG$, X%)) * Y% (Г юG$,1)="W" ▄ 2%a%=≈аG$,2,1)─&DF:Гa%=≈"A":a%=≈"W" 13 ▄ a%=a%-1 рrow%=a% │ 5: column%=a% ┐ 5 э*E=E+(column%-c%)*m% : N=N+(r%-row%)*m% ФА П: ЗNщєA(phip)=a/І(1-esqd*(╣(phip))^2) :Т rad of curv of lat at lat phip KщєB(phip)=v*(1-esqd)/(1-esqd*(╣(phip))^2) :Т same for long at lat phip : АщєM(b,n,phip,phi0)=b*(((1+n+(5/4)*n^2+(5/4)*n^3)*(phip-phi0))-((3*n+3*n^2+(21/8)*n^3)*╣(phip-phi0)*⌡(phip+phi0))+(((15/8)*n^2+(15/8)*n^3)*╣(2*(phip-phi0))*⌡(2*(phip+phi0)))-((35/24)*n^3*╣(3*(phip-phi0))*⌡(3*(phip+phi0)))) "тtmp=b * (((1+n+(5/4)*n^2+(5/4)*n^3)*(phip-phi0))-((3*n+3*n^2+(21/8)*n^3)*╣(phip-phi0)*⌡(phip+phi0))+(((15/8)*n^2+(15/8)*n^3)*╣(2*(phip-phi0))*⌡(2*(phip+phi0)))-((35/24)*n^3*╣(3*(phip-phi0))*⌡(3*(phip+phi0)))) ,=tmp 6: @JТ This takes easting and northing, does iteration for M, returns phip, JТ v, rho, y, and itasqd. T: ^щР_Mvrho(E,N) h y=E-E0 rphip=(N-N0)/a+phi0 |M=єM(b,n,phip,phi0) ├У ░phin=(N-N0-M)/a+phip phip=phin єM=єM(b,n,phip,phi0) ўЩ ■(N-N0-M)<.0015 ╦: бv=єA(phip) лrho=єB(phip) жitasqd=v/rho-1 ЮР_F ЙА Т: ЧGТ This contains equations leading to phi, lambda. For equation IX I DТ need to divide (720*rho*v^5) by v^2 to get within range, which IТ multiplies IX by v^2. So then divide IX by v^2. Similarly for XIIA; IТ v^7 outside number range. In XIIA reduce v^7 to v^3 and divide XIIA &Т by v^4 0: :щР_toll DVII=Ї(phip)/(2*rho*v) NIVIII=(Ї(phip)/(24*rho*v^3))*(5+3*Ї(phip)^2+itasqd-9*Ї(phip)^2*itasqd) X>IX= (Ї(phip)/(720*rho*v^3))*(61+90*Ї(phip)^2+45*Ї(phip)^4) b IX=IX/v^2 l$phi=phip-y^2*VII+y^4*VIII-y^6*IX vX=1/(v*⌡(phip)) ─.XI=1/(6*v^3*⌡(phip))*(v/rho+2*(Ї(phip))^2) ┼?XII=(1/(120*v^5*⌡(phip)))*(5+28*(Ї(phip))^2+24*(Ї(phip))^4) ■UXIIA=(1/(5040*v^3*⌡(phip)))*(61+662*(Ї(phip))^2+1320*(Ї(phip))^4+720*(Ї(phip))^6) ·XIIA=XIIA/v^4 ╗2lambda=lambda0+y*X-y^3*XI+y^5*XII-y^5*XIIA*y^2 ╡P=lambda-lambda0 ╪А ф: пGТ This derives deviation (C) between true north and grid north from зCТ phi and lambda (PROC_toll) and itasqd (PROC_Mvrho). When C is Д-Т negative, true N is East of grid north. НщР_C ЬXIII=╣(phi) 7XIV=((╣(phi)*(⌡(phi))^2)/3)*(1+3*itasqd+2*itasqd^2) .XV=((╣(phi)*(⌡(phi))^4)/15)*(2-(Ї(phi))^2) #C=(P*(XIII)+P^3*(XIV)+P^5*(XV)) А *: 4-Т This derives deviation (C) from E and N >щР_C_EN HXVI=Ї(phip)/v R8XVII=Ї(phip)/(3*v^3)*(1+Ї(phip)^2-itasqd-2*itasqd^2) \6XVIII=Ї(phip)/(15*v^5)*(2+5*Ї(phip)^2+3*Ї(phip)^4) fC=y*XVI-y^3*XVII+y^5*XVIII pА z: └щРoutput(lambda,phi,C) ▌ lambda=²(lambda): phi=²(phi) ≤ew$ = "East" ╒*Г lambda <=0:ew$="West":lambda=-lambda ╛Рdegminsec(phi) І degN=deg: minN=min: secN=sec юРdegminsec(lambda) й degE=deg: minE=min: secE=sec тА ч: ХщРlldisplay РDЯ ┼10) "The Latitude is "; degN; " Degs "; minN; " Mins "; secN; ЭЯ " Secs North"' EЯ ┼10) "The Longitude is "; degE; " Degs "; minE; " Mins "; secE; Я " Secs "; ew$' C=²(C) $%Г C>0 C$="West" ▀ C$="East": C=-C .Рdegminsec(C) 8WЯ ┼10) "True North is ";deg " Degs ";min " Mins ";sec;" Secs ";C$;" of grid North"' B(Я ┼10) "easting= ";E;" northing=";N' LА V: `щР_bd_output jG=²(alpha) tРdegminsec(G) ~Рnumplace(sec," ") ┬sec=num ▓=Я ┼15)"Grid bearing is ";deg" deg ";min" min ";sec" sec"' °)Г C<0 ─ alpha<=■(C) ▄ alpha=alpha+2*╞ іTBg=²(alpha+C-convangle) ╟Рdegminsec(TBg) ╨Рnumplace(sec," ") дsec=num н>Я ┼15)"True bearing is ";deg;" deg ";min" min ";sec" sec"' ьРnumplace(s," ") Б s=num Л*Я ┼15)"Grid distance is ";s;" metres"' ЖРnumplace(dist," ") dist=num -Я ┼15)"True distance is ";dist;" metres"' Я ┼20)"Print vars? Y/N"; !Г єanswer("YN")=1 Р_variables (А 2: <щРerror F!Г ÷=17 ▄ А:Т REPORT:PRINT:END P3Т IF ERL<>1240 THEN REPORT:PRINT" at line ";ERL ZЯ dЯ ┼13)"Press any key"; n У Щ ╔ xА ┌: ▄4Т This draws yellow border, with blue background √щРborder ═1О19,0,4,0,0;19,3,3,0,0; :Т select yellow/blue ╙+Л 0,0:ъ 0,959:ъ 1279,959 :Т draw border Є!ъ 1279,0:ъ 0,0 :Т draw border Ў.Т VDU 28,4,58,78,4 :REM set up text window х8Т VDU 24,10;10;1269;949; :REM set up graphics window рА э: Фщєanswer(alternatives$) П*FX 15,1 ЗУ "A%=їalternatives$, Ґ(╔ ─ &DF)) Щ A%>0 =A% ": ,;Т converts decimal degrees to degrees, minutes, seconds 6щРdegminsec(theta) @6ang=theta*3600+0.0005 : sec=(ang*1000)┐60000 /1000 J+deg=ang │60 : min=deg ┐60 : deg=deg │60 TА ^: h6Т removes last j$ sig figs.Pass space string to j$ rщРnumplace(Y,j$) | num$=ц(Y) ├ бnum$)=j$ ░num=╩(num$) А є: ўOТ Corrects for convergence as in navigation,and checks with OS calc page 18 ╦щР_converge бf=2*y1%+y2% лР_Mvrho(E1,Nm) жXXIII=1/(6*rho*v) Юconvangle=(f*q*XXIII) ЙА Т: ЧIТ Case when E1 and E2 opposite sides central meridian. Nc is point of Т crossing central meridian щРalphasplit ratio=■((y1%)/(E2-E1)) &Nc=N1+ratio*(N2-N1) 0 Nm1=(Nc+N1)/2: Nm2=(Nc+N2)/2 :Nm=Nm1: y2%=0: q=N1-Nc DР_converge Nconvangle1=convangle X"Nm=Nm2:y1%=0:y2%=E2-E0:q=Nc-N2 bР_converge lconvangle2=convangle v#convangle=convangle1+convangle2 ─А ┼: ■щР_F ·XXI=1/(2*rho*v) ╗$XXII=(1-4*itasqd)/(24*rho^2*v^2) ╡F=F0*(1+y^2*XXI+y^4*XXII) ╪А ф: п/Т s is grid distance; dist is True distance зщР_distance Дs=І((E2-E1)^2+(N2-N1)^2) НР_Mvrho(Em,Nm):Fm=F ЬX=(1/6)*(1/F1+4/Fm+1/F2) F=1/X dist=s/F А : *щР_variables 4Я"F= ";F >Я"phip= ";phip HЯ"v= ";v RЯ"rho= ";rho \Я"itasqd= ";itasqd fЯ"M= ";M pЯ"y= ";y zЯ"VII= ";VII └Я"VIII= ";VIII ▌Я"IX= ";IX ≤Я"X= ";X ╒Я"XI= ";XI ╛Я"XII= ";XII ІЯ"XIIA= ";XIIA юЯ"XIII= ";XIII йЯ"XIV= ";XIV тЯ"XV= ";XV чЯ"XXI= ";XXI ХЯ"XXII= ";XXII РЯ"XXIII= ";XXIII ЭА : 0щєconvert(deg,min,sec)=(min*60+sec)/3600+deg : $щРgrid_letters .Г grid$="OSGB" ▄ 8)Рfirst_letter(east$,north$) :A$=Ґ(a%) B)Рsecond_letter(east$,north$):B$=Ґ(a%) LА Vм `Г grid$="OSI" ▄ j0Рfirst_letter(east$,north$) :A$=Ґ(a%):A$="I" t)Рsecond_letter(east$,north$):B$=Ґ(a%) ~А ┬м ▓Г grid$="ITM" ▄ °0Рfirst_letter(east$,north$) :A$=Ґ(a%):A$="I" і)Рsecond_letter(east$,north$):B$=Ґ(a%) ╟А ╨м дГ grid$="UTM" ▄ нA$=Ґ(≈юnorth$,1)+34) ьB$=Ґ(╩юnorth$,2)+32) Б!Г B$>"V":B$=Ґ(╩юnorth$,2)+10) Лnorth$=аnorth$,2) ЖА м A$="?":B$="?" А : (щРfirst_letter(E$,N$) 2column=0:row=0 <Рcolumn_row(E$) FГ E<0 ▄ Pcolumn=init │ 5 + 2 Zл dcolumn=init │ 5 + 3 nм xРcolumn_row(N$) ┌Г N<0 ▄ ▄ row=4 √л ═row=3-init │ 5 ╙м ЄРread_letter ЎА х: рщРsecond_letter(E$,N$) эРcolumn_row(E$) ФГ E<0 ▄ Пcolumn=5+init ┐ 5 Зл column=1+init ┐ 5 м Рcolumn_row(N$) "Г N<0 ▄ ,init=init+.1:N=N-.1 6 Г бц╗N,2)<>"00":init=init-.1 @row= -init ┐ 5 Jл Trow=4-init ┐ 5 ^м hРread_letter rА |: ├CТ "column_row" returns the initial one or two numbers, with "-" ░Т if present, as "init". : єщРcolumn_row(EN$) ўneg$=юEN$,1) ╦ M%=╘(EN$) бГ neg$="-" ▄ л Г M% =8 ▄ жРleft(EN$,3) Юм Й Г M% =7 ▄ ТРleft(EN$,2) Чм л Г M%=7 ▄ Рleft(EN$,2) &м 0Г M%=6 ▄ :Рleft(EN$,1) Dм Nм XА b: lBТ Preserves leading zeros, and repositions "-" when necessary. vщРen_string(EN) ─EN$=ц(╗((EN+0.5)*10/10)) ┼ L%=╘(EN$) ■Г EN<0 ▄ ·L%=L%-1 ╗neg$="-" ╡Г L%>7 ▄ ╪EN$=бEN$,7) фen$=б"0000000"+EN$,7) пл зEN$=бEN$,L%) Дen$=б"000000"+EN$,6) Нм Ьen$=("-"+en$) м Г EN>=0 ▄ Г L%>6 ▄ en$=б"0000000"+EN$,7) *л 4en$=б"000000"+EN$,6) >м Hм RА \: fщРleft(x$,y) pГ ╩(x$)<0 ▄ y=2 zinit=╩(юx$,y)) └А ▌: ≤щРread_letter ╒a%=5*row+column ╛Г a%<9 a%=a%-1 Іa%=a%+65 юА й: тщРconstants_OSGB ч8E0=400000: N0=-100000 : Т grid coords of True origin Х3phi0=╡(49) : Т latitude of True origin Р=lambda0=╡(-2) : Т longitude of True origin (West -) Э7F0=0.9996012717 : Т scale on central meridian Da=6375020.481 : Т semi-major axis in metres scaled by F0 9b=6353722.490 : Т semi-minor similarly scaled /esqd=0.006670539762 : Т e is eccentricity $Dn=(a-b)/(a+b) : Т 0.001673220250 : REM n is (a-b)/(a+b) .east%=0:grid$="OSGB" 8А B: LщРconstants_OSI V8E0=200000: N0=250000 : Т grid coords of True origin `3phi0=╡(53.50) : Т latitude of True origin j=lambda0=╡(-8) : Т longitude of True origin (West -) t7F0=1.000035 : Т scale on central meridian ~Da=6377563.395906615 : Т semi-major axis in metres scaled by F0 ┬9b=6356256.908205645 : Т semi-minor similarly scaled ▓/esqd=0.00667054015 : Т e is eccentricity °.n=(a-b)/(a+b) : Т n is (a-b)/(a+b) іeast%=0:grid$="OSI" ╟А ╨: дщРconstants_ITM н8E0=600000: N0=750000 : Т grid coords of True origin ь3phi0=╡(53.50) : Т latitude of True origin Б=lambda0=╡(-8) : Т longitude of True origin (West -) Л7F0=0.999820 : Т scale on central meridian ЖDa=6376988.93534 : Т semi-major axis in metres scaled by F0 9b=6355608.0987234548 : Т semi-minor similarly scaled /esqd=0.00669438002301 : Т e is eccentricity .n=(a-b)/(a+b) : Т n is (a-b)/(a+b) east%=0:grid$="ITM" (А 2: <щРconstants_UTM F:E0=500000: N0=0 : Т grid coords of True origin P5phi0=╡(0) : Т latitude of True origin Z?lambda0=╡(-3) : Т longitude of True origin (West -) d9F0=0.9996 : Т scale on central meridian nFa=6378388.000 : Т semi-major axis in metres scaled by F0 x;b=6356911.946 : Т semi-minor similarly scaled ┌1esqd=0.0067226700223333 : Т e is eccentricity ▄Fn=(a-b)/(a+b) : Т 0.0016863407 : REM n is (a-b)/(a+b) √east%=0:grid$="UTM" ═А ╙: Є#щРcmd(A$):Т Command line access Ў"Н ┘ Ж:РClose_All:Я·:Рexit(÷):Ю хГ os%=32:РWin_TextIO р,in$=єcl("-i ",1):Г in$="":Рline ▀ Рinput эА Ф: П щРline ЗГ єcl("-?",0):Рsyntax:А Г єcl("-help",0):Рhelp:А Рconstants_OSGB Г юA$,1)>"@" ▄ "G$=єcl("",0):comment$=A$ ,Рngr 6Р_Mvrho(E,N):Р_toll:Р_C_EN @Рoutput(lambda,phi,C) J.Г comment$<>"":Я comment$;┴(16-╘comment$); T"Я ;degN;",";minN;",";secN;" "; ^.Я ;degE;",";minE;",";secE;",";юew$,1);" "; hЯ rл |$Г єcl("-osgb",0):Рconstants_OSGB ├"Г єcl("-osi",0):Рconstants_OSI ░"Г єcl("-itm",0):Рconstants_ITM "Г єcl("-utm",0):Рconstants_UTM є.lat$=єcl("",0):long$=єcl("",0):comment$=A$ ў1degN=╩lat$:A%=їlat$+",",","):lat$=аlat$,A%+1) ╦1minN=╩lat$:A%=їlat$+",",","):lat$=аlat$,A%+1) б1secN=╩lat$:A%=їlat$+",",","):lat$=аlat$,A%+1) л: ж5degE=╩long$:A%=їlong$+",",","):long$=аlong$,A%+1) Ю5minE=╩long$:A%=їlong$+",",","):long$=аlong$,A%+1) Й5secE=╩long$:A%=їlong$+",",","):long$=аlong$,A%+1) Т ew$=long$ Ч: Р_lltoEN:Р_ENtongr(E,N) : GL$=юNGR$,2):NGR$=аNGR$,3) &.Г comment$<>"":Я comment$;┴(16-╘comment$); 0Я GL$;NGR$;" "; :"Я GL$;юNGR$,4);аNGR$,6,4);" "; D"Я GL$;юNGR$,3);аNGR$,6,3);" "; N"Я GL$;юNGR$,2);аNGR$,6,2);" "; X"Я GL$;юNGR$,1);аNGR$,6,1);" "; bЯ GL$;" "; lЯ ;east$;" ";north$;" "; vЯ ─м ┼: ■А ·: ╗щРinput ╡in%=▌(in$):Гin%=0:А ╪ У:У:A$=Ў#in%:Ще#in% └ A$<>"" ф Г A$<>"" ─ юA$,1)<>"#":Рline пЩе#in%:ы#in%:in%=0:А з: ДJщРsyntax:Я"Syntax: NGRCalc [-i infile] [-help]":А Н щРhelp Ь9Я"NGRCalc ";ver$" - Convert National Grid References" 4Я" NGRCalc [-osgb -osi -itm -utm]" LЯ" Converts comma-seperated latitude and longitude to grid reference" ;Я" -osgb - Use OSGB (Great Britain) grid (default)" ,Я" -osi - Use OSI (old Irish) grid" *,Я" -itm - Use ITM (new Irish) grid" 42Я" -utm - Use UTM (Channel Islands) grid" ><Я" Outputs space-seperated string of grid references" HЯ" Examples:" R,Я" Command: NGRCalc 50,0,0 3,0,0,W" \FЯ" Output: SY2834011644 SY28341164 SY283116 SY2811 SY21 SY " f-Я" Command: NGRCalc 49,15,0 2,5,0,W" pFЯ" Output: XD9393627793 XD93932779 XD939277 XD9327 XD92 XD " z2Я" Command: NGRCalc 49,15,0 2,5,0,w -utm" └FЯ" Output: WV6674157944 WV66745794 WV667579 WV6657 WV65 WV " ▌Я ≤Я" NGRCalc " ╒<Я" Converts grid reference to latitude and longitude" ╛Я" Examples:" І"Я" Command: NGRCalc NZ9011" ю.Я" Output: 54,29,10.375 0.36,38.071,W" й"Я" Command: NGRCalc IN0050" т#Я" Output: 53,30,0 8,0,0,W" чА Х: Р+Т > BLib.Generic.ProgEnv 1.06 20Apr2020 Э:щєOS_GetEnv:ЙA$,A%,X%,Y%:X%=1:os%=((╨&FFF4)─&FF00)│256 NГ!░=&D7C1C7C5:run$=ARGV$(0):ЦA%=0╦ARGC:A$=A$+ARGV$(A%)+" ":М:=аA$,╘run$+1) MТ How do I detect SDL *before* filtering out DOS which has no INKEY-256 ? =Т How do I detect that SDL is actually running on Unix-y? $TГos%>31:Г░>&FFFFF:os%=32:чX%Й256:х≥"GetModuleFileName",0,X%,255:run$=$$X%:=@cmd$ .CA%=&600-&7B00*(░>&8000)-&3F00*(░>&C000):Г!(░-&100)=@%:A%=░-&300 8;Г?(╦P-3)=0:A%=&100:Г⌠<&FFFF:A%=░-&300:Г!(⌠+512)=@%:A%=⌠ BcA$=$A%:Г░=&8F00:run$=A$:х≥16╦A$,,A%:х≥72,"",A%:A$=аA$,1+їA$+" "," ",1+їA$," "))):Г╘A$=0:A$=run$ LUЦY%=-1╦0:A$=" "+A$:УA$=аA$,2):Щ≈A$<>32:Г≈A$=34:A%=їA$,"""",2)+1 ▀ A%=їA$+" "," ") V9ГY%:run$=аA$,1-(≈A$=34),A%-1+2*(≈A$=34)):A$=аA$,A%+1) ` М:=A$ j#щРos(A$):Г≈A$=42:ЪA$ ▀ Г╘A$:вA$ tА ~FщРexit(A%):Ъ"FX1,"+ц(A%─255):quit$=quit$:A$=quit$:quit$="":Рos(A$) ┬Г░>&FFFFF:х≤ A% ▀ Ю ▓А °: і#Т > BLib.CmdLine 1.10 27Jul2009 ╟cщєcl(l$,n%):Гl$="":A$=єs(A$):Г≈A$=34:A%=їA$+" "" ",""" ",2):l$=аA$,2,A%-2):A$=єs(аA$,A%+1)):=l$ ╨<Гl$="":A%=їA$+" "," "):l$=юA$,A%-1):A$=єs(аA$,A%+1)):=l$ дKГn%=0:Гl$<>"":A%=їA$,l$):ГA%:A$=єs(юA$,A%-1)+аA$,їA$," ",A%)+1))+" ":=╧ нГn%=0:Гl$<>"":=ё ь%A%=їю" ",≈l$=32)+A$,l$):ГA%=0:="" Б'A$=юA$,A%-1)+єs(аA$,їA$," ",A%)+1)) ЛCГ≈l$=32:l$=аA$,A%):A$=юA$,A%-1):=аl$,1-(≈l$=34),╘l$+2*(≈l$=34)) ЖPГ≈аA$,A%,1)<>34:l$=аA$,A%,їA$+" "," ",A%)-A%):A$=юA$,A%-1)+аA$,A%+╘l$+1):=l$ Jl$=аA$,A%+1,їA$+""" ",""" ",A%+1)-A%-1):A$=юA$,A%-1)+аA$,A%+╘l$+3):=l$ /щєs(A$):ГюA$,1)=" ":УA$=аA$,2):ЩюA$,1)<>" " +ГбA$,1)=" ":УA$=юA$,╘A$-1):ЩбA$,1)<>" " =A$ (: 2Т > TextIO <щРWin_TextIO F0х≥ "GetStdHandle",-10 ╦ @hfile%(1):*INPUT 13 P1х≥ "GetStdHandle",-11 ╦ @hfile%(2):*OUTPUT 14 Z$х≥ "SetConsoleMode",@hfile%(1),0 dА n: x!Т > BLib.Close 1.00 09Aug1998 ┌щРClose_All:*EXEC ▄"in%=in%:Гin%:A%=in%:in%=0:ы#A% √'out%=out%:Гout%:A%=out%:out%=0:ы#A% ═А ╙: ЄТ WINLIB5.BBC Version 1.5 ЎщР_setfocus(H%) хЙ K%,M%,O%,P% р ч P% Й 30 э [OPT 2 Ф*.K% push H% : call "SetFocus" : ret 16 П2.M% cmp dword [esp+8],&500 : jz K% : jmp [^O%] З] 'х≥ "GetWindowLong", @hwnd%, -4 ╦ O% &х≥ "SetWindowLong", @hwnd%, -4, M% (х≥ "SendMessage", @hwnd%, &500, 0, 0 "&х≥ "SetWindowLong", @hwnd%, -4, O% ,х≥ "Sleep", 0 6А Ъ