[ncl-talk] sub: read ascii error

dale zuri dalezuri at gmail.com
Tue Mar 30 11:31:14 MDT 2021


Thanks for your response. I'm able to read the file. But not able to figure
out the dimension coordinates.



On Tue, Mar 30, 2021 at 11:27 AM david brown <david.ian.brown at gmail.com>
wrote:

> I think you should change:
>
> fname = "avhrr_seaice"
> data = readAsciiTable("avhrr_seaice", ncol, "float", nhead)
>
> to
>
> fname = "avhrr_seaice"
> data = readAsciiTable(fname, ncol, "float", nhead)
>
>
>
>
> On Tue, Mar 30, 2021 at 10:24 AM dale zuri via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> 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/9d74ea66/attachment.html>


More information about the ncl-talk mailing list