Re: Issue plotting data on icosahedral grid

From: jim rosinski <James.Rosinski_at_nyahnyahspammersnyahnyah>
Date: Wed Dec 21 2011 - 13:40:43 MST

Excellent! Thanks Dave. I'd be happy to test a fix whenever you'd like
me to try. I have access to bluefire if that's the easiest place to test.

Jim R.

On 12/21/2011 12:46 PM, David Brown wrote:
> Hi Jim,
> I have determined that there is a bug in the code that tests whether a point is inside the projectable area that is causing this drop-out. It does not affect the plot of Dave Randall's data just because of small differences in the coordinate values near the pole. (A complete explanation would be too involved.) I am currently testing a fix to make sure it is reliable and doesn't break something else. I intend to check it in today so it should be available soon if you need it.
> -dave
>
>
> On Dec 16, 2011, at 3:55 PM, Jim Rosinski wrote:
>
>> Hello;
>>
>> I am having an issue using ncl to plot data on an irregular (icosahedral) grid. The model is FIM, a forecast model under development at NOAA. There is a template ncl script and test data to plot icosahedral data from Dave Randall (CSU) provided with the ncl distribution. I tried to follow the format both in creation of the model data, and the ncl script to plot it, as closely as possible.
>>
>> We're 99% of the way there. The attached cgm file shows what appears to be a correct plot, except data at the north and south poles is strangely blacked out. I've spoken with the model developers--they've looked at the data and are convinced that the required values for grid point centers and corner points is correct. Can anyone help with an explanation of what might be happening and/or propose a fix? Many thanks in advance.
>>
>>
>> Regards,
>> Jim Rosinski
>> james.rosinski@noaa.gov
>>
>>
>> FYI, on bluefire:/ptmp/rosinski/fim2nc/testdata, the model input file "testicos.nc" contains the model data we're plotting. In the same directory, ploticos.ncl is the ncl script we're using to plot the data. geodesic.ncgm is the cgm file with the data for the attached plot. geo_1.ncl and geo_2.ncl are copies of the ncl script to plot Dave Randall's data. Mysteriously, they seem to work fine on his data.
>> <Screenshot-geodesic.ncgm-2.png>_______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Dec 21 13:41:04 2011

This archive was generated by hypermail 2.1.8 : Wed Jan 04 2012 - 10:18:09 MST