SPEED
MAKE_ELTENSOR_FOR_CASES.f90
Go to the documentation of this file.
1! Copyright (C) 2012 The SPEED FOUNDATION
2! Author: Ilario Mazzieri
3!
4! This file is part of SPEED.
5!
6! SPEED is free software; you can redistribute it and/or modify it
7! under the terms of the GNU Affero General Public License as
8! published by the Free Software Foundation, either version 3 of the
9! License, or (at your option) any later version.
10!
11! SPEED is distributed in the hope that it will be useful, but
12! WITHOUT ANY WARRANTY; without even the implied warranty of
13! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14! Affero General Public License for more details.
15!
16! You should have received a copy of the GNU Affero General Public License
17! along with SPEED. If not, see <http://www.gnu.org/licenses/>.
18
49
50 subroutine make_eltensor_for_cases(tcase,vcase,&
51 nn,rho_el,lambda_el,mu_el,gamma_el,&
52 nn_loc, zs_elev, zs_all, vs_nodes, thick_nodes, &
53 cs_nnz_loc, cs_loc, ielem, &
54 sub_tag_all, zs, mpi_id, local_n_num, &
55 damping_type, qs, qp, &
56 xs, ys, check_case, label_case)
57
58
59 implicit none
60
61 integer*4 :: tcase, check_case, label_case
62 integer*4 :: vcase, mpi_id
63 integer*4 :: nn
64 integer*4 :: p, q, r, ic
65 integer*4 :: nn_loc
66 integer*4 :: cs_nnz_loc
67 integer*4 :: is,in,ielem , damping_type
68
69 integer*4, dimension(nn_loc) :: local_n_num
70 integer*4, dimension(nn_loc) :: sub_tag_all
71 integer*4, dimension(0:cs_nnz_loc) :: cs_loc
72
73 real*8 :: depth, depth_real, qs_all, qp_all, thickness, vs30, pig
74 real*8 :: vs,vp, rho,lambda,mu,gamma,ni, qs,qp
75 real*8 :: x1,y1,x2,y2,coef_a, coef_b, coef_c, numer, den, distance, f_distance
76
77 real*8, dimension(nn_loc) :: zs_elev
78 real*8, dimension(nn_loc) :: zs_all
79 real*8, dimension(nn_loc) :: vs_nodes, thick_nodes
80
81 real*8, dimension(nn_loc) :: zs, xs, ys
82
83 real*8, dimension(nn,nn,nn) :: rho_el,lambda_el,mu_el,gamma_el, qs_h, qp_h
84
85 real*8 :: stat_id1_x, stat_id1_y, stat_id2_x, stat_id2_y
86
87 character*70 :: filename
88 character*5 :: filesuffix
89
90
91 pig = 4.d0*datan(1.d0);
92
93! STRESS CALCULATION
94
95 if (check_case .eq. 1) then
96
97 filesuffix = '.dat'
98 write(filename, '(A,I5.5,A5)') 'NHCheck', mpi_id, filesuffix
99 !write(*,*) mpi_id, ielem, filename
100 !read(*,*)
101 open(1000 + mpi_id,file=filename,position='APPEND')
102 endif
103
104 do r = 1,nn
105 do q = 1,nn
106 do p = 1,nn
107 is = nn*nn*(r -1) +nn*(q -1) +p
108 in = cs_loc(cs_loc(ielem -1) +is)
109 ic = in
110
111 if (ic .eq. 0 ) write(*,*) 'Error in MAKE_ELTENSOR_FOR_CASES '
112
113 if (tcase.eq.1) then
114 ! CASE 1: GRENOBLE 1
115 call make_mech_prop_case_001(rho,lambda,mu,gamma,qs,qp, & !outputs
116 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
117 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
118
119 elseif (tcase.eq.2) then
120 ! CASE 2: GRENOBLE 2
121 call make_mech_prop_case_002(rho,lambda,mu,gamma,qs,qp, & !outputs
122 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
123 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
124
125 elseif (tcase.eq.3) then
126 ! CASE 3: GUBBIO
127 call make_mech_prop_case_003(rho,lambda,mu,gamma,qs,qp, & !outputs
128 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
129 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
130
131 elseif (tcase.eq.4) then
132 ! CASE 4: SULMONA NOT honoring
133 call make_mech_prop_case_004(rho,lambda,mu,gamma,qs,qp, & !outputs
134 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
135 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
136
137 elseif (tcase.eq.5) then
138 ! CASE 5: VOLVI CASHIMA benchmark - NOT honoring
139 call make_mech_prop_case_005(rho,lambda,mu,gamma,qs,qp, & !outputs
140 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
141 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
142
143
144 elseif (tcase.eq.6) then !NOT FOUND
145 ! CASE 6: FRIULI (Tagliamento river valley)
146 call make_mech_prop_case_006(rho,lambda,mu,gamma,qs,qp, & !outputs
147 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
148 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
149
150 elseif (tcase.eq.7) then
151 ! CASE 7: AQUILA (Evangelista et al. 2017)
152 call make_mech_prop_case_007(rho,lambda,mu,gamma,qs,qp, & !outputs
153 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
154 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
155
156 elseif (tcase.eq.8) then
157 ! CASE 8: SANTIAGO
158 call make_mech_prop_case_008(rho,lambda,mu,gamma,qs,qp, & !outputs
159 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
160 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
161
162 elseif (tcase.eq.10) then
163 ! CASE 10: CHRISTCHURCH INGV - Staircase
164 call make_mech_prop_case_010(rho,lambda,mu,gamma,qs,qp, & !outputs
165 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
166 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
167
168 elseif (tcase.eq.11) then
169 ! CASE 11: CHRISTCHURCH (last model)
170 call make_mech_prop_case_011(rho,lambda,mu,gamma,qs,qp, & !outputs
171 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
172 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
173
174 elseif (tcase.eq.12) then
175 ! CASE 12: PO PLAIN (new model)
176
177 !straight line passing by (x1,y1,0), (x2,y2,0)
178 x1 = 654957.002352; y1 = 4974060.299450;
179 x2 = 688420.525202; y2 = 4957613.600935;
180 !coefficient of the line ax + by + c = 0
181 coef_a = 1.d0/(x2-x1);
182 coef_b = 1.d0/(y1-y2);
183 coef_c = - y1/(y1-y2) + x1/(x1-x2);
184 !distance between (x,y,0) and the line ax + by + c = 0
185 numer = coef_a*xs(ic) + coef_b*ys(ic) + coef_c
186 den = dsqrt(coef_a**2 + coef_b**2)
187 distance = dabs(numer/den)
188 f_distance = 150.d0 + 1850.d0/(1.d0 + dexp(-0.0012*(distance-5000.d0)));
189
190 call make_mech_prop_case_012(rho,lambda,mu,gamma,qs,qp, & !outputs
191 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
192 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic), f_distance)
193
194
195 elseif (tcase.eq.13) then
196 ! CASE 13: PO PLAIN-BEDROCK
197 call make_mech_prop_case_013(rho,lambda,mu,gamma,qs,qp, & !outputs
198 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
199 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
200
201 elseif (tcase.eq.14) then
202 ! CASE 14: WELLINGTON (simplified model, Benites et al. 2005)
203 call make_mech_prop_case_014(rho,lambda,mu,gamma,qs,qp, & !outputs
204 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
205 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
206
207 elseif (tcase.eq.15) then
208 ! CASE 15: MARSICA
209 call make_mech_prop_case_015(rho,lambda,mu,gamma,qs,qp, & !outputs
210 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
211 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
212
213 elseif (tcase.eq.16) then
214 ! CASE 16: ISTANBUL
215 call make_mech_prop_case_016(rho, lambda, mu, gamma, qs, qp, & !outputs
216 xs(ic), ys(ic), zs(ic),&
217 zs_elev(ic), zs_all(ic), &
218 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
219
220 elseif (tcase.eq.18) then
221 ! CASE 18: BEIJING -- SIMPLIFIED MODEL
222 call make_mech_prop_case_018(rho,lambda,mu,gamma,qs,qp, & !outputs
223 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
224 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
225
226 elseif (tcase.eq.19) then
227 ! CASE 19: SALONICCO
228 call make_mech_prop_case_019(rho,lambda,mu,gamma,qs,qp, & !outputs
229 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
230 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
231
232 elseif (tcase.eq.20) then
233 ! CASE 20: ATENE
234 call make_mech_prop_case_020(rho,lambda,mu,gamma,qs,qp, & !outputs
235 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
236 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
237
238 elseif (tcase.eq.21) then
239 ! CASE 21: BEIJING
240 call make_mech_prop_case_021(rho,lambda,mu,gamma,qs,qp, & !outputs
241 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
242 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
243
244 elseif (tcase.eq.22) then
245 !-------------------------------------------------------------------
246 ! CASE 22: NORCIA
247 call make_mech_prop_case_022(rho,lambda,mu,gamma,qs,qp, & !outputs
248 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
249 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
250
251 elseif (tcase.eq.27) then
252 !-------------------------------------------------------------------
253 ! CASE 27: AQUILA-BEDROCK OUTCROP BEDROCK
254 call make_mech_prop_case_027(rho,lambda,mu,gamma,qs,qp, & !outputs
255 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
256 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
257
258 elseif (tcase.eq.28) then
259 !-------------------------------------------------------------------
260 ! CASE 28: NORCIA OUTCROP BEDROCK
261 call make_mech_prop_case_028(rho,lambda,mu,gamma,qs,qp, & !outputs
262 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
263 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
264
265 elseif (tcase.eq.29) then
266 ! CASE 29: THESSALONIKI-OUTCROP BEDROCK
267 call make_mech_prop_case_029(rho,lambda,mu,gamma,qs,qp, & !outputs
268 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
269 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
270
271 elseif (tcase.eq.30) then
272 !-------------------------------------------------------------------
273 ! CASE 30: ATHENS - PARTHENON
274 call make_mech_prop_case_030(rho,lambda,mu,gamma,qs,qp, & !outputs
275 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
276 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
277
278 elseif (tcase.eq.31) then
279 !------------------------------------------------------------------
280 ! CASE 31: GRONINGEN
281 call make_mech_prop_case_031(rho,lambda,mu,gamma,qs,qp, & !outputs
282 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
283 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
284
285 elseif (tcase.eq.32) then
286 !------------------------------------------------------------------
287 ! CASE 32: GRONINGEN - LAYERED MODEL
288 call make_mech_prop_case_032(rho,lambda,mu,gamma,qs,qp, & !outputs
289 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
290 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
291
292 elseif (tcase.eq.33) then
293 !------------------------------------------------------------------
294 ! CASE 33: GRONINGEN - 2ND LAYERED MODEL
295 call make_mech_prop_case_033(rho,lambda,mu,gamma,qs,qp, & !outputs
296 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
297 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
298
299 elseif (tcase.eq.35) then
300 ! CASE 35: THESS+MYGD-FINAL
301 call make_mech_prop_case_035(rho,lambda,mu,gamma,qs,qp, & !outputs
302 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
303 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
304
305 elseif (tcase.eq.38) then
306 ! CASE 38: MONTELIMAR
307 call make_mech_prop_case_038(rho,lambda,mu,gamma,qs,qp, & !outputs
308 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
309 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
310
311 elseif (tcase.eq.40) then
312 !-------------------------------------------------------------------
313 ! CASE 40: KUTCH
314 call make_mech_prop_case_040(rho,lambda,mu,gamma,qs,qp, & !outputs
315 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
316 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
317
318 elseif (tcase.eq.45) then
319 !-------------------------------------------------------------------
320 ! CASE 45: MEXICO CITY DG
321 call make_mech_prop_case_045(rho,lambda,mu,gamma,qs,qp, & !outputs
322 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
323 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
324
325
326 elseif (tcase.eq.46) then
327 !-------------------------------------------------------------------
328 ! CASE 46: KUMAMOTO
329 call make_mech_prop_case_046(rho,lambda,mu,gamma,qs,qp, & !outputs
330 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic))
331
332 elseif (tcase.eq.50) then
333 !-------------------------------------------------------------------
334 ! CASE 50: PLANE-WAVE benchmark - MULTI NOT honoring
335 call make_mech_prop_case_050(rho,lambda,mu,gamma,qs,qp, & !outputs
336 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
337 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
338
339 elseif (tcase.eq.60) then
340 !-------------------------------------------------------------------
341 ! CASE 60: JAKARTA
342 call make_mech_prop_case_060(rho,lambda,mu,gamma,qs,qp, & !outputs
343 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
344 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
345
346 elseif (tcase.eq.70) then
347 !-------------------------------------------------------------------
348 ! CASE 70: L'AQUILA - MULTI BASINS
349 call make_mech_prop_case_070(rho,lambda,mu,gamma,qs,qp, & !outputs
350 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
351 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
352
353 elseif (tcase .eq. 91) then
354 !-------------------------------------------------------------------
355 ! CASE 91: IRPINIA
356 call make_mech_prop_case_091(rho,lambda,mu,gamma,qs,qp, & !outputs
357 xs(ic),ys(ic),zs(ic),zs_elev(ic),zs_all(ic), &
358 vs_nodes(ic), thick_nodes(ic), sub_tag_all(ic))
359
360
361 elseif (tcase.eq.98) then
362 !-------------------------------------------------------------------
363 ! CASE 98: TEST honoring (only TOPOGRAPHY)
364 vs = 100 !VS: S velocity in m/s
365 vp = 200 !VS: S velocity in m/s
366 rho = 2000 !RHO: MASS DENSITY in kg/m^3
367 lambda = rho * (vp**2 - 2*vs**2)
368 mu = rho * vs**2
369 gamma = 0.0d0
370
371 elseif (tcase.eq.99) then
372 !-------------------------------------------------------------------
373 ! CASE 99: TEST honoring (TOPOGRAPHY&ALLUVIAL)
374 depth = zs_elev(ic) !D: depth in m
375 if ((depth .ge. 0.0d0) .and. (zs_all(ic) .ge. 0.0d0)) then
376 vs = 1000.d0 !VS: S velocity in m/s
377 vp = 2081.9942d0 !VP: P velocity in m/s
378 rho = 2000.d0 !rho_el: MASS DENSITY in kg/m^3
379 lambda = rho * (vp**2 - 2*vs**2)
380 mu = rho * vs**2
381 qp = 200;
382 qs = 100;
383 gamma = 4.d0*datan(1.d0)/qs;
384 else
385 vs = 2000.d0 !VS: S velocity in m/s
386 vp = 3463.9976 !VP: P velocity in m/s
387 rho = 2500.d0 !rho_el: MASS DENSITY in kg/m^3
388 lambda = rho * (vp**2 - 2*vs**2)
389 mu = rho * vs**2
390 qp = 200;
391 qs = 100;
392 gamma = 4.d0*datan(1.d0)/qs;
393 endif
394
395 elseif (tcase.eq.100) then
396 !-------------------------------------------------------------------
397 ! CASE 99: TEST plane wave (TOPOGRAPHY&ALLUVIAL)
398 depth = zs_elev(ic) !D: depth in m
399 if ((depth .ge. 0.0d0) .and. (zs_all(ic) .ge. 0.0d0)) then
400 vs = 300.d0 !VS: S velocity in m/s
401 vp = 600.d0 !VP: P velocity in m/s
402 rho = 1800.d0 !rho_el: MASS DENSITY in kg/m^3
403 lambda = rho * (vp**2 - 2*vs**2)
404 mu = rho * vs**2
405 qp = 60;
406 qs = 30;
407 gamma = 4.d0*datan(1.d0)/qs;
408 else
409 vs = 2000.d0 !VS: S velocity in m/s
410 vp = 4000.d0 !VP: P velocity in m/s
411 rho = 2200.d0 !rho_el: MASS DENSITY in kg/m^3
412 lambda = rho * (vp**2 - 2*vs**2)
413 mu = rho * vs**2
414 qp = 400;
415 qs = 200;
416 gamma = 4.d0*datan(1.d0)/qs;
417 endif
418 endif
419
420 if (check_case .eq. 1) &
421 write(1000+mpi_id,*) xs(ic),ys(ic),zs(ic), &
422 dsqrt(mu/rho), &
423 dsqrt((lambda + 2.d0*mu)/rho), &
424 rho, lambda, mu, &
425 qp, qs, gamma, zs_elev(ic), zs_all(ic)
426
427
428 rho_el(p,q,r) = rho
429 lambda_el(p,q,r) = lambda
430 mu_el(p,q,r) = mu
431 gamma_el(p,q,r) = gamma
432 qp_h(p,q,r) = qp !New
433 qs_h(p,q,r) = qs !New
434
435 enddo
436 enddo
437 enddo
438
439 if (check_case .eq. 1) close (1000+mpi_id)
440
441 if (damping_type .eq. 2) then
442 qs = 0; qp = 0;
443 qs_all = 0.d0; qp_all=0.d0; !q*_all instead of v*_all
444 do r = 1,nn
445 do q = 1,nn
446 do p = 1,nn
447 ! Changed vs_all = vs_all + dsqrt(mu_el(p,q,r)/rho_el(p,q,r))
448 ! Changed vp_all = vp_all + dsqrt( (lambda_el(p,q,r) + 2.d0*mu_el(p,q,r))/rho_el(p,q,r) );
449 qs_all = qs_all + qs_h(p,q,r)
450 qp_all = qp_all + qp_h(p,q,r)
451 enddo
452 enddo
453 enddo
454
455 ! Changed qs = 0.1d0*vs_all/nn**3;
456 ! Changed qp = 0.1d0*vp_all/nn**3;
457 qs = qs_all/nn**3;
458 qp = qp_all/nn**3;
459
460 endif
461
462
463 return
464
465 end subroutine make_eltensor_for_cases
466
make_eltensor_for_cases
subroutine make_eltensor_for_cases(tcase, vcase, nn, rho_el, lambda_el, mu_el, gamma_el, nn_loc, zs_elev, zs_all, vs_nodes, thic
Assignes material properties node by node.
Definition MAKE_ELTENSOR_FOR_CASES.f90:53
make_mech_prop_case_001
subroutine make_mech_prop_case_001(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_001.f90:26
make_mech_prop_case_002
subroutine make_mech_prop_case_002(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_002.f90:26
make_mech_prop_case_003
subroutine make_mech_prop_case_003(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_003.f90:26
make_mech_prop_case_004
subroutine make_mech_prop_case_004(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_004.f90:26
make_mech_prop_case_005
subroutine make_mech_prop_case_005(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_005.f90:26
make_mech_prop_case_006
subroutine make_mech_prop_case_006(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_006.f90:26
make_mech_prop_case_007
subroutine make_mech_prop_case_007(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_007.f90:26
make_mech_prop_case_008
subroutine make_mech_prop_case_008(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_008.f90:26
make_mech_prop_case_010
subroutine make_mech_prop_case_010(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_010.f90:26
make_mech_prop_case_011
subroutine make_mech_prop_case_011(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_011.f90:26
make_mech_prop_case_012
subroutine make_mech_prop_case_012(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all, f_distance)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_012.f90:26
make_mech_prop_case_013
subroutine make_mech_prop_case_013(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_013.f90:26
make_mech_prop_case_014
subroutine make_mech_prop_case_014(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_014.f90:26
make_mech_prop_case_015
subroutine make_mech_prop_case_015(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_015.f90:26
make_mech_prop_case_016
subroutine make_mech_prop_case_016(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_016.f90:26
make_mech_prop_case_018
subroutine make_mech_prop_case_018(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_018.f90:26
make_mech_prop_case_019
subroutine make_mech_prop_case_019(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_019.f90:26
make_mech_prop_case_020
subroutine make_mech_prop_case_020(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_020.f90:27
make_mech_prop_case_021
subroutine make_mech_prop_case_021(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_021.f90:27
make_mech_prop_case_022
subroutine make_mech_prop_case_022(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_022.f90:27
make_mech_prop_case_027
subroutine make_mech_prop_case_027(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_027.f90:26
make_mech_prop_case_028
subroutine make_mech_prop_case_028(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_028.f90:27
make_mech_prop_case_029
subroutine make_mech_prop_case_029(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_029.f90:26
make_mech_prop_case_030
subroutine make_mech_prop_case_030(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_030.f90:27
make_mech_prop_case_031
subroutine make_mech_prop_case_031(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_031.f90:27
make_mech_prop_case_032
subroutine make_mech_prop_case_032(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_032.f90:27
make_mech_prop_case_033
subroutine make_mech_prop_case_033(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_033.f90:27
make_mech_prop_case_035
subroutine make_mech_prop_case_035(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_035.f90:26
make_mech_prop_case_038
subroutine make_mech_prop_case_038(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_038.f90:27
make_mech_prop_case_040
subroutine make_mech_prop_case_040(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_040.f90:27
make_mech_prop_case_045
subroutine make_mech_prop_case_045(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_045.f90:27
make_mech_prop_case_046
subroutine make_mech_prop_case_046(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_046.f90:27
make_mech_prop_case_050
subroutine make_mech_prop_case_050(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_050.f90:27
make_mech_prop_case_060
subroutine make_mech_prop_case_060(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_060.f90:26
make_mech_prop_case_070
subroutine make_mech_prop_case_070(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_070.f90:26
make_mech_prop_case_091
subroutine make_mech_prop_case_091(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_091.f90:26