doc/figures/Terrefusee.mp

   1 input mp-geo
   2 input mp-solid;
   3
   4 vardef mappemonde(expr longobs,latobs)=
   5   projection:="non";
   6   Initialisation(5,longobs,latobs,400);
   7   numeric phim,phip,phii;%phi moins -- phi plus - phi intermédiaire
   8   phim=Phi+arcsind(rayon/Rho)-90;
   9   phip=Phi+90-arcsind(rayon/Rho);
  10   color pte[];
  11   pte1=rayon*(cosd(phim)*cosd(Theta),cosd(phim)*sind(Theta),sind(phim));
  12   pte2=rayon*(cosd(phip)*cosd(Theta),cosd(phip)*sind(Theta),sind(phip));
  13   pte3=1/2[pte1,pte2];
  14   pte4-pte3=Normal((0,0,0),pte1,pte2);
  15   if (Phi>90):
  16     phip:=180-phip;
  17     phii:=180-phim;
  18     phim:=phip;
  19     phip:=phii;
  20   fi;
  21   if (Phi<-90):
  22     phip:=-180-phip;
  23     phii:=-180-phim;
  24     phim:=phip;
  25     phip:=phii;
  26   fi;
  27   fill cercles(pte3,pte1,pte3,pte1,pte4) withcolor ciel;
  28   Lecture("Cameriquesud.dat");
  29   Lecture("Ccaraibes.dat");
  30   Lecture("Cameriquecentrale.dat");
  31   Lecture("Cameriquenord.dat");
  32   Lecture("Casie.dat");
  33   Lecture("Ceurope.dat");
  34   Lecture("Cafrique.dat");
  35   if volcans=true:
  36     Lecturevolcans;
  37   fi;
  38   if lacs=true:
  39     Lecturelacs;
  40     Lecturelacssup;
  41   fi;
  42   Lectureiles;
  43   if capitales=true:
  44     Lecturecapitales;
  45   fi;
  46   if fleuves=true:
  47     Lecturerivieres;
  48     %Lecturerivieressup;
  49   fi;
  50   if maillage=true:
  51     drawoptions(withcolor gris);
  52     MaillageSphere;
  53     drawoptions();
  54   fi;
  55   if maille=true:
  56     Maille;
  57   fi;
  58   draw cercles(pte3,pte1,pte3,pte1,pte4);
  59 enddef;
  60
  61 vardef Transform(expr PT)=RotY(RotX(RotZ(PT)))
  62 enddef;
  63
  64 transformation:=true;
  65 %2'
  66 figureespace(-7u,-7u,8u,7u);
  67 picture terre;
  68 fleuves:=false;
  69 lacs:=false;
  70 capitales:=false;
  71 terre=image(
  72   mappemonde(80,20);
  73   );
  74 fill feuillet;
  75 for k=1 upto 100:
  76   fill fullcircle scaled (2*uniformdeviate(1)*mm) shifted((uniformdeviate(1))[coinbg,coinbd]+uniformdeviate(1)*(coinhg-coinbg)) withcolor jaune;
  77 endfor;
  78 draw terre;
  79
  80 Initialisation(500,160,20,50);
  81 unit:=0.2;
  82 traits:=false;
  83 angy:=-60;
  84 outcolor:=rouge;
  85 path rocket[];
  86 rocket1=(0,4.6)--(-0.8,4.6)--(-0.85,5.6)--(-0.95,6.6)--(-1.05,7.6)--(-1.1,8.6)--(-1.2,9.6);
  87 rocket2=(-1.2,9.6)--(-1.25,10.6)--(-1.3,11.6)--(-1.35,12.6)--(-1.4,13.6)--(-1.4,14.6)--(-1.375,15.6)--(-1.27,16.6)--(-1.2,17.6)--(-1.05,18.6);
  88 rocket3=(-1.05,18.6)--(-0.85,19.6)--(-0.65,20.6)--(-0.35,21.6)--(-0.1,22.2)--(-0.1,23.75)--(0,23.75);
  89 subh:=24;
  90 perso1:=true;
  91 Ferme1:=false;
  92 string couleurperso;
  93 couleurperso="if ((tapj mod 48)=0) or ((tapj mod 48)=1) or ((tapj mod 48)=2) or ((tapj mod 48)=6) or ((tapj mod 48)=7) or ((tapj mod 48)=8) or ((tapj mod 48)=12) or ((tapj mod 48)=13) or ((tapj mod 48)=14) or ((tapj mod 48)=18) or ((tapj mod 48)=19) or ((tapj mod 48)=20) or ((tapj mod 48)=27) or ((tapj mod 48)=28) or ((tapj mod 48)=29) or ((tapj mod 48)=33) or ((tapj mod 48)=34) or ((tapj mod 48)=35) or ((tapj mod 48)=39) or ((tapj mod 48)=40) or ((tapj mod 48)=41) or ((tapj mod 48)=45) or ((tapj mod 48)=46) or ((tapj mod 48)=47):red else: white fi;";
  94 TR:=(0,-15,0);
  95 ObjetAnneau1("nbp=10",rocket2);
  96 outcolor:=rouge;
  97 for k=2 upto 6:
  98   Ferme[k]:=true;
  99 endfor;
 100 ObjetAnneau2("nbp=7",rocket1);
 101 ObjetAnneau3("nbp=7",rocket3);
 102 path amortisseur;
 103 amortisseur=for k=-90 step 10 until 0:(cosd(k),sind(k)+1)--endfor for k=0 step 2 until 22:(cosd(k)*11.78-10.78,11.78*sind(k)+1)-- endfor (0,5.75);
 104 TR:=(0,-15,0)+(0,-5,0);
 105 ObjetAnneau4("nbp=22",amortisseur);
 106 angz:=150;
 107 TR:=(0,0,0);
 108 TR:=(0,-15,0)+Image((0,-5,0));
 109 ObjetAnneau5("nbp=22",amortisseur);
 110 angz:=-90;
 111 TR:=(0,0,0);
 112 TR:=(0,-15,0)+Image((0,-5,0));
 113 ObjetAnneau6("nbp=22",amortisseur);
 114 TR:=(0,-15,0)+(0,0,0);
 115 angz:=0;
 116 subh:=1;
 117 ObjetBiface7((0,-0.8,4.6),(0,-0.85,5.6),(0,-0.95,6.6),(0,-1.05,7.6),(0,-1.1,8.6),(0,-1.2,9.6),(0,-2,9.2),(0,-3,8.3),(0,-4,7.2),(0,-5,5.75),(0,-4.5765,4.6402),(0,-4.35,3.8498),(0,-4.1789,3.0456),(0,-4.0645,2.2313),(0,-4.0072,1.4111),(0,-4,1),(0,-3.2,2.2),(0,-2.5,3),(0,-1.6,4));
 118 angz:=150;
 119 ObjetBiface8((0,-0.8,4.6),(0,-0.85,5.6),(0,-0.95,6.6),(0,-1.05,7.6),(0,-1.1,8.6),(0,-1.2,9.6),(0,-2,9.2),(0,-3,8.3),(0,-4,7.2),(0,-5,5.75),(0,-4.5765,4.6402),(0,-4.35,3.8498),(0,-4.1789,3.0456),(0,-4.0645,2.2313),(0,-4.0072,1.4111),(0,-4,1),(0,-3.2,2.2),(0,-2.5,3),(0,-1.6,4));
 120 angz:=-90;
 121 ObjetBiface9((0,-0.8,4.6),(0,-0.85,5.6),(0,-0.95,6.6),(0,-1.05,7.6),(0,-1.1,8.6),(0,-1.2,9.6),(0,-2,9.2),(0,-3,8.3),(0,-4,7.2),(0,-5,5.75),(0,-4.5765,4.6402),(0,-4.35,3.8498),(0,-4.1789,3.0456),(0,-4.0645,2.2313),(0,-4.0072,1.4111),(0,-4,1),(0,-3.2,2.2),(0,-2.5,3),(0,-1.6,4));
 122 nbobj:=9;
 123 DessineFusion;
 124 finespace;
 125 end