Т > 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Т dТ Len's original comments nТ ----------------------- xDТ 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. ,: 62A$=єOS_GetEnv:ver$="1.03":Г os%=32 ▄ *FLOAT 64 @Г A$<>"":Рcmd(A$):Рexit(0) J8Г os%=32:х≥ "ShowWindow",@hwnd%,3:Р_setfocus(@hwnd%) T: ^К 27+128 hРborder rН ┘ Р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"' &Я┼5)"Press a key:"; 0х▌ єanswer("ABCDE") й :и 1 Dlatlongrepeat=╧ Nх∙ latlongrepeat=╧ X ш:Я'' bРlatlong lн vи 2 ─bdrepeat=╧ ┼х∙ bdrepeat=╧ ■ ш:Я'' ·Р_bd ╗н ╡и 3 ╪ lltongr=╧ фх∙ lltongr=╧ п ш:Я'' з Р_lltongr Дн Ни 4 ЬРchoosegrid и 5 Я:quit%=╧:Н ┘ quit%=ё Гquit%:х≤ Ю *к 4А >: HщРchoosegrid RЯ''┼5)"SELECT GRID:"' \0Я┼5)"(A) OSGB - Great Britain & Isle of Man" fЯ┼5)"(B) ITM - Ireland" p$Я┼5)"(C) UTM - Channel Islands" zЯ'┼5)"Press a key:"; └х▌ єanswer("ABCD") й ▌и 1:Рconstants_OSGB ≤и 2:Рconstants_ITM ╒и 3:Рconstants_UTM ╛&и 4:Рconstants_OSI:Т Hidden option Ік юА й: т щРlatlong чРchoice2 Х ш:Я'' РЯ ┼10)"GRID REFERENCE: "G$' ЭР_Mvrho(E,N) Р_toll Р_C_EN $Рoutput(lambda,phi,C):Рlldisplay $<Я '┼10) "Do you want another Lat Long calculation? Y/N"; .&Г єanswer("YN")<>1 latlongrepeat=ё 8А B: L щР_bd VЯ ┼20)"FIRST POSITION"' `Рchoice2 j E1=E:N1=N tР_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 к (Г N2N1 ▄ Falpha=alpha+2*╞ Pм Zм d%Nm=(N1+N2)/2:q=N1-N2:Em=(E1+E2)/2 nГ y1%<0 ┌ y2%<0 ▄ xР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; "Я " Secs North"' ,>Я ┼15) "Longitude "; degE; " Degs "; minE; " Mins "; secE; 6Я " Secs "; ew$' @Р_ENtongr(E,N) JР_ngrdisplay TА ^: h щР_lltoEN r7Т 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) E=E0+P*IV+P^3*V+P^5*VI &А 0: :щР_ENtongr(E,N) DГgrid$="ITM":E=E+100000 NГ E<0:E=E+.1:Т IF E<0 E=E-1 XГ N<0:N=N+.1:Т IF N<0 N=N-1 bРen_string(E) l east$=en$ vР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$' Рcurtail(easting$) *easting$=en$ 4Рcurtail(northing$) >northing$=en$ H!NGR$=A$+B$+easting$+northing$ R5Я ┼15)"6-digit National Grid Reference is ";NGR$' \DЯ ┼15)"4-digit National Grid Reference is ";юNGR$,4);аNGR$,6,2)' f: p:Я '┼15)"Another lat/long to NGR calculation? Key Y/N"; z: └ Г є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 "; $Я "numeral form "' .Х ┼15) G$' 8Рngr Bи 2 L:Я ┼10) "ENTER FULL EASTING AND NORTHING, IN FORM E,N"' VFЯ ┼5 ) "(six or seven figures in each, plus decimals if need be)"' `Х ┼20) E,N' jG$=цE+","+цN tк ~А ┬: ▓-Т This puts NGR into easting and northing ° щРngr іГ (╘ G$ ─ 1)=1:G$="I"+G$ ╟$Г юG$,1)>"`":G$=Ґ(≈G$-32)+аG$,2) ╨Р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" ▄ %a%=≈аG$,2,1)─&DF:Гa%=≈"A":a%=≈"W" !N=N+(a%-32)*100000:E=E+500000 (л 2Р_div(X$,500000,2,4) <Р_div(Y$,100000,0,5) FГ юG$,1)="I" ▄ P(Г grid$="OSI":E=E-500000:N=N-1000000 Z'Г grid$="ITM":E=E-100000:N=N-500000 dм nм x: ┌А ▄: √@Т columns count from 0 to 4 from left, rows 1 to 5 from top. ═щР_div(Q$,m%,c%,r%) ╙0a%=(≈(Q$) ─ &DF)-≈("A")+5: Г a%>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 : "JТ This takes easting and northing, does iteration for M, returns phip, ,Т v, rho, y, and itasqd. 6: @щР_Mvrho(E,N) J y=E-E0 Tphip=(N-N0)/a+phi0 ^M=єM(b,n,phip,phi0) hУ rphin=(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 : щР_toll &VII=Ї(phip)/(2*rho*v) 0IVIII=(Ї(phip)/(24*rho*v^3))*(5+3*Ї(phip)^2+itasqd-9*Ї(phip)^2*itasqd) :>IX= (Ї(phip)/(720*rho*v^3))*(61+90*Ї(phip)^2+45*Ї(phip)^4) D IX=IX/v^2 N$phi=phip-y^2*VII+y^4*VIII-y^6*IX XX=1/(v*⌡(phip)) b.XI=1/(6*v^3*⌡(phip))*(v/rho+2*(Ї(phip))^2) l?XII=(1/(120*v^5*⌡(phip)))*(5+28*(Ї(phip))^2+24*(Ї(phip))^4) vUXIIA=(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)) А : -Т This derives deviation (C) from E and N щР_C_EN *XVI=Ї(phip)/v 48XVII=Ї(phip)/(3*v^3)*(1+Ї(phip)^2-itasqd-2*itasqd^2) >6XVIII=Ї(phip)/(15*v^5)*(2+5*Ї(phip)^2+3*Ї(phip)^4) HC=y*XVI-y^3*XVII+y^5*XVIII RА \: fщРoutput(lambda,phi,C) p lambda=²(lambda): phi=²(phi) zew$ = "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) WЯ ┼10) "True North is ";deg " Degs ";min " Mins ";sec;" Secs ";C$;" of grid North"' $(Я ┼10) "easting= ";E;" northing=";N' .А 8: BщР_bd_output LG=²(alpha) VРdegminsec(G) `Рnumplace(sec," ") jsec=num t=Я ┼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 А : щРerror (!Г ÷=17 ▄ А:Т REPORT:PRINT:END 23Т IF ERL<>1240 THEN REPORT:PRINT" at line ";ERL <Я FЯ ┼13)"Press any key"; P У Щ ╔ ZА d: n4Т This draws yellow border, with blue background xщР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 щРdegminsec(theta) "6ang=theta*3600+0.0005 : sec=(ang*1000)┐60000 /1000 ,+deg=ang │60 : min=deg ┐60 : deg=deg │60 6А @: J6Т removes last j$ sig figs.Pass space string to j$ TщРnumplace(Y,j$) ^ num$=ц(Y) h бnum$)=j$ rnum=╩(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) Nm1=(Nc+N1)/2: Nm2=(Nc+N2)/2 Nm=Nm1: y2%=0: q=N1-Nc &Р_converge 0convangle1=convangle :"Nm=Nm2:y1%=0:y2%=E2-E0:q=Nc-N2 DР_converge Nconvangle2=convangle X#convangle=convangle1+convangle2 bА l: vщР_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 Я"F= ";F Я"phip= ";phip *Я"v= ";v 4Я"rho= ";rho >Я"itasqd= ";itasqd HЯ"M= ";M RЯ"y= ";y \Я"VII= ";VII fЯ"VIII= ";VIII pЯ"IX= ";IX zЯ"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" ▄ )Рfirst_letter(east$,north$) :A$=Ґ(a%) $)Рsecond_letter(east$,north$):B$=Ґ(a%) .А 8м BГ grid$="OSI" ▄ L0Рfirst_letter(east$,north$) :A$=Ґ(a%):A$="I" V)Рsecond_letter(east$,north$):B$=Ґ(a%) `А jм tГ 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$) column=0:row=0 Рcolumn_row(E$) (Г E<0 ▄ 2column=init │ 5 + 2 <л Fcolumn=init │ 5 + 3 Pм ZРcolumn_row(N$) dГ N<0 ▄ n row=4 xл ┌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 Г бц╗N,2)<>"00":init=init-.1 "row= -init ┐ 5 ,л 6row=4-init ┐ 5 @м JРread_letter TА ^: hCТ "column_row" returns the initial one or two numbers, with "-" rТ 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) м Г M%=6 ▄ Рleft(EN$,1) &м 0м :А D: NBТ Preserves leading zeros, and repositions "-" when necessary. XщРen_string(EN) bEN$=ц(╗((EN+0.5)*10/10)) l L%=╘(EN$) vГ 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) л en$=б"000000"+EN$,6) м *м 4А >: HщРleft(x$,y) RГ ╩(x$)<0 ▄ y=2 \init=╩(юx$,y)) fА p: zщР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" А $: .щРconstants_OSI 88E0=200000: N0=250000 : Т grid coords of True origin B3phi0=╡(53.50) : Т latitude of True origin L=lambda0=╡(-8) : Т longitude of True origin (West -) V7F0=1.000035 : Т scale on central meridian `Da=6377563.395906615 : Т semi-major axis in metres scaled by F0 j9b=6356256.908205645 : Т semi-minor similarly scaled t/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" А : щРconstants_UTM (:E0=500000: N0=0 : Т grid coords of True origin 25phi0=╡(0) : Т latitude of True origin "@" ▄ G$=єcl("",0):comment$=A$ Рngr Р_Mvrho(E,N):Р_toll:Р_C_EN "Рoutput(lambda,phi,C) ,.Г comment$<>"":Я comment$;┴(16-╘comment$); 6"Я ;degN;",";minN;",";secN;" "; @.Я ;degE;",";minE;",";secE;",";юew$,1);" "; JЯ Tл ^$Г єcl("-osgb",0):Рconstants_OSGB h"Г єcl("-osi",0):Рconstants_OSI r"Г є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$); Я GL$;NGR$;" "; "Я GL$;юNGR$,4);аNGR$,6,4);" "; &"Я GL$;юNGR$,3);аNGR$,6,3);" "; 0"Я GL$;юNGR$,2);аNGR$,6,2);" "; :"Я GL$;юNGR$,1);аNGR$,6,1);" "; DЯ GL$;" "; NЯ ;east$;" ";north$;" "; XЯ bм l: 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" 2Я" -utm - Use UTM (Channel Islands) grid" <Я" Outputs space-seperated string of grid references" *Я" Examples:" 4,Я" Command: NGRCalc 50,0,0 3,0,0,W" >FЯ" Output: SY2834011644 SY28341164 SY283116 SY2811 SY21 SY " H-Я" Command: NGRCalc 49,15,0 2,5,0,W" RFЯ" Output: XD9393627793 XD93932779 XD939277 XD9327 XD92 XD " \2Я" Command: NGRCalc 49,15,0 2,5,0,w -utm" fFЯ" Output: WV6674157944 WV66745794 WV667579 WV6657 WV65 WV " pЯ zЯ" 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.ProgEnv 1.04 09Jan2007 ч:щєOS_GetEnv:ЙA$,A%:X%=1:os%=((╨&FFF4)─&FF00)│256:чX%-1 ХJГos%=32:Г░>&FFFF:х≥"GetModuleFileName",0,X%,255:A$=$$X%:run$=A$:=@cmd$ РГos%=32:A$=$&100 Э{Г╘A$=0:Г░>&7FFF:run$=$&8100:х≥"OS_GetEnv"╦A$,,A%:х≥"OS_WriteEnv","",A%:A$=аA$,1+їA$+" "," ",1+їA$," "))):Г╘A$=0:A$=run$ 2Г╘A$=0:Г?(╦P-3):A$=$&600 ▀ Г╘A$=0:A$=$(░-&300) 7A%=їA$+" "," "):run$=юA$,A%-1):Гrun$<>"":=аA$,A%+1) bY%=X%│256:A%=9:?X%=0:X%!1=X%+16:X%!16=0:ж&FFD1:A%=X%+16:Г!A%─?A%+A%?2<>8:A%?(A%+1)=13:=$(A%+1) $="" .&щРos(A$):Г≈A$=42:ЪA$ ▀ ГA$<>"":вA$ 8А B@щРexit(A%):Ъ"FX1,"+цA%:quit$=quit$:A$=quit$:quit$="":Рos(A$) LГos%=32:х≤ A% VГos%<6:Ю ▀ *Quit `А j: t#Т > BLib.CmdLine 1.00 09Aug1998 ~Hщєcl(l$,n%):Гl$="":A%=їA$+" "," "):l$=юA$,A%-1):A$=єs(аA$,A%+1)):=l$ ┬Г≈l$=32 ─ A$<>"":A$=" "+A$ ▓bA%=їA$,l$):l$="":ГA%>0─n%>0:l$=аA$,їA$+" "," ",A%+1)+1):ГаA$,A%,1)<>" ":l$=юl$,їl$+" "," ")-1) °qГA%:ГаA$,A%,1)=" ":A$=аA$,2+(≈A$<>32),A%-2-(A%=1)) ▀ ГA%:A$=юA$,A%-1)+аA$,їA$+" "," ",їA$+" "," ",A%)+╘l$)+1) і"A$=єs(A$):Гn%:=єs(l$) ▀ =A%<>0 ╟/щєs(A$):ГюA$,1)=" ":УA$=аA$,2):ЩюA$,1)<>" " ╨+ГбA$,1)=" ":УA$=юA$,╘A$-1):ЩбA$,1)<>" " д=A$ н: ьТ > TextIO БщРWin_TextIO Л0х≥ "GetStdHandle",-10 ╦ @hfile%(1):*INPUT 13 Ж1х≥ "GetStdHandle",-11 ╦ @hfile%(2):*OUTPUT 14 $х≥ "SetConsoleMode",@hfile%(1),0 А : !Т > BLib.Close 1.00 09Aug1998 (щРClose_All:*EXEC 2"in%=in%:Гin%:A%=in%:in%=0:ы#A% <'out%=out%:Гout%:A%=out%:out%=0:ы#A% FА P: ZТ WINLIB5.BBC Version 1.5 dщ Р_setfocus(H%) nЙ K%,M%,O%,P% x ч 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 эА Ъ