plot3D

geofile=fc_simesh.get_geo(3,3,'cylinderkey03');
  meshfile=fc_oogmsh.buildmesh3d(geofile,10);
  Th=fc_simesh.siMesh(meshfile);
  u=Th.eval(@(x,y,z) 3*x.^2-y.^3+z.^2+x.*y);
  figure(1)
  Th.plot(u);
  axis off;axis image;colorbar
  figure(2)
  Th.plot(u,'d',2,'labels',[10,11,31,1000,1020,1021,2000,2020,2021])
  hold on;axis off;axis image;colorbar
  Th.plot(u,'d',2,'labels',1,'FaceColor','interp', 'EdgeColor','none','FaceAlpha',0.4)
  figure(3)
  Th.plot(u,'d',2,'labels',1,'FaceColor','interp', 'EdgeColor','none','FaceAlpha',0.4)
  hold on;axis off;axis image;colorbar
  Th.plot(u,'d',2,'labels',[10,11,1000,2000])
  Th.plot(u,'d',2,'labels',[31,1020,1021,2020,2021], 'FaceColor','interp','EdgeColor','none')
  figure(4)
  Th.plot(u,'d',2,'labels',1,'FaceColor','interp', 'EdgeColor','none','FaceAlpha',0.4)
  hold on;axis off;axis image;colorbar
  Th.plot(u,'d',1,'LineWidth',2)
  fc_tools.graphics.monitors.autoGrid()

[fc-oogmsh] Input file : /home/cuvelier/Travail/Recherch/Matlab/fc-simesh/geodir/3d/cylinderkey03.geo
[fc-oogmsh] Mesh file <fc-oogmsh>/meshes/cylinderkey03-10.msh [version 2.2] already exists.
  -> Use "force" flag to rebuild if needed.
  Mesh cylinderkey03-10.msh is a 3-dimensional mesh
    Force dimension to 3