#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Sun Oct 29 12:44:09 2017

@author: thomas

This file contains data sets for layered media models.

List of list structure

	[[layer1],	Where Layers [C11, C13, C33, C44, density, thickness]	
	[layer2],
	  ....
	[Half Space]]	Half space same as layers except thickness = inf (not needed)	
"""



""" Isotropic model from Dalton et al. Forward Modelling """

Iso1 = [[1.98*10**10, .88*10**10, 2200, 100],
       [1.98*10**10, .88*10**10, 2200, 100],
       [1.98*10**10, .88*10**10, 2200, 100],
       [1.98*10**10, .88*10**10, 2200, 100],
       [1.98*10**10, .88*10**10, 2200, 100],
       [109850000000.0, 26650000000.0, 109850000000.0, 41600000000.0, 2600]]

Iso2 = [[23232000000.0, 14344000000.0, 23232000000.0, 4444000000.0, 2200, 5],
        [45144000000.0, 25564000000.0, 45144000000.0, 9790000000.0, 2200, 5],
        [68508000000.0, 55792000000.0, 68508000000.0, 6358000000.0, 2200, 5],
        [32604000000.0, 21076000000.0, 32604000000.0, 5764000000.0, 2200, 5],
        [70730000000.0, 57882000000.0, 70730000000.0, 6424000000.0, 2200, 5],
        [35200000000.0, 23936000000.0, 35200000000.0, 5632000000.0, 2200, 5],
        [36079999999.99999, 8139999999.999992, 36079999999.99999, 13970000000.0, 2200, 5],
        [39732000000.0, 20680000000.0, 39732000000.0, 9526000000.0, 2200, 5],
        [69234000000.0, 33990000000.0, 69234000000.0, 17622000000.0, 2200, 5],
        [38082000000.0, 21538000000.0, 38082000000.0, 8272000000.0, 2200, 5],
        [109850000000.0, 26650000000.0, 109850000000.0, 41600000000.0, 2600]]

Iso3 = [[23232000000.0, 14344000000.0, 23232000000.0, 4444000000.0, 2200, 10],
        [45144000000.0, 25564000000.0, 45144000000.0, 9790000000.0, 2200, 10],
        [68508000000.0, 55792000000.0, 68508000000.0, 6358000000.0, 2200, 10],
        [32604000000.0, 21076000000.0, 32604000000.0, 5764000000.0, 2200, 10],
        [70730000000.0, 57882000000.0, 70730000000.0, 6424000000.0, 2200, 10],
        [35200000000.0, 23936000000.0, 35200000000.0, 5632000000.0, 2200, 10],
        [36079999999.99999, 8139999999.999992, 36079999999.99999, 13970000000.0, 2200, 10],
        [39732000000.0, 20680000000.0, 39732000000.0, 9526000000.0, 2200, 10],
        [69234000000.0, 33990000000.0, 69234000000.0, 17622000000.0, 2200, 10],
        [38082000000.0, 21538000000.0, 38082000000.0, 8272000000.0, 2200, 10],
        [109850000000.0, 26650000000.0, 109850000000.0, 41600000000.0, 2600]]

""" Transverseley Isotropic model from David new paper """

Ti1 = [[5.89448*10**10,7.66207*10**9,3.34552*10**10,1.408*10**10,2200,100],
      [5.89448*10**10,7.66207*10**9,3.34552*10**10,1.408*10**10,2200,100],
      [5.89448*10**10,7.66207*10**9,3.34552*10**10,1.408*10**10,2200,100],
      [5.89448*10**10,7.66207*10**9,3.34552*10**10,1.408*10**10,2200,100],
      [5.89448*10**10,7.66207*10**9,3.34552*10**10,1.408*10**10,2200,100],
      [109850000000.0, 26650000000.0, 109850000000.0, 41600000000.0, 2600]]

""" Alternating TI model from David new paper """

Ti_Alt = [[30201160000.0, 12658800000.0, 36883000000.0, 12207360000.0, 2200, 5],
         [17741020000.0, 5409360000.0, 15579740000.0, 4097500000.0, 2200, 5],
         [30201160000.0, 12658800000.0, 36883000000.0, 12207360000.0, 2200, 5],
         [17741020000.0, 5409360000.0, 15579740000.0, 4097500000.0, 2200, 5],
         [30201160000.0, 12658800000.0, 36883000000.0, 12207360000.0, 2200, 5],
         [17741020000.0, 5409360000.0, 15579740000.0, 4097500000.0, 2200, 5],
         [30201160000.0, 12658800000.0, 36883000000.0, 12207360000.0, 2200, 5],
         [17741020000.0, 5409360000.0, 15579740000.0, 4097500000.0, 2200, 5],
         [30201160000.0, 12658800000.0, 36883000000.0, 12207360000.0, 2200, 5],
         [17741020000.0, 5409360000.0, 15579740000.0, 4097500000.0, 2200, 5],
         [109850000000.0, 26650000000.0, 109850000000.0, 41600000000.0, 2600]]

Ti_Alt_s = [[17741020000.0, 5409360000.0, 15579740000.0, 4097500000.0, 2200, 5],
         [30201160000.0, 12658800000.0, 36883000000.0, 12207360000.0, 2200, 5],
         [17741020000.0, 5409360000.0, 15579740000.0, 4097500000.0, 2200, 5],
         [30201160000.0, 12658800000.0, 36883000000.0, 12207360000.0, 2200, 5],
         [17741020000.0, 5409360000.0, 15579740000.0, 4097500000.0, 2200, 5],
         [30201160000.0, 12658800000.0, 36883000000.0, 12207360000.0, 2200, 5],
         [17741020000.0, 5409360000.0, 15579740000.0, 4097500000.0, 2200, 5],
         [30201160000.0, 12658800000.0, 36883000000.0, 12207360000.0, 2200, 5],
         [17741020000.0, 5409360000.0, 15579740000.0, 4097500000.0, 2200, 5],
         [30201160000.0, 12658800000.0, 36883000000.0, 12207360000.0, 2200, 5],
         [109850000000.0, 26650000000.0, 109850000000.0, 41600000000.0, 2600]]

Ti_Alt_s50 = [[17741020000.0, 5409360000.0, 15579740000.0, 4097500000.0, 2200, 50],
         [30201160000.0, 12658800000.0, 36883000000.0, 12207360000.0, 2200, 50],
         [17741020000.0, 5409360000.0, 15579740000.0, 4097500000.0, 2200, 50],
         [30201160000.0, 12658800000.0, 36883000000.0, 12207360000.0, 2200, 50],
         [17741020000.0, 5409360000.0, 15579740000.0, 4097500000.0, 2200, 50],
         [30201160000.0, 12658800000.0, 36883000000.0, 12207360000.0, 2200, 50],
         [17741020000.0, 5409360000.0, 15579740000.0, 4097500000.0, 2200, 50],
         [30201160000.0, 12658800000.0, 36883000000.0, 12207360000.0, 2200, 50],
         [17741020000.0, 5409360000.0, 15579740000.0, 4097500000.0, 2200, 50],
         [30201160000.0, 12658800000.0, 36883000000.0, 12207360000.0, 2200, 50],
         [109850000000.0, 26650000000.0, 109850000000.0, 41600000000.0, 2600]]

""" Transverseley Isotropic model calculated from Anderson/Harkrider using
    Anderson equation for c13and Red Book 
