C
C *** SYNTHETIC WAVES FROM THE HASKELL,S UNILATERAL DISLOCATION SOURCE
C     IN THE 3D LAYERED HALF-SPACE USING THE EXACT GREEN FUNCTION 
C     BASED ON THE RECIPROCITY RELATION ***
C
      PROGRAM MHASKU
C
      PARAMETER(MOM=60,ML=15,MOB=5)
      IMPLICIT REAL*8(A-H,O-Z)
      COMPLEX*16 VP(ML),VS(ML)
     *          ,UX(MOM,MOB),UY(MOM,MOB),UZ(MOM,MOB)
      DIMENSION DNS(ML),THK(ML),ROB(MOB)
      COMMON / CN / C1,N1,C2,N2 /FDAT/ OM,NL,RR
      CHARACTER FNAME*15
C
      WRITE(*,1000)
 1000 FORMAT(' MHASKU: WAVES FROM A HASKELL''S UNILATERAL SOURCE'/
     *       '   IN THE 3D LAYERED HALF-SPACE USING THE EXACT GREEN'/
     *       '   FUNCTION BASED ON THE RECIPLOCITY RELATION '/
     *       '  (SP AND OB ARE IN ANY LAYER),')
C
C *** DATA INPUT ***
C
      CALL DATINP(DOM,NOM,ML,NL1,DNS,VP,VS,THK,LS,SDEP,LO,ODEP
     *           ,MOB,NOB,ROB,RAD)
      NL=NL1
C
      CALL DATHKU(FLL,FLW,STR,DIP,FLD,RST,VR,NVR)
C
C *** COMPUTATION ***
C
      CALL HKUXYZ(DOM,MOM,NOM,ML,DNS,VP,VS,THK
     *           ,LS,SDEP,LO,ODEP,MOB,NOB,ROB,RAD
     *           ,FLL,FLW,STR,DIP,FLD,RST,VR,NVR
     *           ,UX,UY,UZ)
C
C *** OUTPUT OF UX, UY AND UZ ***
C
      CALL DATOUT(DOM,MOM,NOM,NOB,MOB,ROB,RAD,UX,UY,UZ)
C
      STOP
      END
