completely filling in a state with data

From: Jeff Van Dorn <vandorn_at_nyahnyahspammersnyahnyah>
Date: Wed, 27 Sep 2006 15:41:50 -0500


I've produced this plot of California:

which is nice and pretty, but I'd ideally like to have the entire
state filled in, instead of having all the jagged boxes showing just
inside the perimeter of the state. I came across what I believe
might be a similar issue in the archives:
(basically trying to fill in the eastern coastline).

So first, is using linmsg (as was suggested in the archives) the
suggested way to go about doing what I want? I'm certainly content
in just doing a simple interpolation of the last data points near the
state boundary, which seems to be what linmsg is doing. However, I'm
having problems using it, or probably more likely, I'm using it
incorrectly. The originally suggested code to implement was:

arr = linmsg(arr,(/0,3/)) ; to fill in Cape Cod Bay area
do gg = 0,dimsizes(arr&lon)-1
     temp = ind(.not.ismissing(arr(gg,:)))
     arr(gg,temp(0)-1)=(/ ((arr(gg,temp(0))*2))-arr(gg,temp(0)+1) /)
end do

However, I believe this is testing for missing data within the data
itself, is that right? If so, I don't actually have any missing data
points in my data, so I'd ideally like to interpolate up to the
state's boundary if that's even possible (thereby implying that the
gap between my data points and the state's boundary is the missing
data). Or do you have any other suggestions in how to go about
getting rid of/reducing the jaggedness in the plot?

I'm using NCL 4.2.0.a033
Script and data at:

some script excerpts:
           res_at_mpOutlineOn = True
           res_at_mpPerimOn = True
           res_at_mpFillOn = True
           res_at_mpOutlineBoundarySets = "NoBoundaries"
           res_at_mpDataBaseVersion = "mediumres" ;
select database
           res_at_mpDataSetName = "Earth..2" ;
select database
           res_at_mpOutlineSpecifiers = "California . " + ca_counties
           res_at_cnFillOn = True
           res_at_cnFillMode = "AreaFill"
           res_at_cnSmoothingOn = True

Thanks for any suggestions!

ncl-talk mailing list
Received on Wed Sep 27 2006 - 14:41:50 MDT

This archive was generated by hypermail 2.2.0 : Mon Oct 09 2006 - 09:24:50 MDT