[ncl-talk] Radar PPI dimension size problem with Conform
Ryan Connelly
rconne01 at gmail.com
Wed Apr 19 14:28:56 MDT 2017
Hi,
I'm looking at the example here:
http://www.ncl.ucar.edu/Applications/Scripts/radar_2.ncl but encountering
the error:
fatal:conform: the dimensions sizes of the second argument do not match
those indicated by the third argument
Here is the code up to that point:
begin
fils = systemfunc ("csh -c 'cd " + diri + " ; ls *.nc'")
do j = 0,dimsizes(fils)-1,1
fili = fils(j) ;+".nc"
rf = addfile(diri+fili,"r")
; print(getfilevarnames(rf))
rd = rf->REF ;short2flt(rf->REF)
dsizes = dimsizes(rd)
angles = rf->azimuth
angles(0:63) = angles(0:63)-360 ; fix to make angles monotonic
DEGTORAD = 0.017453292519943
xcenter = 0.0
ycenter = 0.0
radius = 960 * 0.25 ; this is radius in kilometers
; Create 2D coordinate arrays.
;
printVarSummary(rd(0))
printVarSummary(angles)
inc = radius / (dsizes(0) - 1)
x = inc * ispan(0,dsizes(0)-1,1)
angles2d = conform(rd,angles,0)
x2d = conform(rd,x,1)
xarr = xcenter + x2d * cos(DEGTORAD * angles2d)
yarr = ycenter + x2d * sin(DEGTORAD * angles2d)
It seems that there's a mismatch between rd(0) and angles? rd is
one-dimensional with 6701760 points, and angles is also one-dimensional,
with 5400 points.
Thanks,
Ryan
--
Ryan Connelly
M.S. Student in Atmospheric Sciences, Stony Brook University
B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
University
rconne01 at gmail.com
ryan.connelly at stonybrook.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170419/82a9293e/attachment.html
More information about the ncl-talk
mailing list