赵志远的博客

Transcription

1

2 4 TauP TauP GMT GMT GMT GMT GMT bash bat CAP Cap SAC SAC rms GAMIT GAMIT Linux distaz Fetch Script Hinet CENC SOD Anaconda ObsPy ObsPy

3 2 50 bat bat bat Windows Unix bat ftp C Windows C Java Java Matlab Matlab Matlab cpt Perl Perl gitalk GitHub Pages HTTPS texlive pandoc cpdf pdf convert Atom guake Source Code Pro Centos git clone git apk

4 CDO IP

5 4

6 TauP 1.1 TauP Taup : : Java : Java Taup Java 1.8 Windows Taup Taup zip CMD 1 > jar -xvf Taup zip PATH C:\programs\TauP-2.4.4\bin; CMD 1 > taup Linux Taup 1 $ wget 2 $ tar -xvf TauP tgz 3 $ sudo mv TauP /opt 4 $ echo 'export TAUPHOME=/opt/TauP-2.4.4' >> ~/.bashrc 5 $ echo 'export PATH=${TAUPHOME}/bin:${PATH}' >> ~/.bashrc 6 $ source ~/.bashrc 7 $ taup # 5

7 GMT 2.1 GMT GMT GMT4 GMT4 GMT4 GMT GMT4 GMT4 gmtdefaults GMT JX6c/-3c JX6/-3 -K -O > >> bash GMT 1 # 2 gmtset ANNOT_FONT_PRIMARY 35 ANNOT_FONT_SIZE_PRIMARY 12 3 gmtset HEADER_FONT 35 HEADER_FONT_SIZE 18 4 gmtset PLOT_DEGREE_FORMAT ddd.xxf 5 # 6 R=108/125/33/43 7 J=M6i 8 B=3/5 9 PS=HuaBei.ps 10 # 11 ## 12 psxy -R$R -J$J -B$B -T -K > $PS 13 ## 14 pscoast -R -J -B -W0.5p,black -Sblue -K -O >> $PS 15 psxy fault.dat -R -J -W0.5p,black -m -K -O >> $PS 16 psxy city.dat -R -J -Sc0.2c -W0.5p,black -K -O >> $PS 17 pstext capital.dat -R -J -Wwhite -D0/-0.2c -K -O >> $PS 18 ## 19 psxy -R -J -T -O >> $PS 6

8 GMT 7 20 ## jpg Windows -C-sFONTPATH=C:\ Windows\fonts 21 ## GMT4 bug -A 22 ps2raster $PS -A -P -Tj -E300 -R -W GMT4 2.2 GMT4 GMT GMT4 Linux : ftp://ftp.soest.hawaii.edu/gmt/ : Ubuntu $ sudo apt-get update 2 $ sudo apt-get install gcc g++ make libc6 # 3 $ sudo apt-get install libnetcdf-dev libgdal-dev python-gdal # netcdf 4 $ sudo apt-get install libxaw7-dev # X 5 $ sudo apt-get install libice-dev libsm-dev libx11-dev 6 $ sudo apt-get install libxext-dev libxmu-dev libxt-dev Centos7 1 $ sudo yum install gcc gcc-c++ make glibc # 2 $ sudo yum install netcdf netcdf-devel gdal gdal-devel gdal-python # netcdf 3 $ sudo yum install libxaw-devel # X 4 $ sudo yum install libice-devel libsm-devel libx11-devel 5 $ sudo yum install libxext-devel libxmu-devel libxt-devel gmt gshhg 1 $ wget 2 $ wget 3 $ md5sum gmt src.tar.bz2 gshhg-gmt tar.gz # d3d8a49e5edb4f79c8f2e2c58b gmt src.tar.bz2 5 8ee2653f9daf84d49fefbf990bbfa1e7 gshhg-gmt tar.gz

9 GMT 8 gmt 1 $ tar -xvf gmt src.tar.bz2 2 $ cd gmt $./configure --prefix=/opt/gmt # gmt 4 $ make 5 $ sudo make install-all # install-all install 1 $ cd../ 2 $ tar -xvf gshhg-gmt tar.gz 3 $ sudo mv gshhg-gmt /opt/gmt /share/coast 1 $ echo 'export GMT4HOME=/opt/GMT ' >> ~/.bashrc 2 $ echo 'export PATH=${GMT4HOME}/bin:$PATH'>> ~/.bashrc 3 $ echo 'export LD_LIBRARY_PATH=${LD_LIBRARY_PATH}:${GMT4HOME}/lib64'>> ~/.bashrc # lib64 --> lib 4 $ exec $SHELL -l 1 $ psxy - 2 psxy [64-bit] - Plot lines, polygons, and symbols on maps gs Ubuntu $ sudo apt-get install poppler-data 2 $ sudo apt-get install fonts-arphic-uming fonts-arphic-ukai # gs Linux 3 $ sudo mkdir /usr/share/fonts/winfonts/ # winfonts 4 $ sudo cp simhei.ttf simkai.ttf simsun.ttc simfang.ttf /usr/share/fonts/winfonts # Windows winfonts 5 $ sudo gedit /etc/ghostscript/cidfmap.d/90gs-cjk-resource-gb1.conf # gs 6 7 % STSong-Light 8 /BousungEG-Light-GB <</FileType /TrueType /Path (/usr/share/fonts/truetype/arphic/ uming.ttc) /SubfontId 0 /CSI [(GB1) 4] >> ; 9 /GBZenKai-Medium <</FileType /TrueType /Path (/usr/share/fonts/truetype/arphic/ ukai.ttc) /SubfontId 0 /CSI [(GB1) 4] >> ; 10 /Song-Medium /GBZenKai-Medium ; 11 /Adobe-GB1 /BousungEG-Light-GB ;

10 GMT 9 12 /Adobe-GB1-Bold /GBZenKai-Medium ; % Windows 15 /STSong-Light <</FileType /TrueType /Path (/usr/share/fonts/winfonts/simsun.ttc) / SubfontId 0 /CSI [(GB1) 4] >> ; 16 /STFangsong-Light <</FileType /TrueType /Path (/usr/share/fonts/winfonts/simfang. ttf) /SubfontId 0 /CSI [(GB1) 4] >> ; 17 /STHeiti-Regular <</FileType /TrueType /Path (/usr/share/fonts/winfonts/simhei.ttf ) /SubfontId 0 /CSI [(GB1) 4] >> ; 18 /STKaiti-Regular <</FileType /TrueType /Path (/usr/share/fonts/winfonts/simkai.ttf ) /SubfontId 0 /CSI [(GB1) 4] >> ; $ sudo update-gsfontmap Centos7 1 $ sudo yum install ghostscript-chinese-zh_cn # gs 2 $ sudo mkdir /usr/share/fonts/winfonts/ # winfonts 3 $ sudo cp simhei.ttf simkai.ttf simsun.ttc simfang.ttf /usr/share/fonts/winfonts # Windows winfonts 4 $ sudo gedit /usr/share/ghostscript/conf.d/cidfmap.zh_cn # gs Windows 5 6 % 7 /BousungEG-Light-GB <</FileType /TrueType /Path (/usr/share/fonts/wqy-zenhei/wqyzenhei.ttc) /SubfontId 0 /CSI [(GB1) 4] >> ; 8 /GBZenKai-Medium <</FileType /TrueType /Path (/usr/share/fonts/wqy-zenhei/wqyzenhei.ttc) /SubfontId 0 /CSI [(GB1) 4] >> ; 9 /MSungGBK-Light /BousungEG-Light-GB ; 10 /Adobe-GB1 /BousungEG-Light-GB ; % Windows 13 /STSong-Light <</FileType /TrueType /Path (/usr/share/fonts/winfonts/simsun.ttc) / SubfontId 0 /CSI [(GB1) 4] >> ; 14 /STFangsong-Light <</FileType /TrueType /Path (/usr/share/fonts/winfonts/simfang. ttf) /SubfontId 0 /CSI [(GB1) 4] >> ; 15 /STHeiti-Regular <</FileType /TrueType /Path (/usr/share/fonts/winfonts/simhei.ttf ) /SubfontId 0 /CSI [(GB1) 4] >> ; 16 /STKaiti-Regular <</FileType /TrueType /Path (/usr/share/fonts/winfonts/simkai.ttf ) /SubfontId 0 /CSI [(GB1) 4] >> ; gmt 1 $ gedit /opt/gmt /share/pslib/ps_font_info.d 2 3 # 4 STSong-Light--UniGB-UTF8-H STFangsong-Light--UniGB-UTF8-H

11 GMT 10 6 STHeiti-Regular--UniGB-UTF8-H STKaiti-Regular--UniGB-UTF8-H $ pstext -L # gmt STSong-Light--UniGB-UTF8-H STFangsong-Light--UniGB-UTF8-H STHeiti-Regular--UniGB-UTF8-H STKaiti-Regular--UniGB-UTF8-H 1 $ gedit gmt4_cn.sh 2 3 #!/bin/bash 4 gmtset HEADER_FONT 35 5 pstext -R0/10/0/3 -JX15c/3c -B1/1:."GMT ": -P > cn.ps <<EOF LM GMT LM GMT LM GMT LM GMT 10 EOF 11 rm.gmt* $ bash gmt4_cn.sh # 14 $ gs cn.ps # ps GMT4 Windows C:\programs\ gs919w32.exe gsv50w32.exe gmt win32.exe gmt pdf-win32.exe gshhg win32.exe gsview name : Registered S/N : gmt PS_font_info.d C:\programs\gmt4\share\pslib\

