Re: ce map to polar stereographic

There are numerous examples at
    http://www.ncl.ucar.edu/Applications/polar.shtml

Don't know why it won't plot ... this is done commonly

latLimit = 65.0

  wks = gsn_open_wks("ps","fisel") ; open ps file
   gsn_define_colormap(wks,"WhBlGrYeRe") ; choose colormap
  ;gray = NhlNewColor(wks,0.7,0.7,0.7) ; add gray to colormap

   res = True
   res@gsnPolar = "NH" ; choose hemisphere
   res@gsnMaximize = True

   res@cnFillOn = True ; turn on color
   res@cnLinesOn = False ; no contour lines

   res@mpLandFillColor = "grey" ; color of land
   res@mpMinLatF = latLimit ; adjust map limit

  ;res@cnLevelSelectionMode = "ManualLevels" ; set manual contour levels
  ;res@cnMinLevelValF = 0. ; set min contour level
  ;res@cnMaxLevelValF = 5. ; set max contour level
  ;res@cnLevelSpacingF = 0.5 ; set contour spacing

   res@gsnSpreadColors = True ; use full colormap
  ;res@gsnSpreadColorStart = 2
  ;res@gsnSpreadColorEnd = 97

   res@tiMainString = "FISEL"

   nt = 0 or whatever
   plot = gsn_csm_contour_map_polar(wks,remap(nt,{latLimit:},:),res)

On 2/22/11 4:17 PM, Brandon Fisel wrote:
> Hi,
>
> I have global ECMWF data that have 1D lat/long that I wish to plot using
> a polar stereographic projection. The data itself is 3D (time, lat, lon)
> and I'm only interested in data above 65 N. At the moment, the data will
> only plot correctly if I use a cylindrical equidistant projection.
>
> Is there a method to do this?
>
>
> -------------------------------------------------------------
> Brandon Fisel
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Feb 22 16:29:13 2011