src/figures/cantilever-sim/i-dep.asy

   1 import wtk_graph;
   2
   3 size(15cm,10cm,IgnoreAspect);
   4 scale(Linear, Linear);
   5 real xscale=1;
   6 real yscale=1e12;
   7
   8 graphFile("i-dep.d/i_dep_127_8_1e-6", xscale, yscale, phard, m8,
   9           t=math(units("127","pN/nm"))+", 8 domains");
  10 graphFile("i-dep.d/i_dep_27_8_1e-6", xscale, yscale, pmed, m8,
  11           t=math(units("27","pN/nm"))+", 8 domains");
  12 graphFile("i-dep.d/i_dep_127_30_1e-6", xscale, yscale, phard, m30,
  13           t=math(units("127","pN/nm"))+", 30 domains");
  14 graphFile("i-dep.d/i_dep_27_30_1e-6", xscale, yscale, pmed, m30,
  15           t=math(units("27","pN/nm"))+", 30 domains");
  16 graphFile("i-dep.d/i_dep_0.1_1_1e-6", xscale, yscale, psoft, m1,
  17           t=math(units("0.1","pN/nm"))+", 1 domain");
  18 graphFile("i-dep.d/i_dep_0.1_30_1e-6", xscale, yscale, psoft, m30,
  19           t=math(units("0.1","pN/nm"))+", 30 domains");
  20
  21 ylimits(140, 320);
  22
  23 label(sLabel("Peak index dependence ($"+units("1.0","$\mu$m/s")+"$)"),
  24       point(N), N);
  25 xaxis(sLabel("Peak index"), BottomTop, LeftTicks);
  26 yaxis(sLabel("Unfolding force (pN)"), LeftRight, RightTicks);
  27
  28 add(legend(), point(E), 20E, UnFill);