12 GMT 11 gsview gsview -> options -> advanced configure... -> Ghostscript Include Path ;C:\windows\fonts; Unixtools exe gmt4\bin 2.3 GMT4 gmtdefaults GMT 1 $ gmtdefautls -D # 2 $ gmtdefautls -L # : 1 $ gmtdefaults -L grep PLOT 2 PLOT_CLOCK_FORMAT = hh:mm:ss 3 PLOT_DATE_FORMAT = yyyy-mm-dd 4 PLOT_DEGREE_FORMAT = ddd:mm:ss GMT4 1 $ gmtdefaults -L 2 # 3 # Plot Media Parameters PAGE_COLOR = white # 5 PAGE_ORIENTATION = landscape # portrait 6 PAPER_MEDIA = a4 # 7 # Basemap Annotation Parameters ANNOT_MIN_ANGLE = 20 # 9 ANNOT_MIN_SPACING = 0 # 10 ANNOT_FONT_PRIMARY = Helvetica # 11 ANNOT_FONT_SIZE_PRIMARY = 14p # 12 ANNOT_OFFSET_PRIMARY = 0.2c # 13 ANNOT_FONT_SECONDARY = Helvetica # 14 ANNOT_FONT_SIZE_SECONDARY = 16p # 15 ANNOT_OFFSET_SECONDARY = 0.2c # 16 DEGREE_SYMBOL = ring # 17 HEADER_FONT = Helvetica # 18 HEADER_FONT_SIZE = 36p # 19 HEADER_OFFSET = 0.5c # 20 LABEL_FONT = Helvetica # 21 LABEL_FONT_SIZE = 24p # 22 LABEL_OFFSET = 0.3c # 23 OBLIQUE_ANNOTATION = 1 # 24 PLOT_CLOCK_FORMAT = hh:mm:ss # 25 PLOT_DATE_FORMAT = yyyy-mm-dd #

13 GMT PLOT_DEGREE_FORMAT = ddd:mm:ss # ddd.xf 27 Y_AXIS_TYPE = hor_text # Y ver_text 28 # Basemap Layout Parameters BASEMAP_AXES = WESN # WesN 30 BASEMAP_FRAME_RGB = black # 31 BASEMAP_TYPE = fancy # inside graph fancy plain 32 FRAME_PEN = 1.25p # plain 33 FRAME_WIDTH = 0.2c # fancy 34 GRID_CROSS_SIZE_PRIMARY = 0c # 0 35 GRID_PEN_PRIMARY = 0.25p # 36 GRID_CROSS_SIZE_SECONDARY = 0c # 0 37 GRID_PEN_SECONDARY = 0.5p # 38 MAP_SCALE_HEIGHT = 0.2c # 39 POLAR_CAP = 85/90 # 40 TICK_LENGTH = 0.2c # -JX6i/-3i Y 41 TICK_PEN = 0.5p # 42 X_AXIS_LENGTH = 25c # x 43 Y_AXIS_LENGTH = 15c # y 44 X_ORIGIN = 2.5c # 45 Y_ORIGIN = 2.5c # 46 UNIX_TIME = FALSE # -U 47 UNIX_TIME_POS = BL/-2c/-2c # 48 UNIX_TIME_FORMAT = %Y %b %d %H:%M:%S # 49 # Color System Parameters COLOR_BACKGROUND = black # Z cpt 51 COLOR_FOREGROUND = white # Z cpt 52 COLOR_NAN = 128 # Z NAN 53 COLOR_IMAGE = adobe # PS 54 COLOR_MODEL = rgb # cpt 55 HSV_MIN_SATURATION = 1 # S 56 HSV_MAX_SATURATION = 0.1 # S 57 HSV_MIN_VALUE = 0.3 # V 58 HSV_MAX_VALUE = 1 # V 59 # PostScript Parameters CHAR_ENCODING = ISOLatin1+ # 61 DOTS_PR_INCH = 300 # 62 GLOBAL_X_SCALE = 1 # x 63 GLOBAL_Y_SCALE = 1 # y 64 N_COPIES = 1 # 65 PS_COLOR = rgb # PS 66 PS_IMAGE_COMPRESS = lzw # PS 67 PS_IMAGE_FORMAT = ascii # PS 68 PS_LINE_CAP = butt # butt round square 69 PS_LINE_JOIN = miter # 70 PS_MITER_LIMIT = 35 # miter 71 PS_VERBOSE = FALSE # PS

14 GMT TRANSPARENCY = 0 # 73 # I/O Format Parameters D_FORMAT = %.12lg # 75 FIELD_DELIMITER = tab # GMT 76 GRIDFILE_FORMAT = nf # 77 GRIDFILE_SHORTHAND = FALSE # 78 INPUT_CLOCK_FORMAT = hh:mm:ss # 79 INPUT_DATE_FORMAT = yyyy-mm-dd # 80 IO_HEADER = FALSE # 81 N_HEADER_RECS = 1 # 1 82 NAN_RECORDS = pass # NAN 83 OUTPUT_CLOCK_FORMAT = hh:mm:ss # 84 OUTPUT_DATE_FORMAT = yyyy-mm-dd # 85 OUTPUT_DEGREE_FORMAT = D # 86 XY_TOGGLE = FALSE # 87 # Projection Parameters ELLIPSOID = WGS-84 # 89 MAP_SCALE_FACTOR = default # 90 MEASURE_UNIT = cm # 91 # Calendar/Time Parameters TIME_FORMAT_PRIMARY = full # full January abbrevaite Jan character J 93 TIME_FORMAT_SECONDARY = full # 94 TIME_EPOCH = T12:00:00 # 95 TIME_IS_INTERVAL = OFF # 96 TIME_INTERVAL_FRACTION = 0.5 # 97 TIME_LANGUAGE = us # 98 TIME_UNIT = d # 99 TIME_WEEK_START = Sunday # 100 Y2K_OFFSET_YEAR = 1950 # # Miscellaneous Parameters HISTORY = TRUE # 103 INTERPOLANT = akima # 104 LINE_STEP = 0.025c # 105 VECTOR_SHAPE = 0 # v 106 VERBOSE = FALSE # -V 2.4 GMT4 1. D:\MapSIS\A13 _ \ A _ _.ORA 2. 3.

15 GMT bat 1 rem 2 gawk -F ^ "{if (NR>2) print $16,$17,$6}" A _ _.ORA gawk "!a[$0]++" >.txt 3 rem 4 gawk "{if ($1>=109 && $1<=121 && $2>=33 && $2 <=43) print $1,$2,$3}".txt >.txt txt Windows gmt gmt gmt set MAP_TICK_LENGTH_PRIMARY 2p 2 gmt set FONT_LABEL 10p 3 set PS= -colortopo.ps 4 set R=109/121/33/43 5 gmt grdcut earth_relief_01m.grd -R%R% -GcutTopo.grd 6 gmt grd2cpt cuttopo.grd -Cglobe -S-6800/8200/10 -D > colortopo.cpt 7 gmt psbasemap -R%R% -JM5i -B2 -BWesN -K > %PS% 8 gmt grdimage cuttopo.grd -R -J -CcolorTopo.cpt -O -K >>%PS% 9 gmt psxy -.txt -R -J -St0.2c -Gblack -K -O >> %PS% 10 gmt pstext -R -J -F+f9p,39+jBC -Dj0/0.2c -.txt -K -O >> %PS% 11 gmt psscale -Bx3000+u"m" -By+l"Topo" -CcolorTopo.cpt -Dx1.3c/0+w4i/0.4c+o0/-1c+h - O --FONT_ANNOT_PRIMARY=10p >> %PS% 12 gmt psconvert %PS% -P -A -Tg -Z -C-sFONTPATH=C:/windows/fonts 13 del gmt* 1: 850m.cpt 0 193/228/ /228/ /228/ /241/ /241/ /109/ /108/ /226/143 2: m

16 GMT : 1

17 GMT : 2

18 GMT bash bat GMT GMT Linux bash bash bash Windows bat bash bat 1. ## => rem 2. PS=MAP.PS => set PS=MAP.PS 3. $PS => %PS% 4. rm gmt.* => del gmt.* Notepad++ Ctrl+F ## rem $PS $R $J %PS% %R% %J% bat bash 1 T1=$(gmt grdinfo world6m.grd -T1000) 2 gmt makecpt -Crainbow $T1 > 1.cpt bat for 1 for /f %%i in ('gmt grdinfo world6m.grd -T1000') do set T1=%%i 2 gmt makecpt -Crainbow %T1% > 1.cpt bash 1 # EOF gmt psxy 2 gmt psxy -R -J -Sc0.5c -Gred -W0.5p,black -K -O << EOF >> $PS EOF bat 1 echo > tmp 2 echo >> tmp 3 echo >> tmp 4 gmt psxy tmp -R -J -Sc0.5c -Gred -W0.5p,black -O -K >> %PS% bash \ 1 awk '$4>=0 && $4<70 {print $5,$4,$3*0.03}' $nearloc \ 2 gmt psxy -R -J -Sc -Gpurple -W0.1p,black -O -K >> $PS bat ^ 1 gawk "$4>=0 && $4<70 {print $5,$4,$3*0.03}" %nearloc% gmt psxy -R -J -Sc - Gpurple -W0.1p,black -O -K >> %PS% bat ^ > ^

19 GMT 18 1 rem ^ 2 echo S 0.1i c 0.20 red 0.1p,black 0.25i ^>300 km > tmp grep cut wc cat UnixTools GMT bin bash awk gawk bash 1 #!bin/bash 2 PS=example.ps 3 # 4 R=-R108/125/33/43 5 # 6 J=-JM6i 7 # 8 gmt psbasemap $R $J -B3 -K > $PS 9 T1=$(gmt grdinfo huabei.grd -T100) 10 gmt makecpt -Ctopo $T1 > huabei.cpt 11 # 12 gmt grdimage -R -J huabei.grd -Chuabei.cpt -K -O >> $PS 13 # 14 gmt psxy -R -J -St0.5c -Gblack -K -O >> $PS << EOF EOF # event.dat awk '{if (NR>1) print $4,$3}' event.dat gmt psxy -R -J -Sc0.4c -Gblack -O >> $PS 22 # 23 rm gmt.* huabei.cpt bat 1 set PS=example.ps 2 rem 3 set R=-R108/125/33/43 4 rem 5 set J=-JM6i 6 rem 7 gmt psbasemap %R% %J% -B3 -K > %PS% 8 for /f %%i in ('gmt grdinfo huabei.grd -T100') do set T1=%%i 9 gmt makecpt -Ctopo %T1% > huabei.cpt 10 rem 11 gmt grdimage -R -J huabei.grd -Chuabei.cpt -K -O >> %PS% 12 rem 13 echo > tmp 14 echo >> tmp 15 echo >> tmp 16 gmt psxy tmp -R -J -St0.5c -Gblack -K -O >> %PS%

