【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