gremlin_help.f 8.55 KB
Newer Older
1
2
c this is <gremlin_help.f>
c------------------------------------------------------------------------------
3
c $Id: gremlin_help.f,v 1.3 2002-04-19 15:20:34 forbrig Exp $
4
5
6
7
8
9
10
11
12
c
c 25/03/98 by Thomas Forbriger (IfG Stuttgart)
c
c Here is some extra help on gremlin commands
c
c REVISIONS and CHANGES
c    25/03/98   V1.0   Thomas Forbriger
c    14/01/99   V1.1   removed link to old gremlin code
c    04/03/99   V1.2   file/mtt explanation
13
14
c    24/05/00   V1.3   - introduced reso option wgpd
c                      - introduced reso option wense
15
c    19/04/02   V1.4   introduced command swense
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
c
c==============================================================================
c
      subroutine gremhelp_main
c
c give help for main menu
c
      call gremhelp_tit('main')
      call gremhelp_subtit('display data')
      print *,' dgr            display real green data in different modes'
      print *,' dtt            display travel time data'
      print *,' dgi index      plot green data from data array index'
      print *,' dti index      plot traveltime data from data array index'
      call gremhelp_subtit('display/edit model')
      print *,' dre            display reference model in a table'
      print *,' dwo            display working copy of model in a table'
      print *,' med            edit actual model (model will be written to'
      print *,'                file edit.p.mod and will be read again)'
      print *,' dmo file       display basic control screen for model in "file"'
      print *,' dmi index      display named model of index'
      call gremhelp_subtit('set/display free parameters')
      print *,' spc            set free model parameters'
      print *,' imo mode       preset parameters for inversion modes'
      print *,'                (what you would do by "spc" normally)'
      print *,'     sgrad        shear wave gradient'
      print *,'     pvel         p-wave model'
      print *,' weight para    call mod_weight with "para" (modifies'
      print *,'                "spc" settings)'
      print *,' sano para      call par_sano with "para" (modifies'
      print *,'                "spc" settings)'
      print *,' dpc            display free parameter settings in a table'
      print *,'                similar to the "spc" mask'
      print *,' tpc            display free parameter settings in a table'
      call gremhelp_subtit('set/display inversion parameters')
      print *,' spa            set various parameters'
      print *,' dpa            display all parameter settings'
      call gremhelp_subtit('control inversion')
      print *,' dda            display basic crontrol screen'
      print *,' opt lim,step   optimize model with extra break condition'
      print *,'     lim          stop when X2 reaches "lim"'
      print *,'     step         stop after a maximum os "step" iterations'
      print *,' ofi nu,lim,step,mode'
      print *,'                optimize at fixed "nu"'
      print *,'     lim          stop when X2 reaches "lim"'
      print *,'     step         stop after a maximum os "step" iterations'
      print *,'     mode         1: plot improvement'
      print *,'                  2: plot improvement, model and data'
      print *,' x2c numin,numax,npts'
      print *,'                plot X2 development'
      print *,'     numin        start with this nu'
      print *,'     numax        end with this nu'
      print *,'     npts         number of points to plot'
      print *,' lx2 numin,numax,npts'
      print *,'                plot X2 development for linearized problem'
      print *,'     numin        start with this nu'
      print *,'     numax        end with this nu'
      print *,'     npts         number of points to plot'
      print *,' tpa para       fetch new model for given "nu"'
      print *,'     para         parameter to calculate "nu" from'
      call gremhelp_subtit('misc')
      print *,' verb           set verbosity modes'
      print *,' mon            switch monitor devices on/off'
      print *,' pgpar          set plot parameters'
      print *,' dpg            display pgplot parameters'
      print *,' dev device     change pgplot device'
      print *,' file           read/write various files'
      print *,' reso           linear resolution analysis'
c      print *,' old            enter old ancient gremlin code'
c      print *,'                (on your own risk)'
      print *,' term           use "term" to exit the program'
      return
      end
c
c----------------------------------------------------------------------
c 
      subroutine gremhelp_resan
c
c help on resolution analysis menu
c
      call gremhelp_tit('resolution analysis')
      print *,' '
      print *,' parder         rate partial derivatives'
      print *,' tpd ival       plot travel time partial derivatives for parameter ival'
      print *,' tpda           plot travel time partial derivatives for all parameters'
      print *,' gpd ival       plot green partial derivatives for parameter ival'
      print *,' gpda           plot green partial derivatives for all parameters'
102
103
      print *,' wgpd name      write green partial derivatives to files with'
      print *,'                base ''name'' '
104
105
106
107
108
109
      print *,' ense f,t,nu,rms   rate model parameters with index from ''f'' to'
      print *,'                ''t'' when optimizing all other free parameters'
      print *,'                stabilized with ''nu'' and accepting an rms-error'
      print *,'                of ''rms'' in a linear approximation'
      print *,' sense f,t,nu,rms  just like ''ense'' but puts rating for each'
      print *,'                parameter on a different page'
110
111
112
      print *,' wense f,t,nu,rms,name'
      print *,'                just like ''sense'' but puts rating for each'
      print *,'                parameter to a file with base ''name'' '
113
114
115
116
117
118
      print *,' swense f,t,nu,err,name'
      print *,'                just like ''wense'' but checks square'
      print *,'                error rather than rms error'
      print *,'                ''err'' is the allowed increase of the' 
      print *,'                square error relative to the misfit'
      print *,'                of the reference synthetics'
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
      call gremhelp_subtit('tests')
      print *,' orth           calculate scalar products and normalized scalar'
      print *,'                products of vectors of partial derivatives'
      call gremhelp_subtit('parameter control')
      print *,' spc            set free model parameters'
      print *,' dpc            display free parameter settings in a table'
      print *,'                similar to the "spc" mask'
      print *,' tpc            display free parameter settings in a table'
      print *,' dgr            display real green data in different modes'
      print *,' dtt            display travel time data'
      print *,' spa            set various parameters'
      print *,' dpa            display all parameter settings'
      print *,' dda            display basic crontrol screen'
      call gremhelp_subtit('misc')
      print *,' exit           leave submenu'
      print *,' quit           leave submenu'
      return
      end
c 
c----------------------------------------------------------------------
c 
      subroutine gremhelp_files
c
c help on files menu
c
      call gremhelp_tit('files')
      print *,' '
      print *,'green file     read green data from "file"'
      print *,'gwread file    read green weights from "file"'
      print *,'tt file        read travel-time data from "file"'
      print *,'mtt file       merge travel-time data from "file"'
      print *,'               this will be used for sections below asphalt'
      print *,'model file     read polynomial model from "file"'
      print *,'para file      read gremlin parameters from "file"'
      print *,'save file      write actual work model to "file"'
      print *,'resp file      write complex response file'
      print *,'wsave file     save model parameter weights to file'
      print *,'wread file     read model parameter weights from file'
      print *,'exit           leave submenu'
      print *,'quit           leave submenu'
      return
      end
c 
c----------------------------------------------------------------------
c
      subroutine gremhelp_subtit(title)
c
c print help subtitle
c
      character title*(*)
c
      print *,' '
      print *,title
      return
      end
c
c----------------------------------------------------------------------
c
      subroutine gremhelp_tit(title)
c
c print help title
c
      character title*(*)
c
      print *,' '
      print *,'GREMLIN (help): ',title
      return
      end
c
c ----- END OF gremlin_help.f -----