[ncl-talk] Radar PPI dimension size problem with Conform
Barry Lynn
barry.h.lynn at gmail.com
Wed Apr 19 20:45:15 MDT 2017
Hi Ryan:
Is it possible that the format of the radar data has changed?
I had to write several Fortran programs to read the radar data from
different years (as processed by the NOAA Weather and Climate Tool). I
write out the data and then map it objectively using a wrf output program
as the gridded background.
I can send you these programs if you like.
Barry
On Wed, Apr 19, 2017 at 11:28 PM, Ryan Connelly <rconne01 at gmail.com> wrote:
> 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
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
--
Barry H. Lynn, Ph.D
Senior Lecturer,
The Institute of the Earth Science,
The Hebrew University of Jerusalem,
Givat Ram, Jerusalem 91904, Israel
Tel: 972 547 231 170
Fax: (972)-25662581
C.E.O, Weather It Is, LTD
Weather and Climate Focus
http://weather-it-is.com
Jerusalem, Israel
Local: 02 930 9525
Cell: 054 7 231 170
Int-IS: x972 2 930 9525
US 914 432 3108
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170420/0ef14936/attachment.html
More information about the ncl-talk
mailing list