Changeset 27232 for issm/trunk/test/NightlyRun/test2006.py
- Timestamp:
- 08/25/22 16:50:29 (3 years ago)
- Location:
- issm/trunk
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk
-
issm/trunk/test
- Property svn:mergeinfo changed
/issm/trunk-jpl/test merged: 27036,27115-27118,27133,27142,27160,27167,27183,27186,27196,27202,27209,27229
- Property svn:mergeinfo changed
-
issm/trunk/test/NightlyRun/test2006.py
r27035 r27232 18 18 md = model() 19 19 md.cluster = generic('name', gethostname(), 'np', 5) 20 md.mesh = gmshplanet('radius', 6.371012 * 1e3, 'resolution', 700.) #700 km resolution mesh 20 21 21 22 # Load precomputed mesh … … 163 164 for i in range(len(locations)): 164 165 if i == 0: 165 md.outputdefinition.definitions = [ 166 nodalvalue( 167 'name', 'SNode', 168 'definitionstring', 'Outputdefinition1', 169 'model_string', 'Sealevel', 170 'node', locations[i] 171 ) 172 ] 166 md.outputdefinition.definitions = [nodalvalue('name', 'SNode', 167 'definitionstring', 'Outputdefinition1', 168 'model_string', 'Sealevel', 169 'node', locations[i])] 173 170 else: 174 171 md.outputdefinition.definitions.append( 175 nodalvalue( 176 'name', 'SNode', 177 'definitionstring', 'Outputdefinition' + str(i + 1), 178 'model_string', 'Sealevel', 179 'node', locations[i] 180 ) 181 ) 172 nodalvalue('name', 'SNode', 173 'definitionstring', 'Outputdefinition' + str(i + 1), 174 'model_string', 'Sealevel', 175 'node', locations[i])) 182 176 #}}} 183 177 184 178 # Algorithm #{{{ 185 179 md.qmu.method = dakota_method.dakota_method('nond_samp') 186 md.qmu.method = dmeth_params_set( 187 md.qmu.method, 188 'seed', 1234, 189 'samples', 10, 190 'sample_type', 'random' 191 ) 180 md.qmu.method = dmeth_params_set(md.qmu.method, 181 'seed', 1234, 182 'samples', 10, 183 'sample_type', 'random') 192 184 md.qmu.output = 1 193 185 #}}} … … 235 227 236 228 # Compare statistics with our own here 237 svalues = mds.results.StatisticsSolution[-1].SealevelSamples # all values at locations229 svalues = mds.results.StatisticsSolution[-1].SealevelSamples # all values at locations 238 230 239 231 dvalues = np.zeros((md.qmu.method.params.samples, len(locations)))
Note:
See TracChangeset
for help on using the changeset viewer.