echo off
% ...............рисуем карту в северной пол?рной проекции
map('pn',0,360,70,90);
pack
grid
%
load list_n        % список станций
load list_n.dat    % координаты станций (lat,lon)
list_n(:,2) = lon360(list_n(:,2));
%k=list_n(:,2)<0;
%list_n(k,2)=list_n(k,2)+360;
%
que = 'y';
while que=='y'
   % ............... графический ввод региона станций
   l_mm=[2;1];
   while l_mm(1)>=l_mm(2);     % контроль ввода по долготе
      [l_mm,f_mm]=ginput(2);
   end
   l_mm = r2d(l_mm);          % долгота в градусах
   l_mm = lon360(l_mm);
   f_mm = 90-f_mm;            % широта в градусах
   k_st = list_n(:,2)>=l_mm(1) & list_n(:,2)<=l_mm(2); % станции из региона
   c_color=input('Задайте код цвета и метки [напр. - c12o] ... ','s');
   % ............... перебор по станци?м
   list1=list(k_st,:);         % список выбранных станций
   list_n1 = list_n(k_st,:);   % и их географ. координат
   [m n]=size(list1);
   for i_s = 1:m
      tmp = ['load ' list1(i_s,:) '.dat']; % строка загрузки файла данных
      eval (tmp)
      lon = eval([list1(i_s,:) '(:,2)']);  % вектор долгот ВМП
      lat = eval([list1(i_s,:) '(:,3)']);  % вектор широт ВМП
      polar_n(lon,lat,c_color);
      eval(['clear ' list1(i_s,:)])
   end
   que = input ('Повторим ? [y/n] ... ','s');
end
