Changeset 24179
- Timestamp:
- 09/29/19 14:08:28 (5 years ago)
- Location:
- issm/trunk-jpl/src/m/partition
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/partition/AreaAverageOntoPartition.m
r17686 r24179 29 29 %finally, project vector: 30 30 vector=project2d(md3d,vector,layer); 31 md.qmu. partition=project2d(md3d,md3d.qmu.partition,layer);31 md.qmu.vpartition=project2d(md3d,md3d.qmu.vpartition,layer); 32 32 end 33 33 34 34 %ok, first check that part is Matlab indexed 35 part=md.qmu. partition+1;35 part=md.qmu.vpartition+1; 36 36 37 37 %some check: 38 38 if md.qmu.numberofpartitions~=max(part), 39 error('AreaAverageOntoPartition error message: ''npart'' should be equal to max(md.qmu. partition)');39 error('AreaAverageOntoPartition error message: ''npart'' should be equal to max(md.qmu.vpartition)'); 40 40 end 41 41 -
issm/trunk-jpl/src/m/partition/AreaAverageOntoPartition.py
r23465 r24179 32 32 #finally, project vector: 33 33 vector = project2d(md3d,vector,layer) 34 md.qmu. partition = project2d(md3d,md3d.qmu.partition,layer)34 md.qmu.vpartition = project2d(md3d,md3d.qmu.vpartition,layer) 35 35 36 36 37 37 #ok, first check that part is Matlab indexed 38 part = (md.qmu. partition).copy()38 part = (md.qmu.vpartition).copy() 39 39 part = part.flatten() + 1 40 40 41 41 #some check: 42 42 if md.qmu.numberofpartitions != max(part): 43 raise RuntimeError('AreaAverageOntoPartition error message: ''npart'' should be equal to max(md.qmu. partition)')43 raise RuntimeError('AreaAverageOntoPartition error message: ''npart'' should be equal to max(md.qmu.vpartition)') 44 44 45 45
Note:
See TracChangeset
for help on using the changeset viewer.