****************************************************************************** BSW Electronic Mail 20-Nov-1997 23:05:05 CET Message Number 0021 ****************************************************************************** Author: Yuki Hatanaka Subject: Bug in UPDPOL There is a bug in POLUPD.f. The problem is that when a pole file is updated, the last entry in the original pole file is lost. I attach a part of the source code of $FG/POLUPD.f (line 73-). The problem in reading the IGS ERP file is solved by modifying just one line indicated in the list which is supposed to be a bug. But I'm not sure if this change cause no problem in reading ERP files of other formats. - Yuki >C >C READ INFO OF ALL INPUT POLE FILES >C --------------------------------- > NVAL=1 > NFLCOL=1 > CALL GTFILE('INFIL ',NFLCOL,MAXFIL,NFIL,INFIL) > DO 30 I1=1,NFIL >C WRITE(*,*)'nfil,i1,filnam= ',NFIL,I1,' ',INFIL(I1) > IFORM=0 > CALL OPNFIL(LFN001,INFIL(I1),'OLD','FORMATTED', > 1 'READONLY',' ',IOSTAT) > CALL OPNERR(LFNERR,LFN001,IOSTAT,INFIL(I1),'POLUPD') > DO 40 I2=NVAL,MAXVAL > CALL RDPOL(LFN001,IERB1D,IFORM,CODTIT,NUTCOD,POLTIM(I2), > 1 POLCOO(1,I2),GPSUTC(I2),REM(I2),RMSPOL(1,I2), > 2 IEND) > IF (IFORM.EQ.5.OR.IFORM.EQ.6) THEN > IF (REM(I2).EQ.'DEF')REM(I2)='B5F' > IF (REM(I2).EQ.'EXT')REM(I2)='B5E' > ENDIF > IF (IEND.GT.0) GOTO 50 >40 CONTINUE > GOTO 910 >50 CONTINUE >c NVAL=I2-1 <---- I modified this line as the next line. > NVAL=I2 > CLOSE(LFN001) >30 CONTINUE > > [Mailed from: jjohnson@unavco.ucar.edu (James Johnson)]