这个项目是仿真全钙钛矿太阳能电池IV特性的
代码如下：
go atlas

mesh
x.m loc=0.0 spac=0.25
x.m loc=2.0 spac=0.25
x.m loc=5.0 spac=0.25
y.m loc=0.0 spac=0.02
y.m loc=1.0 spac=0.02

region num=1 material=Silver x.min=0.0 x.max=5.0 y.max=0.2
region num=2 material=TiO2 x.min=0.0 x.max=5.0 y.min=0.2 y.max=0.3 donors=3.0e17
region num=3 user.material=MaPbI3 x.min=0.0 x.max=5.0 y.min=0.3 y.max=0.75 accept=1e13
region num=4 user.material=Spiro x.min=0.0 x.max=5.0 y.min=0.75 y.max=0.9 accept=3e17
region num=5 material=ITO x.min=0.0 x.max=5.0 y.min=0.9 y.max=1.0
region num=6 material=TiO2 x.min=0.0 x.max=5.0 y.min=0.8 y.max=0.9 donors=3e17
region num=7 user.material=MaPbI3 x.min=0.0 x.max=5.0 y.min=0.9 y.max=1.35 accept=1e13
region num=8 user.material=Spiro x.min=0.0 x.max=5.0 y.min=1.35 y.max=1.5 accept=3e17
region num=9 material=Silver x.min=0.0 x.max=5.0 y.min=1.5 y.max=1.7

electrode name=anode top
electrode name=cathode bottom

material material=Spiro user.group=semiconductor user.default=Organic
permittivity=3 eg300=3 affinity=2.45 nc300=2.5e20 nv300=2.5e20 mun=2.0e-4 mup=2.0e-4 taun0=1e-7 taup0=1e-7

material material=MaPbI3 user.group=semiconductor user.default=organic
permittivity=6.5 eg300=1.5 affinity=3.9 nc300=2.5e20 nv300=2.5e20 mun=50 mup=50 taun0=6e-8 taup0=6e-8
copt=3.27e-11 augn=0.88e-29 augp=0.88e-29

material material=TiO2 permittivity=9 eg300=3.2 affinity=4 nc300=2.2e18 nv300=1.8e19 mun=20 mup=10 taun0=6e-10 taup0=6e-10
#
models conmob fermi srh optr auger bgn
method itlimit=40 maxtraps=10
#
beam num=1 x.o=2.5 y.o=-0.5 angle=90
wavel.start=0.305 wavel.end=3.455 wavel.num=300
power.file=solarex03_15.spec out.power=solarex03_15.log
output opt.int

solve b1=0.1
log outf=solar_iv.log
save outfile =my_solar.str
tonyplot solar_iv.log

问题是：运行到solve语句的时候出现can not trap这个问题，但是我真不知道哪里下手修改。。