$RandomizedR

$RandomizedR is a special command to setup the random alloy structures. The format is

$RandomizedR
Region_Number_used Random_type location_N_used seeding_num
location_1 location_2 ...... location_N
Region_1 method P1 P2 P3 P4 P5  .... 
Region_2 method P1 P2 P3 P4 P5  .... 
...
Region_N method  P1 P2 P3 P4 P5  .... 

Method=0 fully random number
Method=1 P1=dx(um), P2=average_composition, P3=sigma, P4=indium_low P5=delta- In%
Method=2 P1=dx(um), P2=average_composition, P3=sigma, P4=indium_low P5=delta- In%
Method=3 P1=dx(um), P2=average_composition, P3=sigma, P4=indium_low P5=delta- In%

If method ==11 12 13, the format would be

$RandomizedR
Region_Number_used Random_type location_N_used seeding_num
location_1 location_2 ...... location_N
Region_1 method Region_ini  Region_fin P1 P2 P3 P4 P5  .... 
Region_2 method Region_ini  Region_fin P1 P2 P3 P4 P5  .... 
...
Region_N method Region_ini  Region_fin  P1 P2 P3 P4 P5  .... 

Method=11 use_region_min use_region_max  P1=dx(um), P2=average_composition, P3=sigma, P4=indium_low P5=delta- In%
Method=12 use_region_min use_region_max  P1=dx(um), P2=average_composition, P3=sigma, P4=indium_low P5=delta- In%
Method=13 use_region_min use_region_max  P1=dx(um), P2=average_composition, P3=sigma, P4=indium_low P5=delta- In%

The difference between 1(11) and 2(12) is that Method=2 (12) will not replace the trap-related coefficients, nonradiative lifetime, doping, impurity, and generation of the original region. 
Method==2(12)  use the original trap-related coefficients, nonradiative lifetime, doping, impurity, and generation of the original region
Method==3(13)  use the original polarization, trap-related coefficients, nonradiative lifetime, doping, impurity, and generation of the original region

The idea is to generate virtual regions. The region number is "location_N_used" 
The second line defines the ID of virtual regions (location_N_used). For example location_N_used=10
13 14 15 16 17 18 19 20 21 22 23    means these regions are virtual region
P1 is the lattice size. typically 2.83e-4 um. 
P2=average_com of this QW . For example 0.15 is 15% average indium composition
P3=sigma  -> broaden factor of averaging local indium composition. Typically 2-4 lattice size. Unit is um. 
P4=indium_low  define indium composition of virtual region 1. 
P5=delta-In%   defines the composition difference of each virtual region. It should be linear.

Note that no matter which method is used, the mobility, B, and C will use the original region setting.

Example:

$RandomizedR
6 1 20 -10       #total 20 regions as virtual region
13 14 15 16 17 18 19 20 21 22 23 .... 32     # The region id of 20 virtual region 
3 0                                       （QW location)
4 0
5 0
6 0
7  1  2.83e-4  0.15  9.8e-4 0.01 0.02
8  2  2.83e-4  0.15  9.8e-4 0.01 0.02         # use the original lifetime, doping, , impurity, traps, and generation of the original region.
9  3  2.83e-4  0.15  9.8e-4 0.01 0.02         # use the original polarization, lifetime, doping, , impurity, traps, and generation of the original region.
10 11  1  10   2.83e-4  0.15  9.8e-4 0.01 0.02   # use virtual region 1-10 as 1% to 19%  if composition smaller than 1%, it would be 1%. If > 19% it would be 19%
10 11  11 20   2.83e-4  0.40  9.8e-4 0.21 0.02    # use virtual region 11-20 as 21% to 39%  if composition smaller than 21%, it would be 21%. If > 39% it would be 39%