V1.3 release

Fix some bugs in jointly invert both phase and group velocity

example/.gmtcommands4

@@ -1,10 +0,0 @@
-# GMT common arguments shelf
--Ba0.2f0.1:S velocity (km/s):
--JM2i
--R121.35/121.605/24.975/25.2
--x2.3i
--Y5i
--y-2.7i
--jM2i
-EOF
-

example/.gmtdefaults4

@@ -1,107 +0,0 @@
-#
-#	GMT-SYSTEM 4.5.16 [64-bit] Defaults file
-#
-#-------- Plot Media Parameters -------------
-PAGE_COLOR		= white
-PAGE_ORIENTATION	= landscape
-PAPER_MEDIA		= a4
-#-------- Basemap Annotation Parameters ------
-ANNOT_MIN_ANGLE		= 20
-ANNOT_MIN_SPACING	= 0
-ANNOT_FONT_PRIMARY	= Helvetica
-ANNOT_FONT_SIZE_PRIMARY	= 6p
-ANNOT_OFFSET_PRIMARY	= 0.2c
-ANNOT_FONT_SECONDARY	= Helvetica
-ANNOT_FONT_SIZE_SECONDARY	= 16p
-ANNOT_OFFSET_SECONDARY	= 0.2c
-DEGREE_SYMBOL		= ring
-HEADER_FONT		= Helvetica
-HEADER_FONT_SIZE	= 36p
-HEADER_OFFSET		= 0.5c
-LABEL_FONT		= Helvetica
-LABEL_FONT_SIZE		= 5p
-LABEL_OFFSET		= 0.1c
-OBLIQUE_ANNOTATION	= 1
-PLOT_CLOCK_FORMAT	= hh:mm:ss
-PLOT_DATE_FORMAT	= yyyy-mm-dd
-PLOT_DEGREE_FORMAT	= ddd:mm:ss
-Y_AXIS_TYPE		= hor_text
-#-------- Basemap Layout Parameters ---------
-BASEMAP_AXES		= WESN
-BASEMAP_FRAME_RGB	= black
-BASEMAP_TYPE		= fancy
-FRAME_PEN		= 1.25p
-FRAME_WIDTH		= 0.1c
-GRID_CROSS_SIZE_PRIMARY	= 0c
-GRID_PEN_PRIMARY	= 0.25p
-GRID_CROSS_SIZE_SECONDARY	= 0c
-GRID_PEN_SECONDARY	= 0.5p
-MAP_SCALE_HEIGHT	= 0.2c
-POLAR_CAP		= 85/90
-TICK_LENGTH		= 0.1c
-TICK_PEN		= 0.3p
-X_AXIS_LENGTH		= 25c
-Y_AXIS_LENGTH		= 15c
-X_ORIGIN		= 2.5c
-Y_ORIGIN		= 2.5c
-UNIX_TIME		= FALSE
-UNIX_TIME_POS		= BL/-2c/-2c
-UNIX_TIME_FORMAT	= %Y %b %d %H:%M:%S
-#-------- Color System Parameters -----------
-COLOR_BACKGROUND	= black
-COLOR_FOREGROUND	= white
-COLOR_NAN		= 128
-COLOR_IMAGE		= adobe
-COLOR_MODEL		= rgb
-HSV_MIN_SATURATION	= 1
-HSV_MAX_SATURATION	= 0.1
-HSV_MIN_VALUE		= 0.3
-HSV_MAX_VALUE		= 1
-#-------- PostScript Parameters -------------
-CHAR_ENCODING		= ISOLatin1+
-DOTS_PR_INCH		= 300
-GLOBAL_X_SCALE		= 1
-GLOBAL_Y_SCALE		= 1
-N_COPIES		= 1
-PS_COLOR		= rgb
-PS_IMAGE_COMPRESS	= lzw
-PS_IMAGE_FORMAT		= ascii
-PS_LINE_CAP		= butt
-PS_LINE_JOIN		= miter
-PS_MITER_LIMIT		= 35
-PS_VERBOSE		= FALSE
-TRANSPARENCY		= 0
-#-------- I/O Format Parameters -------------
-D_FORMAT		= %.12lg
-FIELD_DELIMITER		= tab
-GRIDFILE_FORMAT		= nf
-GRIDFILE_SHORTHAND	= FALSE
-INPUT_CLOCK_FORMAT	= hh:mm:ss
-INPUT_DATE_FORMAT	= yyyy-mm-dd
-IO_HEADER		= FALSE
-N_HEADER_RECS		= 1
-NAN_RECORDS		= pass
-OUTPUT_CLOCK_FORMAT	= hh:mm:ss
-OUTPUT_DATE_FORMAT	= yyyy-mm-dd
-OUTPUT_DEGREE_FORMAT	= D
-XY_TOGGLE		= FALSE
-#-------- Projection Parameters -------------
-ELLIPSOID		= WGS-84
-MAP_SCALE_FACTOR	= default
-MEASURE_UNIT		= cm
-#-------- Calendar/Time Parameters ----------
-TIME_FORMAT_PRIMARY	= full
-TIME_FORMAT_SECONDARY	= full
-TIME_EPOCH		= 2000-01-01T12:00:00
-TIME_IS_INTERVAL	= OFF
-TIME_INTERVAL_FRACTION	= 0.5
-TIME_LANGUAGE		= us
-TIME_UNIT		= d
-TIME_WEEK_START		= Sunday
-Y2K_OFFSET_YEAR		= 1950
-#-------- Miscellaneous Parameters ----------
-HISTORY			= TRUE
-INTERPOLANT		= akima
-LINE_STEP		= 0.025c
-VECTOR_SHAPE		= 0
-VERBOSE			= FALSE

