[ARPSSUPPORT] ARPSPLT terrain outline bug

[Xue, Ming]

 I found the bug:
          DO i = 2,nx-2  in the following (two places) was mis-typed as DO i=2,ny-2.
And I changed line 235 to (added +1 to kend).
        235 kend = MIN(kend+1, nz-1)

Yunheng, please update the official version.
Ming

 IF(zbgn /= zend) THEN
      kend = 2
      DO k = 2,nz-1
        DO j = 2,ny-2
          DO i = 2,nx-2
            IF(zpc(i,j,k) < zend) THEN
              kend=k
              GO TO 225
            END IF
          END DO
        END DO
        GO TO 235
  225   CONTINUE
      END DO
  235 kend = MIN(kend+1, nz-1)

    print*,'zend, kend =', zend, kend

      kbgn= nz-1
      DO k = nz-1,2,-1
        DO j = 2,ny-2
          DO i = 2,nx-2
            IF(zpc(i,j,k) > zbgn) THEN
              kbgn=k
              GO TO 250
            END IF
          END DO
        END DO
        GO TO 245
  250   CONTINUE
      END DO
  245 kbgn = MAX(kbgn,2)
Ming

> -----Original Message-----
> From: arpssupport-bounces at caps.ou.edu [mailto:arpssupport-bounces at caps.ou.edu] On Behalf Of Michael
> Kiefer
> Sent: 2009年7月22日 12:01
> To: arpssupport at caps.ou.edu
> Subject: [ARPSSUPPORT] ARPSPLT terrain outline bug
> 
> Hello,
> 
> I am having an issue with ARPSPLT when I plot vertical cross-sections in uneven terrain.  The plots have a
> strange "terrain shadow" effect wherein the upper part of the cross-section is not displayed; the shape of the
> blacked-out region is the same shape as the terrain profile.  This is seen in the figure I have attached.  I have
> zeroed in on the place in arpspltlib.f90 where it calls the zxplot subroutines to construct the terrain profile, but I
> am unsure of where to go from there.  Any help is much appreciated.
> 
> Mike
> 
> 
> Michael Kiefer
> Research Associate
> Department of Geography
> Michigan State University
> East Lansing, MI 48823