pro plot_site_locations,LOCALE=locale

;   PSYM -  The following integer values of PSYM will create the
;        corresponding plot symbols
;        0 - circle
;        1 - downward arrow (upper limit), base of arrow begins at plot
;             value
;        2 - upward arrow (lower limt)
;        3 - 5 pointed star
;        4 - triangle
;        5 - upside down triangle
;        6 - left pointing arrow
;        7 - right pointing arrow
;        8 - square
;
; *** Add 2DVD locations to map
;
;  SN25  -92.464,  42.239
;  SN35  -92.366,  42.183
;  SN36  -92.282,  42.126
;  SN37  -92.072,  41.992
;  SN38  -91.874,  41.860
;  SN70  -91.42,   41.641
;
    TRUE=1 & FALSE=0
    if(NOT KEYWORD_SET(LOCALE)) then locale = 'WALLOPS'
    circle=0 & down_arrow=1 & up_arrow=2 & star=3 & triangle=4
    down_triangle=5 & left_arrow=6 & right_arrow=7
    color = 0
    sym_size=0.75

;    print,'Printing instrument locations for ' + LOCALE

    CASE LOCALE of 
        'WALLOPS': BEGIN
;
; *** Draw box for Pocomoke grid
; 
;            grid_lons = [-75.603D, -75.603D, -75.545D, -75.545D]
;            grid_lats = [ 38.101D,  38.055D,  38.055D,  38.101D]
            lon1 = -75.603D & lon2 = -75.545D
            lat1 =  38.055D & lat2 =  38.101D
            plots,[lon1,lon2],[lat2,lat2],thick=3,color=1
            plots,[lon2,lon2],[lat2,lat1],thick=3,color=1
            plots,[lon2,lon1],[lat1,lat2],thick=3,color=1
            plots,[lon1,lon1],[lat1,lat2],thick=3,color=1
;
; *** Plot disdrometer locations
; 
            Disd_ID = ['SN25','SN35','SN36','SN37','SN38','SN70']
            Disd_Lons = [ 37.944320,  37.944320,  37.934302,  37.937610, $
                          37.929408,  37.934320]
            Disd_Lats = [-75.480618, -75.481138, -75.471147, -75.456412, $
                         -75.473177, -75.471275]

            plotsym2, star, sym_size, /fill
            plots,disd_lons,disd_lats,psym=8,color=color
        END

        'IFLOODS': BEGIN
            disd_lons = [-92.464, -92.366, -92.282, -92.072, -91.874, -91.542]
            disd_lats = [ 42.239,  42.183,  42.126,  41.992,  41.860,  41.641]
            plotsym2, star, sym_size, /fill
            plots,disd_lons,disd_lats,psym=8,color=color
;
; *** Plot U. Iowa XPOL locations
;
            xpol_names = ['XPOL1','XPOL2','XPOL3','XPOL4']
            xpol_lons  = [ -91.858537, -91.409198, -91.734063,-91.623129]
            xpol_lats  = [  43.178575,  42.922222,  41.887048, 41.644025]
            n_xpol = n_elements(xpol_names)
            
            sym_size=0.75
            plotsym2, down_triangle, sym_size, /fill
            plots,xpol_lons,xpol_lats,psym=8,color=color
            for i=0,n_xpol-1 do begin
                xyouts,xpol_lons[i]-0.2,xpol_lats[i]+0.025,xpol_names[i], $
                  color=color,charsize=0.5
            endfor
;
; *** Plot WSR-88D Locations
;
            wsr88d_names = ['KDVN','KDMX']
            wsr88d_lons  = [-90.580848, -93.723000]
            wsr88d_lats  = [ 41.611539,  41.731000]
            n_88d = n_elements(wsr88d_names)
            
            sym_size=0.75
            plotsym2, circle, sym_size, /fill
            plots,wsr88d_lons,wsr88d_lats,psym=8,color=color
            
            for i=0,n_88d-1 do begin
                xyouts,wsr88d_lons[i]-0.2,wsr88d_lats[i]+0.025,wsr88d_names[i],$
                  color=color,charsize=0.75
            endfor
        END
    ENDCASE

    return
    end