pro plot_map, win, minLat, maxLat, minLon, maxLon
wset,win
; Tzipi Sidel May 2004
latCenter = 0.5*(minLat+maxLat)
lonCenter = 0.5*(minLon+maxLon)
if (minLon gt 0 and maxLon lt 0) then begin
        lonCenter  = minLon+0.5*((180-minLon)+(180+maxLon))
        if (lonCenter gt 180) then lonCenter = maxLon+(lonCenter -180)
endif
if (minLon gt 0 and maxLon gt 0 and maxLon lt minLon) then begin
        lonCenter = 0.5*(minLon+maxLon)-180.0
endif
difLat=abs(minLat-maxLat)
if (abs(minLat) gt 50 and abs(maxLat) gt 50 and difLat le 40) then begin
    projection='lambert' 
endif else begin
    projection='cylindrical'
endelse
green=getColor('gray',1)
blue=getColor('blue',2)
red=getColor('red',3)
backgroundcolor= getcolor('blue', 4)
black=getColor('black',0)
white=getColor('white',255)
!P.Background = backgroundcolor
case projection of
   'lambert': begin

      map_set, $
         label=2, /grid, /continent, /hires, /lambert, /isotropic, $
         latCenter, lonCenter, $
         limit=[minLat,minLon,maxLat,maxLon], $
         e_continents={fill:1, color:green}
   end
   
   'cylindrical': begin
      map_set, $
         label=2, /grid, /continent, /hires, /isotropic, /cylindrical, $
         latCenter, lonCenter, $
         limit=[minLat,minLon,maxLat,maxLon], $
         e_continents={fill:1, color:green}
   end
   else:
endcase
end