input mp-geo-v0-59;
arborescence := "../../DATA/";
maillage:=true;
figure(-10u,-10u,10u,10.5u);
show xa;
Bonne(2,46.5,echelle(2,46.5,50));
Lecturevillesp("France");
draw fullcircle scaled(2*abs(bonnec(2931/60,139/60)-bonnec(50.6333,3.05))) shifted bonnec(2931/60,139/60);
draw bonnec(2931/60,139/60)--bonnec(50.6333,3.05);
%pour échelle;
message("distance dessin (cm) de la partie du meridien -central- concerne");
show abs(bonnec(46.5,2)-bonnec(46.5,7))/cm;%distance en cm sur le dessin
message("distance reelle (km) de la partie du meridien -central- concerne");
show pi*cosd(46.5)*6340/36;%longueur en km d'un arc de méridien.
message("Longueur reelle (km) pour 1 cm sur le dessin");
show (pi*cosd(46.5)*6340/36)/(abs(bonnec(46.5,2)-bonnec(46.5,7))/cm);%longueur en km pour 1 cm sur le dessin
%marque_p:="creux";
%pointe(bonnec(46.5,2),bonnec(46.5,7));
zoom(1);
message("sans zoom, distance dessin (cm) de la partie du meridien -central- concerne");
show 31.7*abs(bonnec(46.5,2)-bonnec(46.5,7))/cm;
message("zoom a choisir pour un échelle 1:5 000 000");
show (cm*pi*cosd(46.5)*6340/36)/(50*abs(bonnec(46.5,2)-bonnec(46.5,7)));
draw ((Xa,Ya)+u*(1,1))--((Xa,Ya)+u*(2,1));
label.top(btex 0 etex,(Xa,Ya)+u*(1,1));
label.top(TEX(""&decimal(Long)&"~km"),(Xa,Ya)+u*(2,1));
fin;
end