anchor.diff:  !----
         if (model.eq.'WELLSCOPP') then
             call wellscopp(mom,mag)
  +      elseif (model.eq.'MAIBEROZA') then
  +          call maiberoza(mom,mag)
  +      elseif (model.eq.'ALLENHAYES') then
  +          call allenhayes(mom,mag)
         else
             call usage('!! Error: no relationship named '//trim(model))
         endif
  
positive.code:       b(2) = 0.08d0   ! 0.21
        b(3) = 0.63d0   ! 0.24
        b(4) = 0.69d0   ! 0.08
        write(*,1008)
        do 106 i = 1,4
            avgslp = 10**(a(i)+b(i)*mag)
            write(*,1002) typ(i),avgslp
    106 continue
   1001 format('Surface rupture length:')
   1002 format(4X,A6,':',F12.3,X,'m')
   1003 format('Total subsurface rupture length:')
   1004 format('Downdip rupture width:')
   1005 format('Total rupture area:')
   1006 format(4X,A6,':',F12.3,X,'m^2')
   1007 format('Maximum slip:')
   1008 format('Average slip:')
        RETURN
        END
  
  
negative.code: !      write(0,*)
  !     1 '                       HANKSBAKUN (Hanks and Bakun, 2008)'
        write(0,*)
       1 '                       *BLASERETAL (Blaser et al., 2010): ',
       2                         'global subduction zone events (Mw 4.8-',
       3                         '9.5)'
        write(0,*)
       1 '                       *STRASSERETAL (Strasser et al., 2010): ',
       2                         'global subduction zone events (Mw 6.0-',
       3                         '9.5)'
        write(0,*)
       1 '                       *YENMA (Yen and Ma, 2011): collision ',
       2                         'zone (Taiwan) events (Mw 4.6-8.9)'
  !      write(0,*)
  !     1 '                       STIRLINGETAL (Stirling et al., 2013)'
        write(0,*)
       1 '                       ALLENHAYES (Allen and Hayes, 2017): ',
       2                         'global finite fault models of Mw 7.1-',
       3                         '9.2 subduction interface events'
        write(0,*)
       1 '                       * not yet implemented'
        STOP
        END
  