gremlin_help.f 11.1 KB
Newer Older
1
2
c this is <gremlin_help.f>
c------------------------------------------------------------------------------
thomas.forbriger's avatar
thomas.forbriger committed
3
c $Id$
4
c
5
c Copyright 1998, 2010 by Thomas Forbriger (IfG Stuttgart)
6
7
8
c
c Here is some extra help on gremlin commands
c
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
c ----
c This program is free software; you can redistribute it and/or modify
c it under the terms of the GNU General Public License as published by
c the Free Software Foundation; either version 2 of the License, or
c (at your option) any later version. 
c 
c This program is distributed in the hope that it will be useful,
c but WITHOUT ANY WARRANTY; without even the implied warranty of
c MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
c GNU General Public License for more details.
c 
c You should have received a copy of the GNU General Public License
c along with this program; if not, write to the Free Software
c Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA
c ----
c
c
26
27
28
29
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
30
31
c    24/05/00   V1.3   - introduced reso option wgpd
c                      - introduced reso option wense
32
c    19/04/02   V1.4   introduced command swense
33
c    06/05/02   V1.5   introduced sqense and ssense
thomas.forbriger's avatar
thomas.forbriger committed
34
c    07/05/02   V1.6   new data weight routines
thomas.forbriger's avatar
thomas.forbriger committed
35
36
c    07/08/02   V1.7   - found old command npa
c                      - explain rms and err parameter in resolution analysis
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
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
53
      print *,' med            edit current model (model will be written to'
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
      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'
73
      print *,' weights th1,th2,f1,f2'
thomas.forbriger's avatar
thomas.forbriger committed
74
75
76
77
      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 ;-)'
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
      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
101
102
103
104
      print *,'                  ''para'' is the value taken from'
      print *,'                  the x2c-abscissa.'
      print *,' npa para       like tpa, but calculate partial'
      print *,'                  derivatives first'
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
      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
129
      print *,' qtpd ival      like tpd, but take existing values'
130
131
      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
132
      print *,' qgpd ival      like gpd, but take existing values'
133
      print *,' gpda           plot green partial derivatives for all parameters'
134
135
      print *,' wgpd name      write green partial derivatives to files with'
      print *,'                base ''name'' '
136
137
      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
138
139
140
141
142
      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'
143
144
      print *,' sense f,t,nu,rms  just like ''ense'' but puts rating for each'
      print *,'                parameter on a different page'
145
146
147
      print *,' wense f,t,nu,rms,name'
      print *,'                just like ''sense'' but puts rating for each'
      print *,'                parameter to a file with base ''name'' '
148
149
150
151
      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
152
      print *,'                square error given as a fraction of misfit'
153
      print *,'                of the reference synthetics'
154
155
156
157
158
159
      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
160
161
162
163
      print *,' qense, qsense, qssense, qsqense'
      print *,'                these are versions of ''ense'','
      print *,'                ''sense'', ''ssense'', and ''sqense'' '
      print *,'                that use preexisting partial derivatives'
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
189
190
191
192
193
194
195
196
197
      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
198
      print *,'save file      write current work model to "file"'
199
200
201
      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
202
203
204
205
      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'
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
      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 -----