[ncl-talk] plot gc_latlon

[Adam Herrington]

In regards to the only transect example (
https://www.ncl.ucar.edu/Applications/transect.shtml), is there any reason
we should doubht that the transect computed from gc_latlon:

  dist     = gc_latlon(leftlat,leftlon,rightlat,rightlon,npts,2)

is accurately depicted by:

  map = gsn_csm_map_ce(wks,mres)         ; create map

  gsn_polyline(wks,map,(/leftlon,rightlon/),(/leftlat,rightlat/),pres)

? I am not sure that gc_latlon and gsn_polyline (as portrayed on
gsn_csm_map_ce, e.g. a global gridded dataset) use the same procedure to
compute a transect, given identical input (rightlon, leftlon, etc...)

If the two procedures are not consistent, does anyone know a way to
(a) actually plot the transect from gc_latlon on a map (directly plot
variable 'trans', in the example?), or
(b) compute a transect the same way that gsn_polyline and gsn_csm_map_ce do.

Thanks in advance!

Adam
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160529/1f49efc2/attachment.html