anchor.diff:        write(*,*)
        1 '-w WVLNTH  Wavelength of radar signal (cm)'
         write(*,*)
  -     1 '-f IFILE   Input LOS displacement file (lo la uLOS)'
  +     1 '-f IFILE   Input LOS displacement file (stlo stla stdp uLOS)'
         write(*,*)
  -     1 '-o OFILE   Output phase file (lo la phase(rad))'
  +     1 '-o OFILE   Output phase file (stlo stla stdp phase(rad))'
         write (*,*)
        1 '-h         Online help (this message)'
         write (*,*)
  
positive.code: c----
  c Convert line of sight displacement to phase difference from satellite
  c----
     13 read (10,*,end=14) stlo,stla,ulos
            ulos = 100.0d0*ulos
            phase = 2.0d0*tpi*ulos/wvl
            phase = dmod(phase,tpi)
            if (phase.lt.-pi) phase = phase + tpi
            if (phase.gt.pi) phase = phase - tpi
            write (20,*) stlo,stla, phase
            goto 13
     14 continue
  
        END
  
  
negative.code:       IMPLICIT NONE
        REAL*8 pi,d2r,tpi
        PARAMETER (pi=4.0d0*datan(1.0d0),d2r=pi/180.0d0,tpi=2.0d0*pi)
        REAL*8 stlo,stla,stdp,ulos,wvl,phase
        INTEGER narg,i
        CHARACTER*40 ifile,ofile
        LOGICAL ex
  