[ncl-talk] Radar PPI dimension size problem with Conform

[Ryan Connelly]

Hi Barry,

No idea about the format, though it's possible since DBZ is not a variable
in the netCDF (it's REF now).  I use radx to convert the Level II files
from their native format to CfRadial netCDF's, nothing involving the WCT.

Thanks,
Ryan

On Wed, Apr 19, 2017 at 10:45 PM, Barry Lynn <barry.h.lynn at gmail.com> wrote:

> 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/Applic
>> ations/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 <(914)%20432-3108>
>



-- 
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/20170420/0ddbb05c/attachment.html