【NCL】获取plot图经纬度的页面坐标-datatondc

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
begin
        infile = addfile("fnl_20101013_00_00_c.nc", "r")
        var = infile->tmp(12,:,:)
;       printVarSummary(var)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
        wks = gsn_open_wks("x11","name")
        res = True
        res@gsnFrame = False
        plot = gsn_csm_contour_map_ce(wks, var, res)
;
        lons = (/120, -120/)
        lats = (/10, -20/)
        x_out = new(dimsizes(lons),float)
        y_out = new(dimsizes(lats),float)
        datatondc(plot,lons,lats,x_out,y_out)
;       print(x_out)
;       print(y_out)
;
        txres = True
        txres@txFontHeightF = 0.015
        txres@txJust = "TopLeft"
        gsn_text_ndc(wks,"Hello",x_out(0),y_out(0),txres)
        gsn_text_ndc(wks,"World !",x_out(1),y_out(1),txres)
        frame(wks)
end