'Sample input data len-mod.dat: data for the model. 1999, March 21.'
0 1 1 / (Cartesian coordinates, velocities, loss factors)
0 60 -20 20 -25 10 (boundaries of the model)
4 surfaces
5 simple blocks:
-1 2 -3 /
-1 3 /
1 2 -4 /
4 /
-2 / (free space)
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 5 I I NCB I 0 3 I
I-----------I-----------I I-----------I-----------I
I S.Block 1 I 1 8 I I C.Block 1 I 1 5 I
I S.Block 2 I 2 10 I I C.Block 2 I 2 6 I
I S.Block 3 I 3 13 I I C.Block 3 I 3 7 I
I S.Block 4 I 4 14 I I-----------I-----------I
I S.Block 4 I 5 15 I I C.Block 1 I 4 1 I
I-----------I-----------I I I 5 3 I
I I 6 -1 I I-----------I-----------I
I S.Block 1 I 7 2 I I C.Block 2 I 6 2 I
I I 8 -3 I I-----------I-----------I
I-----------I-----------I I C.Block 3 I 7 4 I
I S.Block 2 I 9 -1 I I-----------I-----------I
I I 10 3 I
I-----------I-----------I
I I 11 1 I
I S.Block 3 I 12 2 I
I I 13 -4 I
I-----------I-----------I
I S.Block 4 I 14 4 I
I-----------I-----------I
I S.Block 5 I 15 -2 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 60 I
I I Class 2 I C.Block 2 I 8 107 I
I I I C.Block 3 I 9 154 I
I------------I------------I------------I------------I
I I Surface 1 I I 10 207 I
I I------------I------------I------------I
I I Surface 2 I I 11 242 I
I I------------I------------I------------I
I I Surface 3 I I 12 295 I
I I------------I------------I------------I
I I Surface 4 I I 13 337 I
I I------------I------------I------------I
I I I P velocity I 14 358 I
I I I Density I 15 379 I
I I C.Block 1 I No data I 16 379 I
I (D) I I No data I 17 379 I
I Functions I I No data I 18 379 I
I I I . I . . I
I I I : I : : I
I I I No data I 60 379 I
I I------------I------------I------------I
I I I P velocity I 61 386 I
I I I Density I 62 393 I
I I C.Block 2 I No data I 63 393 I
I I I No data I 64 393 I
I I I No data I 65 393 I
I I I . I . . I
I I I : I : : I
I I I No data I107 393 I
I I------------I------------I------------I
I I I P velocity I108 400 I
I I I Density I109 407 I
I I C.Block 3 I No data I110 407 I
I I I No data I111 407 I
I I I No data I112 407 I
I I I . I . . I
I I I : I : : I
I I I No data I154 407 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 155 1 I I 208 1 I I 243 1 I I 296 1 I
I 156 1.000 I I 209 1.000 I I 244 1.000 I I 297 1.000 I
I 157 1 I I 210 1 I I 245 1 I I 298 1 I
I 158 2 I I 211 -3 I I 246 2 I I 299 -3 I
I 159 -3 I I 212 0 I I 247 -3 I I 300 0 I
I 160 0.000 I I 213 0.000 I I 248 0.000 I I 301 0.000 I
I 161 3 I I 214 4 I I 249 3 I I 302 5 I
I 162 3 I I 215 0.000 I I 250 3 I I 303 0.000 I
I 163 0.000 I I 216 20.000 I I 251 0.000 I I 304 10.000 I
I 164 20.000 I I 217 40.000 I I 252 20.000 I I 305 25.000 I
I 165 40.000 I I 218 60.000 I I 253 30.000 I I 306 40.000 I
I 166 -20.000 I I 219 3.000 I I 254 -20.000 I I 307 60.000 I
I 167 0.000 I I 220 2.822 I I 255 0.000 I I 308 -16.000 I
I 168 20.000 I I 221 -1.289 I I 256 20.000 I I 309 -10.691 I
I 169 0.000 I I 222 3.500 I I 257 -12.000 I I 310 -14.537 I
I 170 0.000 I I 223 .004 I I 258 -7.200 I I 311 -10.231 I
I 171 0.000 I I 224 .000 I I 259 -12.000 I I 312 -15.000 I
I 172 0.000 I I 225 .002 I I 260 -8.000 I I 313 .005 I
I 173 -12.000 I I 226 .250 I I 261 2.200 I I 314 .000 I
I 174 0.000 I I 227 .167 I I 262 -8.000 I I 315 .007 I
I 175 0.000 I I 228 .000 I I 263 -12.000 I I 316 .345 I
I 176 0.000 I I 229 -.007 I I 264 -7.200 I I 317 .257 I
I 177 0.000 I I 230 .004 I I 265 -12.000 I I 318 .000 I
I 178 .002 I I 231 .000 I I 266 .004 I I 319 -.020 I
I 179 .000 I I 232 .250 I I 267 .000 I I 320 .008 I
I 180 .002 I I 233 .004 I I 268 .004 I I 321 .000 I
I 181 .167 I I 234 -.007 I I 269 .067 I I 322 .292 I
I 182 .167 I I 235 .000 I I 270 .067 I I 323 .009 I
I 183 .000 I I 236 .250 I I 271 .000 I I 324 -.014 I
I 184 -.007 I I 237 .000 I I 272 -.015 I I 325 .005 I
I 185 .000 I I 238 .002 I I 273 .000 I I 326 .151 I
I 186 .000 I I 239 .000 I I 274 .000 I I 327 .345 I
I 187 .000 I I 240 .004 I I 275 .000 I I 328 .006 I
I 188 .002 I I 241 .167 I I 276 .005 I I 329 -.010 I
I 189 .000 I I 242 .250 I I 277 .000 I I 330 .000 I
I 190 .002 I I-------------I I 278 .005 I I 331 .229 I
I 191 .167 I I 279 .333 I I 332 .000 I
I 192 .167 I I 280 .333 I I 333 .003 I
I 193 .002 I I 281 .002 I I 334 .000 I
I 194 .000 I I 282 .000 I I 335 .009 I
I 195 .002 I I 283 .002 I I 336 .117 I
I 196 .167 I I 284 .167 I I 337 .151 I
I 197 .167 I I 285 .167 I I-------------I
I 198 .000 I I 286 .000 I
I 199 -.007 I I 287 -.007 I
I 200 .000 I I 288 .000 I
I 201 .000 I I 289 .000 I
I 202 .000 I I 290 .000 I
I 203 .002 I I 291 .002 I
I 204 .000 I I 292 .000 I
I 205 .002 I I 293 .002 I
I 206 .167 I I 294 .167 I
I 207 .167 I I 295 .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 338 1 I I 359 3 I I 380 1 I I 394 1 I
I 339 1.000 I I 360 1.000 I I 381 1.000 I I 395 1.000 I
I 340 3 I I 361 4 I I 382 0 I I 396 0 I
I 341 0 I I 362 0 I I 383 0 I I 397 0 I
I 342 0 I I 363 0 I I 384 0 I I 398 0 I
I 343 0.000 I I 364 0.000 I I 385 0.000 I I 399 0.000 I
I 344 2 I I 365 2 I I 386 6.000 I I 400 7.000 I
I 345 -20.000 I I 366 0.000 I I-------------I I-------------I
I 346 0.000 I I 367 10.000 I
I 347 6.000 I I 368 1.700 I C.B.2: Density: C.B.3: Density:
I 348 4.000 I I 369 3.700 I I-------------I I-------------I
I 349 .002 I I 370 1.000 I I K IPAR(K) I I K IPAR(K) I
I 350 .000 I I 371 .000 I I RPAR(K) I I RPAR(K) I
I 351 .000 I I 372 .000 I I-------------I I-------------I
I 352 .167 I I 373 .167 I I 387 3 I I 401 3 I
I 353 .000 I I 374 .000 I I 388 1.000 I I 402 1.000 I
I 354 .000 I I 375 .000 I I 389 0 I I 403 0 I
I 355 .000 I I 376 .000 I I 390 0 I I 404 0 I
I 356 .002 I I 377 1.000 I I 391 0 I I 405 0 I
I 357 .000 I I 378 .000 I I 392 0.000 I I 406 0.000 I
I 358 .167 I I 379 .167 I I 393 2.900 I I 407 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,154) are increased by 9 and shifted to
(IPAR(K),K=16,156).
(E):
Parameters (*PAR(K),K=155,337) are shifted to
(*PAR(K),K=157,339).
Parameters (*PAR(K),K=338,407) are shifted to
(*PAR(K),K=347,416).
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 346 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 340 1 I
I 341 1.000 I
I 342 -1 I
I 343 0 I
I 344 0 I
I 345 0.000 I
I 346 0.000 I
I-------------I
========================================================================