anchor.diff:            end if
   
             numberOfStarsInSample = stars_count
  -
  -          write(6,*) 'DEBUG: cone.f min l: ', min_gen_l
  -          write(6,*) 'DEBUG: cone.f max l: ', max_gen_l
  -          write(6,*) 'DEBUG: cone.f min b: ', min_gen_b
  -          write(6,*) 'DEBUG: cone.f max b: ', max_gen_b
         end subroutine
   
   
  
positive.code: 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
  
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)
  