;; Plot box on map

;;DIAG_NORCPM; LONLATBOX: 120,130,-5,5
;;DIAG_NORCPM; OUTFILE: fig
;;DIAG_NORCPM; PLOTRES: ;; no other figure set

begin
    lonlatbox = (/-150,-90,-5,5/) 
    lonbe = (/lonlatbox(0),lonlatbox(1)/)
    latbe = (/lonlatbox(2),lonlatbox(3)/)
    outfile = "box_nino3"
    xpoints = (/lonbe(0), lonbe(1), lonbe(1), lonbe(0), lonbe(0)/)
    ypoints = (/latbe(0), latbe(0), latbe(1), latbe(1), latbe(0)/)
    
    res = True
        res@mpMaxLatF =  90.
        res@mpMinLatF = -90.
        res@mpCenterLonF = 210.
        res@mpLandFillColor = "darkgreen"

        res@gsnDraw  = False
        res@gsnFrame = False

    polyres = True
        polyres@gsFillColor = "red"

    

    wks = gsn_open_wks("ps",outfile)
    plot = gsn_csm_map(wks,res)
    gsn_polygon(wks,plot,xpoints,ypoints,polyres)

    if ( .not.res@gsnDraw )then
        draw(plot)
    end if 
    if ( .not.res@gsnFrame )then
        frame(wks)
    end if 
end