20 GMT rem event.dat gawk "{if (NR>1) print $4,$3}" event.dat gmt psxy -R -J -Sc0.4c -Gblack -O >> % PS% 20 rem 21 del gmt.* huabei.cpt tmp gmt4 gmt5 1 gmt set MAP_FRAME_TYPE plain 2 gmt set MAP_TICK_LENGTH -0.1c 3 gmt set FONT_ANNOT_PRIMARY=8p 4 gmt set FORMAT_GEO_MAP ddd 5 gmt set MAP_FRAME_PEN 0.03c 6 set PS=seismicray.ps 7 set file=earthquake.txt 8 rem 9 gawk "{if (NR>1) print $9,$8}" resultsn.dat > %file% 10 rem 10cm 11 gmt psbasemap -R115.5/120/38/41.5 -JM10c -Ba1 -BWS -K > %PS% 12 gmt psbasemap -R -J -B0 -Bne -K -O >> %PS% 13 rem gmt psxy -Fr 14 rem 15 for /f "tokens=1,2" %%i in (tstz.txt) do ( 16 gmt psxy %file% -R -J -W0.3p -Fr%%i/%%j -O -K >> %PS% 17 ) 18 gmt psxy %file% -R -J -Sc0.1c -Gblue -Wblack -K -O >> %PS% 19 gmt psxy tstz.txt -R -J -St0.25c -Gred -Wblack -K -O >> %PS% 20 echo H > legend.txt 21 echo S 0.2c t 0.25c red >> legend.txt 22 echo S 0.2c c 0.15c blue >> legend.txt 23 gmt pslegend legend.txt -R -J -F+p+gwhite -DjBL+w1.5c -L1.5 -O >> %PS% 24 gmt psconvert %PS% -Tg -A -P -C-sFONTPATH=C:\Windows\fonts -Z 25 del.gmt* gmt.* %file% legend.txt

21 GMT :

22 CAP 3.1 Cap Oh-My-Cap perl cpanm AXI.SC BHE SAC RESP.SC.AXI.00.BHE SAC oh-my-cap/data qzgy.fk qzgy event.conf event.info RunMe.sh RunMe.sh 1 # DATA RunMe.sh 2 $ cd oh-my-cap/data 3 $ ls # # 6 Catlog_cap.txt # 7 chosedata.sh # 8 draw_st.sh # 9 earth_relief_02m.grd # 10 event.conf # cap 11 event.info # 12 gmt4depth.pl # 13 GMT_globe.cpt # cpt 14 mod.txt # 15 my-process.pl # 16 rename.pl # 17 resp # RESP 18 RunMe.sh # 19 sn.sh # 20 $ bash RunMe.sh 21

23 CAP $ cat Glib/qzgy.fk 2 # 5.1( 4.9 ) 3 # P Vs=Vp/sqrt(3) 4 # FK configuration 5 DIST: 15/500/5 6 DEPTH: 2/36/1 7 NT: DT: FLAT: NO # NO or YES. if undefined, scripts will use YES $ cd Glib 19 $ perl run_parallel.pl qzgy.fk event.info SAC 1 SAC> lh kzdate kztime b delta npts e o a f evla evlo evel evdp mag dist az event.info Cap Taup cut cut 1 # 8h 2 $ echo T11:26: > event.info event.conf example event.conf 1 $ cat event.conf 2 # CAP paramete# CAP parameters 3 MODEL: qzgy 4 DEPTH: 3/36/3 5 MAG: 4.0 6

24 CAP 23 7 # CAP options 8 -J: 0/0.1/0/0.1 process.pl rename.pl 1 #!/usr/bin/env perl 2 use strict; 3 use warnings; 4 $ENV{SAC_DISPLAY_COPYRIGHT} = 0; 5 == 1 or die "Usage: perl $0 dir\n"; 7 my ($dir) 8 9 chdir $dir; # 12 foreach my $file (glob "*20*.BH?") { 13 my ($sta, $net, $time, $loc, $chn) = split /\./, $file; 14 rename $file, "$net.$sta.$loc.$chn.sac"; 15 } chdir ".."; process.pl my-process.pl ( rdseed.pl rename.pl) 1 $ perl my-process.pl chosedata 1 # 2 $ perl weight.pl # 4 $ perl inversion.pl # 6 $ perl get_depth.pl Warning: flag=21 => the minimum 0.0/57.8/ 6.4 is at boundary P no minimum found SAC RESP PZ

25 CAP 24 event.info ±5 1 #!/bin/bash 2 awk '{print $3,$2}' event.info awk '{ 3 lon1=int($1-5);lon2=int($1+5); 4 lat1=int($2-5);lat2=int($2+5); 5 print "grdcut *.grd -R"lon1"/"lon2"/"lat1"/"lat2" -G1.grd"; 6 print "grdgradient 1.grd -Gint.grd -Nt1 -A30/60"; 7 print "grdimage -R -JM5i -CGMT_globe.cpt -Iint.grd 1.grd -K -P -B1 > st.ps"; 8 print "psxy -Sa0.3i -R -J -Gred -K -O >> st.ps <<EOF"; 9 print $1,$2; 10 print "EOF";}' sh 11 pscoast -R -J -W1p -K -O >> st.ps 12 awk '{print $2,$3}' st_info psxy -St0.2i -R -J -G55/0/250 -K -O >> st.ps 13 awk '{print $2,$3+0.1," BC",$6}' st_info pstext -R -J -O >> st.ps 14 ps2raster st.ps -A -Tf 15 rm.gmt* *.grd st.ps Z Pnl (t0) 100s (e) 100s 1 #!/bin/bash 2 echo \# stlo stla az dist s/n kstnm > sn_info 3 ls -1 *.z awk '{print "r",$1; 4 print "setbb var1 (&1,t0&)"; 5 print "setbb var2 (&1,t0&+100)"; 6 print "setbb var3 (&1,e&-100)"; 7 print "setbb var4 (&1,e&)"; 8 print "mtw %var1% %var2%"; 9 print "rms to user5"; 10 print "mtw %var3% %var4%" 11 print "rms to user6"; 12 print "setbb sn (&1,user5&/&1,user6&)"; 13 print "ch user7 %sn%"; 14 print "lh stlo stla az dist user7 kstnm"; 15 }END{print "q";}' sac \ 16 awk '/ stlo /{printf "%-8.4f ",$3} 17 / stla /{printf "%-8.4f ",$3} 18 / az /{printf "%-8.4f ",$3} 19 / dist /{printf "%-8.4f ",$3} 20 / user7 /{printf "%-8.4f ",$3} 21 / kstnm/{printf "%s\n",$3}' sort -k 5 -rn >> sn_info

26 CAP Cap Hz f1 f2 f3 f trans cap.pl -C km 3. ISO CLVD event.conf -J: 0/0.1/0/0.1 h RunMe.sh 1 #!/bin/bash 2 begin_time=$(date +%s) 3 # ( ) 4 dir= if [! -d $dir ]; then 7 echo ===== $dir 8 exit 1 9 fi echo ===== $dir 12 rm -rf../$dir 13 # 14 echo ===== event.info event.conf $dir 15 cp event.conf event.info $dir 16 echo ===== $dir 17 cp resp/* $dir 18 cd $dir 19 rm -rf chosedata 20 mkdir chosedata 21 echo ===== saclst km 22 saclst stlo stla az dist kstnm f *20*.BHZ awk '$5>80 && $5<300 {print $0}' sort -k 4 -rn > st_info 23 echo ===== chosedata 24 awk -F. '{print "cp event.info event.conf st_info ""*"$1"* chosedata"}' st_info sh 25 cd.. 26 echo ===== chosedata example $dir 27 cp -r $dir/chosedata../$dir 28 echo =====

27 CAP cp my-process.pl rename.pl../ 30 cd.. 31 perl my-process.pl $dir 32 echo ===== 33 perl weight.pl $dir 34 echo ===== 35 perl inversion.pl $dir 36 echo ===== 37 perl get_depth.pl $dir 38 rm my-process.pl rename.pl 39 echo ===== 40 cp DATA/draw_st.sh./$dir 41 data=data/earth_relief_02m.grd 42 cpt=data/gmt_globe.cpt 43 cp $data $cpt./$dir 44 cp DATA/sn.sh./$dir 45 cd $dir 46 bash draw_st.sh 47 echo ===== sn_info 48 bash sn.sh 49 end_time=$(date +%s) 50 declare -i cost_time=$end_time-$begin_time 51 echo ===== Done! ===== 52 echo ===== $cost_time =====

28 SAC 4.1 SAC rms SAC rms 1 SAC> r SC_HMS.z 2 SAC> lh b t0 3 4 FILE: SC_HMS.z b = e-02 7 t0 = e SAC> setbb var1 (&1,t0&) 10 SAC> setbb var2 (&1,t0&+100) 11 SAC> mtw %var1% %var2% 12 SAC> rms to user1 13 SAC> mtw b t0 14 SAC> rms to user2 15 SAC> mtw SAC> rms to user3 17 SAC> mtw b SAC> rms to user4 19 SAC> mtw b %var1% 20 SAC> rms to user5 21 SAC> mtw &1,b& %var1% 22 SAC> rms to user6 23 SAC> lh user1 user2 user3 user4 user5 user FILE: SC_HMS.z user1 = e user2 = e user3 = e user4 = e user5 = e user6 = e-06 27

29 SAC 28 bug mtw SAC

