anchor.diff:       &                                      / thick_disk_sfr_param)
                         sfr_try = thick_disk_max_sfr * ran(iseed)
                         if (sfr_try <= time_try_sfr) then
  -                          starBirthTime(stars_count) = time_try
  +                          starBirthTime(stars_count) = (
  +     &                        time_try + thick_disk_birth_init_time)
                             exit
                         end if
                     end do
  +              else if (random_value > thick_disk_stars_fraction 
  +     &                 .and. random_value <= thick_disk_stars_fraction 
  +     &                                       + halo_stars_fraction) then
  +                  disk_belonging(stars_count) = 3
  +                  starBirthTime(stars_count) = halo_birth_init_time
  +     &                + halo_stars_formation_time * ran(iseed)
                 else
                     disk_belonging(stars_count) = 1
                     starBirthTime(stars_count) = thin_disk_birth_init_time 
  
positive.code:      &                total_generated_mass_in_bin + star_mass
                end if
  
                if (disk_belonging(stars_count) == 1) then
                    scale_height = THIN_DISK_SCALE_HEIGHT_KPC
                else if (disk_belonging(stars_count) == 2) then
  
negative.code:       mrep = mrep * massReductionFactor
        normal_mrep = mrep
        burst_mrep = mrep * BURST_FORMATION_FACTOR
  
        max_age = max(thin_disk_age, thick_disk_age, halo_age)
        burst_init_time = max_age - burst_age
        thin_disk_birth_init_time = max_age - thin_disk_age
        thick_disk_birth_init_time = max_age - thick_disk_age
        halo_birth_init_time = max_age - halo_age
  C     This can be easily proved by taking derivative from
  C     y = t * exp(-t / tau)
        thick_disk_max_sfr_relative_time = thick_disk_sfr_param
        thick_disk_max_sfr = (thick_disk_max_sfr_relative_time 
       &                      * exp(-thick_disk_max_sfr_relative_time 
       &                              / thick_disk_sfr_param))
  
        write(6,*) '      Mass reduction factor = ', massReductionFactor
  