Bugs in HYPODD_1.0 :	
====================

in src/ph2dt/Makefile:
----------------------
The Makefile for ph2dt refers to a missplaced include file. The correct
place for it is the directory named "include". To fix:
- Remove ph2dt.inc from the directory HYPODD/src/ph2dt/
- Change last line in HYPODD/src/ph2dt/Makefile to read:
        'ph2dt.o : ../../include/ph2dt.inc'
- Recompile ph2dt
Otherwise, you might change one file, re-compile, and have no effect
because of the computer using the other file.

______________________________________________________________________

in util/eqplot.m:
----------------
---line 128:
        for k=1:length(id); if(id(k)==cusp(i));
should read:
        for k=1:length(id); if(id(k)==cusp0(i));

---line 162:
        for k=1:length(id); if(id(k)==cusp(i));
should read
        for k=1:length(id); if(id(k)==cusp0(i));

---in section STATION PLOT:
replace:
  axis([min([slon lon']) max([slon lon']) min([slat lat']) max([slat lat'])]);
  title('STATION MAP'); xlabel('longitude'); ylabel('latitude');
  axis('equal');box('on');
with 
  axis('equal');box('on');
  axis([min([slon lon']) max([slon lon']) min([slat lat']) max([slat lat'])]);
  title('STATION MAP'); xlabel('longitude'); ylabel('latitude');

_____________________________________________________________________________

in src/hypoDD/hypoDD.f:
----------------------

---replace line 466:
      wt_ccs= awt_ccp(i)
with
      wt_ccs= awt_ccs(i)

---replace line 872:
         sc= (mod(real(ev_time(i)),10000)/100) + src_t(i)
with 
         sc= (mod(real(ev_time(i)),10000)/100) - src_t(i)

---replace lines 861-863
            write(*,'("FATAL ERROR (src_t). Please report to ",
     &                "felix@andreas.wr.usgs.gov")')
            stop
with
            write(*,*)'WARNING: org time diff > 5s for ',src_cusp(i)

---replace line 953 
         write(fu3,'(a5,1x,f9.4,1x,f9.4,1x,f9.4,1x,f9.4,1x,i7,1x,
with
         write(fu3,'(a7,1x,f9.4,1x,f9.4,1x,f9.4,1x,f9.4,1x,i7,1x,

____________________________________________________________________________

in src/hypoDD/trialsrc.f:
----------------------
---replace line 52:
            src_cusp(i) = 0
with 
            src_cusp(i) = ev_cusp(i)

____________________________________________________________________________

in src/ph2dt/ph2dt.f:
----------------------
Problem: if last line in dt.ct is a pair with no phases, the ph2dt 
does not go out of the data reading loop.
---replace line 393:
130   read (2,'(a)',end=160) line               ! read header line
with
130   read (2,'(a)',end=159) line               ! read header line

---replace line 415:
160   nobs_ct(i)=k-1
with
159   if(k.eq.1) goto 200
160   nobs_ct(i)=k-1

____________________________________________________________________________

in src/hypoDD/getdata.f:
------------------------
Problem: Uniqueness of event ID in event files (event.dat/event.sel) is not 
checked. If event ID is non-unique in these files, a fatal error may result 
during matrix scaling.  
To fix
---replace line 194 
     &      'or missing in event file!")')
with
     &      'or missing/repeated in event file!")')

---remove line 195,
         write (*,*) 'Missing events:'

---replace line 199,
               if (icusp(i).eq.ev_cusp(j)) k = 1
with
               if (icusp(i).eq.ev_cusp(j)) k = k+1

---and replace line 201,  
            if (k.eq.0) write(*,*) icusp (i)
with the 5 lines
            if (k.eq.0) write(*,*) icusp (i),' is missing.'
            if (k.ge.2) then
                write(*,*) icusp (i),' is non-unique.' 
                stop'Event ID must be unique!'
            endif

____________________________________________________________________________

in src/hypoDD/partials.f:   (04/06/28)
------------------------
Array size definition problem, for layered S-velocity models with more than 
20 layers.

---replace line 42,
        real            vs(20)
with
        real            vs(MAXLAY)

____________________________________________________________________________

in src/hypoDD/lsfit_svd.f:   (04/07/30)
-------------------------
The design matrix g is NOT reset to zero before filling it up with the 
partial derivatives. Since g is a local variable, this shouldn't matter under 
unix, but I am not sure about other platforms. So to be safe, add after 
line 94 "c     Set up full G matrix") the following 6 lines:

c initiate G matrix:
      do i=1,ndt
         do j=1,nev*4
            g(i,j)= 0
         enddo
      enddo

____________________________________________________________________________

10/2004: delaz.f delaz2.f replaced.
____________________________________________________________________________