support for gx1v6 grid

[l_vankampenhout@uu_nl]

Hi Dave,

The ice diagnostic script "ice_diag_20120323" seems to be lacking support for gx1v6 grids:

fatal:ContourPlotDraw: coordinates are out of range for drawing over a map: standard AreaFill rendering method will not work;
consider setting the resource trGridType to "TriangularMesh" if coordinates contain missing values
fatal:ContourPlotDraw: coordinates are out of range for drawing over a map: standard AreaFill rendering method will not work;
consider setting the resource trGridType to "TriangularMesh" if coordinates contain missing values

It took me a while to figure out what was wrong. It turns out NCL doesn't like the missing values in the TLAT and TLON arrays of this grid. Although I can create a dirty hack myself, perhaps you do have a proper fix for this?

Thanks,
Leo

 

[dbailey]

Hi Leo,Glad to hear you are using the diagnostic package. One quick and dirty fix is to use NCL version 6.1 or earlier. There was something in 6.2 that changed how missing values are used in the coordinate arrays. I do have a more up to date version, but I don't have it on the release repository yet. If you have access to yellowstone, you can have a look at:

/glade/u/home/dbailey/ice_diag_trunkor/glade/u/home/dbailey/ice_diag_swift

The latter is a new Swift (parallelized) version. We have it working on geyser/caldera, but I haven't tried it on many other systems. Otherwise I will try to put out a new release version in the coming months. I am not currently the sea ice liaison, so I'm not being paid to work on this.Dave

 

[l_vankampenhout@uu_nl]

Hi Dave,the version located at /glade/u/home/dbailey/ice_diag_trunkworks well for me. I already hoped I would not be the first person discovering this. Thanks!Leo