Download CUSVM version 1.3 (272MB)
GeodataQ source code, surfaces,velocity grids, vs30 grids and control file

Summary and General Requirements

The velocity model extension is:

Minimum longitude= -92.9
Maximum longitude= -85.1
Minimum latitude= 33.1
Maximum latitude= 39.9
Minimum depth= 0
Maximum depth= 60e3

However, If you request a value outside this region it will return a value.

The tarball includes all the directories with the files you will need:

How to Un-Tar and Compile

  1. Un-tar the files in DDDIR (or wherever you put the tar file):
     % gunzip cusvm1.3.tar.gz
              % tar -xvf cusvm1.3.tar
        
  2. Make executable the file compile_copy and run it:
     % chmod +x compile_copy
              % ./compile_copy
        
  3. Geodata must be created in GDDIR.

Example, Input, and Output File Format (Points)

There are several ways you can query the velocity model (lines,grids, request depth to surfaces). The most useful for wave propagation simulation grids is to request points. In order to inspect the model we recommend to use the GUI.

The command line for geodataq is:

Usage:

 % geodataq typeofquery velorveldepthorbdrortopo db in rout

where

                  typeofquery:  (0) points or (1) lines or planes
     velorveldepthorbdrortopo:  (0) velocity & density (asl)
                                (1) velocity & density (depth from surface)
                                (2) velocity & density (squeeze to msl != 1)
                                (3) elevation topo
                                (4) elevation bedrock
                                (5) depth to bedrock
                           db:  path to a database
                           in:  input file, query specifications
                         rout:  output path (directory),
                                output file named output.out

The typeofquery for points is (0). The input file format is:

NumberOfStations
                lon_point_1  lat_point_1 elevation_or_depth_point_1
                lon_point_2  lat_point_2 elevation_or_depth_point_2

                For example, the query of two point would be:

                2
                -90.5 35.5 10
                -88.3 37.4 100

In EXAMPLES you should find a file called input_example.in.

The command line to run input_example.in if geodataq or geodataq.exe is in GDDIR:

        % ./geodataq 0 0 GDDIR/DATABASE/ GDDIR/EXAMPLE/input_example.in  GDDIR/

A file named output.out (at GDDIR) is created. The output has the same information as the input and the values of vp,vs and rho in m/s and kg/m^3. The last column indicates the unit where the point is located. The output.out file for input_example.in is:

	-90.500000 35.500000 10.000000 1550.292660 382.796149 1667.000000 3
	-88.300000 37.400000 100.000000 2657.534211 409.354537 2139.454134 7

The damping is left as a parameter to be chosen by the user.

How to Set Up the Matlab GUI for Linux

The GUI provided in this distribution (CUSVM Version 1.3) calls the code "geodataq" and takes advantage of the built in visualization capabilities of Matlab. It is intended to serve as a help to inspect the model.

  1. From GDDIR (directory where you untar and unzip cusvm1.3.tar.gz) run:
     % chmod +x compile_copy
              % ./compile_copy
        
  2. Open matlab and run:

     % cd GDDIR/GUI
              % geodataqgui
        

    gui

  3. Choose the DATABASE path (Init Database).

    gui_dbpath

  4. Add references.

    gui_gpath

  5. Choose the type of query and the property to be displayed.

    gui_queryexample

References


Use of this model should reference:

Contact: Leonardo Ramirez-Guzman