Re: NARR subset plots

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Fri, 15 Aug 2008 15:58:51 -0600

Try something like

   wks = gsn_open_wks ("x11", "narr_test") ; open workstation
   gsn_define_colormap (wks,"amwg")

   res = True ; plot mods
   res_at_cnFillOn = True ; color fill
   res_at_cnLinesOn = False ; no contour lines
   res_at_gsnMaximize = True ; make large
   res_at_gsnSpreadColors = True ; use total colormap
  ;res_at_gsnSpreadColorStart = 4
  ;res_at_gsnSpreadColorEnd = -1
   res_at_lbLabelAutoStride = True
   res_at_mpGridAndLimbOn = True
   res_at_mpFillOn = False
   res_at_pmTickMarkDisplayMode = "Always" ; turn on tickmarks
   res_at_tmXTOn = False
   res_at_gsnAddCyclic = False ; regional data

   res_at_mpMinLatF = min(lat2d) - 1 ; extra space
   res_at_mpMaxLatF = max(lat2d) + 1
   res_at_mpMinLonF = min(lon2d) - 1
   res_at_mpMaxLonF = max(lon2d) + 1

   res_at_cnLevelSelectionMode = "ManualLevels" ; set manual contour levels
   res_at_cnMinLevelValF = 230. ; set min contour level
   res_at_cnMaxLevelValF = 315. ; set max contour level
   res_at_cnLevelSpacingF = 5. ; set contour spacing

   res_at_tiMainString = "Original grid: "+date(0)
   plot = gsn_csm_contour_map_ce(wks,x(0,:,:),res) ; Draw original
grid on map

   res_at_mpMinLatF = 35.
   res_at_mpMaxLatF = 45.
   res_at_mpMinLonF = -120.
   res_at_mpMaxLonF = -105.
   res_at_tiMainString = "subset grid: "+date(0)
   plot = gsn_csm_contour_map_ce(wks,x(0,:,:),res)
   end

prabhakar shrestha wrote:
> I have been trying to plot NARR datasets however the best I can do is
> plot the whole NARR data in its own projection.
> Is there a way to plot only a subset of the data (specified boundary lat
> long) in geographic co-ordinates?
>
> Regards,
> Prabhakar Shrestha
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Aug 15 2008 - 15:58:51 MDT

This archive was generated by hypermail 2.2.0 : Sun Aug 17 2008 - 23:24:33 MDT