[MOVED] TACD simulated solarcell error message

go atlas 
mesh auto

x.m loc=0.0 s=0.25
x.m loc=1.0 s=0.25

# Here we specify an auxiliary quantum tunneling mesh used to
# help model the tunneling interface between the diodes.
#
qtx.m loc=0.0 s=0.25
qtx.m loc=1.0 s=0.25
qty.m l=0.66    s=0.002
qty.m l=0.675   s=0.001
qty.m l=0.69    s=0.002
qty.m l=4.34    s=0.002
qty.m l=4.40    s=0.001
qty.m l=4.54    s=0.002

# The tandem cell is formed by three p-n photodiodes stacked
# on top of one another.  The three cells are separated by a GaAs tunnel junction. 

#Top Cell  
region material=InAlP   bot thick=0.0300  accept=9e19 ny=10  
region material=InGaP   bot thick=0.1000  accept=5e17 ny=10  
region material=InGaP   bot thick=0.5000  donors=1e17 ny=10  
region material=InAlGaP bot thick=0.0300  donors=1e18 ny=10  
#
#Tunnel Junction 1 t=0.05
region material=GaAs bot thick=0.10 donors=2e20 ny=100
region material=GaAs bot thick=0.10 accept=1e20 ny=100 
#
# Middle Cell  t=2.84
region material=InGaAs   bot thick=0.100  accept=5e18   ny=10 
region material=InGaAs   bot thick=3.450  donors=1.2e17 ny=10 
region material=GaInP    bot thick=0.100  donors=1e18   ny=10
#
#Tunnel Junction 2 t=0.05
region material=GaAs   bot thick=0.10 donors=1e20 ny=10 
region material=GaAs   bot thick=0.10 accept=6e19 ny=10
#
#Bottom Cell  t=2.84
region material=InGaAs  bot thick=0.10   accept=1.5e18 ny=10   
region material=Ge      bot thick=0.05   accept=1.5e18 ny=10  
region material=Ge      bot thick=180.00 donors=5e18   ny=10  
# Electrodes on top and bottom
electrode name=anode   top
electrode name=cathode bottom
#
Tonyplot 
# The complex index data will mostly use the Sopra database.
# The index data for InGaP is supplied in a file.
#
material mat=GaAs             sopra=Gaas.nk
material mat=InAlGaP          sopra=Againp.all
material mat=AlGaAs           sopra=Algaas.all
material material=InGaP      sopra=againp0.nk            
material material=InAlP       sopra=Againp10.nk
material material=InGaP       sopra=againp0.nk                          
material material=InGaAs      sopra=InGaAs.nk
material material=Ge          sopra=ge.nk

#
# Here we define the affinity for InGaP and AlGaAs to
# set up the band alignment
#
material material=InGaP   affinity=4.08
material material=AlGaAs  affinity=3.54
#
# Enable physical models 
#
models srh fermi conmob optr auger bgn
#
 method itlimit=25 maxtraps=10
#

#
  beam num=1 x.o=0.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         
  
  tonyplot input solar spectrum
#
tonyplot solarex03_15.log -set solarex03_15.set
#
# Here is an alternate simplified am0 solar spectrum
#
     power.file=solarex03_15.spec out.power=solarex03_15.log
#
# Output optical intensity.
#
output band.temp opt.int
#
# save the initial structure
#
   save outf=solarex03_15.str


solve vanode=0.01
solve vanode=0.5
solve name=anode vstep=0.1 vfinal=2.7
#
# Set NLBBT .
#
models srh fermi ni.fermi conmob optr auger bgn bbt.nonlocal bbt.nlderivs  
#
solve b1=0.1
solve b1=1
#
log outf=solarex03_15.log
#
# Ramp the bias back.
#
solve name=anode vstep=-0.01   vfinal=2.2
solve name=anode vstep=-0.1    vfinal=0
#
# Extract important figures of merit.
#
extract init infile="solarex03_1.log"
extract name="Jsc" max(curve(v."anode", i."cathode"))
extract name="JscmAcm2" $Jsc*1e08*1e03
extract name="Voc" x.val from curve(v."anode", i."cathode") where y.val=0.0 
extract name="Pm" max(curve(v."anode", (v."anode" * i."cathode"))) 
extract name="Vm" x.val from curve(v."anode", (v."anode"*i."cathode") ) \
	where y.val=$"Pm"
extract name="Im" $"Pm"/$"Vm"
extract name="FF" ($"Pm"/($"Jsc"*$"Voc"))*100
extract name="Opt_int" max(beam."1")
extract name="Eff" (1e8*$Pm/$Opt_int)*100
#
# Plot the results
#
tonyplot solarex03_15.log -set solarex03_15.set

quit