R/defaultVEL.R

In Rquake: Seismic Hypocenter Determination

defaultVEL<-function(kind=1)
  {
    if(kind==1)
      {
         vel=list()
         vel$'zp'=c(0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,
               20,22,24,26,32,42,50,150)
         vel$'vp'=c(5.47,5.57,5.68,5.76,5.81,5.87,5.92,5.97,6.03,6.08,
               6.13,6.18,6.24,6.29,6.35,
                   6.4,6.45,6.51,6.56,6.62,6.67,6.78,6.89,7.06,7.8,
               7.84,7.87,8.87)
         vel$'ep'=c(0.07,0.07,0.07,0.07,0.07,0.07,0.07,0.07,0.07,
               0.07,0.07,0.07,0.07,0.07,0.07,
                   0.07,0.07,0.07,0.07,0.07,0.07,0.07,0.07,0.07,
               0.07,0.07,0.07,0.07)
         vel$'zs'=c(0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,
               18,19,20,22,24,26,32,42,50,150)
         vel$'vs'=c(3.24,3.29,3.34,3.38,3.41,3.43,3.46,3.48,3.51,
               3.53,3.56,3.59,3.61,3.64,3.66,
                   3.69,3.72,3.74,3.77,3.79,3.82,3.87,3.92,4.01,
               4.41,4.43,4.44,5.14)
         vel$'es'=c(0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,
               0.1,0.1,0.1,0.1,0.1,0.1,
                   0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1)
         vel$'descriptor'=c("##MODELFU  JAPAN step VELOCITY MODEL (lees , 1990)   (LEES, 1988   BO  ",
                   "# P DEPTH   P VEL      PERR      S DEPTH    S VEL    SERR ")
         vel$'name'=c("fuj1.vel.txt")
      }
    if(kind==2)
      {

        vel=list()
        vel$'zp'=c(0,5,10,20,30,50)
        vel$'vp'=c(5.06,5.8,5.95,6.07,6.13,7.61)
        vel$'ep'=c(0,0,0,0,0,0)
        vel$'zs'=c(0,5,10,20,30,50)
        vel$'vs'=c(3.01,3.45,3.54,3.61,3.65,4.53)
        vel$'es'=c(0,0,0,0,0,0)
        vel$'descriptor'=c("# MODEL LITHOS  Seismic Velocity Model for Ecuador Guiller et al, 1996",
              "#  P DEPTH   P VEL      PERR      S DEPTH    S VEL    SERR")
        vel$'name'=c("LITHOS.vel.txt")

      }
    if(kind==3)
      {
        vel=list()
        vel$'zp'=c(0,2.5,5,10,15,20,25,32,42,50,150)
        vel$'vp'=c(5.47,5.73,5.87,6.13,6.4,6.67,6.94,7.8,7.84,7.87,8)
        vel$'ep'=c(0.07,0.07,0.07,0.07,0.07,0.07,0.07,0.07,0.07,0.07,0.07)
        vel$'zs'=c(0,2.5,5,10,15,20,25,32,42,50,150)
        vel$'vs'=c(3.24,3.37,3.43,3.56,3.69,3.82,3.94,4.41,4.43,4.44,6.44)
        vel$'es'=c(0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1)
        vel$'descriptor'=c("##MODEL J2  JAPAN step VELOCITY MODEL (lees , 1990)",
              "# P DEPTH   P VEL      PERR      S DEPTH    S VEL    SERR        ")
        vel$'name'=c("jap1.vel.txt")
      }
    if(kind==4)
      {
        StHelens4.vel=list()
        vel$'zp'=c(0,0.25,1.5,2,2.5,3,3.5,4.75,6,7.5,9,12,16,20,25,32,41)
        vel$'vp'=c(4.6,4.66,4.94,5.05,5.3,5.55,5.8,6.11,6.2,6.35,6.5,6.65,
              6.75,6.84,6.93,7.06,7.64)
        vel$'ep'=c(0.06,0.05,0.03,0.02,0.02,0.02,0.02,0.02,0.01,0.01,0.01,
              0.02,0.04,0.09,0.17,0.3,0.09)
        vel$'zs'=c(0,0.2,1.5,2,2.5,3,3.5,4.8,6,7.5,9,12,16,20,25,32,41)
        vel$'vs'=c(2.66,2.69,2.85,2.91,3.06,3.2,3.35,3.52,3.58,3.66,3.75,3.84,
              3.9,3.95,4,4.08,4.41)
        vel$'es'=c(0.03,0.03,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.02,0.05,0.1,0.17,0.05)
        vel$'descriptor'=c("##MODEL S3HEL:  ST. HELENS VELOCITY MODEL     Steve Mal(LEES, 1988   L3  ",
              "# P DEPTH   P VEL      PERR     S DEPTH    SVEL     SERR                 ")
        vel$'name'=c("StHelens4.vel.txt")
      }
    if(kind==5)
      {
        vel=list()
        vel$'zp'=c(-1,0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20)
        vel$'vp'=c(3.55,3.73,4.47,5.17,5.51,5.6,5.64,5.72,5.81,5.91,6.02,6.14,6.26,
              6.39,6.51,6.63,6.73,6.83,6.91,6.96,6.99,7)
        vel$'ep'=c(0.07,0.07,0.07,0.07,0.07,0.07,0.07,0.07,0.07,0.07,0.07,0.07,0.07,
              0.07,0.07,0.07,0.07,0.07,0.07,0.07,0.07,0.07)
        vel$'zs'=c(-1,0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20)
        vel$'vs'=c(2.05,2.15,2.58,2.99,3.19,3.23,3.26,3.3,3.36,3.41,3.48,3.55,3.62,
              3.69,3.76,3.83,3.89,3.95,3.99,4.02,4.04,4.05)
        vel$'es'=c(0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,
              0.1,0.1,0.1,0.1,0.1,0.1,0.1)
        vel$'descriptor'=c("##", "#")
        vel$'name'=c("montserat.vel.txt")
      }

  if(kind==6)
    {
      vel=list()
      vel$'zp'=c(0,0.5,1,1.5,2,2.5,3,3.5,4,5.5,12,20)
      vel$'vp'=c(4.5,4.51,4.92,4.92,5.46,5.46,5.54,5.54,5.58,5.58,6.05,7.2)
      vel$'ep'=c(0,0,0,0,0,0,0,0,0,0,0,0)
      vel$'zs'=c(0,0.5,1,1.5,2,2.5,3,3.5,4,5.5,12,20)
      vel$'vs'=c(2.43,2.59,2.97,2.97,3.15,3.15,3.27,3.27,3.42,3.42,3.49,4.15)
      vel$'es'=c(0,0,0,0,0,0,0,0,0,0,0,0)
      vel$'descriptor'=c("#MODEL WU COSO REGINAL FINE LAYERS REGIONAL VELOCITY MODEL",
            "#P DEPTH   P VEL      PERR      S DEPTH    S VEL    SERR        ")
      vel$'name'=c("wu_coso.vel.txt")
      
    }
    
    return(vel)
  }