[ncl-talk] Terrain and SST

[Rick Brownrigg]

Hi,

Not really sure what you mean by "...nothing works" -- what result do you
get?

On a general note, in your first call to the overlay() function, usually
the map-based plot is the first argument, i.e, the basis of the overlay.
Finally, the call to maximize_output() is either called in lieu of
draw()/frame(), or must be called before draw()/frame();  see the docs at:

    http://ncl.ucar.edu/Document/Graphics/Interfaces/maximize_output.shtml

Rick

On Fri, Feb 8, 2019 at 10:20 AM Andrea Perez <andreaperez151 at yahoo.es>
wrote:

> Hello,
> I am using MM5 ouput and I want to plot terrain, wind vector and SST in a
> map.
> I have tried with a mask, res_base but nothing works. Can any one help me??
>
> code:
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>
> begin
> f1=addfile("data.nc","r")
>
> lat=f1->latitcrs
> lon=f1->longicrs
> dims=dimsizes(lat)
> lat2d=(/lat/)
> lat2d at units = "degreeN"
> lon2d=(/lon/)
> lon2d at units = "degreeE"
>
> ter=f1->terrain
> ter at lat2d=lat2d
> ter at lon2d=lon2d
>
> u=f1->u10
> v=f1->v10
> u at lat2d=lat2d
> u at lon2d=lon2d
> v at lat2d=lat2d
> v at lon2d=lon2d
> med=ter
>
>        ts = f1->tseasfc(24,:,:)
>        ts = (ts*1)
>        ts at lat2d=lat2d
>        ts at lon2d=lon2d
>
>
> ;******************************************************************************
> wks = gsn_open_wks("eps","terrain_sst")          ; open a workstation
>
> ;******************************************************************************
>
>   res           = True
>   res at gsnDraw   = False
>   res at gsnFrame  = False
>   res at mpProjection      = "Mercator"
>   res at cnLinesOn       = False
>   res at cnFillOn        = True               ; color plot desired
>   res at mpOutlineDrawOrder     = "PostDraw"
>   res at mpFillDrawOrder        = "Predraw"
>   res at gsnSpreadColors     = True
>   res at cnLevelSelectionMode = "ManualLevels"
>   res at lbLabelAutoStride =   True
>   res at cnLineLabelsOn  = False              ; turn off contour lines
>   res at cnFillDrawOrder      = "Predraw"
>   res at mpDataBaseVersion     = "HighRes"    ;
>
>   res at mpProjection      = "Mercator"
>   res at mpDataBaseVersion     = "HighRes"
>   res at mpLimitMode       = "Corners"            ; choose range of map
>   res at mpLeftCornerLatF  = lat2d(0,0)
>   res at mpLeftCornerLonF  = lon2d(0,0)
>   res at mpRightCornerLatF = lat2d(dims(0)-1,dims(1)-1)
>   res at mpRightCornerLonF = lon2d(dims(0)-1,dims(1)-1)
>
> ;  res at tfDoNDCOverlay = True
>
>   res at mpGeophysicalLineColor = "black"                 ; color of
> continental outlines
>   res at mpPerimOn              = True                    ; draw box around
> map
>   res at mpGridLineDashPattern  = 2                       ; lat/lon lines as
> dashed
>   res at mpOutlineBoundarySets = "GeophysicalAndUSStates" ; add state
> boundaries
>   res at mpUSStateLineColor    = "black"                     ; make them red
>   res at mpOutlineOn            = True
>   res at mpDataBaseVersion     = "MediumRes"
>   res at mpDataSetName         = "Earth..4"
>   res at mpOutlineBoundarySets = "National"
>   res at mpGeophysicalLineColor      = "black"       ; color of cont.
> outlines
>   res at mpGeophysicalLineThicknessF = 3          ; thickness of outlines
>   res at mpNationalLineThicknessF = 3 ; interior boundaries
>
> ;**********************************************************
>   ter_res                      = res
> ;***********************Terrain****************************
>   ter_res at cnMinLevelValF = 0
>   ter_res at cnMaxLevelValF = 3000
>   ter_res at cnLevelSpacingF   =   250
>   ter_res at cnFillPalette = "OceanLakeLandSnow"   ; choose colormap
>   ter_res at mpFillOn              = False
> ;**************************SST*****************************
>   res_base                        = res
>   res_base at cnFillOn             = True     ; turn on color fill
>   res_base at cnLinesOn            = False    ; turn off contour lines
>   res_base at lbLabelStride        = 1        ; every 4th label
>
>   res_base at lbLabelAutoStride =   True
>   res_base at cnLevelSelectionMode = "ManualLevels"
>   res_base at cnMinLevelValF = 296
>   res_base at cnMaxLevelValF = 303
>   res_base at cnLevelSpacingF   =   1
>   res_base at cnLineLabelsOn  = False              ; turn off contour lines
>
>   res_base at tiMainString         = "Draw land ON TOP of contours"      ;
> title
>   res_base at cnFillDrawOrder      = "Predraw"  ; draw contours firs
>
> ;**************Create the two plots************************
>
>   ter_plot = gsn_csm_contour_map(wks,ter,ter_res)
>   sst_plot = gsn_csm_contour(wks,ts,res_base)
>
> ;*************Overlay the sst plot on the terrain plot*****
>   overlay(sst_plot,ter_plot)
>
> ;---Drawing the terrain plot will also draw dbz plot
>   draw(ter_plot)
>   frame(wks)
>
>
> ;---Remove the dbz plot from terrain plot so we can re-overlay it on a new
> plot
>   NhlRemoveOverlay(ter_plot,sst_plot,False)
>
> ;---Recreate the terrain plot with a fully opaque labelbar
>   ter_res at lbOverrideFillOpacity = True   ; don't apply opacity to
> labelbar colors
>                                       ; (NCL V6.4.0 and later)
>   ter_plot = gsn_csm_contour_map(wks,ts,ter_res)
>
> ;--Re-overlay the dbz plot on the new terrain plot.
>   overlay(ter_plot,sst_plot)
>
> ;---Drawing the terrain plot will also draw dbz plot
>   draw(ter_plot)
>   frame(wks)
>
>   maximize_output(wks,True)
>
> end
>
> _______________________________________________
> 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/20190211/76f4a274/attachment.html>