gremlin_help.f 10.4 KB
Newer Older
1
2
c this is <gremlin_help.f>
c------------------------------------------------------------------------------
thomas.forbriger's avatar
thomas.forbriger committed
3
c $Id$
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
c    06/05/02   V1.5   introduced sqense and ssense
thomas.forbriger's avatar
thomas.forbriger committed
17
c    07/05/02   V1.6   new data weight routines
thomas.forbriger's avatar
thomas.forbriger committed
18
19
c    07/08/02   V1.7   - found old command npa
c                      - explain rms and err parameter in resolution analysis
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
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'
thomas.forbriger's avatar
thomas.forbriger committed
36
      print *,' med            edit current model (model will be written to'
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
      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'
56
      print *,' weights th1,th2,f1,f2'
thomas.forbriger's avatar
thomas.forbriger committed
57
58
59
60
      print *,'                calculate misfit-dependend data weights'
      print *,'                you are invited to have a look in'
      print *,'                subroutine dat_mmweights to learn about'
      print *,'                the calculation rule ;-)'
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
      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'
thomas.forbriger's avatar
thomas.forbriger committed
84
85
86
87
      print *,'                  ''para'' is the value taken from'
      print *,'                  the x2c-abscissa.'
      print *,' npa para       like tpa, but calculate partial'
      print *,'                  derivatives first'
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
      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'
thomas.forbriger's avatar
thomas.forbriger committed
112
      print *,' qtpd ival      like tpd, but take existing values'
113
114
      print *,' tpda           plot travel time partial derivatives for all parameters'
      print *,' gpd ival       plot green partial derivatives for parameter ival'
thomas.forbriger's avatar
thomas.forbriger committed
115
      print *,' qgpd ival      like gpd, but take existing values'
116
      print *,' gpda           plot green partial derivatives for all parameters'
117
118
      print *,' wgpd name      write green partial derivatives to files with'
      print *,'                base ''name'' '
119
120
      print *,' ense f,t,nu,rms   rate model parameters with index from ''f'' to'
      print *,'                ''t'' when optimizing all other free parameters'
thomas.forbriger's avatar
thomas.forbriger committed
121
122
123
124
125
      print *,'                stabilized with ''nu'' and accepting an'
      print *,'                increase in the rms-error'
      print *,'                by ''rms'' in a linear approximation'
      print *,'                ''rms'' is given as a fraction of the'
      print *,'                rms misfit of the reference model'
126
127
      print *,' sense f,t,nu,rms  just like ''ense'' but puts rating for each'
      print *,'                parameter on a different page'
128
129
130
      print *,' wense f,t,nu,rms,name'
      print *,'                just like ''sense'' but puts rating for each'
      print *,'                parameter to a file with base ''name'' '
131
132
133
134
      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' 
thomas.forbriger's avatar
thomas.forbriger committed
135
      print *,'                square error given as a fraction of misfit'
136
      print *,'                of the reference synthetics'
137
138
139
140
141
142
      print *,' ssense f,t,nu,err'
      print *,'                just like ''sense'' but in misfit,'
      print *,'                not in rms mode'
      print *,' sqense f,t,nu,err'
      print *,'                just like ''ense'' but in misfit,'
      print *,'                not in rms mode'
thomas.forbriger's avatar
thomas.forbriger committed
143
144
145
146
      print *,' qense, qsense, qssense, qsqense'
      print *,'                these are versions of ''ense'','
      print *,'                ''sense'', ''ssense'', and ''sqense'' '
      print *,'                that use preexisting partial derivatives'
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
      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"'
thomas.forbriger's avatar
thomas.forbriger committed
181
      print *,'save file      write current work model to "file"'
182
183
184
      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'
thomas.forbriger's avatar
thomas.forbriger committed
185
186
187
188
      print *,'wwrt file      write original traveltime weights'
      print *,'wwct file      write calculated traveltime weights'
      print *,'wwrg file      write original green weights'
      print *,'wwcg file      write calculated green weights'
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
      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 -----