source: SH_SHM/trunk/command/READHORIZ_LNX.SHC @ 51

Revision 51, 1.1 KB checked in by marcus, 14 years ago (diff)

r36 | svn | 2008-08-28 08:47:03 +0200 (Do, 28 Aug 2008) | 1 line

updated station info; phasmap reads from db; new command procedures

Line 
1! file READHORIZ.SHC
2!      =============
3!
4! version 1, 13-Mar-2003
5!
6! Reads horizontal component (R or T) of a station set
7! K. Stammler, 13-Mar-2003
8
9default 1 db:      path
10default 2 ;;       read time
11default 3 ;;       read length
12default 4 ;;       station list
13default 5 ;;       component
14default 6 ;;       channel
15default 7 ;;       azimuth
16
17sdef ncnt 1           ! n counter
18sdef ecnt 2           ! e counter
19sdef maxtrc           ! number of traces
20
21nr
22
23!reads/noswap/invhdr #1 #2 #3 #4 ne #6
24! changed for use in beijing
25reads/noswap #1 #2 #3 #4 ne #6
26demean all
27
28calc i &maxtrc = $dsptrcs
29
30loop_start:
31
32   if  "ecnt gti "maxtrc  goto/forward loop_exit:
33   if  ^comp("ncnt) eqs n  goto/forward n_ok:
34      del "ncnt
35      goto loop_start:
36   n_ok:
37   if  ^comp("ecnt) eqs e  goto/forward e_ok:
38      del "ecnt
39      goto loop_start:
40   e_ok:
41
42   rotaz |"ncnt|,|"ecnt| #7
43
44   calc i &ncnt = "ncnt + 2
45   calc i &ecnt = "ecnt + 2
46
47goto loop_start:
48loop_exit:
49
50del |1-|"maxtrc|
51if  #5 eqs t  goto/forward del_r:
52   del _comp(t)
53   goto/forward del_ok:
54del_r:
55   del _comp(r)
56del_ok:
57
58rd
59
60return
Note: See TracBrowser for help on using the repository browser.