[ncl-talk] sub: read ascii error
Alessandra Giannini
alesall at iri.columbia.edu
Tue Mar 30 12:00:56 MDT 2021
Hi there,
can you print out:
printVarSummary(data)
?
are lat and lon really included in the matrix called “data”?
The way you define “templat” and “templon”, they should have the same dimension. (Is said dimension nlat x nlon?)
At any rate, lon and lat should be 1-dimensional.
Rule #0 of ncl-talk is that you should have an idea of the structure of your data…
Maybe you do, but you are not telling us!
warm regards, alessandra
> On Mar 30, 2021, at 1:45 PM, dale zuri via ncl-talk <ncl-talk at mailman.ucar.edu> wrote:
>
> Hi Adam,
> I would appreciate any help to plot this ASCII data.
>
> Thanks
> DZ
>
>
> ---------- Forwarded message ---------
> From: dale zuri <dalezuri at gmail.com <mailto:dalezuri at gmail.com>>
> Date: Tue, Mar 30, 2021 at 10:24 AM
> Subject: sub: read ascii error
> To: Ncl-talk <ncl-talk at ucar.edu <mailto:ncl-talk at ucar.edu>>
>
>
> Dear NCL user,
> I would appreciate any help to plot this data. I'm getting this error.
>
> Thanks in advance
> DZ
>
> fatal:Variable (lat) is undefined
> fatal:["Execute.c":8637]:Execute: Error occurred at or near line 18 in file test.ncl
>
>
>
>
> begin
> nlat=720
> nlon=1440
>
> ncol = 6
> nhead = 1
>
> fname = "avhrr_seaice"
> data = readAsciiTable("avhrr_seaice", ncol, "float", nhead)
> temp12D = onedtond(data(2,:),(/nlat,nlon/)) ; convert 1D array to a 2D array
> temp22D = onedtond(data(3,:),(/nlat,nlon/)) ; convert 1D array to a 2D array
> templat = data(0,:) ; convert 1D array to a 2D array
> templon = data(1,:) ; convert 1D array to a 2D array
> ; Assign named dimensions
> temp12D!0 = "lat"
> temp12D!1 = "lon"
> temp12D&lat = templat ; assign coord variables
> temp12D&lon = templon
> temp22D!0 = "lat"
> temp22D!1 = "lon"
> temp22D&lat = templat ; assign coord variables
> temp22D&lon = templon
> temp12D&lat at units = "degrees_north"
> temp12D&lon at units = "degrees_east"
> temp22D&lat at units = "degrees_north"
> temp22D&lon at units = "degrees_east"
> ; Assign coordinate variables
> printVarSummary(temp12D)
> printVarSummary(temp22D)
>
> wks = gsn_open_wks("x11","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 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,temp12D(:,0),res)
> res at gsnCenterString = "JJAS"
> plot(1) = gsn_csm_contour_map_polar(wks,temp22D(:,1),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
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210330/66a6c016/attachment.html>
More information about the ncl-talk
mailing list