function[]=line_style(axis_handle,lev,colors)
%disp('_____________________________________________________________________')
%disp('Usage: line_style(axis_handle,lev,colors)')
%disp('Progam changes linewidth to 1 at z=lev, linestyle to ''--'' at z<lev')
%disp('arguments are optional, default: axis_handle=gca, lev=0,colors=''kkk''')
%disp('colors: 1-below lev,2-lev,3-above lev')
%disp('_____________________________________________________________________')
%disp('senya@atmos.umd.edu, v0.1 11/15/01, v0.2 12/4/02')
%disp('_____________________________________________________________________')

if(nargin==0)
    axis_handle=gca;
    lev=0;
    colors='kkk';
end
if(nargin==1)
    lev=0;
    colors='kkk'
end
%switch axes layer to <top>
set(axis_handle,'layer','top');


%get handles to patch (or line) objects in current axes
hh = findobj('type','line','parent',axis_handle);
color='color';
if (length(hh)==0)
 hh = findobj('type','patch','parent',axis_handle);
 color='EdgeColor';
end   

count =0;
for i = 1:length(hh);
   if (length((get(hh(i),'UserData'))) == 1 & ~isnan(get(hh(i),'UserData')))
      count = count + 1;
%      h(count) = hh(i);
      cdata(count) = get(hh(i),'UserData');
   end
end
%contour_levels = unique(cdata);
ind=find(cdata==lev);
set(hh(ind),'linewidth',1,color,colors(2),'linestyle','-')
ind=find(cdata>lev);
set(hh(ind),'linewidth',0.5,color,colors(3),'linestyle','-')
ind=find(cdata<lev);
set(hh(ind),'linewidth',0.5,color,colors(1),'linestyle','--')
