/com	ANSYS 11.0
! Next100_pmt screen, sipm at y=0, triangular solid model 
/filname,el_sipm,1
vp=-11340 	! voltage at top ( arbitrary, to give 0.3 V/cm drift)
/com	vc=000
vg=-11250
va=0
v0=0
rw=.00015	!wire radius
xw=.00025	!wire location, x
y0=0.000	!wire locations, y
ywa=.002	!anode
ywg=.007	!gate
yp=.01		!buffer loc cathode

dy=0.0005
efy0=0		!initialization value
ymin=0

nx=20

nysa0=20
nya=50
nyw=nx
ny2=30

/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





/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,,y0
k,2,xw,y0
k,3,xw,y0,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


/com	VSBV, NV1, NV2, SEPO, KEEP1, KEEP2
VSBV,1,all,,delete,delete

asel,s,area,,16,19
aesize,all,.000008
alls
esiz,.00005
smrtsize,4
/com	mopt,vmesh,alte
vmesh,4

da,1,volt,v0
da,5,volt,vp
da,16,volt,va
da,17,volt,va
da,18,volt,vg
da,19,volt,vg

fini

/solu
anty,stat
solve
fini

/post1
/gropts,view,1
set,last

path,axx,2,6,400
ppath,1,,0,0.0
/com	ppath,2,,0,ywg-(0.5*(ywg-ywa))
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