Ignore:
Timestamp:
08/25/22 16:50:29 (3 years ago)
Author:
Mathieu Morlighem
Message:

merged trunk-jpl and trunk for revision 27230

Location:
issm/trunk
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk

  • issm/trunk/test

  • issm/trunk/test/NightlyRun/test2006.py

    r27035 r27232  
    1818md = model()
    1919md.cluster = generic('name', gethostname(), 'np', 5)
     20md.mesh = gmshplanet('radius', 6.371012 * 1e3, 'resolution', 700.)  #700 km resolution mesh
    2021
    2122# Load precomputed mesh
     
    163164for i in range(len(locations)):
    164165    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])]
    173170    else:
    174171        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]))
    182176#}}}
    183177
    184178# Algorithm #{{{
    185179md.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 )
     180md.qmu.method = dmeth_params_set(md.qmu.method,
     181                                 'seed', 1234,
     182                                 'samples', 10,
     183                                 'sample_type', 'random')
    192184md.qmu.output = 1
    193185#}}}
     
    235227
    236228# Compare statistics with our own here
    237 svalues = mds.results.StatisticsSolution[-1].SealevelSamples # all values at locations
     229svalues = mds.results.StatisticsSolution[-1].SealevelSamples  # all values at locations
    238230
    239231dvalues = np.zeros((md.qmu.method.params.samples, len(locations)))
Note: See TracChangeset for help on using the changeset viewer.