Hi All,
This is my first time to use NCL for plotting my research data. Thus, I have not yet had any experience to do it. Hope your sharing for that.
I want to make a contour plot from the attached file (ASCII file). I use the following code:
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
begin
; read in data
ncols = 9 ; -10:2.5:10 Latitude
nrows = 33 ; 80:2.5:160E Longitude
data = asciiread ("djfOLR.txt", (/nrows,ncols/), "float")
wks = gsn_open_wks("ps" ,"ce") ; open a ps file
gsn_define_colormap(wks,"BlAqGrYeOrRe") ; choose colormap
res = True ; plot mods desired
res@cnFillOn = True ; turn on color fill
res@cnLinesOn = False ; turn of contour lines
res@cnLevelSpacingF = 0.5 ; contour spacing
res@gsnSpreadColors = True ; use full range of color map
res@lbLabelStride = 4
res@pmTickMarkDisplayMode = "Always"; use NCL default lat/lon labels
res@gsnAddCyclic = False ; data already has cyclic point
; this must also be set for any zoom
; note that the gsn_csm_*map_ce templates automatically set
; res@mpLimitMode="LatLon" for you. If you are plotting a different projection,
; you may have to set this resource.
res@mpMinLatF = -10 ; range to zoom in on
res@mpMaxLatF = 10.
res@mpMinLonF = 80.
res@mpMaxLonF = 160.
plot = gsn_csm_contour_map_ce(wks,data, res)
;************************************************
end
I have the result, but only the plot of Indonesian map without my data. Hope your sharing for this problem..
Regards,
JuKY
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2.1.8 : Wed Dec 21 2011 - 10:44:06 MST