/com	ANSYS 11.0
! Next100_pmt screen, sipm at y=0, triangular solid model 
vp=000 	! voltage at  PMT screen
vc=000
vg=12000
va=0
v0=0
rw=.000015	!wire radius
xw=.00025	!wire location, x
y0=0.000	!wire locations, y
ywa=.002	!anode
ywg=.007	!gate
ywc=.08		!cathode	
yp=.13		!buffer loc cathode

dy=0.0005

efy0=0
ymin=0

nx=20

nysa0=140
nya=50
nyw=nx
ny2=100

/PNUM,KP,1  
/PNUM,LINE,1
/PNUM,AREA,1
/PNUM,VOLU,0
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0 
/psym,cs,1
/PSF,PRES,NORM,2,0,1
/SHRINK,0   
/ESHAPE,1   
/EFACET,1   
/RATIO,1,1,1


y=y0


/prep7
et,1,122	!3d 20-node  estat solid
msha,1,3d
mp,perx,1,1.
mp,perx,2,1.
mp,perx,3,1.

/com	 
k,1,,y
k,2,xw,y
k,3,xw,y,xw
k,4,,yp
l,1,2 
l,2,3 
l,3,1 
al,1,2,3
l,1,4
vdrag,1,,,,,,4

/com	CYL4, XCENTER, YCENTER, RAD1, THETA1, RAD2, THETA2, DEPTH
CYL4,xw,ywa,rw,,,,xw
CYL4,xw,ywg,rw,,,,xw
CYL4,xw,ywc,rw,,,,xw
/com	VSBV, NV1, NV2, SEPO, KEEP1, KEEP2
VSBV,1,all,,delete,delete

asel,s,area,,20,25
aesize,all,.000003
alls
esiz,.00005
smrtsize,4
/com	mopt,vmesh,alte
vmesh,5

da,1,volt,0
da,5,volt,vp
da,20,volt,va
da,21,volt,va
da,22,volt,vg
da,23,volt,vg
da,24,volt,vc
da,25,volt,vc

fini

/solu
anty,stat
solve
fini

/post1
/gropts,view,1
set,last

path,axx,2,6,2000
ppath,1,,0,0.0
ppath,2,,0,yp
pdef,efy,ef,y
plpath,efy
/eof

*GET,ymin,NODE,1,loc,y
*GET,efy0,NODE,1,EF,y
*cfopen,el_buf_3d.out,txt,,append
*vwrite,ymin,efy0
(f10.3,f10.7)
*CFCLOS
fini
parsav,all,pars,txt
/cle
parres,new,pars,txt
*enddo
/eof