[ncl-talk] sub: plot ASCII data
dale zuri
dalezuri at gmail.com
Tue Mar 30 13:51:27 MDT 2021
Dear ncl users,
I managed to get a map. But still some error.
I would appreciate any help to fix this error.
Thanks
Error:
warning:ScalarFieldSetValues: irregular coordinate array sfXArray
non-monotonic: defaulting sfXArray
warning:ScalarFieldSetValues: irregular coordinate array sfYArray
non-monotonic: defaulting sfYArray
warning:ScalarFieldSetValues: irregular coordinate array sfXArray
non-monotonic: defaulting sfXArray
warning:ScalarFieldSetValues: irregular coordinate array sfYArray
non-monotonic: defaulting sfYArray
;----------------------------------------------------------------------
; dataonmap_5.ncl
;
begin
nlat=720
nlon=1440
ncol = 6
nhead = 1
fname = "avhrr_seaice"
nrow = numAsciiRow (fname)
data = asciiread (fname,(/nrow,ncol/), "float")
data at FillValue = -999.9
printVarSummary(data)
lat = data (::nlat,0)
mlat = dimsizes(lat)
lon = data (::nlon,1)
mlon = dimsizes(lon)
precip = data (:,0)
printVarSummary(precip)
pre = onedtond(precip,(/mlat,mlon/))
printVarSummary(pre)
lat!0 = "lat"
lon!0 = "lon"
pre!0 = "lat"
pre!1 = "lon"
;-----------------------------------------------------------------------
; Assign coordinate variables
;-----------------------------------------------------------------------
lat&lat = lat
lon&lon = lon
pre&lat = lat
pre&lon = lon
pre&lat at units="degrees_north"
pre&lon at units="degrees_east"
;-----------------
printVarSummary(pre)
wks = gsn_open_wks("pdf","test") ; send graphics to PNG file
plot = new(2,graphic)
cmap = read_colormap_file("UKM_hadcrut"); choose colormap
res = True ; plot mods desired
res at cnFillOn = True ; turn on color
;res at mpShapeMode = "FreeAspect"
res at vpWidthF = 0.3
res at vpHeightF = 0.2
; res at cnFillPalette 2 = "grads_rainbow" ; choose colormap
res at cnFillPalette = cmap ; choose colormap
res at gsnPolar = "NH" ; specify the hemisphere
res at gsnDraw = False
res at gsnFrame = False
res at gsnTickMarksOn = False
res at mpMinLatF = 50 ; specify min lat
res at gsnAddCyclic = False
; res at mpMinLatF = minlatt ; specify min lat
;res at mpMaxLatF = maxlatt
res at trGridType = "TriangularMesh"
res at cnLinesOn = False ; True is default
;res at mpFillOn = False ; turns off continent gray
res at cnLineLabelsOn = False ; True is default
res at mpGridAndLimbOn = False
res at mpLabelFontHeightF = "22"
res at mpPerimLineThicknessF = 2
res at mpOutlineOn = True
; res at mpLandFillColor = "tan"
;res at tiMainString = ""+type+"-"+mon+""
res at lbLabelBarOn = False ; turn off individual lb's
res at gsnLeftString = "" ; no titles
res at gsnRightString = ""
res at tiXAxisString = ""
res at tiYAxisString = ""
res at cnLevelSelectionMode = "ManualLevels" ; set manual contour levels
res at cnMinLevelValF = 0;0 ; set min contour level
res at cnMaxLevelValF = 100 ;6 ; set max contour level
res at cnLevelSpacingF = 10;0.5 ; set contour spacing
;res at lbOrientation = "Horizontal"
;res at lbLabelFontHeightF = 0.014
;res at lbLabelFont = "22"
; res at lbTitleString = ""+type+""
; res at lbTitleString = " (deg C)" ; bar title
; res at lbTitlePosition = "Right" ; title
location
; res at lbTitleDirection = "Across" ; letter
angle
; res at lbTitleAngleF = 90. ; title
angle
; res at lbTitleFontHeightF = 0.02 ; font
height
res at gsnStringFont = "helvetica-bold" ; will be used for
gsnCenterString later
res at gsnCenterStringFontHeightF = 0.014
res at gsnCenterString = "NDJFM"
plot(0) = gsn_csm_contour_map_polar(wks,pre,res)
res at gsnCenterString = "JJAS"
plot(1) = gsn_csm_contour_map_polar(wks,pre,res)
;********************************
;;;;;Contour plot;;;;;;;;;;;;;;;
;********************************
; panel plot only resources
;********************************
; panel plot only resources
; panel plot only resources
resP = True ; modify the panel plot
resP at gsnMaximize = True ; large format
resP at gsnPanelLabelBar = True ; add common colorbar
resP at lbLabelFontHeightF = 0.012
resP at gsnPanelBottom = 0.2;;;increases figure size ;
shrink panel plot by setting bottom edge of plot
resP at gsnPanelTop = 0.9 ; shrink panel plot by
setting top edge of plot
resP at pmLabelBarWidthF = 0.50 ; make thinner/small number
; resP at pmLabelBarHeightF = 0.05
resP at pmLabelBarHeightF = 0.08
resP at pmLabelBarOrthogonalPosF = -0.016
; resP at gsnPanelYWhiteSpacePercent = 5. ; increase spacing along
Y-axis between panel plots
resP at gsnPanelYWhiteSpacePercent = 2;;increases the figure size ;
increase spacing along Y-axis between panel plots
res at lbLabelFontHeightF = 0.018
resP at lbLabelFont = "22"
; resP at lbTitleString = ""+tic+""
resP at cnLinesOn = False ; turn off contour lines
resP at lbLabelFont = "helvetica-bold" ; add common colorbar
resP at lbTitleFont = "helvetica-bold" ; add common colorbar
resP at lbTitleFontHeightF= .015 ; make title smaller
resP at cnLineLabelsOn = False
resP at txFontHeightF = 0.03
resP at txFont = "helvetica-bold"
gsn_panel(wks,plot,(/1,2/),resP)
;---This will maximize the size of all the paneled stuff.
;end do
end
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210330/04e937cb/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test.pdf
Type: application/pdf
Size: 568402 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210330/04e937cb/attachment-0001.pdf>
More information about the ncl-talk
mailing list