<div dir="ltr"><div><div>Please do not write directly to me:   <a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>  only<br><br>---<br></div>A possible reason you have had no response is that the question(s) posed are not clear ... at least to me.<br><br>--<br>Question 1:<br><br>What do you mean: "Say I have WRF output on standard lat/lon grid".<br>To what 'standard lat/lon grid' are you referring?<br><br>Typically:<br>WRF lat/lon are (Time, south_north, west_east).<br>WRF Variables are: (Time, bottom_top, south_north, west_east)  and (Time, south_north, west_east)<br>WRF grids are regional not global<br><br>Climate models are often global (time,level,lat,lon), (time,lat,lon)<br>Grids are rectilinear: lat(lat), lon(lon)<br>Some are regional.<br><br>---<br>Cylindrical coordinates; I have never used cylindrical  x,y,z coordinates.<br>Do you mean (time,z,lat,lon)? ... That is *not* a math cylindrical coord system.<br><br>This is a possible generic approach.<br>---<br><br>   xp = ... ; (time,Np,lat2d,lon2d)<br>   xz <=== xp               ; <=== *you* must interpolate to 'z' (Nz)<br><br><br>   clat = ....            ; center lat at each time<br>   clon = ....            ;        lon<br>   crad = (/ 1,3,5,.../)  ; radii (degrees)<br><br><br>   Ntim = dimsizes(clat)  ; # time steps<br>   Nrad = 5               ; # of radii<br>   Npts = 100             ; # of points at each radii<br> <br>   Rlat = new( Npts, "float")   ; lat at each radius and time<br>   Rlon = new( Npts, "float")<br><br>   Drad = new( (/Ntim,Nrad,Nz,Npts/), "float")  ; data at each time/radii<br>  <br>   do nt=0,Ntim-1<br>     do nr=0,Nrad-1<br>        nggcog(clat(nt), clon(nt), crad(nr), Rlat, Rlon)<br>       do nz=0,KZ-1<br>            ... interpolate to Z ... some WRF function<br>            ; curvilinear grid<br>            Drad(nt,nr,nz,:) = rcm2points_Wrap (lat2d, lon2d, x(nt,nz,:,:), Rlat, Rlon, 0)<br>            ; linint2 for gloable rectilinear grids<br>     end do<br>   end do<br><br>; Likely very slow ... *save to netCDF once it works correctly<br><br>;----------<br>Question 2: "cylindrical equidistant projection"  is a map projection.<br><br>Is your WRF grid global or limited area (regional)<br><br>Using an FFT on a regional grids makes no sense to me; only on global grids.<br></div>You can not get appropriate wave numbers on a regional grid.<br><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Oct 5, 2017 at 6:26 PM, Prashanth Bhalachandran <span dir="ltr"><<a href="mailto:prashanth.bhalachandran@gmail.com" target="_blank">prashanth.bhalachandran@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div style="font-family:PalatinoLinotype-Roman"><font size="3" face="Palatino">Hello and Warm Greetings !!</font></div><div style="font-family:PalatinoLinotype-Roman"><font size="3" face="Palatino"><br></font></div><div style="font-family:PalatinoLinotype-Roman"><font size="3" face="Palatino">I am Prashanth, a PhD student at Purdue and I primarily study hurricanes. I had a quick question regarding NCL that I was hoping you’d help me with. </font></div><div style="font-family:PalatinoLinotype-Roman"><font size="3" face="Palatino"><br></font></div><div style="font-family:PalatinoLinotype-Roman"><font size="3" face="Palatino">Say I have WRF output on standard lat/lon grid. I am interested in converting all the data to storm-centric (I know the coordinates of the moving storm center), cylindrical coordinates. That is lat-lon (x,y,z) to r, theta and z. Can you please tell me as to what the best way to do this is using NCL? </font></div><div style="font-family:PalatinoLinotype-Roman"><font size="3" face="Palatino"><br></font></div><div style="font-family:PalatinoLinotype-Roman"><font size="3" face="Palatino">I also noticed that NCL had a cylindrical equidistant projection? However, I understand that this is only for plotting? I am guessing I need to do a coordinate transform compulsorily because I intend to take a Fourier transform in the theta direction. </font></div><div style="font-family:PalatinoLinotype-Roman"><font size="3" face="Palatino"><br></font></div><div style="font-family:PalatinoLinotype-Roman"><font size="3" face="Palatino">Any direction from your end will be greatly appreciated. </font></div><div style="font-family:PalatinoLinotype-Roman"><font size="3" face="Palatino"><br></font></div><div style="font-family:PalatinoLinotype-Roman"><font size="3" face="Palatino">Please help me out. </font></div><span class="HOEnZb"><font color="#888888"><div style="font-family:PalatinoLinotype-Roman"><font size="3" face="Palatino"><br></font></div><div style="font-family:PalatinoLinotype-Roman"><font size="3" face="Palatino">Prashanth </font></div></font></span></div><br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>