Changeset 26907
- Timestamp:
- 02/25/22 00:34:01 (3 years ago)
- Location:
- issm/trunk-jpl/test/NightlyRun
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test257.py
r26810 r26907 43 43 md.smb = SMBautoregression() 44 44 md.smb.num_basins = 3 # number of basins 45 md.smb.basin_id = idbasin - 1 # prescribe basin ID number to elements; # NOTE: offset because of 1-based vertex indexing45 md.smb.basin_id = idbasin # prescribe basin ID number to elements; 46 46 md.smb.beta0 = np.array([[0.5, 1.2, 1.5]]) # intercept values of SMB in basins [m ice eq./yr] 47 47 md.smb.beta1 = np.array([[0.0, 0.01, -0.01]]) # trend values of SMB in basins [m ice eq./yr^2] -
issm/trunk-jpl/test/NightlyRun/test542.py
r26638 r26907 1 1 #Test Name: PigTranFrontalforcingsrignot 2 2 import numpy as np 3 4 from calvingvonmises import * 5 from frontalforcingsrignot import * 3 from calvingvonmises import calvingvonmises 4 from frontalforcingsrignot import frontalforcingsrignot 6 5 from socket import gethostname 7 6 from model import * 8 from parameterize import *9 from setflowequation import *10 from setmask import *11 from solve import *12 from triangle import *7 from parameterize import parameterize 8 from setflowequation import setflowequation 9 from setmask import setmask 10 from solve import solve 11 from triangle import triangle 13 12 14 13 … … 35 34 md.frontalforcings = frontalforcingsrignot() 36 35 md.frontalforcings.num_basins = 2 37 md.frontalforcings.basin_id = idbasin - 1 # NOTE: offset because of 1-based vertex indexing36 md.frontalforcings.basin_id = idbasin 38 37 md.frontalforcings.subglacial_discharge = 0.1 * np.ones((md.mesh.numberofvertices,)) 39 38 md.frontalforcings.thermalforcing = 0.5 * np.ones((md.mesh.numberofvertices,)) 40 39 for elem in range(md.mesh.numberofelements): 41 40 if idbasin[elem] == 2: 42 md.frontalforcings.thermalforcing[md.mesh.elements[elem, 0:3] - 1] = 1.5 # 41 md.frontalforcings.thermalforcing[md.mesh.elements[elem, 0:3] - 1] = 1.5 #NOTE: offset because of 1-based vertex indexing 43 42 44 43 md.transient.ismovingfront = 1 … … 51 50 52 51 # Fields and tolerances to track changes 53 field_names = [ 54 'Vx1', 'Vy1', 'Vel1', 'Pressure1', 'Bed1', 'Surface1', 'Thickness1', 'MaskIceLevelset1', 'CalvingMetlingRate1', 55 'Vx2', 'Vy2', 'Vel2', 'Pressure2', 'Bed2', 'Surface2', 'Thickness2', 'MaskIceLevelset2', 'CalvingMetlingRate2', 56 'Vx10', 'Vy10', 'Vel10', 'Pressure10', 'Bed10', 'Surface10', 'Thickness10', 'MaskIceLevelset10', 'CalvingMetlingRate10' 57 ] 52 field_names = ['Vx1', 'Vy1', 'Vel1', 'Pressure1', 'Bed1', 'Surface1', 'Thickness1', 'MaskIceLevelset1', 'CalvingMetlingRate1', 53 'Vx2', 'Vy2', 'Vel2', 'Pressure2', 'Bed2', 'Surface2', 'Thickness2', 'MaskIceLevelset2', 'CalvingMetlingRate2', 54 'Vx10', 'Vy10', 'Vel10', 'Pressure10', 'Bed10', 'Surface10', 'Thickness10', 'MaskIceLevelset10', 'CalvingMetlingRate10'] 58 55 field_tolerances = [ 59 56 1e-11, 2e-11, 2e-11, 1e-12, 2e-11, 6e-12, 9e-12, 1e-12, 1e-9, 60 57 2e-11, 1e-11, 1e-11, 9e-12, 2e-1, 2e-11, 2e-11, 1e-11, 1e-9, 61 2e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-9, 1e-9 62 ] 58 2e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-9, 1e-9] 63 59 field_values = [ 64 60 md.results.TransientSolution[0].Vx, -
issm/trunk-jpl/test/NightlyRun/test543.py
r26674 r26907 1 1 #Test Name: PigTranStochasticforcing 2 2 import numpy as np 3 4 from frontalforcingsrignotautoregression import * 3 from frontalforcingsrignotautoregression import frontalforcingsrignotautoregression 5 4 from socket import gethostname 6 5 from model import * 7 from parameterize import *8 from setflowequation import *9 from setmask import *10 from solve import *11 from triangle import *6 from parameterize import parameterize 7 from setflowequation import setflowequation 8 from setmask import setmask 9 from solve import solve 10 from triangle import triangle 12 11 13 12 … … 41 40 idb_df[ii] = 2 42 41 #Dimensionalities 43 nb_tf 44 nb_clv 42 nb_tf = 2 43 nb_clv = 2 45 44 nb_flmlt = 2 46 45 47 46 # Calving parameters 48 47 md.mask.ice_levelset = 1e4 * (md.mask.ice_levelset + 0.5) 49 md.calving.calvingrate = 0.3 *np.ones((md.mesh.numberofvertices,))48 md.calving.calvingrate = 0.3 * np.ones((md.mesh.numberofvertices,)) 50 49 md.levelset.spclevelset = np.full((md.mesh.numberofvertices,), np.nan) 51 50 md.levelset.migration_max = 10.0 … … 53 52 md.frontalforcings = frontalforcingsrignotautoregression() 54 53 md.frontalforcings.num_basins = nb_tf 55 md.frontalforcings.basin_id = idb_tf - 1 # NOTE: offset because of 1-based vertex indexing54 md.frontalforcings.basin_id = idb_tf 56 55 md.frontalforcings.subglacial_discharge = 0.1 * np.ones((md.mesh.numberofvertices,)) 57 56 md.frontalforcings.beta0 = np.array([[0.05, 0.01]]) # intercept values of TF in basins [C] … … 62 61 md.frontalforcings.phi = np.array([[0.1, -0.1, 0.01, -0.01], [0.2, -0.2, 0.1, 0.0]]) # autoregressive parameters 63 62 #Floating Ice Melt parameters 64 md.basalforcings.floatingice_melting_rate = 0.1 *np.ones((md.mesh.numberofvertices,))63 md.basalforcings.floatingice_melting_rate = 0.1 * np.ones((md.mesh.numberofvertices,)) 65 64 66 65 #Covariance matrix 67 covtf = 1e-4*np.identity(nb_tf)68 covclv = 1e-1*np.identity(nb_clv)69 covclv[0, 0] = 1/10*covclv[0,0]70 covflmlt = 0.05*np.identity(nb_flmlt)66 covtf = 1e-4 * np.identity(nb_tf) 67 covclv = 1e-1 * np.identity(nb_clv) 68 covclv[0, 0] = 1 / 10 * covclv[0, 0] 69 covflmlt = 0.05 * np.identity(nb_flmlt) 71 70 #covglob = np.zeros([6,6]) 72 71 #covglob[0:2,0:2] = covtf … … 75 74 76 75 #Hard-coding covariance matrix because python is complaining 77 covglob = np.array([[1e-4,0.,0.,0.,0.,0.],[0.,1e-4,0.,0.,0.,0.],[0.,0.,1e-2,0.,0.,0.],[0.,0.,0.,1e-1,0.,0.],[0.,0.,0.,0.,0.05,0.],[0.,0.,0.,0.,0.,0.05]]) 76 covglob = np.array([[1e-4, 0., 0., 0., 0., 0.], 77 [0., 1e-4, 0., 0., 0., 0.], 78 [0., 0., 1e-2, 0., 0., 0.], 79 [0., 0., 0., 1e-1, 0., 0.], 80 [0., 0., 0., 0., 0.05, 0.], 81 [0., 0., 0., 0., 0., 0.05]]) 78 82 #testchol = np.linalg.cholesky(covglob) 79 83 #print(testchol) … … 81 85 # Stochastic forcing 82 86 md.stochasticforcing.isstochasticforcing = 1 83 md.stochasticforcing.fields = ['FrontalForcingsRignotAutoregression', 'DefaultCalving','FloatingMeltRate']87 md.stochasticforcing.fields = ['FrontalForcingsRignotAutoregression', 'DefaultCalving', 'FloatingMeltRate'] 84 88 md.stochasticforcing.defaultdimension = 2 85 md.stochasticforcing.default_id = idb_df-1 #NOTE: offset because of 1-based vertex indexing86 md.stochasticforcing.covariance = covglob# global covariance among- and between-fields87 md.stochasticforcing.randomflag = 0# determines true/false randomness89 md.stochasticforcing.default_id = idb_df 90 md.stochasticforcing.covariance = covglob # global covariance among- and between-fields 91 md.stochasticforcing.randomflag = 0 # determines true/false randomness 88 92 89 md.transient.ismovingfront 93 md.transient.ismovingfront = 1 90 94 md.transient.isgroundingline = 1 91 95 92 md.transient.requested_outputs = ['default', 'CalvingCalvingrate', 'CalvingMeltingrate','BasalforcingsFloatingiceMeltingRate']96 md.transient.requested_outputs = ['default', 'CalvingCalvingrate', 'CalvingMeltingrate', 'BasalforcingsFloatingiceMeltingRate'] 93 97 md.cluster = generic('name', gethostname(), 'np', 2) 94 98 md = solve(md, 'Transient') … … 96 100 # Fields and tolerances to track changes 97 101 field_names = [ 98 'Vx1' ,'Vy1' ,'Vel1' ,'Thickness1' ,'MaskIceLevelset1' ,'CalvingCalvingrate1' ,'CalvingMeltingrate1' ,'BasalforcingsFloatingiceMeltingRate1', 99 'Vx2' ,'Vy2' ,'Vel2' ,'Thickness2' ,'MaskIceLevelset2' ,'CalvingCalvingrate1' ,'CalvingMeltingrate2' ,'BasalforcingsFloatingiceMeltingRate2', 100 'Vx10','Vy10','Vel10','Thickness10','MaskIceLevelset10','CalvingCalvingrate10','CalvingMeltingrate10','BasalforcingsFloatingiceMeltingRate10' 101 ] 102 'Vx1', 'Vy1', 'Vel1', 'Thickness1', 'MaskIceLevelset1', 'CalvingCalvingrate1', 'CalvingMeltingrate1', 'BasalforcingsFloatingiceMeltingRate1', 103 'Vx2', 'Vy2', 'Vel2', 'Thickness2', 'MaskIceLevelset2', 'CalvingCalvingrate1', 'CalvingMeltingrate2', 'BasalforcingsFloatingiceMeltingRate2', 104 'Vx10', 'Vy10', 'Vel10', 'Thickness10', 'MaskIceLevelset10', 'CalvingCalvingrate10', 'CalvingMeltingrate10', 'BasalforcingsFloatingiceMeltingRate10'] 102 105 103 106 field_tolerances = [ 104 1e-11,2e-11,2e-11,1e-11,1e-9,1e-10,1e-10,1e-10, 105 2e-11,1e-11,1e-11,9e-11,2e-9,1e-10,1e-10,1e-10, 106 2e-10,1e-10,1e-10,1e-10,5e-9,1e-10,1e-10,1e-10 107 ] 107 1e-11, 2e-11, 2e-11, 1e-11, 1e-9, 1e-10, 1e-10, 1e-10, 108 2e-11, 1e-11, 1e-11, 9e-11, 2e-9, 1e-10, 1e-10, 1e-10, 109 2e-10, 1e-10, 1e-10, 1e-10, 5e-9, 1e-10, 1e-10, 1e-10] 108 110 field_values = [ 109 111 md.results.TransientSolution[0].Vx, … … 130 132 md.results.TransientSolution[9].CalvingCalvingrate, 131 133 md.results.TransientSolution[9].CalvingMeltingrate, 132 md.results.TransientSolution[9].BasalforcingsFloatingiceMeltingRate 133 ] 134 md.results.TransientSolution[9].BasalforcingsFloatingiceMeltingRate] -
issm/trunk-jpl/test/NightlyRun/test544.py
r26837 r26907 1 1 #Test Name: PigTranAutoregressionAndStochasticforcings 2 2 import numpy as np 3 4 from autoregressionlinearbasalforcings import * 5 from SMBautoregression import * 6 from stochasticforcing import * 3 from autoregressionlinearbasalforcings import autoregressionlinearbasalforcings 4 from SMBautoregression import SMBautoregression 5 from stochasticforcing import stochasticforcing 7 6 from socket import gethostname 8 7 from model import * 9 from parameterize import *10 from setflowequation import *11 from setmask import *12 from solve import *13 from triangle import *8 from parameterize import parameterize 9 from setflowequation import setflowequation 10 from setmask import setmask 11 from solve import solve 12 from triangle import triangle 14 13 15 14 … … 37 36 md.smb = SMBautoregression() 38 37 md.smb.num_basins = nb_bas # number of basins 39 md.smb.basin_id = idb - 1 # prescribe basin ID number to elements; # NOTE: offset because of 1-based vertex indexing38 md.smb.basin_id = idb # prescribe basin ID number to elements; 40 39 md.smb.beta0 = np.array([[0.5, 1.2]]) # intercept values of SMB in basins [m ice eq./yr] 41 40 md.smb.beta1 = np.array([[0.0, 0.01]]) # trend values of SMB in basins [m ice eq./yr^2] … … 49 48 md.calving.calvingrate = 0.1*np.ones((md.mesh.numberofvertices,)) 50 49 md.levelset.spclevelset = np.full((md.mesh.numberofvertices,), np.nan) 51 md.levelset.migration_max = 10.0 50 md.levelset.migration_max = 10.0 52 51 md.frontalforcings.meltingrate = np.zeros((md.mesh.numberofvertices,)) 53 52 54 53 #Basal forcing implementation 55 54 md.basalforcings = autoregressionlinearbasalforcings() 56 md.basalforcings.num_basins 57 md.basalforcings.basin_id = idb - 1 # NOTE: offset because of 1-based vertex indexing58 md.basalforcings.beta0 59 md.basalforcings.beta1 55 md.basalforcings.num_basins = nb_bas 56 md.basalforcings.basin_id = idb 57 md.basalforcings.beta0 = np.array([[1.0, 2.50]]) # intercept values of DeepwaterMelt in basins [m/yr] 58 md.basalforcings.beta1 = np.array([[0.2, 0.01]]) # trend values of DeepwaterMelt in basins [m/yr^2] 60 59 md.basalforcings.ar_initialtime = md.timestepping.start_time # initial time in the AR model parameterization [yr] 61 md.basalforcings.ar_order 62 md.basalforcings.ar_timestep 63 md.basalforcings.phi 64 md.basalforcings.deepwater_elevation 65 md.basalforcings.upperwater_elevation 66 md.basalforcings.upperwater_melting_rate 60 md.basalforcings.ar_order = 1 61 md.basalforcings.ar_timestep = 1.0 # timestep of the autoregressive model [yr] 62 md.basalforcings.phi = np.array([[0.0], [0.1]]) # autoregressive parameters 63 md.basalforcings.deepwater_elevation = np.array([[-1000, -1520]]) 64 md.basalforcings.upperwater_elevation = np.array([[0, -50]]) 65 md.basalforcings.upperwater_melting_rate = np.array([[0,0]]) 67 66 md.basalforcings.groundedice_melting_rate = np.zeros((md.mesh.numberofvertices,)) 68 67 69 68 #Covariance matrix (hard-coding) 70 covglob = np.array([[3.0,0.,0.,0.,0.,0.],[0.,3.0,0.,0.,0.,0.],[0.,0.,0.01,0.,0.,0.],[0.,0.,0.,0.1,0.,0.],[0.,0.,0.,0.,400,0.],[0.,0.,0.,0.,0.,400]]) 69 covglob = np.array([[3.0, 0., 0., 0., 0., 0.], 70 [0., 3.0, 0., 0., 0., 0.], 71 [0., 0., 0.01, 0., 0., 0.], 72 [0., 0., 0., 0.1, 0., 0.], 73 [0., 0., 0., 0., 400, 0.], 74 [0., 0., 0., 0., 0., 400]]) 71 75 72 76 #Stochastic forcing 73 77 md.stochasticforcing.isstochasticforcing = 1 74 md.stochasticforcing.fields = ['SMBautoregression', 'DefaultCalving','BasalforcingsDeepwaterMeltingRateAutoregression']78 md.stochasticforcing.fields = ['SMBautoregression', 'DefaultCalving', 'BasalforcingsDeepwaterMeltingRateAutoregression'] 75 79 md.stochasticforcing.defaultdimension = 2 76 md.stochasticforcing.default_id = idb-1 #NOTE: offset because of 1-based vertex indexing77 md.stochasticforcing.covariance = covglob# global covariance among- and between-fields78 md.stochasticforcing.randomflag = 0# determines true/false randomness80 md.stochasticforcing.default_id = idb 81 md.stochasticforcing.covariance = covglob # global covariance among- and between-fields 82 md.stochasticforcing.randomflag = 0 # determines true/false randomness 79 83 80 md.transient.ismovingfront 81 md.transient.requested_outputs = ['default', 'SmbMassBalance','BasalforcingsFloatingiceMeltingRate','BasalforcingsSpatialDeepwaterMeltingRate']82 md.transient.isstressbalance 83 md.transient.ismasstransport 84 md.transient.issmb 85 md.transient.isthermal 86 md.transient.isgroundingline 84 md.transient.ismovingfront = 1 85 md.transient.requested_outputs = ['default', 'SmbMassBalance', 'BasalforcingsFloatingiceMeltingRate', 'BasalforcingsSpatialDeepwaterMeltingRate'] 86 md.transient.isstressbalance = 1 87 md.transient.ismasstransport = 1 88 md.transient.issmb = 1 89 md.transient.isthermal = 0 90 md.transient.isgroundingline = 1 87 91 88 92 md.cluster = generic('name', gethostname(), 'np', 2) … … 91 95 # Fields and tolerances to track changes 92 96 field_names = [ 93 'Vx1','Vy1','Vel1','Thickness1','SmbMassBalance1','BasalforcingsFloatingiceMeltingRate1','BasalforcingsSpatialDeepwaterMeltingRate1', 94 'Vx5','Vy5','Vel5','Thickness5','SmbMassBalance5','BasalforcingsFloatingiceMeltingRate5','BasalforcingsSpatialDeepwaterMeltingRate5', 95 'Vx10','Vy10','Vel10','Thickness10','SmbMassBalance10','BasalforcingsFloatingiceMeltingRate10','BasalforcingsSpatialDeepwaterMeltingRate10' 96 ] 97 'Vx1', 'Vy1', 'Vel1', 'Thickness1', 'SmbMassBalance1', 'BasalforcingsFloatingiceMeltingRate1', 'BasalforcingsSpatialDeepwaterMeltingRate1', 98 'Vx5', 'Vy5', 'Vel5', 'Thickness5', 'SmbMassBalance5', 'BasalforcingsFloatingiceMeltingRate5', 'BasalforcingsSpatialDeepwaterMeltingRate5', 99 'Vx10', 'Vy10', 'Vel10', 'Thickness10', 'SmbMassBalance10', 'BasalforcingsFloatingiceMeltingRate10', 'BasalforcingsSpatialDeepwaterMeltingRate10'] 97 100 98 101 field_tolerances = [ 99 1e-11,1e-11,2e-11,1e-11,1e10,1e-9,1e-10, 100 1e-11,1e-11,2e-11,9e-11,1e10,1e-9,1e-10, 101 2e-10,2e-10,2e-10,1e-10,1e10,1e-9,1e-10 102 ] 102 1e-11, 1e-11, 2e-11, 1e-11, 1e10, 1e-9, 1e-10, 103 1e-11, 1e-11, 2e-11, 9e-11, 1e10, 1e-9, 1e-10, 104 2e-10, 2e-10, 2e-10, 1e-10, 1e10, 1e-9, 1e-10] 103 105 field_values = [ 104 106 md.results.TransientSolution[0].Vx, … … 122 124 md.results.TransientSolution[9].SmbMassBalance, 123 125 md.results.TransientSolution[9].BasalforcingsFloatingiceMeltingRate, 124 md.results.TransientSolution[9].BasalforcingsSpatialDeepwaterMeltingRate 125 ] 126 md.results.TransientSolution[9].BasalforcingsSpatialDeepwaterMeltingRate] -
issm/trunk-jpl/test/NightlyRun/test621.py
r26679 r26907 12 12 13 13 md = triangle(model(), '../Exp/79North.exp', 6000) 14 md = setmask(md, '../Exp/79NorthShelf.exp', '')14 md = setmask(md, '../Exp/79NorthShelf.exp', '') 15 15 md = parameterize(md, '../Par/79North.py') 16 16 md = setflowequation(md, 'SSA', 'all') 17 17 18 18 #Default friction 19 md.friction 20 md.friction.coefficient = 30 *np.ones(md.mesh.numberofvertices)21 md.friction.p = 1*np.ones((md.mesh.numberofelements))22 md.friction.q = 1*np.ones((md.mesh.numberofelements))19 md.friction = friction() 20 md.friction.coefficient = 30 * np.ones(md.mesh.numberofvertices) 21 md.friction.p = np.ones((md.mesh.numberofelements)) 22 md.friction.q = np.ones((md.mesh.numberofelements)) 23 23 24 24 # Basin separation default 25 idb_df = np.zeros((md.mesh.numberofelements ,))26 iid1 = np.where(md.mesh.y <=-1.08e6)[0]25 idb_df = np.zeros((md.mesh.numberofelements)) 26 iid1 = np.where(md.mesh.y <= -1.08e6)[0] 27 27 for ii in range(md.mesh.numberofelements): 28 28 for vertex in range(3): … … 33 33 idb_df[ii] = 2 34 34 #Covariance matrix 35 covPw = np.array([[0.75e10, 0.0],[0.0,0.5e10]])35 covPw = np.array([[0.75e10, 0.0], [0.0, 0.5e10]]) 36 36 37 37 # Stochastic forcing 38 38 md.stochasticforcing.isstochasticforcing = 1 39 md.stochasticforcing.fields 40 md.stochasticforcing.defaultdimension 41 md.stochasticforcing.default_id = idb_df-1 #NOTE: offset because of 1-based vertex indexing42 md.stochasticforcing.covariance = covPw# global covariance43 md.stochasticforcing.randomflag = 0# determines true/false randomness39 md.stochasticforcing.fields = ['FrictionWaterPressure'] 40 md.stochasticforcing.defaultdimension = 2 41 md.stochasticforcing.default_id = idb_df 42 md.stochasticforcing.covariance = covPw # global covariance 43 md.stochasticforcing.randomflag = 0 # determines true/false randomness 44 44 45 md.transient.issmb = 0;46 md.transient.ismasstransport = 1;47 md.transient.isstressbalance = 1;48 md.transient.isthermal = 0;49 md.transient.isgroundingline = 0;45 md.transient.issmb = 0 46 md.transient.ismasstransport = 1 47 md.transient.isstressbalance = 1 48 md.transient.isthermal = 0 49 md.transient.isgroundingline = 0 50 50 51 51 md.transient.requested_outputs = ['default', 'FrictionWaterPressure'] 52 52 md.timestepping.start_time = 0 53 md.timestepping.time_step 53 md.timestepping.time_step = 1 54 54 md.timestepping.final_time = 5 55 md.cluster = generic('name', gethostname(),'np',3)55 md.cluster = generic('name', gethostname(), 'np', 3) 56 56 md = solve(md, 'Transient') 57 57 58 58 # Fields and tolerances to track changes 59 field_names = [ 60 'Vx1','Vy1','Vel1','Thickness1','FrictionWaterPressure1', 61 'Vx2','Vy2','Vel2','Thickness2','FrictionWaterPressure2', 62 'Vx10','Vy10','Vel10','Thickness10','FrictionWaterPressure10' 63 ] 59 field_names = ['Vx1', 'Vy1', 'Vel1', 'Thickness1', 'FrictionWaterPressure1', 60 'Vx2', 'Vy2', 'Vel2', 'Thickness2', 'FrictionWaterPressure2', 61 'Vx10', 'Vy10', 'Vel10', 'Thickness10', 'FrictionWaterPressure10'] 64 62 65 field_tolerances = [ 66 2e-10,2e-10,2e-10,2e-10,2e-10, 67 4e-10,4e-10,4e-10,4e-10,4e-10, 68 8e-10,8e-10,8e-10,8e-10,8e-10 69 ] 63 field_tolerances = [2e-10, 2e-10, 2e-10, 2e-10, 2e-10, 64 4e-10, 4e-10, 4e-10, 4e-10, 4e-10, 65 8e-10, 8e-10, 8e-10, 8e-10, 8e-10] 70 66 71 field_values = [ 72 md.results.TransientSolution[0].Vx, 73 md.results.TransientSolution[0].Vy, 74 md.results.TransientSolution[0].Vel, 75 md.results.TransientSolution[0].Thickness, 76 md.results.TransientSolution[0].FrictionWaterPressure, 77 md.results.TransientSolution[1].Vx, 78 md.results.TransientSolution[1].Vy, 79 md.results.TransientSolution[1].Vel, 80 md.results.TransientSolution[1].Thickness, 81 md.results.TransientSolution[1].FrictionWaterPressure, 82 md.results.TransientSolution[4].Vx, 83 md.results.TransientSolution[4].Vy, 84 md.results.TransientSolution[4].Vel, 85 md.results.TransientSolution[4].Thickness, 86 md.results.TransientSolution[4].FrictionWaterPressure 87 ] 88 89 90 91 67 field_values = [md.results.TransientSolution[0].Vx, 68 md.results.TransientSolution[0].Vy, 69 md.results.TransientSolution[0].Vel, 70 md.results.TransientSolution[0].Thickness, 71 md.results.TransientSolution[0].FrictionWaterPressure, 72 md.results.TransientSolution[1].Vx, 73 md.results.TransientSolution[1].Vy, 74 md.results.TransientSolution[1].Vel, 75 md.results.TransientSolution[1].Thickness, 76 md.results.TransientSolution[1].FrictionWaterPressure, 77 md.results.TransientSolution[4].Vx, 78 md.results.TransientSolution[4].Vy, 79 md.results.TransientSolution[4].Vel, 80 md.results.TransientSolution[4].Thickness, 81 md.results.TransientSolution[4].FrictionWaterPressure]
Note:
See TracChangeset
for help on using the changeset viewer.