wangliang1989

#!/usr/bin/env perl
use strict;
use warnings;

my $x = 70;
my $y = 50;
my $d = 2;
my $j = "H22c";
my $r = "g";
my $meca = "1,2,3,-4,-5,-10";# gmt model surport NED only M11,M22,M33,M12,M13,M23 1为北向（N），2为东向（E），3向下（D）

wangliang1989

#!/usr/bin/env perl
use strict;
use warnings;

@ARGV == 2 or die "Usage: perl $0 latitude longitude\n";
my ($lat, $lon) = @ARGV;
print "make sure latitude first!\n\n";
my $filename = 'jkdwedd';
open (CRUST, "| ./getCN1point > $filename");
print CRUST "$lat $lon\n";

wangliang1989

#Network | Station | Latitude | Longitude | Elevation | SiteName | StartTime | EndTime 
#台湾地区为例
EM|TW001|23.625|119.506|28.9|PengHu, Taiwan|2006-10-11T08:49:23|2007-06-04T04:37:19
EM|TW005|24.772|120.931|89.0|Hsinchu, Taiwan|2007-05-27T02:44:21|2007-06-02T04:55:02
EM|TW010|24.728|121.072|226.7|Chutung, Taiwan|2007-04-10T08:54:59|2007-04-28T09:41:35
EM|TW020|24.691|121.148|383.8|Hohsing, Taiwan|2007-04-02T02:34:28|2007-04-28T08:34:51
EM|TW030|24.696|121.223|562.7|Chiale, Taiwan|2007-04-02T06:10:12|2007-04-28T07:24:38
EM|TW040|24.673|121.294|1330.0|Yufeng, Taiwan|2007-04-02T09:01:00|2007-04-28T05:57:31
EM|TW050|24.683|121.359|984.4|Baling, Taiwan|2007-04-01T08:06:33|2007-05-03T08:04:17
EM|TW060|24.645|121.435|1140.9|Mingchih, Taiwan|2007-04-01T05:23:02|2007-05-03T06:54:05

wangliang1989

#!/usr/bin/env perl
use strict;
use warnings;
use Time::Local;

#将如下格式的地震目录转为dmc要求的格式
#5152004|2015-05-25T10:46:03|28.5718|87.2411|42.44|us|NEIC PDE|us|us10002c08|mb|4.3|us|XIZANG
# ============= 需要自己改的 ===============
my $infile = "cat.txt";
my $outfile = "20161208.txt";