30 GAMIT 5.1 GAMIT Linux gamit10.5 Ubuntu14.04LTS source testdata root ~/ gamit10.5 Ubuntu14.04LTS 1 $ cd ~ 2 $ sudo apt-get update 3 $ sudo apt-get install gcc gfortran csh tcsh libx11-dev 4 $ sudo mv source /opt/gamit $ cd /opt/gamit $ csh install_software csh chomod 1 $ chmod +x install_software 2 $ csh install_software x11 no Makefile.config 1 $ gedit /opt/gamit10.5/libraries/makefile.config 2 3 # Makefile.config X11LIBPATH /usr/lib 5 X11INCPATH /usr/include 6 7 # Makefile.config 107 MAXATM 25 8 GAMIT size dependent variables 9 MAXSIT MAXSAT MAXATM MAXEPC # Makefile.config 284 linux 29

31 GAMIT 30 arc_lib.a 1 $ cd /opt/gamit10.5/gamit/arc 2 $ gfortran -c generalrel.f90 3 $ make 4 $ ls arc make_gamit $ gedit /opt/gamit10.5/com/make_gamit 2 3 # # set MODELES=c$MODELES arc cversn.f 1 $ gedit /opt/gamit10.5/gamit/clean/cversn.f 2 3 #. 4 return 5 end Makefile.generic 32 -m m64 1 $ gedit /opt/gamit10.5/gamit/solve/makefile.generic 2 3 # m64 -m32 4 FFLAGS = -03 -Wuninitialized -Wunused -fno-f2c -ffast-math -m32 #gftn 1 $ gedit ~/.bashrc 2 3 #.bashrc 4 export PATH="$PATH:/opt/gamit10.5/gamit/bin:/opt/gamit10.5/com:/opt/gamit10.5/ kf/bin" 5 export HELP_DIR=/opt/gamit10.5/help/ 1 $ source ~/.bashrc

32 GAMIT 31 1 $ cd /opt/gamit $ csh install_software doy sh_get_rinex 1 $ doy 2 $ sh_get_rinex

33 6.1 distaz distaz distaz C Matlab distaz.exe distaz C 1 #include <stdio.h> 2 #include <stdlib.h> 3 #include <math.h> 4 5 /* 6 c 7 c compile on Linux use : 8 c $ cc distaz.c -o distaz -lm 9 c compile on Windows use : 10 c > gcc distaz.c -o distaz.exe -lm 11 c 12 */ int main(int argc, char *argv[]) 15 { 16 double stalat, stalon, evtlat, evtlon; 17 double delta, az, baz; 18 double scolat, slon, ecolat, elon; 19 double a,b,c,d,e,aa,bb,cc,dd,ee,g,gg,h,hh,k,kk; 20 double rhs1,rhs2,sph,rad,del,daz,dbaz,pi,piby2; if (argc!= 5) { 23 printf("usage: distaz sta_lat sta_lon evt_lat evt_lon\n"); 24 printf(" Returns: Delta Baz Az\n"); 25 exit(1); 26 } else { 27 stalat = atof(argv[1]); 28 stalon = atof(argv[2]); 32

34 33 29 evtlat = atof(argv[3]); 30 evtlon = atof(argv[4]); 31 } 32 if ((stalat == evtlat)&&(stalon == evtlon)) { 33 printf("%6.3f %6.3f %6.3f\n", 0.0, 0.0, 0.0); 34 exit(0); 35 } pi= ; 38 piby2=pi/2.0; 39 rad=2.*pi/360.0; sph=1.0/ ; scolat=piby2 - atan((1.-sph)*(1.-sph)*tan(stalat*rad)); 44 ecolat=piby2 - atan((1.-sph)*(1.-sph)*tan(evtlat*rad)); 45 slon=stalon*rad; 46 elon=evtlon*rad; a=sin(scolat)*cos(slon); 49 b=sin(scolat)*sin(slon); 50 c=cos(scolat); 51 d=sin(slon); 52 e=-cos(slon); 53 g=-c*e; 54 h=c*d; 55 k=-sin(scolat); aa=sin(ecolat)*cos(elon); 58 bb=sin(ecolat)*sin(elon); 59 cc=cos(ecolat); 60 dd=sin(elon); 61 ee=-cos(elon); 62 gg=-cc*ee; 63 hh=cc*dd; 64 kk=-sin(ecolat); del=acos(a*aa + b*bb + c*cc); 67 delta=del/rad; rhs1=(aa-d)*(aa-d)+(bb-e)*(bb-e)+cc*cc - 2.; 70 rhs2=(aa-g)*(aa-g)+(bb-h)*(bb-h)+(cc-k)*(cc-k) - 2.; 71 dbaz=atan2(rhs1,rhs2); 72 if (dbaz<0.0) { 73 dbaz=dbaz+2*pi; 74 } 75 baz=dbaz/rad; 76

35 34 77 rhs1=(a-dd)*(a-dd)+(b-ee)*(b-ee)+c*c - 2.; 78 rhs2=(a-gg)*(a-gg)+(b-hh)*(b-hh)+(c-kk)*(c-kk) - 2.; 79 daz=atan2(rhs1,rhs2); 80 if(daz<0.0) { 81 daz=daz+2*pi; 82 } 83 az=daz/rad; 84 /* 85 c 86 c Make sure 0.0 is always 0.0, not c 88 */ 89 if(abs(baz-360.) <.00001) baz=0.0; 90 if(abs(az-360.) <.00001) az=0.0; printf("%6.3f %6.3f %6.3f\n", delta, baz, az); 93 exit(0); 94 } distaz Matlab Visualization/VTK_ParaView/matlab/distaz.m 1 function [dk,dd,daze,dazs]=distaz (sta,sto,epa,epo) 2 rad=pi/180.0d0; 3 4 sa = atan( *tan(sta*rad)); 5 ea = atan( *tan(epa*rad)); 6 ssa = sin(sa); 7 csa = cos(sa); 8 so = sto*rad; 9 eo = epo*rad; 10 sea = sin(ea); 11 cea = cos(ea); 12 ces = cos(eo-so); 13 ses = sin(eo-so); if (sa==ea) 16 if (sto==epo) 17 dk =0.00; 18 dd =0.00; 19 daze=0.0; 20 dazs=0.0; 21 return 22 end 23 end 24

36 35 25 if sta== if epa== dk =0.00; 28 dd =0.00; 29 daze=0.00; 30 dazs=0.00; 31 return 32 end 33 end if sta== if epa== dk =0.00; 40 dd =0.00; 41 daze=0.00; 42 dazs=0.00; 43 return 44 end 45 end 46 dd = ssa*sea+csa*cea*ces; 47 if dd~=0., dd=atan(sqrt(1.0-dd*dd)/dd); end 48 if dd==0., dd=pi/2.0; end 49 if dd<0.0, dd=dd+pi; end 50 dd = dd/rad; 51 dk = dd*111.19; dazs = atan2(-ses,(ssa/csa*cea-sea*ces)); 54 daze = atan2(ses,(sea*csa/cea-ssa*ces)); 55 dazs = dazs/rad; 56 daze = daze/rad; 57 if (dazs<0.00), dazs=dazs+360.0; end 58 if (daze<0.00), daze=daze+360.0; end 59 return

37 36

38 7.1 Fetch Script plenv perl Fetch Scripts plenv plenv perl perl 1 # plenv 2 $ git clone git://github.com/tokuhirom/plenv.git ~/.plenv 3 $ echo 'export PATH="$HOME/.plenv/bin:$PATH"' >> ~/.bashrc 4 $ echo 'eval "$(plenv init -)"' >> ~/.bashrc 5 $ exec $SHELL -l 6 # plenv build 7 $ git clone git://github.com/tokuhirom/perl-build.git ~/.plenv/plugins/perl-build/ 8 # perl 9 $ plenv install --list 10 # 11 $ plenv install # 13 $ plenv rehash 14 # perl 15 $ plenv versions 16 $ plenv global cpanm perl cpan cpanm cpan perl perl 1 # cpanm 2 $ plenv install-cpanm 3 # 4 $ cpanm install XML::SAX 5 $ cpanm install Bundle::LWP 6 $ cpanm install XML::SAX::ExpatXS 37

39 perl No patch utility found patch 1 $ sudo yum install patch Centos7 Bundle::LWP libwww-perl 1 $ yum search libwww-perl 2 axelget, fastestmirror, langpacks 3 No metadata available for base 4 No metadata available for epel 5 No metadata available for extras 6 No metadata available for google-chrome-mirrors 7 No metadata available for nux-dextop 8 No metadata available for updates 9 Loading mirror speeds from cached hostfile 10 * base: mirrors.aliyun.com 11 * epel: mirror01.idc.hinet.net 12 * extras: mirrors.aliyun.com 13 * nux-dextop: li.nux.ro 14 * updates: mirrors.sohu.com 15 =========================== N/S matched: libwww-perl =========================== 16 perl-libwww-perl.noarch : A Perl interface to the World-Wide Web 17 $ sudo yum install perl-libwww-perl.noarch XML::SAX::ExpatXS expat.h 1 $ sudo apt-get install libexpat1-dev # Ubuntu 2 $ sudo yum install expat-devel # Centos7 XML::SAX::ExpatXS Web Service Fetch Script perl python perl FetchData FetchEvent FetchMetadata python FetchSyn 1 $ mv FetchData ~/bin/fetchdata Hinet HinetPy Github NIED Hi-net Python HinetPy Python 3.4 Python Anaconda

40 Hinet : HinetPy 1 $ pip install HinetPy # 2 $ pip install HinetPy --upgrade # 3 $ pip install git+ # 4 $ pip uninstall HinetPy # HinetPy Win32tools 1 $ tar -xvf win32tools.tar.gz # 2 $ cd win32tools/ 3 $ make 4 $ cp catwin32.src/catwin32 win2sac.src/win2sac_32 ~/bin/ # Hi-net T00:00 GMT from HinetPy import Client, win32 3 # Hi-net username, password 4 client = Client("username", "password") 5 # 6 # '0101' Hi-net 7 data, ctable = client.get_waveform('0101', ' ', 20) 8 # 9 # cnt channle table 10 # 0101_ _20.cnt 0101_ ch 11 # cnt SAC 12 win32.extract_sac(data, ctable) 13 # channel table SAC PZ 14 win32.extract_pz(ctable) 7.3 CENC