"""

Ti_And = [[2899232000, 2628832000.00000, 2899232000.00000, 135200000, 2000, 5],
         [4428288000, 3621884000.00000, 4428288000.00000, 403202000, 2000, 5],
         [3880898000, 1737602000.00000, 3880898000.00000, 1071648000, 2000, 5],
         [5803281000, 1303281000.00000, 5803281000.00000, 2250000000, 2250, 5], 
         [6502500000, 837233003.311483, 5891264999.99999, 2653641000, 2250, 5],
         [6502500000, 1060066724.65262, 5891264999.99999, 2542430250, 2250, 5], 
         [7249556250, 164204573.338301, 6568097962.49999, 2542430250, 2250, 5],
         [7249556250, 1233965981.06737, 6568097962.49999, 2009306250, 2250, 5],
         [7834401000, 2044076178.95201, 4920003828.00000, 2250000000, 2250, 5],
         [7834401000, 1441600669.38677, 4920003828.00000, 2542430250, 2250, 5],
         [9000000000, -4269893000.00504, 5652000000.00000, 4517750250, 2250, 5],
         [11329956000, -1978370008.75475, 7115212368.00000, 4517750250, 2250, 5],
         [37966522500, 20687597381.9136, 36068196375.0000, 6448090000, 2500, 5],
         [73116160000, 51016884341.6222, 69460351999.9999, 6831022500, 2500, 5],
         [73116160000, 54835023309.5214, 62148736000.0000, 6642250000, 2500, 5],
         [73116160000, 54835023309.5214, 62148736000.0000, 6642250000, 2500, 5],
         [80997160000, 49942258638.5152, 68847586000.0000, 10312402500, 2500, 5],
         [78316022500, 52559389945.8997, 66568619125.0000, 7841102500, 2500, 5],
         [75707522500, 51810313140.8765, 68136770250.0000, 8482410000, 2500, 5],
         [93544222500, 70123451119.1532, 84189800250.0000, 7430440000, 2500, 5],
         [101075585000,71022091873.1579, 90968026500.0000,10400000000,2600]]
    
    
""" Anderson Ti model only layer 6 over half space depth adjusted to total 100 """

Ti_Andy_L6 = [[6502500000, 1047079799.58546, 5891264999.99999, 2542430250, 2250, 20],
             [6502500000, 1047079799.58546, 5891264999.99999, 2542430250, 2250, 20],
             [6502500000, 1047079799.58546, 5891264999.99999, 2542430250, 2250, 20],
             [6502500000, 1047079799.58546, 5891264999.99999, 2542430250, 2250, 20],
             [6502500000, 1047079799.58546, 5891264999.99999, 2542430250, 2250, 20],
             [101075585000, 70865099275.7238, 90968026500.0, 10400000000, 2600]]

""" model from keetel and chen papers """

chen = [[100800000000., 32200000000., 100800000000., 34300000000., 2800., 18000.],
        [115101000000., 37830500000., 115101000000., 38635250000., 2900., 6000.],
        [139159000000., 44857000000., 139159000000., 47151000000., 3100., 6000.],
        [221892000000., 76098000000., 221892000000., 72897000000., 3300., ]]

haskell = [[101788920000.0, 39731580000.0, 101788920000.0, 31028670000.0, 2700, 13600.0],
           [81675000000.0, 27068040000.0, 81675000000.0, 27303480000.0, 2700, 11850.0],
           [209668800000.0, 79933800000.0, 209668800000.0, 64867500000.0, 3000]]
