'SAMPLE INPUT DATA MODEL.DAT: DATA FOR THE MODEL. 1996, JULY 29.'
0 1 1 / (Cartesian coordinates, velocities, loss factors)
0 60 -20 20 -25 10 (boundaries of the model)
4 SURFACES
4 SIMPLE BLOCKS:
-1 2 -3 /
-1 3 /
1 2 -4 /
4 /
3 COMPLEX BLOCKS:
1 3 /
2 /
4 /
'SURFACE' 1
1 2 -3 0 / (i.e. W(X1,X2)-X3=0, tension=0)
3 3 (numbers of grid points)
0 20 40 (X1 grid coordinates)
-20 0 20 (X2 grid coordinates)
0 0 0
0 -12 0
0 0 0 (X3 coordinates at grid points)
'SURFACE' 2
1 -3 0 0 / (i.e. W(X1)-X3=0, tension=0)
4 (numbers of grid points)
0 20 40 60 (X1 grid coordinates)
3 3 0 3.5 (X3 coordinates at grid points)
'SURFACE' 3
1 2 -3 0 / (i.e. W(X1,X2)-X3=0, tension=0)
3 3 (numbers of grid points)
0 20 30 (X1 grid coordinates)
-20 0 20 (X2 grid coordinates)
-12 -12 -12
-12 -5 -12
-12 -12 -12 (X3 coordinates at grid points)
'SURFACE' 4
1 -3 0 0 / (i.e. W(X1)-X3=0, tension=0)
5 (numbers of grid points)
0 10 25 40 60 (X1 grid coordinates)
-16 -17 -20 -17 -15 (X3 coordinates at grid points)
'END OF SURFACES' /
'COMPLEX BLOCK' 1
'VP ' 1 (P wave velocity)
3 0 0 0 / (i.e. VP=W(X3), tension=0)
2 (numbers of grid points)
0 -20 (X3 grid coordinates)
4 6 (velocities at grid points)
'DEN' 1 (density)
4 0 0 0 / (i.e. density=W(VP), tension=0)
2 (numbers of grid points)
0 10 (VP grid coordinates)
1.7 3.7 (densities at grid points, density=1.7+0.2*VP)
'COMPLEX BLOCK' 2
'VP ' 1 (P wave velocity)
0 0 0 0 / (i.e. VP=constant, tension=0)
6 (the value of velocity)
'DEN' 1 (density)
0 0 0 0 / (i.e. density=constant, tension=0)
2.9 (the value of density)
'COMPLEX BLOCK' 3
'VP ' 1 (P wave velocity)
0 0 0 0 / (i.e. VP=constant, tension=0)
7 (the value of velocity)
'DEN' 1 (density)
0 0 0 0 / (i.e. density=constant, tension=0)
3.1 (the value of density)
'END OF COMPLEX BLOCKS, END OF THE INPUT DATA FOR THE MODEL' /
------------------------------------------------------------------------
This 3-D model consists of two layers and of the lenticular inclusion
with edges, situated in the upper layer, see the schematic Figure 5 in
V.Cerveny, L.Klimes, I.Psencik: Complete Seismic-Ray Tracing in 3-D
Structures (In: D.Doornbos, ed.: Seismological Algorithms, Academic
Press 1988). For the general description of the input data format see
the subroutine packages 'model', 'srfc' and 'parm'. P wave velocity
in the upper layer (Complex Block 1) increases linearly with depth and
has no lateral variations. The lenticular inclusion (Complex Block 2)
and the bottom layer (Complex Block 3) are homogeneous. S wave
velocities and loss factors are not specified by these input data, the
default values (VS=0.57735*VP, zero loss factors) are assumed. The
comments in brackets (starting from the column 23) in these input data
have no influence on reading the data.
------------------------------------------------------------------------
Debugging information:
The input data are stored in the common blocks /MODELT/, /MODELC/ and
/VALC/. For the description of the common blocks /MODELT/ and
/MODELC/ refer to the subroutine package 'model', for the common block
/VALC/ refer to the subroutine package 'val'. The arrays KSB and KCB
of the common block /MODELC/ and the arrays IPAR and RPAR of the
common block /VALC/ are filled by the parameters of the model in the
following way:
Arrays KSB and KCB of the common block /MODELC/:
I-----------I-----------I I-----------I-----------I
I I K KSB(K) I I I K KCB(K) I
I-----------I-----------I I-----------I-----------I
I NSB I 0 4 I I NCB I 0 3 I
I-----------I-----------I I-----------I-----------I
I S.Block 1 I 1 7 I I C.Block 1 I 1 5 I
I S.Block 2 I 2 9 I I C.Block 2 I 2 6 I
I S.Block 3 I 3 12 I I C.Block 3 I 3 7 I
I S.Block 4 I 4 13 I I-----------I-----------I
I-----------I-----------I I C.Block 1 I 4 1 I
I I 5 -1 I I I 5 3 I
I S.Block 1 I 6 2 I I-----------I-----------I
I I 7 -3 I I C.Block 2 I 6 2 I
I-----------I-----------I I-----------I-----------I
I S.Block 2 I 8 -1 I I C.Block 3 I 7 4 I
I I 9 3 I I-----------I-----------I
I-----------I-----------I
I I 10 1 I
I S.Block 3 I 11 2 I
I I 12 -4 I
I-----------I-----------I
I S.Block 4 I 13 4 I
I-----------I-----------I
Arrays IPAR and RPAR of the common block /VALC/:
I--------------------------------------I------------I
I Description I K IPAR(K) I
I--------------------------------------I------------I
I (A) Number of the classes I 0 2 I
I-------------------------I------------I------------I
I (B) Classes I Class 1 I 1 6 I
I I Class 2 I 2 9 I
I------------I------------I------------I------------I
I I I Surface 1 I 3 10 I
I I Class 1 I Surface 2 I 4 11 I
I I I Surface 3 I 5 12 I
I (C) Groups I I Surface 4 I 6 13 I
I I------------I------------I------------I
I I I C.Block 1 I 7 18 I
I I Class 2 I C.Block 2 I 8 23 I
I I I C.Block 3 I 9 28 I
I------------I------------I------------I------------I
I I Surface 1 I I 10 81 I
I I------------I------------I------------I
I I Surface 2 I I 11 116 I
I I------------I------------I------------I
I I Surface 3 I I 12 169 I
I I------------I------------I------------I
I I Surface 4 I I 13 211 I
I I------------I------------I------------I
I I I P velocity I 14 232 I
I I I Density I 15 253 I
I I C.Block 1 I No data I 16 253 I
I (D) I I No data I 17 253 I
I Functions I I No data I 18 253 I
I I------------I------------I------------I
I I I P velocity I 19 260 I
I I I Density I 20 267 I
I I C.Block 2 I No data I 21 267 I
I I I No data I 22 267 I
I I I No data I 23 267 I
I I------------I------------I------------I
I I I P velocity I 24 274 I
I I I Density I 25 281 I
I I C.Block 3 I No data I 26 281 I
I I I No data I 27 281 I
I I I No data I 28 281 I
I------------I------------I------------I------------I
(E) Parameters of individual functions:
Surface 1: Surface 2: Surface 3: Surface 4:
I-------------I I-------------I I-------------I I-------------I
I K IPAR(K) I I K IPAR(K) I I K IPAR(K) I I K IPAR(K) I
I RPAR(K) I I RPAR(K) I I RPAR(K) I I RPAR(K) I
I-------------I I-------------I I-------------I I-------------I
I 29 1 I I 82 1 I I 117 1 I I 170 1 I
I 30 1.000 I I 83 1.000 I I 118 1.000 I I 171 1.000 I
I 31 1 I I 84 1 I I 119 1 I I 172 1 I
I 32 2 I I 85 -3 I I 120 2 I I 173 -3 I
I 33 -3 I I 86 0 I I 121 -3 I I 174 0 I
I 34 0.000 I I 87 0.000 I I 122 0.000 I I 175 0.000 I
I 35 3 I I 88 4 I I 123 3 I I 176 5 I
I 36 3 I I 89 0.000 I I 124 3 I I 177 0.000 I
I 37 0.000 I I 90 20.000 I I 125 0.000 I I 178 10.000 I
I 38 20.000 I I 91 40.000 I I 126 20.000 I I 179 25.000 I
I 39 40.000 I I 92 60.000 I I 127 30.000 I I 180 40.000 I
I 40 -20.000 I I 93 3.000 I I 128 -20.000 I I 181 60.000 I
I 41 0.000 I I 94 2.822 I I 129 0.000 I I 182 -16.000 I
I 42 20.000 I I 95 -1.289 I I 130 20.000 I I 183 -10.691 I
I 43 0.000 I I 96 3.500 I I 131 -12.000 I I 184 -14.537 I
I 44 0.000 I I 97 .004 I I 132 -7.200 I I 185 -10.231 I
I 45 0.000 I I 98 .000 I I 133 -12.000 I I 186 -15.000 I
I 46 0.000 I I 99 .002 I I 134 -8.000 I I 187 .005 I
I 47 -12.000 I I 100 .250 I I 135 2.200 I I 188 .000 I
I 48 0.000 I I 101 .167 I I 136 -8.000 I I 189 .007 I
I 49 0.000 I I 102 .000 I I 137 -12.000 I I 190 .345 I
I 50 0.000 I I 103 -.007 I I 138 -7.200 I I 191 .257 I
I 51 0.000 I I 104 .004 I I 139 -12.000 I I 192 .000 I
I 52 .002 I I 105 .000 I I 140 .004 I I 193 -.020 I
I 53 .000 I I 106 .250 I I 141 .000 I I 194 .008 I
I 54 .002 I I 107 .004 I I 142 .004 I I 195 .000 I
I 55 .167 I I 108 -.007 I I 143 .067 I I 196 .292 I
I 56 .167 I I 109 .000 I I 144 .067 I I 197 .009 I
I 57 .000 I I 110 .250 I I 145 .000 I I 198 -.014 I
I 58 -.007 I I 111 .000 I I 146 -.015 I I 199 .005 I
I 59 .000 I I 112 .002 I I 147 .000 I I 200 .151 I
I 60 .000 I I 113 .000 I I 148 .000 I I 201 .345 I
I 61 .000 I I 114 .004 I I 149 .000 I I 202 .006 I
I 62 .002 I I 115 .167 I I 150 .005 I I 203 -.010 I
I 63 .000 I I 116 .250 I I 151 .000 I I 204 .000 I
I 64 .002 I I-------------I I 152 .005 I I 205 .229 I
I 65 .167 I I 153 .333 I I 206 .000 I
I 66 .167 I I 154 .333 I I 207 .003 I
I 67 .002 I I 155 .002 I I 208 .000 I
I 68 .000 I I 156 .000 I I 209 .009 I
I 69 .002 I I 157 .002 I I 210 .117 I
I 70 .167 I I 158 .167 I I 211 .151 I
I 71 .167 I I 159 .167 I I-------------I
I 72 .000 I I 160 .000 I
I 73 -.007 I I 161 -.007 I
I 74 .000 I I 162 .000 I
I 75 .000 I I 163 .000 I
I 76 .000 I I 164 .000 I
I 77 .002 I I 165 .002 I
I 78 .000 I I 166 .000 I
I 79 .002 I I 167 .002 I
I 80 .167 I I 168 .167 I
I 81 .167 I I 169 .167 I
I-------------I I-------------I
C.Block 1: VP: C.B.1: Density: C.Block 2: VP: C.Block 3: VP:
I-------------I I-------------I I-------------I I-------------I
I K IPAR(K) I I K IPAR(K) I I K IPAR(K) I I K IPAR(K) I
I RPAR(K) I I RPAR(K) I I RPAR(K) I I RPAR(K) I
I-------------I I-------------I I-------------I I-------------I
I 212 1 I I 233 3 I I 254 1 I I 268 1 I
I 213 1.000 I I 234 1.000 I I 255 1.000 I I 269 1.000 I
I 214 3 I I 235 4 I I 256 0 I I 270 0 I
I 215 0 I I 236 0 I I 257 0 I I 271 0 I
I 216 0 I I 237 0 I I 258 0 I I 272 0 I
I 217 0.000 I I 238 0.000 I I 259 0.000 I I 273 0.000 I
I 218 2 I I 239 2 I I 260 6.000 I I 274 7.000 I
I 219 -20.000 I I 240 0.000 I I-------------I I-------------I
I 220 0.000 I I 241 1.000 I
I 221 6.000 I I 242 1.700 I C.B.2: Density: C.B.3: Density:
I 222 4.000 I I 243 1.900 I I-------------I I-------------I
I 223 .002 I I 244 1.000 I I K IPAR(K) I I K IPAR(K) I
I 224 .000 I I 245 .000 I I RPAR(K) I I RPAR(K) I
I 225 .000 I I 246 .000 I I-------------I I-------------I
I 226 .167 I I 247 .167 I I 261 3 I I 275 3 I
I 227 .000 I I 248 .000 I I 262 1.000 I I 276 1.000 I
I 228 .000 I I 249 .000 I I 263 0 I I 277 0 I
I 229 .000 I I 250 .000 I I 264 0 I I 278 0 I
I 230 .002 I I 251 1.000 I I 265 0 I I 279 0 I
I 231 .000 I I 252 .000 I I 266 0.000 I I 280 0.000 I
I 232 .167 I I 253 .167 I I 267 2.900 I I 281 3.100 I
I-------------I I-------------I I-------------I I-------------I
Attention: The storage of the above model coefficients in the common
block /VALC/ may be changed when reading additional input data,
especially when defining auxiliary surfaces in the model.
Example:
If calling subroutine SRFC1 once again to read the data
'AUXILIARY SURFACE' 5
-1 0 0 0 (I.E. W-X1=0, TENSION=0)
0. (W=0)
'END OF AUXILIARY SURFACES' /
in order to define additional surface 5, the above contents of arrays
IPAR and RPAR of the common block /VALC/ would become modified in the
following way:
(A),(B),(C):
Parameters (IPAR(K),K=1,6) are increased by 1.
(C),(D):
Parameters (IPAR(K),K=7,13) are increased by 2 and shifted to
(IPAR(K),K=8,14).
(D):
Parameters (IPAR(K),K=14,28) are increased by 9 and shifted to
(IPAR(K),K=16,30).
(E):
Parameters (*PAR(K),K=29,211) are shifted to
(*PAR(K),K=31,213).
Parameters (*PAR(K),K=214,281) are shifted to
(*PAR(K),K=223,290).
New data added to the common block /VALC/:
I--------------------------------------I------------I
I Description I K IPAR(K) I
I--------------------------------------I------------I
I (C) Groups I Class 1 I Surface 5 I 7 15 I
I------------I------------I------------I------------I
I (D) Funct. I Surface 5 I I 15 220 I
I------------I------------I------------I------------I
(E) Parameters of individual functions:
Surface 5:
I-------------I
I K IPAR(K) I
I RPAR(K) I
I-------------I
I 214 1 I
I 215 1.000 I
I 216 -1 I
I 217 0 I
I 218 0 I
I 219 0.000 I
I 220 0.000 I
I-------------I
========================================================================