example/DSurfTomo.in

@@ -1,16 +1,16 @@
 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
 c INPUT PARAMETERS
 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
-surfdata.dat			 c: data file 
+surfdataTB.dat                   c: data file
 18 18 9                          c: nx ny nz (grid number in lat lon and depth direction)
 25.2  121.35                     c: goxd gozd (upper left point,[lat,lon])
 0.015 0.017                      c: dvxd dvzd (grid interval in lat and lon direction)
-20                              c: nsrc*maxf 
+20                               c: max(sources, receivers)
-4.0  0.0                         c: weight damp
+4.0  0.1                         c: weight damp
-3                                c: nsublayer (numbers of sublayers for each grid interval:grid --> layer)
+3                                c: sablayers (for computing depth kernel, 2~5)
-0.5 2.8                          c: minimum velocity, maximum velocity
+0.5 2.8                          c: minimum velocity, maximum velocity (a priori information)
-10                               c: maxiter (iteration number)
+10                               c: maximum iteration
-0.1                              c: sparsity fraction
+0.2                              c: sparsity fraction 
 26                               c: kmaxRc (followed by periods)
 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3.0
 0                                c: kmaxRg

example/plotslicenewVs.gmt

@@ -1,49 +0,0 @@
-#!/bin/csh
-
-#-----------------------------------------------------------
-# 2015-06-08 Hongjian Fang
-# example:csh plotslicenewVs.gmt DSurfTomo.inMeasure.dat 0.2 0.4 0.8 1.1
-#-----------------------------------------------------------
-
-gmtset ANOT_FONT_SIZE 6
-gmtset LABEL_FONT_SIZE 6
-gmtset FRAME_WIDTH=0.1c
-gmtset LABEL_OFFSET=0.1c
-gmtset LABEL_FONT_SIZE=5p
-gmtset TICK_LENGTH=0.1c
-gmtset TICK_PEN=0.3p
-set inp3D  = $1
-set ps   = horizVs.ps
-set J      = -JM2i     #size for plot
-set cpt = slice.cpt
-
-# start
-makecpt -Cseis -T0.6/1.5/0.1   > $cpt #velocity boundary
-set RMAP = `minmax -C $inp3D | awk '{print "-R"$1"/"$2"/"$3"/"$4}'`
-psbasemap $RMAP $J -Ba0.1f0.05WseN -P -K -Y5i >$ps
-awk '{if($3==depth1) print $1,$2,$4}' depth1=$2 $inp3D|xyz2grd -R -I0.017/0.015 -Gtmp.grd 
-grdimage tmp.grd $J $RMAP -BSenW -E100 -C$cpt -O -K >> $ps
-rm -rf tmp.grd
-psscale -Cslice.cpt -Ba0.1f0.05:'S velocity (km/s)': -D1.0i/-0.15i/3.00/0.2h  -O -K -P >> $ps
-
-makecpt -Cseis -T0.8/1.8/0.1   > $cpt #velocity boundary
-psbasemap $RMAP $J -Ba0.1f0.05NwsE -P -K -O -X2.3i >>$ps
-awk '{if($3==depth2) print $1,$2,$4}' depth2=$3 $inp3D|xyz2grd -R -I0.017/0.015 -Gtmp.grd 
-grdimage tmp.grd $J $RMAP -BSenW -E100 -C$cpt -O -K >> $ps
-psscale -Cslice.cpt -Ba0.2f0.1:'S velocity (km/s)': -D1.0i/-0.15i/3.00/0.2h  -O -K -P >> $ps
-
-makecpt -Cseis -T1.1/2.0/0.1  > $cpt #velocity boundary
-psbasemap $RMAP $J -Ba0.1f0.05WneS -P -K -O -Y-2.7i -X-2.3i >>$ps
-awk '{if($3==depth2) print $1,$2,$4}' depth2=$4 $inp3D|xyz2grd -R -I0.017/0.015 -Gtmp.grd 
-grdimage tmp.grd $J $RMAP -BSenW -E100 -C$cpt -O -K >> $ps
-psscale -Cslice.cpt -Ba0.2f0.1:'S velocity (km/s)': -D1.0i/-0.25i/3.00/0.2h  -O -K -P >> $ps
-
-makecpt -Cseis -T1.3/2.4/0.1   > $cpt #velocity boundary
-psbasemap $RMAP $J -Ba0.1f0.05SwnE -P -K -O -X2.3i >>$ps
-awk '{if($3==depth2) print $1,$2,$4}' depth2=$5 $inp3D|xyz2grd -R -I0.017/0.015 -Gtmp.grd 
-grdimage tmp.grd $J $RMAP -BSenW -E100 -C$cpt -O -K >> $ps
-psscale -Cslice.cpt -Ba0.2f0.1:'S velocity (km/s)': -D1.0i/-0.25i/3.00/0.2h  -O  -P >> $ps
-
-rm *cpt *grd
-ps2eps -f $ps
-evince *eps &

src/main.f90

@@ -121,7 +121,7 @@ program SurfTomo
 
   ! OUTPUT PROGRAM INFOMATION            
   write(*,*)
-  write(*,*),'                         S U R F  T O M O'
+  write(*,*),'                             DSurfTomo (v1.3)'
   write(*,*),'PLEASE contact Hongjain Fang &
     (fanghj@mail.ustc.edu.cn) if you find any bug'
   write(*,*)