[ncl-talk] Shading significant grid points

Anahita Amiri Farahani aamir003 at ucr.edu
Sat Jan 16 11:02:59 MST 2016


Hi Karin,

Thanks for your answer. I tried it, but didn't work.

Cheers,
Ana

On Sat, Jan 16, 2016 at 1:49 AM, Karin Meier-Fleischer <
meier-fleischer at dkrz.de> wrote:

> Hi Ana,
>
> the warning message is telling you exactly what to do. Set
>
>    res at gsnAddCyclic = False
>
> because your grid is not global (longitudes are not 0-360 degrees), as you
> already said
>
> lat ranges from 0.5 to 29.5 and lon ranges from -45.5 to 14.5 and the
> resolution is 1*1 degree.
>
>
> Bye,
> Karin
>
>
> Am 16.01.2016 um 04:19 schrieb Anahita Amiri Farahani <aamir003 at ucr.edu>:
>
> Dear All,
>
> I tried to shade where only significant grid points and wrote this script:
>
>
>
> sig2D!0="lat"
> sig2D!1="lon"
> ;Coordinate
> sig2D&lat=lat
> sig2D&lon=lon
>
> alb_clr_der_summer2D!0="lat"
> alb_clr_der_summer2D!1="lon"
> ;Coordinate
> alb_clr_der_summer2D&lat=lat
> alb_clr_der_summer2D&lon=lon
>
> a = addfile("$NCARG_ROOT/lib/ncarg/data/cdf/landsea.nc","r")
> lsdata = a->LSMASK
> lsm  =
> landsea_mask(lsdata,alb_clr_der_summer2D&lat,alb_clr_der_summer2D&lon)
> alb_clr_der_summer2D = mask(alb_clr_der_summer2D,lsm.eq.1,False)
> ;alb_clr_der_summer2D = mask(alb_clr_der_summer2D,lsm.eq.3,False)
>
>
>
> lsm  = landsea_mask(lsdata,sig2D&lat,sig2D&lon)
> sig2D = mask(sig2D,lsm.eq.1,False)
> ;sig2D = mask(sig2D,lsm.eq.3,False)
>
>
> print("Plot")
>
> wks = gsn_open_wks("eps","1st_der_test")
> gsn_define_colormap(wks,"BlWhRe")   ; choose color
> ;gsn_define_colormap(wks,"BlAqGrYeOrRe")
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>   res                     = True               ; plot mods desired
>   res at cnFillOn            = True              ; turn on color
>   res at gsnSpreadColors     = True               ; use full color map
>   res at cnLinesOn           = False              ; no contour lines
>   res at cnLineLabelsOn      = False              ; no line labels
>
>   res at mpMaxLatF = 29.5                         ; specify the plot domain
>   res at mpMinLatF = 0.5                         ;
>   res at mpMinLonF = -44.5                        ;
>   res at mpMaxLonF = 15.5
>
>   res at gsnAddCyclic        = False
>
>  ; res at cnLevelSelectionMode =  "AutomaticLevels"   ; AutomaticLevels
>
>
>   res at cnLevelSelectionMode =  "ManualLevels"   ; AutomaticLevels
>   res at cnMinLevelValF       = -1.              ; min level
>   res at cnMaxLevelValF       =  1.               ; max level
>   res at cnLevelSpacingF      =  0.1               ; interval
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
> res at tiMainString        = "1st der test"
>
>
> res at gsnDraw              = False           ; Do not draw plot
> res at gsnFrame             = False           ; Do not advance frame
>
> plot= gsn_csm_contour_map_ce(wks,alb_clr_der_summer2D(:,:),res)
>
> ;delete(wks)
>
> ; ========================= PLOT 2 ==============================
>  res2 = True                            ; res2 probability plots
>
>
>   res2 at gsnDraw              = False           ; Do not draw plot
>   res2 at gsnFrame             = False           ; Do not advance frame
>   res2 at cnLevelSelectionMode = "ManualLevels" ; set manual contour levels
>   res2 at cnMinLevelValF      = 0.00        ; set min contour level
>   res2 at cnMaxLevelValF      = 1.05        ; set max contour level
>   res2 at cnLevelSpacingF     = 0.05        ; set contour spacing
>
>   res2 at cnInfoLabelOn       = False       ; turn off info label
>
>   res2 at cnLinesOn           = False       ; do not draw contour lines
>   res2 at cnLineLabelsOn      = False       ; do not draw contour labels
>
>   res2 at cnFillScaleF        = 0.6         ; add extra density
> ;  delete(alpha at long_name)
> ;  delete(alpha at units)
>
>   plot2   = gsn_csm_contour(wks,gsn_add_cyclic_point(sig2D(:,:)), res2)
>   plot2   = ShadeGtContour(plot2, 0.95, 17)  ; shade all areas less than
> the
>                                              ; 0.05 contour level
>   overlay (plot, plot2)
>
>   draw (plot)
>   frame(wks)
>
>
>
>
> but I got this warning:
>
>
> (0)     gsn_add_cyclic: Warning: The range of your longitude data is not
> 360.
>
> (0)     You may want to set gsnAddCyclic to False to avoid a warning
>
> (0)     message from the spline function.
>
> warning:_NhlCreateSplineCoordApprox: Attempt to create spline
> approximation for X axis failed: consider adjusting trXTensionF value
>
> warning:IrTransInitialize: error creating spline approximation for
> trXCoordPoints; defaulting to linear
>
>
> lat ranges from 0.5 to 29.5 and lon ranges from -45.5 to 14.5 and the
> resolution is 1*1 degree.
>
>
> Thanks,
>
> Ana
>
>
>
>
>
>
>
> _______________________________________________
> 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/20160116/fd3b72a3/attachment.html 


More information about the ncl-talk mailing list