Changeset 24187


Ignore:
Timestamp:
09/30/19 13:47:14 (6 years ago)
Author:
schlegel
Message:

CHG: element qmu partition needs to be a column

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/test/NightlyRun/test244.py

    r24174 r24187  
    7070md.qmu.numberofpartitions = md.mesh.numberofelements
    7171md = partitioner(md, 'package', 'linear','type','element')
    72 md.qmu.epartition = (md.qmu.epartition - 1).reshape(-1, 1).T
     72md.qmu.epartition = (md.qmu.epartition - 1)
    7373
    7474#variables
     
    7979for pa in range(np.size(mint_on_partition)):
    8080    vi = np.where(md.qmu.epartition == pa)
    81     mint = telms[vi] * 1.05
     81    mint = telms[0,vi] * 1.05
    8282    pos = np.where(mint < Tmin)
    8383    mint[pos] = Tmin
    84     mint_on_partition[pa] = max(mint / telms[vi])
     84    mint_on_partition[pa] = max(mint / telms[0,vi])
    8585
    8686mint_on_partition[np.where(np.isnan(mint_on_partition))] = 10**-10
Note: See TracChangeset for help on using the changeset viewer.