diff --git a/example/.gmtcommands4 b/example/.gmtcommands4 deleted file mode 100644 index e3d01fa..0000000 --- a/example/.gmtcommands4 +++ /dev/null @@ -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 - diff --git a/example/.gmtdefaults4 b/example/.gmtdefaults4 deleted file mode 100644 index 2ecdffe..0000000 --- a/example/.gmtdefaults4 +++ /dev/null @@ -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 diff --git a/example/DSurfTomo.in b/example/DSurfTomo.in index a7b6c89..bce45ef 100644 --- a/example/DSurfTomo.in +++ b/example/DSurfTomo.in @@ -1,21 +1,21 @@ 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 -4.0 0.0 c: weight damp -3 c: nsublayer (numbers of sublayers for each grid interval:grid --> layer) -0.5 2.8 c: minimum velocity, maximum velocity -10 c: maxiter (iteration number) -0.1 c: sparsity fraction +20 c: max(sources, receivers) +4.0 0.1 c: weight damp +3 c: sablayers (for computing depth kernel, 2~5) +0.5 2.8 c: minimum velocity, maximum velocity (a priori information) +10 c: maximum iteration +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 0 c: kmaxLc -0 c: kmaxLg +0 c: kmaxLg 0 c: synthetic flag(0:real data,1:synthetic) -0.02 c: noiselevel -0.5 c: threshold +0.02 c: noiselevel +0.5 c: threshold diff --git a/example/plotslicenewVs.gmt b/example/plotslicenewVs.gmt deleted file mode 100644 index 8f4ff15..0000000 --- a/example/plotslicenewVs.gmt +++ /dev/null @@ -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 & diff --git a/src/DSurfTomo b/src/DSurfTomo deleted file mode 100755 index dd8cd53..0000000 Binary files a/src/DSurfTomo and /dev/null differ diff --git a/src/main.f90 b/src/main.f90 index 8ab9d53..57c3b4c 100644 --- a/src/main.f90 +++ b/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(*,*)