anchor.diff:  C                 assuming uniform distribution
                     longitude = min_longitude + delta_longitude*ran(iseed)
                     latitude = min_latitude + DELTA_LATITUDE * ran(iseed)
  -
  -                  if (longitude < min_gen_l) then
  -                    min_gen_l = longitude
  -                  end if
  -                  if (longitude > max_gen_l) then
  -                    max_gen_l = longitude
  -                  end if
  -                  if (latitude < min_gen_b) then
  -                    min_gen_b = latitude
  -                  end if
  -                  if (latitude > max_gen_b) then
  -                    max_gen_b = latitude
  -                  end if
  -
  -                  if (stars_count < 6) then
  -                    write(6,*) 'DEBUG: random generated l: ', longitude
  -                  end if
  -                  if (stars_count < 6) then
  -                    write(6,*) 'DEBUG: random generated b: ', latitude
  -                  end if
         
                     inner_do2: do
                         ! Distance between Sun and random point in the cone
  
positive.code:      &                        ttdisk = 12.0,
       &                        ft,
       &                        fz,
       &                        opposite_triangle_side,
       &                        min_gen_l = 1E6,
       &                        max_gen_l = -(1E6),
       &                        min_gen_b = 1E6,
       &                        max_gen_b = -(1E6)
            integer :: stars_count
  
            common /tm/ starBirthTime,
  
negative.code:      &                     coordinate_Zcylindr
            common /patron/ heightPattern
            common /index/ flagOfWD,numberOfWDs,disk_belonging
  
            ! NOTE: this can be infinity
            delta_longitude = DELTA_LATITUDE / dcos(cone_height_latitude)
  