Dave:

I have boiled a problem down to the program listed below.
The data (read from unit 8) is on mss: /FOSTER/test/plot.dat.
The problem is that it crashes on frame 3, as follows:

shavano> a.out
frame   1
frame   2
AR004 - BAD SCALAR ARGUMENT TO ARLIB MATH ROUTINE
TB001 - BEGINNING OF TRACEBACK
	- trbk    	 WAS CALLED BY ARERP%  	AT 312176a
	- ARERP%  	 WAS CALLED BY ATAN2   	AT 312276d
	- ATAN2   	 WAS CALLED BY CPPLPS	AT 102217d (LINE NUMBER 166)
	- CPPLPS	 WAS CALLED BY CPPKLP	AT 60433d (LINE NUMBER 158)
	- CPPKLP	 WAS CALLED BY CPLBAM	AT 54254d (LINE NUMBER 103)
	- CPLBAM	 WAS CALLED BY TST	AT 616c (LINE NUMBER 40)
	- TST	 WAS CALLED BY $START$ 	AT 301c
TB002 - END OF TRACEBACK
Abort
TB001 - BEGINNING OF TRACEBACK
      - $TRBK    WAS CALLED BY f_sig    AT  352514b (LINE NUMBER      102)
      - f_sig    WAS CALLED BY __handlr AT  251620b 
      - __handlr WAS CALLED BY killm    AT  247540a 
      - killm    WAS CALLED BY raise    AT  247615a (LINE NUMBER       15)
      - raise    WAS CALLED BY abort    AT  217035a (LINE NUMBER       20)
      - abort    WAS CALLED BY ARERP%   AT  312177d 
      - ARERP%   WAS CALLED BY ATAN2    AT  312276d 
      - ATAN2    WAS CALLED BY CPPLPS   AT  102217d (LINE NUMBER      166)
      - CPPLPS   WAS CALLED BY CPPKLP   AT   60433d (LINE NUMBER      158)
      - CPPKLP   WAS CALLED BY CPLBAM   AT   54254d (LINE NUMBER      103)
      - CPLBAM   WAS CALLED BY TST      AT     616c (LINE NUMBER       40)
      - TST      WAS CALLED BY $START$  AT     301c 
TB002 - END OF TRACEBACK
Abort (core dumped)

The first 2 frames are ok, and the 3rd has all ezmap stuff ok, but no
contours.

If I change LLP from 3 to 2, it does not crash. Also, if I eliminate
plon, and call maproj('CE',0.,0.,0.) it does not crash. Any ideas??
(I am using the default ncarg libs on shavano, i.e., NOT 3.1.3)
(I think it also runs ok if LLP=3 and LLO=0).

--Ben


      program tst
c
      parameter(imx=73,jmx=36)
      dimension plot(imx,jmx),viewport(4),iasf(13)
      data viewport /.12,.90,.2,.70/
      data lu/8/, nframes/6/
      data iasf/13*1/
c
      parameter(lrwrk=2000,liwrk=1000,liama=300000,lcra=10000)
      common/work/ rwrk(lrwrk),iwrk(liwrk),iama(liama),iiasf(13),
     +  iaia(10),igia(10),xcra(lcra),ycra(lcra)
      external drawcl
c
      call opngks
      call gsclip(0)
      call gsasf(iasf)
      call gsfais(1)
c
      call cpsetr('XC1',-180.)
      call cpsetr('XCM',180.)
      call cpsetr('YC1',-87.5)
      call cpsetr('YCN',87.5)
      call cpseti('SET',0)
      call cpseti('MAP',1)
      call cpseti('LLP -- line label positioning',3)
      call cpseti('LLO -- line label orientation',1)
c
      call mapsti('EL',0)
      call mappos(viewport(1),viewport(2),viewport(3),viewport(4))
      plon = 0.
      do i=1,nframes
        read(lu) plot
        call maproj('CE',0.,plon,0.)
        call mapset('MA',0.,0.,0.,0.)
        call mapint
        call maplot
        call cprect(plot,imx,imx,jmx,rwrk,lrwrk,iwrk,liwrk)
        call arinam(iama,liama)
        call cplbam(plot,rwrk,iwrk,iama)
        call cpcldm(plot,rwrk,iwrk,iama,drawcl)
        call cplbdr(plot,rwrk,iwrk)
        call frame 
        write(6,"('frame ',i3)") i
        plon = plon - 15.
      enddo
      call clsgks
      stop
      end
c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
      subroutine drawcl(xcs,ycs,ncs,iai,iag,nai)
      dimension xcs(*),ycs(*),iai(*),iag(*)
      idr=1
      do i=1,nai
        if (iai(i).lt.0) idr=0
      enddo
      if (idr.ne.0) call curved(xcs,ycs,ncs)
      return
      end