Please help for SIlvaco PNP transistor model

[DZ8205]

My silvaco code as follow:

go atlas
TITLE Bipolar Gummel plot and  IC/VCE with constant IB

mesh
x.m l=0   spacing=0.05
x.m l=0.8 spacing=0.05
x.m l=1.0 spacing=0.05
x.m l=1.5 spacing=0.05
x.m l=5.0 spacing=0.05

y.m l=-0.2  spacing=0.05
y.m l=0.04  spacing=0.05
y.m l=0.06  spacing=0.05
y.m l=0.15  spacing=0.05
y.m l=0.30  spacing=0.05
y.m l=1.0   spacing=0.05

region number=1 x.min=0.0 x.max=5   y.min=0     y.max=1   material=silicon
region number=2 x.min=0.0 x.max=5   y.min=-0.2  y.max=0   material=oxide
region number=3 x.min=2.0 x.max=5   y.min=-0.1  y.max=0   material=SiGe x.composition=0.2
region number=4 x.min=0.0 x.max=0.5 y.min=-0.1  y.max=-0  material=poly

electrode name=emitter   x.min=2.5 x.max=2.5 y.min=-0.2 y.max=-0.1
electrode name=base      x.min=1.0 x.max=1.5 y.min=-0.2 y.max=0.0 
electrode name=collector x.min=0.0 x.max=0.5 y.min=-0.2 y.max=-0.1

doping uniform conc=1e18 n.type region=1
doping uniform conc=1e19 p.type region=3
doping uniform conc=1e19 p.type region=4

# set bipolar models
models conmob fldmob consrh auger print 
contact name=emitter n.poly surf.rec

solve init
save outf=bjtex04_0.str
tonyplot  bjtex04_0.str -set bjtex04_0.set

# Gummel plot

method  newton autonr trap
solve vcollector=-0.025
solve vcollector=-0.1
solve vcollector=-0.25 vstep=-0.25 vfinal=-2 name=collector

solve vbase=-0.025
solve vbase=-0.1
solve vbase=-0.2

log outf=bjtex04_0.log
solve vbase=-0.3 vstep=-0.05 vfinal=-1 name=base

tonyplot  bjtex04_0.log -set bjtex04_0_log.set

#IC/VCE with constant IB

#ramp Vb
log off
solve init
solve vbase=-0.025
solve vbase=-0.05

solve  vbase=-0.1 vstep=-0.1 vfinal=-0.7 name=base

# switch to current boundary conditions
contact name=base current
# ramp IB and save solutions
solve ibase=1.e-6
save outf=bjtex04_1.str master
solve ibase=2.e-6
save outf=bjtex04_2.str master
solve ibase=3.e-6
save outf=bjtex04_3.str master
solve ibase=4.e-6
save outf=bjtex04_4.str master
solve ibase=5.e-6
save outf=bjtex04_5.str master


# load in each initial guess file and ramp VCE
load inf=bjtex04_1.str master
log outf=bjtex04_1.log
solve vcollector=0.0 vstep=-0.25 vfinal=-5.0 name=collector

load inf=bjtex04_2.str master
log outf=bjtex04_2.log
solve vcollector=0.0 vstep=-0.25 vfinal=-5.0 name=collector

load inf=bjtex04_3.str master
log outf=bjtex04_3.log
solve vcollector=0.0 vstep=-0.25 vfinal=-5.0 name=collector

load inf=bjtex04_4.str master
log outf=bjtex04_4.log
solve vcollector=0.0 vstep=-0.25 vfinal=-5.0 name=collector

load inf=bjtex04_5.str master
log outf=bjtex04_5.log
solve vcollector=0.0 vstep=-0.25 vfinal=-5.0 name=collector

# plot results
tonyplot -overlay  bjtex04_1.log bjtex04_2.log bjtex04_3.log bjtex04_4.log bjtex04_5.log  -set bjtex04_1_log.set
quit

Please let me know any possible reasons why this PNP code can not work.

Many thanks.