[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