[ncl-talk] Convert polar streographic coordinates to lat/lon grid

Dennis Shea shea at ucar.edu
Mon Sep 12 07:50:41 MDT 2016


You have sent the same question before, The reason for no reply is that the
question is rather vague.

[1] You can read the variable and use the two-dimensional coordinates
latitude(y,x) ; longitude (y,x) to plot vis:


    f = addfile("foo.nc", "r")
    x = f->X                          ; :(time,lev,y,x)
    x at lat2d = f->latitude      ; (y,x)
    x at lon2d = f->longitude

  wks = gsn_open_wks("png","polar")              ; send graphics to PNG file
  res            = True                          ; plot mods desired
  res at gsnPolar   = "NH"                          ; specify the hemisphere

  nt  = 0
  lev = 10
  plot = gsn_csm_contour_map_polar(wks,x(nt,klev,:,:),res)

or, maybe

  plev = 500   ; natural coordinate
  plot = gsn_csm_contour_map_polar(wks,x(nt,{plev},:,:),res)

===

[2]
If you want the polar data regridded then see:
    http://www.ncl.ucar.edu/Applications/ESMF.shtml

Your source grid would be classified as 'curvilinear'.

Your resulting file would ve 'rectilinear'

Good luck





On Mon, Sep 12, 2016 at 4:00 AM, Sourav Chatterjee <srvsxc at outlook.com>
wrote:

> Dear All
>
>
> How to convert the stereographic projection x-y coordinates to lat lon
> grids? I have the 2d lat lon coordinates but the variable is in x-y
> coordinates e.g
>
>
> temp(time,depth,y,x).
>
>
> and latitude(y,x) ; longitude (y,x)
>
>
> Is it poosible to get temp(time,depth,latitude,longitude) from this
> information?
>
>
> Thanks
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160912/f6ccaf3b/attachment.html 


More information about the ncl-talk mailing list