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:           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
  
  
  
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)
  