File: dim.m

package info (click to toggle)
abinit 9.10.4-3
  • links: PTS, VCS
  • area: main
  • in suites: sid
  • size: 518,712 kB
  • sloc: xml: 877,568; f90: 577,240; python: 80,760; perl: 7,019; ansic: 4,585; sh: 1,925; javascript: 601; fortran: 557; cpp: 454; objc: 323; makefile: 77; csh: 42; pascal: 31
file content (17 lines) | stat: -rw-r--r-- 399 bytes parent folder | download | duplicates (7) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
load file.out
ngx=60
ngy=50
ngz=50
isodens=0.005
x=[1:ngx]
y=[1:ngy]
z=[1:ngz]
v=file(:);
v1=reshape(v,ngx,ngy,ngz);
p=patch(isosurface(y,x,z,v1,isodens));
axis([1 ngy 1 ngx 1 ngz])
set(p,'FaceColor','blue','EdgeColor','none');
light('Position',[0 1 0],'Style','infinite')
light('Position',[1 0 0],'Style','infinite')
light('Position',[0 0 -1],'Style','infinite')
set(gca,'DataAspectRatio',[1 1 1])