Приведены результаты моделирования запятненности 13 активных карликовых звезд спектральных классов G-K на основе многолетних фотометрических наблюдений. Для анализа использовались как результаты ведущихся в КрАО UBVRI наблюдений восьми звезд, так и данные литературы. С помощью усовершенствованной зональной модели запятненности звездных фотосфер, которая позволяет теперь рассматривать присутствие на звезде двух активных долгот, были переопределены параметры запятненности 13 избранных активных красных карликов типа BY Dra. Выполнен анализ изменения характеристик пятенной активности этих систем с течением времени с целью выявления возможной цикличности. У всех звезд, кроме OU Gem и BE Cet, отмечается довольно сильная зависимость изменения широты пятен от их площади с коэффициентом корреляции (по модулю) R( , S) от 0.38 до 0.92. У пяти звезд отмечается антикорреляция средней широты пятен и их площади (R( , S) от –0.24 до –0.73). Мы считаем это проявлением дрейфа пятен к экватору по мере их развития. Восемь звезд выборки показывают положительную корреляцию, то есть дрейф пятен к полюсу с ростом их площади. У 9 звезд наблюдаются циклы активности, проявляемые как в фотометрической переменности, так и в изменениях площадей и средних широт пятен. Для 5 звезд отмечаются периоды дрейфа пятен по широте, причем скорости широтного дрейфа пятен меньше по модулю аналогичного значения для солнечных пятен в 1.5–3 раза.