function ii=inlon(lonmin,lonmax,lon)
% gibt ii=1 zurck, falls lon zwischen lonmin und lonmax
% liegt, wobei lonmin die westliche grenze und lonmax
% die oestliche ist.
% Alle Werte koennen dauch gemischt pos/negativ sein.

% alle Werte auf den Bereich 0-360 abbilden
lonmin=mod(lonmin,360);
lonmax=mod(lonmax,360);
lon=mod(lon,360);

if lonmin>lonmax 
    ii= lon>=lonmin | lon<=lonmax;
else
    ii= lon>=lonmin & lon<=lonmax;
end

% alte Version
% ii=zeros(size(lon));
% if lonmin>lonmax & lonmax<0
%     lonmax=lonmax+360
% end
% 
% if lonmin<lonmax
%     ii1=lonmin<=lon & lonmax>=lon;
%     ii=ii | ii1;
%     ii1=lonmin<=lon-360 & lonmax>=lon-360;
%     ii=ii | ii1;
% else
%     warning('wrong longitude region')
% end
return

