File {
  * input files:
  Grid=   "mesh/3D_real_msh.tdr"
  Doping=   "mesh/3D_real_msh.tdr"
  Parameter = "parameters/models.par"
  Current="mirrors_voltage2/m__FILE_.plt"
  Output= "logs/m__FILE_"
}

Electrode {
  { Name="n_contact"    Voltage=0.0 	Resist=0.00001	AreaFactor=840 } 
  { Name="p_contact" 	Voltage=1.35 	Resist=0.00001	AreaFactor=1000 } # 
}

Physics {

	AreaFactor=30

	Recombination(
		SRH 
			(
				DopingDependence 
				ExpTempDependence
				Tunneling                	# important for Reverce junction
			)
	)

	EffectiveIntrinsicDensity(BandGapNarrowing (OldSlotboom))

	Fermi

	Temperature = 300

	Mobility (DopingDep)
	HeteroInterfaces
	Thermionic

	Laser (

		Optics ( # use this way of defining optics for lasing action, not the way the line above!!!
			FEScalar (
				EquationType = Waveguide
				Symmetry = Symmetric	          		# or NonSymmetric or Periodic or Symmetric
																# if Symmetric is used, reflective boundary is at x=0
				Coordinates = Cartesian					# Specifies the coordinate system 
																# (can also be set to Cylindrical)
				LasingWavelength = 800            	# [nm]

				TargetEffectiveIndex = (3.4 3.34 3.5) 	
																# initial guess. Should be the same number of 
																# array elements as ModeNumber below

				Polarization = TE		            	# for scalar solver only (TE TM)
				Boundary = ("Type1") 					# ("Type2" "Type1") choose boundary conditions
				ModeNumber = 1               			# up to 10 modes
			)
		)


		ModeNumber = 1

		LongitudinalModes				# TransverseModes for edge emitter simulation; could be also LongitudinalModes
		CavityLength = _CavityLength_			# [micrometer]
		lFacetReflectivity = 0._REFLECTIVITY_L_	# Left facet power reflectivity
		rFacetReflectivity = 1.0				# Right facet power reflectivity # has little or NO effect on ModeGain
		OpticalLoss = 0.0				# [1/cm] # has little or NO effect on ModeGain and on SpontEmission
		WaveguideLoss					# activate feedback of loss from Optics. 
											# Without that OpticalLoss is however NOT neglected. Hence, what that
											# WaveguideLoss means really?

		Broadening = 0.02
		CosHyper

		qwTransport 

		qwExtension = 12e-9	# [m], AutoDetect - auto read of QW widths; This is however default
		qwScatmodel 			# carrier scattering is modeled, 
									# and only a fraction of the carriers enter the quantum well

		QWeScatTime = 8e-13   # [s]
		QWhScatTime = 4e-13   # [s]

		eQWMobility = 9200    # [cm^2/Vs]
		hQWMobility = 400     # [cm^2/Vs]

		RefractiveIndex (TemperatureDep)

	) # End of Laser section #####################################################

}

# xFraction is fraction of Al ###################################################
Physics (region="n-emitter") {
	MoleFraction(RegionName = "n-emitter" xFraction=0.50)
}
Physics (region="a_o") {
	MoleFraction(RegionName = "a_o" xFraction=0.08) 
	Active	# to indicate QW region
}
Physics (region="waveguide_1") {
	MoleFraction(RegionName = "waveguide_1" xFraction=0.33)
}
Physics (region="waveguide_2") {
	MoleFraction(RegionName = "waveguide_2" xFraction=0.33)
}
Physics (region="p-emitter") {
	MoleFraction(RegionName = "p-emitter" xFraction=0.50)
}
##################################################################################



Insert= "PlotSection_des.cmd"
Insert= "MathSection_des.cmd"


Solve {
	Quasistationary ( InitialStep=0.002
	MaxStep=0.01  Minstep=0.01
  	Goal {
		name="p_contact" voltage=1.9}) {	
		Coupled {  Electron Hole Poisson QWeScatter QWhScatter PhotonRate }
	} 	
}