Drawing a Boundary Around a Model Domain

[Matt Georgescu]

Thanks to everyone that responded!

The final image looks like this:
http://cep.rutgers.edu/~matt1/boundary.gif

The brilliant idea, thanks to (another) Matt, was to create another data
file, consisting of only values that are on the model domain, and values
that are not, using a simple Fortran program:

      do i=1,num_x
         do j=1,num_y
           IF(i .eq. 1 .or. i .eq. num_x)THEN
               edge(i,j)=1.
             ELSEIF(j .eq. 1 .or. j .eq. num_y)THEN
               edge(i,j)=1.
             ELSEIF(i .eq. 215)THEN
               edge(i,j)=1.
             ELSE
               edge(i,j)=0.
             ENDIF
         enddo
      enddo

I then open up both control files (the one from my model output and the
one corresponding from the data created above) and plot ... the issue here
was that i was not just plotting a box and needed to take curvature, due
to the size of the domain, into account.

matt