41 :8080 # IP fdsnws/event # FDSN Web 4 1 # 5 query? # 6 minmag # 7 maxmag # 8 minlatitude # 9 maxlatitude # 10 minlongitude # 11 maxlongitude # 12 mindepth # 13 starttime # 14 endtime # 15 format # QuakeML text fdsnws/station # FDSN 18 network # 19 station # 20 level # network station channel response 21 format # StationXML text :00:00&endtime= T00:00:00&format=text :00:00&minlatitude=35&maxlatitude=44&minlongitude=114&maxlongitude=120&format= text :00:00&minlatitude=20&maxlatitude=50&minlongitude=74&maxlongitude=136&format= text

42 41 GS GTA metadata 1 &format=text BH 1 minlatitude=35&minlongitude=114&maxlongitude=120&level=channel&format=text GS HE level 1 text GS HE level 1 text HE CLI 1 response minlongitude=114&maxlongitude=120&level=station&format=text minlongitude=114&maxlongitude=120&level=channel&format=text 1 starttime= t15:20:00&endtime= t15:40:00 HTTP 414

43 JWEED JWEED : washington.edu/projects/jweed JWEED Select Servers Service Provider URL Service Provider : CENC-FDSN service.. URL : :8080/fdsnws,1 IRIS FetchData FetchData IRIS Perl SERVICEBASE TIMESERIESWS METADATAWS bash shell 1 $ echo `export SERVICEBASE= >> ~/.bashrc 2 $ echo `export TIMESERIESWS= >> ~/. bashrc 3 $ echo `export METADATAWS= >> ~/.bashrc 4 $ source ~/.bashrc 5 # $ FetchData -N HE -S TAH,TLK -C?HZ -s T10:00:00 -e T10:30:00 -o mseed #!/usr/bin/python 2 # ObsPy FDSN 3 from obspy.core import UTCDateTime 4 from obspy.fdsn import Client 5 # obspy.clients.fdsn obspy.fdsn 6 client=client(" 7 # 8 starttime=utcdatetime(" t10:10:00.000") 9 # 00 Z st=client.get_waveforms("he", "TAH,TLK", "00", "?HZ", starttime, starttime + 600) 11 # 12 st.plot() 13 st.write(' mseed ',format="mseed") 14 inventory=client.get_stations(network="he", station="*") 15 inventory.plot(projection= local )

44 SOD SOD Standing Order for Data Java SOD Java Linux SOD 1 $ sudo yum install java # Centos7 java java 2 $ wget # 3 $ tar -xvf sod tgz # 4 $ sudo mv sod /opt # 5 $ echo 'export PATH=${PATH}:/opt/sod-3.2.8/bin'>> ~/.bashrc # 6 $ source ~/.bashrc # 7 $ sod -h # SOD Windows SOD SOD sod sod-3.2.8/bin PATH CMD 1 > sod.bat -h SOD 1 # linux Windows sod.bat sod 2 $ sod -r --demo > demo.xml # 3 $ sod -f demo.xml # 4 # recipe 5 Congratulations, valid recipe. 6 Fiji Islands Region (-20.6, ) m 2003/01/04 05:15:03 GMT 6.5 MW 7 Luzon, Philippines (15.7, 119.7) m 2003/01/06 23:43:50 GMT 6.0 MW 8 Chile-Argentina Border Region (-33.8, -70.1) m 2003/01/07 00:54:51 GMT 6.0 MB 9 Fiji Islands Region (-19.7, ) m 2003/01/09 02:50:45 GMT 6.0 MW 10 New Ireland Region, P.N.G. (-5.3, 153.7) m 2003/01/10 13:11:56 GMT 6.7 MW 11 Channel: II.AAK.00.BHE 12 Channel: II.AAK.00.BHN 13 Channel: II.AAK.00.BHZ 14 Channel: II.ABKT.00.BHE 15 Channel: II.ABKT.00.BHN 16 Channel: II.ABKT.00.BHZ 17 Channel: II.ALE.00.BHE 18 Channel: II.ALE.00.BHN 19 Channel: II.ALE.00.BHZ 20 Channel: II.ARU.00.BHE 21 Channel: II.ARU.00.BHN

45 44 22 Channel: II.ARU.00.BHZ Anaconda ObsPy ObsPy Python Anaconda ObsPy pyenv Python2 Python3 Python pyenv Python 1 # pyenv 2 $ curl -L bash 3 $ gedit ~/.bashrc 4 # pyenv 5 export PYENV_ROOT="$HOME/.pyenv" 6 export PATH="$PYENV_ROOT/bin:$PATH" 7 eval "$(pyenv init -)" 8 9 $ source ~/.bashrc 10 # 11 $ pyenv rehash 12 # 13 $ pyenv uninstall 14 # pyenv 15 $ pyenv update Anaconda Anaconda Python Python2 Python3 Python Anaconda 1 $ pyenv install --list # python # python 3

46 45 5 anaconda # python 2 6 anaconda # python # python 9 $ sudo yum install readline readline-devel readline-static 10 $ sudo yum install openssl openssl-devel openssl-static 11 $ sudo yum install sqlite-devel 12 $ sudo yum install bzip2-devel bzip2-libs 13 # anaconda $ pyenv install anaconda v 15 # 16 $ pyenv rehash 17 # anaconda3 18 $ pyenv global anaconda # 20 $ pyenv versions system 23 * anaconda (set by /home/zzy/.pyenv/version) # python 26 $ python Python Anaconda (64-bit) (default, May , 13:09:58) 29 [GCC (Red Hat )] on linux 30 Type "help", "copyright", "credits" or "license" for more information. 31 >>> Obspy conda-forge Anaconda 1 $ conda config --add channels conda-forge 2 $ conda config --add channels obspy # 32 obspy 1 $ conda create -n obspy python=3.6 # Anaconda Python Proceed ([y]/n)? y 6 Extracting packages... 7 [ COMPLETE ] ###################################### 100% 8 Linking packages... 9 [ COMPLETE ] ###################################### 100% 10 #

47 46 11 # To activate this environment, use: 12 # $ source activate obspy 13 # 14 # To deactivate this environment, use: 15 # $ source deactivate 16 # obspy 1 $ source activate obspy 2 (obspy)$ which python # Python 3 ~/.pyenv/versions/anaconda /envs/obspy/bin/python obspy 1 (obspy)$ conda install obspy 2 Fetching package metadata:... 3 Solving package specifications: (obspy)$ obspy-runtests # basemap : (obspy) $ conda install basemap ObsPy ObsPy 1 # obs_test.py, 2 # 3 # $ source activate obspy # obspy 4 # $ python obs_test.py # 5 # $ source deactivate # obspy 6 # 7 #, 8 import matplotlib.pyplot as plt 9 plt.style.use("ggplot") #, obspy fdsn IRIS

48 47 12 import obspy 13 from obspy.fdsn import Client 14 # obspy.clients.fdsn obspy.fdsn 15 c = Client("IRIS") #, ( ) 18 cat = c.get_events(starttime=obspy.utcdatetime(2010, 5, 24, 16), 19 endtime=obspy.utcdatetime(2010, 5, 24, 17), 20 minmagnitude=6.4) 21 print("event depth in km:", cat[0].origins[0].depth / ) 22 cat.plot(projection="ortho"); 7.1: 1 #, 2 st = c.get_waveforms(network="iu", station="anmo", location="00", channel="bhz", 3 starttime=obspy.utcdatetime(2010, 5, 24, 16, 18, 28), 4 endtime=obspy.utcdatetime(2010, 5, 24, 17, 18, 28), 5 attach_response=true) 6 st.plot() 1 #,, 2 st.remove_response(output="disp", pre_filt=(0.005, 0.01, 4, 8)) 3 st.plot()

49 48 7.2: 7.3: 1 #, 2 st_synth = obspy.read(" 3 "network=iu&station=anmo&components=z&" 4 "eventid=gcmt:c a&dt=0.05&" 5 "units=displacement&" 6 "model=iasp91_2s&label=tutorial&format=miniseed") 7 st_synth.plot() 7.4: 1 #, 2 from obspy.taup import TauPyModel 3

50 49 4 m = TauPyModel("IASP91") 5 arr = m.get_travel_times( 6 source_depth_in_km=582.1, distance_in_degree=54.16, 7 phase_list=["p"]) 8 print(arr) 9 10 #, P 11 # Combine into a single Stream object. 12 st_all = st.copy() + st_synth.copy() # Filter to a period band from 1 to 10 seconds to highlight 15 # the P phase and force the same band limit on both. 16 st_all.detrend("demean") 17 st_all.taper(0.05) 18 st_all.filter("bandpass", freqmin=0.1, freqmax=1.0, corners=3) # Cut window around the P arrival 21 starttime = st_synth[0].stats.starttime endtime = st_synth[0].stats.starttime st_all.trim(starttime, endtime) st_all.plot() 7.5:

51 50

52 bat 8.1 bat aaa.txt line 1 line 2 line 3 line1 line2 line3 bat off 2 for /f "eol=l tokens=*" %%i in (aaa.txt) do echo %%i >> tmp.txt 3 copy tmp.txt bbb.txt 4 del tmp.txt bbb.txt bat off 2 title 3 color 2f 4 echo. 51

53 BAT 52 5 date /t 6 echo ====================================================================== 7 echo ===== ===== 8 echo ====================================================================== 9 echo. 10 echo. 11 :begin 12 set /p name= ( txt dat ) 13 echo. 14 for /f "delims=" %%a in ('dir/b/a-d *.%name%')do (echo. &echo %%a & set/a n+=1 & echo. 15 (for /f "usebackq delims=" %%i in ("%%a") do for %%k in (%%i) do echo.%%k) > temp 16 move temp "%%a" ) 17 set /a c=0+n 18 echo. 19 echo ===== %cd% %c%!====== 20 echo. 21 echo. 22 goto begin yl txt yl txt yl txt yl TX bat off 2 setlocal enabledelayedexpansion 3 for %%i in (yl*) do ( 4 set "str=%%i" 5 ren "%%i" "!str:yl=tx!") off 2 rem 3 set YE=%date:~0,4% 4 set MO=%date:~5,2%

54 BAT 53 5 set DA=%date:~8,2% 6 set DG=1 7 set/a vy1=%ye% %% set/a vy2=%ye% %% 4 9 set/a vy3=%ye% %% if %vy1%==0 (set var=true) else (if %vy2%==0 (if %vy3%==0 (set var=false) else ( set var=true)) else (set var=false)) 11 set LY=%YE% 12 set LM=%MO% 13 if %MO:~0,1%==0 (set MO=%MO:~1,1%) 14 if %DA:~0,1%==0 (set DA=%DA:~1,1%) 15 if %DA% GTR %DG% (set/a LD=%DA%-%DG%) else ( 16 if %MO%==1 (set/a LY=%YE%-1) & (set/a LM=12+%MO%-1) & (set/a LD=31+%DA%-%DG%) else ( 17 set/a LM=%MO%-1 18 if %MO%==3 (if %var%==false (set/a LD=28+%DA%-%DG%) else (set/a LD=29+%DA%-%DG%)) 19 for %%a in ( ) do (if "%MO%"=="%%a" (set/a LD=31+%DA%-%DG%)) 20 for %%b in ( ) do (if "%MO%"=="%%b" (set/a LD=30+%DA%-%DG%)))) 21 if %LM% LSS 10 set LM=0%LM:~-1% 22 if %LD% LSS 10 set LD=0%LD:~-1% 23 set Han= 24 set Han=%LY%%LM%%LD% 25 echo %Han% off 2 set mydate= :a 4 set YE=%mydate:~0,4% 5 set MO=%mydate:~4,2% 6 set DA=%mydate:~6,2% 7 echo %mydate% 8 set/a vy1=%ye% %% set/a vy2=%ye% %% 4 10 set/a vy3=%ye% %% if %vy1%==0 (set var=true) else (if %vy2%==0 (if %vy3%==0 (set var=false) else ( set var=true)) else (set var=false)) 12 set LY=%YE% 13 set LM=%MO% 14 if %MO:~0,1%==0 (set MO=%MO:~1,1%) 15 if %DA:~0,1%==0 (set DA=%DA:~1,1%) if %MO%==2 ( 18 if %DA%==28 (if %var%==false (set LD=1) &(set LM=3) else (set/a LD=1+%DA%)) 19 if %DA%==29 (set LD=1) &(set LM=3) 20 if %DA% LSS 28 (set/a LD=1+%DA%) )

55 BAT if %MO%==12 (if %DA%==31 (set LD=1) & (set LM=1) & (set /a LY=%YE%+1) else (set/a LD=1+%DA%)) for %%a in ( ) do (if "%MO%"=="%%a" (if "%DA%"==31 (set /a LM=%MO%+1) & (set LD=1) else (set /a LD=%DA%+1))) 25 for %%b in ( ) do (if "%MO%"=="%%b" (if "%DA%"==30 (set /a LM=%MO%+1) & ( set LD=1) else (set /a LD=%DA%+1))) 26 if %LM% LSS 10 set LM=0%LM:~-1% 27 if %LD% LSS 10 set LD=0%LD:~-1% 28 set mydate=%ly%%lm%%ld% 29 echo %mydate% off 2 :: 3 set /p date1= ( :11:3): 4 set /p date2= ( :6:31): 5 for /f "tokens=1,2,3,4,5,6,7 delims=-/:. " %%i in ("%date1%") do ((set Y1=%%i) && (set M1=%%j) && (set D1=%%k) && (set H1=%%l) && (set F1=%%m) && (set S1=%%n) && (set MS1=%%o)) 6 for /f "tokens=1,2,3,4,5,6,7 delims=-/:. " %%i in ("%date2%") do ((set Y2=%%i) && (set M2=%%j) && (set D2=%%k) && (set H2=%%l) && (set F2=%%m) && (set S2=%%n) && (set MS2=%%o)) 7 set /a secs=((d *(y (m2-14)/12)/4+367*(m2-2-(m2-14)/12*12) /12-3*((y (m2-14)/12)/100)/4)*86400+H2*3600+F2*60+S2)-((d *( y (m1-14)/12)/4+367*(m1-2-(m1-14)/12*12)/12-3*((y (m1-14)/12)/100) /4)*86400+H1*3600+F1*60+S1) 8 set /a D=secs/86400,H=(secs%%86400)/3600,M=(secs%%3600)/60,S=secs%%60 9 echo.&echo.%date1% %date2% :%D% %H% %M% %S% bat off 2 reg add HKCR\.bat\ShellNew /v nullfile /f >nul 3 reg add HKCR\batfile /ve /d bat /f >nul 8.2 bat Windows bat 100 txt

56 BAT 55 txt type *.txt > result.dat type off 2 for %%i in (*.txt) do type %%i >>result.dat && echo. >> result.dat OK L13056H00.txt TEC setlocal enabledelayedexpansion 2 for /f %%i in ('dir /s/b "..\*.txt"') do ( 3 set R=-R70/140/15/55 4 set J=-JM8i 5 gmtset BASEMAP_TYPE plain 6 gmtset TICK_LENGTH 0 7 gmtset HEADER_FONT 0 HEADER_FONT_SIZE 18 ANNOT_FONT_SIZE 16 8 set file=%%i 9 echo!file! 10 set PS=!file:~-13,9!.ps 11 gawk "{if (NR>4) print $1,$2,$6}"!file! > data.txt 12 surface!r! -I0.25 data.txt -Gdata.grd 13 rem makecpt -Crainbow -T-18/18/5 -Z > cpt 14 grdimage data.grd -Ccpt!R!!J! -K >!PS! 15 psscale -Ccpt -D8.7/2.9/5.7/0.2 -B5 -K -O >>!PS! 16 psbasemap -R -J -B10/5:."@~D@~VTEC Map (!file:~-13,9! )": -K -O >>!PS! 17 psxy CHINA.xyz -R -J -W0.5p,black -M -K -O >>!PS! 18 psxy.dat -R -J -Sa0.3 -Gblack -W0.5p,black -O >>!PS! 19 ps2raster!ps! -Tj -P -C-sFONTPATH=C:/Windows/fonts 20 del.gmt* data.txt inter.xyz data.grd 21 ) convert jpg

57 BAT :

58 BAT 57 1 convert -delay 20 L13056*.jpg -loop 0 L13056.gif 8.2: 8.3 Windows Unix GMT 15 Windows Unix 1. cat.exe 2. cut.exe 3. gawk.exe 4. grep.exe 5. head.exe n 6. join.exe 7. md5sum.exe MD5 8. paste.exe 9. printf.exe C 10. sed.exe 11. sort.exe 12. split.exe

59 BAT tail.exe n 14. uniq.exe 15. wc.exe 5 1. cpdf.exe PDF 2. curl.exe URL distaz.exe 4. MIF2GMT.exe 5. wget.exe Unix Windows Unix BAT %% 8.4 bat gmt Linux 1 T1=$(grdinfo world6m.grd -T1000) 2 makecpt -Crainbow $T1 > 1.cpt grdinfo world6m.grd -T1000 T1 Windows bat for 1 for /f %%i in ('grdinfo world6m.grd -T1000') do set T1=%%i 2 makecpt -Crainbow %T1% > 1.cpt T1=%%i 8.5 ftp ftp cmd 1 # ftp 2 C:\Users\zhao > ftp # ftp 3 user: # 4 password: # 5 ftp > help # ftp 6 ftp > ls # ftp 7 ftp > dir # ftp 8 ftp >!dir # 9 ftp > pwd # ftp 10 ftp > cd A # ftp A

60 BAT ftp > lcd B # B 12 ftp > ascii # ASCII 13 ftp > bin # 14 ftp > get C.txt # ftp A C.txt B 15 ftp > put D.dat # B D.dat ftp A 16 ftp > mget *.txt # ftp txt 17 ftp > mput *.txt # txt ftp 18 ftp > prompt # mget mput yes no 19 ftp > quit # ftp 20 ftp > bye # bat ftp ftp off 2 color 0a 3 rem 4 :yesterday 5 set YE=%date:~0,4% 6 set MO=%date:~5,2% 7 set DA=%date:~8,2% 8 set DG=1 9 set/a vy1=%ye% %% set/a vy2=%ye% %% 4 11 set/a vy3=%ye% %% if %vy1%==0 (set var=true) else (if %vy2%==0 (if %vy3%==0 (set var=false) else ( set var=true)) else (set var=false)) 13 set LY=%YE% 14 set LM=%MO% 15 if %MO:~0,1%==0 (set MO=%MO:~1,1%) 16 if %DA:~0,1%==0 (set DA=%DA:~1,1%) 17 if %DA% GTR %DG% (set/a LD=%DA%-%DG%) else ( 18 if %MO%==1 (set/a LY=%YE%-1) & (set/a LM=12+%MO%-1) & (set/a LD=31+%DA%-%DG%) else ( 19 set/a LM=%MO%-1 20 if %MO%==3 (if %var%==false (set/a LD=28+%DA%-%DG%) else (set/a LD=29+%DA%-%DG%)) 21 for %%a in ( ) do (if "%MO%"=="%%a" (set/a LD=31+%DA%-%DG%)) 22 for %%b in ( ) do (if "%MO%"=="%%b" (set/a LD=30+%DA%-%DG%)))) 23 if %LM% LSS 10 set LM=0%LM:~-1% 24 if %LD% LSS 10 set LD=0%LD:~-1% 25 set Han= 26 set Han=%LY%%LM%%LD% 27 echo %LY% %LM% %LD% 28 call : 29 call : 30 call :

61 BAT call : rem 34 if %LM:~0,1%==0 (set LM=%LM:~1,1%) 35 if %LD:~0,1%==0 (set LD=%LD:~1,1%) 36 set HAN=%LY%.%LM%.1-%LM%.%LD% 37 ECHO %HAN% 38 call : 39 PAUSE 40 exit : 43 ECHO. 44 echo %LY% %LM% %LD% 45 type.ftp > temp.ftp 46 call :aaa 47 ftp -n -s:"d:\ \temp.ftp" 48 goto :eof : 51 echo %LY% %LM% %LD% 52 type.ftp > temp.ftp 53 call :aaa 54 ftp -n -s:"d:\ \temp.ftp" 55 ECHO. 56 goto :EOF : 59 ECHO. 60 echo %LY% %LM% %LD% 61 type.ftp > temp.ftp 62 call :aaa 63 ftp -n -s:"d:\ \temp.ftp" 64 ECHO. 65 goto :eof : 68 echo %LY% %LM% %LD% 69 type.ftp > temp.ftp 70 call :aaa 71 ftp -n -s:"d:\ \temp.ftp" 72 ECHO. 73 goto :eof : 76 echo %LY% %LM% %LD% 77 type.ftp > temp.ftp 78 call :aaa

62 BAT ftp -n -s:"d:\ \temp.ftp" 80 ECHO. 81 goto :eof rem dst temp.ftp abc 84 :aaa off&setlocal 86 set f="temp.ftp" 87 set src=abc 88 set dst=%han% 89 for /f "usebackq delims=" %%a in (%f%) do (if not defined flag cd.>%f%&set flag=1 90 set v=%%a 91 setlocal enabledelayedexpansion 92 set v=!v:%src%=%dst%! 93 echo.!v!>>%f% 94 endlocal 95 ) 96 goto :eof ftp -n -s:".ftp file".ftp.ftp 1 #.ftp 2 open your IP # your IP IP 3 user username # username 4 password # password 5 lcd d:\ \ 6 bin 7 prompt 8 get abc.dat 9 quit 10 exit 11 # 68.ftp 12 open your IP 13 user username 14 password 15 get 13142X411ZKGD9042abc.EPMS 16 quit 17 exit

63 C 9.1 Windows C Windows gcc MinGW > > > > PATH C:\MinGW\bin PATH : lib : C:\MinGW\lib : include : C:\MinGW\include cmd mingw-get MinGw installation manager cmd gcc g++ gdb MinGw installation manager 1 > mingw-get install gcc 2 > mingw-get install g++ 3 > mingw-get install gdb C test.c 1 # include <stdio.h> 2 int main() 3 { 4 printf("%s\n","hello world"); 5 return 0; 6 } cmd gcc test.c a.exe gdb a.exe 62

64 C gcc a.out -o cc distaz.c -o distaz -lm -lm C math.h gcc -lm undefined reference to 'xxxx' -l Linux cc gcc Linux cc gcc

65 Java 10.1 Java Java Java Ubuntu 1 $ tar -zxvf jdk-8u101-linux-x64.tar.gz 2 $ sudo mv jdk-1.8.0_101 /opt/java 3 $ export "JAVA_HOME=/opt/java" >> ~/.bashrc 4 $ export "JRE_HOME=/opt/java/jre" >> ~/.bashrc 5 $ export "CLASSPATH=.:$CLASSPATH:$JAVA_HOME/lib:$JRE_HOME/lib" >> ~/.bashrc 6 $ export "PATH=$PATH:$JAVA_HOME/bin:$JRE_HOME/bin" >> ~/.bashrc 7 $ source ~/.bashrc 8 $ java -version 9 java version "1.8.0_101" Java Centos 1 $ sudo yum install java Java Windows oracle Windows Java JAVA_HOME : C:\Program Files (x86)\java\jdk1.8.0_91; CLASSPATH.;%JAVA_HOME%\lib\dt.jar;%JAVA_HOME%\lib\tools.jar; Path %JAVA_HOME%\bin;%JAVA_HOME%\jre\bin; CMD java -version Java

66 Matlab 11.1 Matlab Matlab Linux Ubuntu14.04LTS Ubuntu 1 $ sudo mkdir ~/temp_matlab 2 $ sudo mount -o loop path_to_matlab/mathworks_r2014a.iso ~/temp_matlab 3 $ cd ~/temp_matlab 4 $ sudo./install Crack license_405329_r2014a.lic ## Matlab 1 $ sudo chmod 755 /usr/local/matlab 2 $ sudo cp path_to_matlab/crack/linux/libmwservices.so /usr/local/matlab/r2014a/bin /glnxa64 3 $ sudo umount ~/temp_matlab 4 $ matlab matlab 1 $ matlab 2 matlab command not found 1 $ cd /usr/local/bin/ 2 $ sudo ln -s /usr/local/matlab/r2011b/bin/matlab matlab 3 $ matlab

67 MATLAB Matlab cpt Matlab Kelly Kearney Matlab cptcmap cpt cptcmap Kelly Kearney cptcmap-pkg 4 m 1 cptfiles cptcmap.m cpt Matlab cptfiles 22 cpt cptcmap.m cptcmap.m cptcbar.m colorbar parsepv.m cptcbar.m cptcbar.m minmax.m cptcmap cptfiles cpt cpt cpt cptfiles_personal % see more details in cptcmap.m 2 % basic usage 3 cptcmap('name'); % name: cpt file name without.cpt 4 cptcmap('showall'); % display all colormaps in cptfiles 5 % optional input variables 6 cptcmap('name', 'mapping', 'direct'); % or 'scaled' [default] 7 cptcmap('name', 'ncol', num); % number of colors 8 cptcmap('name', 'flip', true); % reverse the colormap order, default is false a = 0:0.01:1; 2 b = a; 3 c = rand(101); 4 pcolor(a,b,c); 5 shading flat; 6 cptcmap('gmt_globe'); 7 cb = colorbar; 8 set(cb, 'fontsize', 12);

68 MATLAB 67 : 11.1:

69 Perl 12.1 Perl SAC yyyy/yyyymm/yyyymmddhhmmss /*.SAC RESP resp SAC SAC 3 1. evalresp RESP transfer SAC SAC ; 2. evalresp fname RESP ; 3. RESP RESP RESP Perl SAC 1 #!/usr/bin/env perl 2 use strict; 3 use warnings; 4 $ENV{SAC_DISPLAY_COPYRIGHT} = 0; 5 my ($f1) = 0.005; 6 my ($f2) = 0.01; 7 my ($f3) = 1.0; 8 my ($f4) = 1.3; 9 my $dir; 10 foreach $dir (<*/*/*>) { 11 if (-d $dir) { 12 system "cp resp/* $dir" 13 chdir $dir; 14 open(sac, " sac") or die "Error in opening sac\n"; 15 print SAC "wild echo off\n"; 16 print SAC "r *.SAC \n"; 17 print SAC "rglitches; rmean; rtrend; taper \n"; 18 print SAC "trans from evalresp to vel freq $f1 $f2 $f3 $f4\n"; 19 print SAC "mul 1e-7\n"; 20 print SAC "w over\n"; 21 print SAC "q\n"; 68

70 PERL close(sac); 23 unlink glob "RESP*"; 24 chdir "../../.."; 25 } 26 }

71 70

72 13.1 Centos7 Hugo Even GitHub Pages zhaozhiyuan.org Travis-ci Hugo Hugo Go Hugo Release hugo hugo.exe $ mv hugo $HOME/bin # hugo 2 $ hugo -h # hugo 3 hugo is the main command, used to build your Hugo site $ hugo new site hugo_even # hugo_even 6 $ cd hugo_even # hugo 7 $ hugo new post/first.md # content/post/ 8 $ git clone themes/even # even 9 $ gedit config.toml # config.toml 10 $ hugo server # more 1 <!- -more- -> # A

73 72 GitHub zhaozhiyuan.org gh-pages setting GitHub Pages Source gh-pages branch Custom domain zhaozhiyuan.org Github Pages zhaozhiyuan.org master.travis.yml setting - Personal access tokens - Generate new token Github Travis-ci zhaozhiyuan.org Environment Variable GH_TOKEN value master gh-pages myproject.git GitHub Pages gh-pages setting custom domain A project CNAME zhaozhiyuan1989.github.io 13.2 gitalk Gitalk Gitalk Gitalk Github Issue Preact en zh-cn zh-tw es-es fr Github

74 Gitalk Github APP 1 - Application name : 2 blog gitalk # 3 - Homepage URL : 4 # 5 - Application description : 6 # 7 - Authorization callback URL : 8 # 9 - Registered application 10 Client ID:192cc5c033afc1e788a9 # gitalk 11 Client Secret:6663a17a2802a6c87c9b4eb fceb73deea # gitalk gitalk /layouts/post/single.html 1 <!-- gitalk begin --> 2 <div id="gitalk-container"></div> 3 </div> 4 <link rel="stylesheet" href=" 5 <script src=" 6 <script type="text/javascript"> 7 var gitalk = new Gitalk({ 8 id: '{{.Date }}', 9 title: '{{.Title }}', 10 clientid: '192cc5c033afc1e788a9', 11 clientsecret: '6663a17a2802a6c87c9b4eb fceb73deea', 12 repo: 'blog-comment', 13 owner: 'zhaozhiyuan1989', 14 admin: ['zhaozhiyuan1989'], 15 body: decodeuri(location.href), 16 }); 17 gitalk.render('gitalk-container'); 18 </script> 19 <!-- gitalk end --> : clientid clientsecret repo owner repo admin repo Github Gitalk

75 74 Github MarkDown Github repo issue GMT PR GMT PR GMT PR hugo-theme-even PR 13.3 GitHub Pages HTTPS GitHub Pages Custom Domain HTTPS Github zhaozhiyuan.org Settings GitHub Pages Enforce HTTPS Unavailable for your site because your domain is not properly configured to support HTTPS A GitHub Pages A GitHub Not yet available for your site because the certificate has not finished being issued Enforce HTTPS HTTPS HTTPS HTTP HTTPS

76 texlive 14.1 texlive texlive2017.iso Windows texlive texlive2017.iso install-tl-advanced.bat cmd 1 > tex -v 2 > latex -v 3 > texworks Windows texlive texworks Linux texlive 1 # 2 $ sudo apt-get install libdigest-perl-md5-perl perl-tk # Ubuntu 3 $ sudo yum install perl-digest-md5 perl-tk # Centos 4 $ sudo mount -o loop texlive2017.iso /mnt/ # 5 $ cd /mnt 6 $ sudo./install-tl # I 7 $ gedit ~/.bashrc 8 9 # TeX Live export MANPATH=${MANPATH}:/usr/local/texlive/2017/texmf-dist/doc/man 11 export INFOPATH=${INFOPATH}:/usr/local/texlive/2017/texmf-dist/doc/info 12 export PATH=${PATH}:/usr/local/texlive/2017/bin/x86_64-linux $ source ~/.bashrc 15 $ sudo umount /mnt/ 16 # texlive 17 $ sudo tlmgr update --self --repository texlive/tlnet/ 75

77 76 18 $ sudo tlmgr update --all --repository texlive/tlnet/ sudo apt-get install texlive Linux texlive 14.2 pandoc pandoc markdown pandoc markdown pdf html docx UTF-8 notepad++ UTF pandoc Windows pandoc pandoc windows.msi pandoc windows.zip PATH MacOS pandoc brew 1 $ brew install pandoc Linux pandoc pandoc linux.tar.gz 1 # pandoc $ sudo tar xvzf pandoc linux.tar.gz --strip-components 1 -C /usr 3 # pandoc rm /usr/bin/pandoc 4 $ which pandoc 5 $ pip install pandocfilters # 6 # 7 $ pandoc -h 8 # 9 $ pandoc --list-input-formats

78 77 10 # 11 $ pandoc --list-input-formats 12 # 13 $ pandoc --list-highlight-languages 14 # 15 $ pandoc --list-highlight-styles f -t -o -s -i --toc --listings listings --mathjax html --pdf-engine=xelatex xelate pdf -V mainfont=simsun theme=default --slide-level= # md2html 2 $ pandoc -f markdown -t html test.md -o test.html -s --mathjax 3 $ pandoc test.md -o test.html 4 5 # md2pdf pdf latex texlive 6 # pandoc pdf documentclass ctex ctexrep ctexbook ctexart 7 $ pandoc test.md --pdf-engine=xelatex -V mainfont=simsun -o test.pdf 8 $ pandoc test.md -o test.pdf --template=zzy.tex --listings --pdf-engine=xelatex 9 10 # md2rst 11 $ pandoc test.md -o test.rst # md2docx 14 $ pandoc test.md -o test.docx markdown pandoc markdown % markdown

79 78 pandoc html DZslides slidy s5 slidous reveal.js reveal.js default beige sky night serif simple solarized black blood league moon white 1 # DZslides 2 $ pandoc test.md -o slides.html -t dzslides -s --mathjax -i 3 # Slidy 4 $ pandoc test.md -o slides.html -t slidy -s --mathjax -i 5 # revealjs 6 $ git clone 7 $ pandoc test.md -o slides.html -t revealjs -s -V theme=beige --mathjax -i 8 # beamer pdf 9 $ pandoc test.md -o slides.pdf -s -i -t beamer --pdf-engine=xelatex -V mainfont= SimSun 10 # beamer documentclass ctexrep 11 $ pandoc -D beamer > test.tex 12 $ pandoc -o test.pdf -t beamer --template test.tex --pdf-engine=xelatex test.md 14.3 cpdf pdf cpdf Coherent PDF cpdf PDF zip PDF PDF PDF PDF PDF

80 79 logo PDF cpdf cpdf-binaries-master.zip exe PATH PDF 1 # in.pdf out.pdf 2 $ cpdf in.pdf 1-3,12-end -o out.pdf 3 # in.pdf $ cpdf in.pdf ~5-~3 -o out.pdf 5 # in.pdf 6 $ cpdf in.pdf even -o out.pdf 7 # in.pdf 8 $ cpdf in.pdf odd -o out.pdf 9 # in.pdf page001.pdf page002.pdf 10 $ cpdf -split in.pdf -o page%%%.pdf # bat page%%%%%%.pdf 11 # in.pdf 10 1 pdf 12 $ cpdf -split in.pdf -o Chunk%%%.pdf -chunk 10 # bat Chunk%%%%%%.pdf 13 # in.pdf 14 $ cpdf -split-bookmarks 0 in.pdf 15 # one.pdf two.pdf three.pdf merge.pdf 16 $ cpdf one.pdf two.pdf three.pdf -o merged.pdf 17 # 2017 pdf 18 $ cpdf 2017*.pdf -o out.pdf 19 # $ cpdf -pad-after 1,3,4 in.pdf -o out.pdf PDF 1 # 2 $ cpdf -info file.pdf 3 # 4 $ cpdf -page-info file.pdf 5 # PDF 6 $ cpdf -blacktext in.pdf -o out.pdf 7 # x y in.pdf 2 8 $ cpdf -scale-page "2 2" in.pdf -o out.pdf

81 80 9 # A4 10 $ cpdf -scale-to-fit a4portrait in.pdf -o out.pdf 11 # in.pdf $ cpdf -twoup-stack in.pdf -o out.pdf 13 # in.pdf x 26pt y 18mm 14 $ cpdf -shift "26pt 18mm" in.pdf -o out.pdf 15 # in.pdf $ cpdf -rotate-contents 90 in.pdf -o out.pdf 17 # in.pdf 600pt x 400pt 18 $ cpdf -crop " pt 400pt" in.pdf -o out.pdf 19 # 20 $ cpdf -crop "20mm 20mm 300mm 300mm" in.pdf -o out.pdf PDF 1 # in.pdf zzy user 2 $ cpdf -encrypt 128bit zzy user in.pdf -o encrypt.pdf 3 # encrypt.pdf 4 $ cpdf -decrypt encrypt.pdf owner=zzy -o decrypt.pdf PDF 1 # 2 $ cpdf -compress in.pdf -o out.pdf 3 # 4 $ cpdf -decompress in.pdf -o out.pdf PDF 1 # in.pdf 2 $ cpdf -list-bookmarks in.pdf 3 0 "Part 1" 1 open 4 1 "Part 1A" "Part 1B" "Part 2" "Part 2a" 5 8 # bookmarks.txt in.pdf 9 $ cpdf -add-bookmarks bookmarks.txt in.pdf -o out.pdf 10 # in.pdf 11 $ cpdf -list-annotations in.pdf Annotation text content Annotation text content 2

82 # from.pdf in.pdf 20 $ cpdf -copy-annotations from.pdf in.pdf -o out.pdf PDF 1 # Split in.pdf 10s 2 $ cpdf -presentation in.pdf 2-end -trans Split -duration 10 -o out.pdf 3 # pdf 4 $ cpdf -set-page-mode FullScreen in.pdf -o out.pdf PDF logo 1 # 2 $ cpdf -add-text "Page %Page of %EndPage" -top 100pt -font "Times-Roman" -fontsize 20 in.pdf -o out.pdf 3 # logo 4 $ cpdf -stamp-on logo.pdf in.pdf -o out.pdf 5 # AND in.pdf in2.pdf 6 $ cpdf -merge in.pdf in2.pdf AND -add-text "Copyright 2017" -o out.pdf 7 # pdf 5 8 $ cpdf -attach-file sheet.xls -to-page 5 in.pdf -o out.pdf 9 # 10 $ cpdf -remove-files in.pdf -o out.pdf % Windows bat %% LICEcap LICEcap Windows OSX.gif LICEcap 230kb LICEcap v1.26 for Windows LICEcap v1.25 for OSX

83 82.gif.lcf shift+space CPU ttygif ttygif Linux ttyrec gif Linux Centos7 1 # 2 $ sudo yum install ImageMagick gcc 3 # ttyrec 4 $ git clone 5 $ cd ttyrec 6 $ make 7 # 3 PATH 8 $ cp ttyrec ttyplay ttytime ~/bin 9 $ cd.. 10 # ttygif 11 $ git clone 12 $ cd ttygif 13 $ make 14 $ cp ttygif ~/bin ttygif 1 C99 for for (int i=0; i<argc; i++) 4 ^ c c99 ttygif.c make 1 # ttygif.c int i; 4 for (i=0; i<argc; i++) ttyrec ttyplay ttytime ttygif ttyrec gif 1 $ ttygif myrecording 2 failed to execute: xwd -id out /tmp/ttygif.ln1nqw/0.xwd

84 83 xwd 1 $ sudo yum install xwd 2 $ ttygif myrecording ScreenToGif gif ScreenToGif gif gif ScreenToGif gif 14.5 convert ImageMagick ImageMagick 200 PNG JPEG GIF TIFF DPX PDF SVG ImageMagick ImageMagick ImageMagick Centos7 ImageMagick7 yum ImageMagick6 1 $ sudo yum install ImageMagick Windows7 ImageMagick7

85 bash bat convert bash 1 # la jpg my.gif 30 ms 0 2 $ convert -delay 30 -loop 0 la*.jpg my.gif 3 # gif 4 $ convert my.gif -coalesce temporary.gif 5 $ convert -size <original size> temporary.gif -resize 24x24 smaller.gif GMT5 1 #!/bin/bash 2 R=70/150/10/60 3 J=M25c 4 B=10 5 PS=0china.ps 6 gmt psxy -J$J -R$R -T -K -Xc -Yc > $PS 7 for code in ; 8 do 9 declare -i red=$random*256/ declare -i green=$random*256/ declare -i blue=$random*256/ gmt pscoast -R$R -J$J -ECN.$code+p1p+g$red/$green/$blue -K -O >> $PS 13 cp $PS $code.ps 14 gmt psxy -R$R -J$J -T -O >> $code.ps 15 done 16 gmt psxy -R$R -J$J -T -O >> $PS 17 gmt psconvert -P -E100 *.ps 18 convert -delay 30 -loop 0 *.jpg china.gif 19 rm gmt.* *.ps *.jpg Atom guake Source Code Pro Atom github markdown guake-terminal F12

86 : source-code-pro 0 o O Atom guake source-code-pro Atom 1 $ sudo yum localinstall atom.x86_64.rpm # centos atom 2 $ sudo dpkg -i atom.x86_64.deb # ubuntu atom guake 1 $ sudo yum install guake # centos guake 2 $ sudo add-apt-repository ppa:webupd8team/unstable # ubuntu guake 3 $ sudo apt-get update 4 $ sudo apt-get install guake