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